얌생 방목장

새루벵

얌생 — Thu, 12 Nov 2009 13:12:14 GMT

하하. 오랜만입니다.
그리고 이 녀석도 오랜만이네요.

색칠까지 한 건 정말, 기억도 안 나네요 -_-;;
덕분에 색칠을 어떻게 해야할지 이것저것 실험해봐야 했습니다.

이번 주말에는 좀 쉬어야겠습니다 ^^;;

tag : 마젠시아, 루벵

Regression Assignment by SAS

얌생 — Wed, 11 Nov 2009 17:21:57 GMT

data hw;
input Y X1 X2 X3 X4 X5 X6;
cards;
443 49 79 76 8 15 205
290 27 70 31 6 6 129
676 115 92 130 0 9 339
536 92 62 92 5 8 247
481 67 42 94 16 3 202
296 31 54 34 14 11 119
453 105 60 47 5 10 212
617 114 85 84 17 20 285
514 98 72 71 12 -1 242
400 15 59 99 15 11 174
473 62 62 81 9 1 207
157 25 11 7 9 9 45
440 45 65 84 19 13 195
480 92 75 63 9 20 232
316 27 26 82 4 17 134
530 111 52 93 11 13 256
610 78 102 84 5 7 266
617 106 87 82 18 7 276
600 97 98 71 12 8 266
480 67 65 62 13 12 196
279 38 26 44 10 8 110
446 56 32 99 16 8 188
450 54 100 50 11 15 205
335 53 55 60 8 0 170
459 61 53 79 6 5 193
630 60 108 104 17 8 273
483 83 78 71 11 8 233
617 74 125 66 16 4 265
605 89 121 71 8 8 283
388 64 30 81 10 10 176
351 34 44 65 7 9 143
366 71 34 56 8 9 162
493 88 30 87 13 0 207
648 112 105 123 5 12 340
449 57 69 72 5 4 200
340 61 35 55 13 0 152
292 29 45 47 13 13 123
688 82 105 81 20 9 268
408 80 55 61 11 1 197
461 82 88 54 14 7 225
;
run;

Proc reg data = hw;
Model Y = X1 X2 X3 X4 X5 X6;
Plot Y * X1; symbol i = rl value = dot;
Plot (student. h.) * obs. ;
Plot (cookd. dffits.) * obs.;
Output out = rout r = r h = h cookd = cookd dffits = dfits;
Run;

Data rout1;
Set rout;
di2 = r ** 2 / 26813; *SSE=26813
PFUN = h / (1-h);
RFUN = 2 / (1-h) * di2 / (1-di2);
HIM = PFUN + RFUN;
obs = _N_;
Run;

Proc print;
Var r cookd dfits HIM;
Run;

Proc gplot;
Plot HIM * obs PFUN * RFUN = obs;
Symbol i = none value = dot;
Run;

Number of Observations Read	40
Number of Observations Used	40

Analysis of Variance
Source	DF	Sum ofSquares	MeanSquare	F Value	Pr > F
Model	6	581506	96918	119.28	<.0001
Error	33	26813	812.51346
Corrected Total	39	608319

Root MSE	28.50462	R-Square	0.9559
Dependent Mean	464.92500	Adj R-Sq	0.9479
Coeff Var	6.13102

Parameter Estimates
Variable	DF	ParameterEstimate	StandardError	t Value	Pr > \|t\|
Intercept	1	35.17718	22.07552	1.59	0.1206
X1	1	2.85471	5.33718	0.53	0.5963
X2	1	3.27531	5.33304	0.61	0.5433
X3	1	3.18628	5.28872	0.60	0.5510
X4	1	3.18778	0.99180	3.21	0.0029
X5	1	-0.66767	0.89337	-0.75	0.4601
X6	1	-1.16584	5.32170	-0.22	0.8279

OBS	r	cookd	dfits	HIM	OBS	r	cookd	dfits	HIM
1	-9.4554	0.00214	-0.12082	0.12902	21	11.6951	0.00385	0.16216	0.15213
2	-15.0285	0.01243	-0.29196	0.27151	22	4.2217	0.00083	0.07519	0.21964
3	-1.7860	0.00047	-0.05676	0.54729	23	-12.2306	0.00580	-0.19900	0.19913
4	19.3467	0.00740	0.22581	0.13320	24	-50.1057	0.08219	-0.79112	0.40021
5	3.9815	0.00128	0.09324	0.34372	25	33.5962	0.02798	0.44651	0.22420
6	-11.4235	0.00538	-0.19157	0.20760	26	7.8561	0.00345	0.15318	0.25919
7	9.6997	0.00828	0.23770	0.37610	27	-28.9024	0.00570	-0.19995	0.10413
8	1.7616	0.00043	0.05389	0.51715	28	11.4789	0.01053	0.26826	0.35263
9	-19.7755	0.01369	-0.30743	0.20475	29	2.9861	0.00040	0.05192	0.21019
10	-24.2996	0.05459	-0.61823	0.44303	30	-6.2403	0.00074	-0.07114	0.10212
11	12.9788	0.00498	0.18460	0.16125	31	17.9501	0.00900	0.24882	0.16709
12	21.9044	0.04371	0.55153	0.42659	32	7.7186	0.00111	0.08701	0.10173
13	-28.7317	0.03443	-0.49262	0.27413	33	31.0298	0.06892	0.70066	0.40826
14	-9.0569	0.00576	-0.19820	0.31403	34	-54.2670	0.20722	-1.28147	0.62871
15	12.1341	0.02595	0.42155	0.63695	35	15.5949	0.00704	0.21979	0.16492
16	-16.6217	0.01476	-0.31859	0.27014	36	-23.4297	0.01683	-0.34203	0.19956
17	49.2736	0.08979	0.82601	0.41316	37	-12.4714	0.00545	-0.19306	0.18406
18	2.0612	0.00014	0.03109	0.16415	38	71.4380	0.24755	1.48031	0.80120
19	17.9102	0.00871	0.24476	0.16340	39	-34.7872	0.02553	-0.42713	0.21305
20	38.1876	0.02640	0.43653	0.21998	40	-46.1918	0.04033	-0.54767	0.28775

scatter plot

cook's distance

studentized residual

DFITS

hadi's influence measure

이건 뭐다냐...
게다가 leverage values는 어디로 갔냐 -_-;

ㅎㅇ

얌생 — Sat, 31 Oct 2009 13:46:34 GMT

오랜만에.

긴 생머리를 그린 게 오랜만이라는 말이기도 하지만 그림 자체도 오랜만입니다.

2005년 4월 14일 개설.

얌생 — Mon, 26 Oct 2009 09:35:39 GMT

600포스트 기념 포스팅.
601포스트가 되었다.

무의미하다.

숫자가 의미를 갖는 것은
숫자에 의미를 부여하기 때문이다.
다른 말로 숫자는 의미를 가지지 않는다.
이 이글루가 나에게 무슨 의미를 가지냐는 의문은
숫자가 나에게 무슨 의미를 가지냐는 의문의 답과 동일하다.
현상에 의미를 부여하지 않으면 의미가 없다.
의미는 현상에 대한 해석이다.

해석은 권력이다.

타인의 해석에 기대는 것은
앎을 타자화하여 권력에 굴복하게 한다.
생각하는대로 살지 않으면 사는대로 생각하게 된다.
어디선가 들은 말은 위와 같은 맥락일 것이다.
내 권력이 아닌 이글루는 무용하다.

이글루를 그만둘까?

-하고 고민하게 되는 것은,
아마 이런 생각들 때문일 것이다.

결국 4개.

얌생 — Wed, 21 Oct 2009 14:46:12 GMT

아뇨 그냥 그렇다구요.
하하.

오랜만에 사진을 올리네요.
괜히 폰으로 찍어놓고 방치하다가 시험기간이 돼서야 방출.

시험기간이니까요.

오랜만에 다시 다니게 된(사진 찍을 당시에는) 학교.

롯데리아 쟁반에 깔려있던 종이.
'기대이상'으로 깨끗하다는 말은 깨끗하다는 말이 아니잖아요.
라는 생각으로 찍었습니다. 오지랖도 넓지.

추석 때 찍은 보름달. 새벽 2시쯤?
이노믜 폰 카메라는 확대 기능도 없다냐.

가을 하늘이 참 맑고 높구나, 하고 집으로 들어가기 전에 찍은 사진.
오른쪽에 어두컴컴하고 북실북실한 건 아마 저...일걸요?

형이 받은 케익. 네모나고 길쭉한 상자라 롤케익이 들어있을 줄 알았는데 믿음을 저버리다니.
꽃모양으로 장식돼 있어서 이뻤습니다. 맛있었고요.
시험공부 하던 도중이라 밑에 깔린 건 강의 노트.

아... 이건 방금 찍긴 했는데 당최 왜 찍었는지 모르겠네요.
생각컨데 볼펜으로 시험공부한 걸 꿰뚫어 찢어발기고 싶은 마음에 하지 않았을까요.
하하하.

얌생님의 말:
제
젠장!
제시카, 환상적 몸매 과시
라고 해서 들어가봤더니
제시카 고메즈야
ㅠㅠ

tag : 꿰뚫어찢어발기리, 가을이구나, 우리시카

Sampling에 의한 모평균 추정 (결과 수정)

얌생 — Mon, 19 Oct 2009 09:26:01 GMT

우선 모집단은 제 최근의 이글루 포스팅에 달린 덧글의 개수입니다.
저의 답글이나 답글에 대한 답글, 링크신고 답글 등의 무용한 덧글은 제외했습니다.
샘플은 '운수 좋은 날'을 포함한 그 이전 10개의 포스팅을 추출하였습니다.
제가 sample을 뽑기로 한 지금 이 시점은 의도한 것이 아니므로 random하게 추출했다고 간주할 수 있습니다.

sample size는 10, 각 sample의 값은 아래와 같습니다.
3 6 6 4 3 6 3 7 5 4

이와 같을 때, 제 최근의 이글루 포스팅에 달린 평균 댓글이 4개라고 할 수 있는가에 대한 가설을 검정해보겠습니다.

계산해보면 sample mean(평균)은 4.7, sample SD(standard deviation)은 1.494434118, 약 1.49입니다.
(sample의 SD이기 때문에 n인 10이 아니라 n-1인 9로 나누었습니다.)
그러므로 sampling distribution의 SE(standard error)는 sample SD를 sample size n의 제곱근으로 나눈 값, 즉 0.472581562(0.47)입니다.

위 값들을 토대로 95% CI(confidence interval)을 추정해보겠습니다.
모집단의 SD를 모르고, sample size가 충분히 크지(n>30) 않으므로 Normal distribution이 아닌 T distribution을 사용합니다.
degrees of freedom(자유도)는 n-1로 9이고 95%를 양측으로 추정하니 t(9,0.025)를 찾으면 2.262라는 값을 찾을 수 있습니다.
95% CI:(mean-t(9,0.025)*se , mean+t(9,0.025)*se) = (4.7-2.262*0.47 , 4.7+2.262*0.47) = (3.63 , 5.77)

그 다음은 드디어 가설 검정의 단계입니다.
앞의 굵은 글씨의 문장에서 H0(null hypothesis:귀무가설):population mean(모평균)=4 라는 걸 알 수 있습니다.
H1(alternative hypothesis:대립가설):population mean≠4. Rejection의 기준값은 t(9,0.025)=2.262
sample mean은 4.7이므로 4.7-4(대립가설의 모평균)/se=1.481225795(1.481)
1.481 < t(9,0.025)=2.262
보시는 바와 같이 null hypothesis에 대해
sample mean이 상당한(significant) 반증이 될 수 없으므로.
95% significance level로 null hypothesis를 기각할 수 없습니다(cannot reject H0).

결론은 제 최근의 평균 덧글은 4개를 벗어나지 못합니다!라고 95%로 추정할 수 있습니다.
아..... 좀 슬프다.

이상 '사회과학도를 위한 통계적 분석(영강ㅠㅠ)'이라는 강의의 복습이었습니다.

tag : 통계적분석

규다.

얌생 — Sun, 18 Oct 2009 02:53:22 GMT

미안.

tag : 규, ㅈㅅ

운수 좋은 날

얌생 — Tue, 13 Oct 2009 08:00:38 GMT

학교 길 건너 버스정류장에서 버스를 타고 집에 오려고 횡단보도를 건너려는데,
얼레-. 내가 도착하니 신호가 바뀌네. 후후.
횡단보도 건너서 버스정류장에서 기다리니,
바로 집에가는 버스가 오네? 와우.
정류장에 정차하려 슬금슬금 다가오는 버스를 보고,
웬 일로 저 버스에 사람이 빈 자리가 있구나. 올레!

아.....

교통카드 잔액이 남아있었다면 더없이 좋았겠구나...
게다가 이노믜 학생증 교통카드는 지하철이 아니면 충전도 안 됐지, 하하.
그리고 내일까지 할 과제가 두개나 있네.

tag : 버스를사왔는데왜먹지를못해

생활의 달인

얌생 — Sat, 10 Oct 2009 17:43:23 GMT

오늘은 명함을 날려서 지나가는 여성분의 팔에 맞추었습니다.
하늘이 놀라고 땅이 놀라고 그분도 놀라고 나도 놀랐습니다.
참으로 진귀한 경험이었습니다.

...........뭐야이거.

내가 말이야,

얌생 — Thu, 08 Oct 2009 03:45:40 GMT

지금까지 한 과제가 10개쯤 되거든?
근데 지금 해야하는 과제가 5개쯤 있네.

ㅎㅅㅎㅅㅎㅅㅎㅅㅎㅅㅎㅅㅎㅅㅎ

이번 주말은 과제만 하자.
...그리고 다음주는 시험공부만 하자.