본문내용
1 50 1 1 49 1 1 50 1 2 50 1 2 48 1 2 51
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc glm data=ex4;
class part operator;
model y=part operator part*operator;
random part operatoer part*operaotr;
test h=part operator e=part*operator;
run;
<축소모형>
data ex4;
input part operator y @@;
cards;
1 1 50 1 1 49 1 1 50 1 2 50 1 2 48 1 2 51
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc glm data=ex4;
class part operator;
model y=part operator;
random part operator;
run;
<분산성분 검정>
data ex4;
input part operator y @@;
cards;
1 1 50 1 1 49 1 1 50 1 2 50 1 2 48 1 2 51
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc varcomp data=ex4 method=type1;
class part operator;
model y=part operator;
run;
(1) 변량효과모형 :
분산성분모형 :
(2) ① 가설 :
② 유의수준 : 0.05
③ 검정통계량 : 7.85
④ 기각역
⑤ 결과 및 해석 : 페인트 브랜드는 F=7.85이고 p=0.0169이다. 유의수준 0.05에서 면적에 대한 페인트 브랜드의 효과는 유의하다.
(3)
3.5
(1)
반복이 있는 이원배치법 모형식
:
가변수를 이용한 회귀 모형식
:
(2)
data ex5;
input fer type y @@;
cards;
1 1 12 1 1 18 1 2 4 1 2 6 1 3 20 1 3 14
2 1 14 2 1 10 2 2 13 2 2 13 2 3 22 2 3 22
3 1 12 3 1 4 3 2 6 3 2 12 3 3 15 3 3 9
4 1 20 4 1 14 4 2 17 4 2 9 4 3 17 4 3 25
;
run;
proc anova data=ex5;
class fer type;
model y=fer type fer*type;
run;
(3)
y가 인자 a의 수준 1에 있을 때
-1 y가 인자 a의 수준 3에 있을 때
0 기타
y가 인자 a의 수준 2에 있을 때
-1 y가 인자 a의 수준 3에 있을 때
0 기타
y가 인자 a의 수준 1에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
y가 인자 a의 수준 2에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
y가 인자 a의 수준 3에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
data ex5_3;
input x1 x2 z1 z2 z3 y;
x1z1=x1*z1;
x1z2=x1*z2;
x1z3=x1*z3;
x2z1=x2*z1;
x2z2=x2*z2;
x2z3=x2*z3;
cards;
1 0 1 0 0 12
1 0 1 0 0 18
0 1 1 0 0 4
0 1 1 0 0 6
-1 -1 1 0 0 20
-1 -1 1 0 0 14
1 0 0 1 0 14
1 0 0 1 0 10
0 1 0 1 0 13
0 1 0 1 0 13
-1 -1 0 1 0 22
-1 -1 0 1 0 22
1 0 0 0 1 12
1 0 0 0 1 4
0 1 0 0 1 6
0 1 0 0 1 12
-1 -1 0 0 1 15
-1 -1 0 0 1 9
1 0 -1 -1 -1 20
1 0 -1 -1 -1 14
0 1 -1 -1 -1 17
0 1 -1 -1 -1 9
-1 -1 -1 -1 -1 17
-1 -1 -1 -1 -1 25
;
run;
proc reg data=ex5_3;
model y=x1 x2 z1 z2 z3 x1z1 x1z2 x1z3 x2z1 x2z2 x2z3;
run;
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc glm data=ex4;
class part operator;
model y=part operator part*operator;
random part operatoer part*operaotr;
test h=part operator e=part*operator;
run;
<축소모형>
data ex4;
input part operator y @@;
cards;
1 1 50 1 1 49 1 1 50 1 2 50 1 2 48 1 2 51
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc glm data=ex4;
class part operator;
model y=part operator;
random part operator;
run;
<분산성분 검정>
data ex4;
input part operator y @@;
cards;
1 1 50 1 1 49 1 1 50 1 2 50 1 2 48 1 2 51
2 1 52 2 1 52 2 1 51 2 2 51 2 2 51 2 2 51
3 1 53 3 1 50 3 1 50 3 2 54 3 2 52 3 2 51
4 1 49 4 1 51 4 1 50 4 2 48 4 2 50 4 2 51
5 1 48 5 1 49 5 1 48 5 2 48 5 2 49 5 2 48
6 1 52 6 1 50 6 1 50 6 2 52 6 2 50 6 2 50
7 1 51 7 1 51 7 1 51 7 2 51 7 2 50 7 2 50
8 1 52 8 1 50 8 1 49 8 2 53 8 2 48 8 2 50
9 1 50 9 1 51 9 1 50 9 2 51 9 2 48 9 2 49
10 1 47 10 1 56 10 1 49 10 2 46 10 2 47 10 2 48
run;
proc varcomp data=ex4 method=type1;
class part operator;
model y=part operator;
run;
(1) 변량효과모형 :
분산성분모형 :
(2) ① 가설 :
② 유의수준 : 0.05
③ 검정통계량 : 7.85
④ 기각역
⑤ 결과 및 해석 : 페인트 브랜드는 F=7.85이고 p=0.0169이다. 유의수준 0.05에서 면적에 대한 페인트 브랜드의 효과는 유의하다.
(3)
3.5
(1)
반복이 있는 이원배치법 모형식
:
가변수를 이용한 회귀 모형식
:
(2)
data ex5;
input fer type y @@;
cards;
1 1 12 1 1 18 1 2 4 1 2 6 1 3 20 1 3 14
2 1 14 2 1 10 2 2 13 2 2 13 2 3 22 2 3 22
3 1 12 3 1 4 3 2 6 3 2 12 3 3 15 3 3 9
4 1 20 4 1 14 4 2 17 4 2 9 4 3 17 4 3 25
;
run;
proc anova data=ex5;
class fer type;
model y=fer type fer*type;
run;
(3)
y가 인자 a의 수준 1에 있을 때
-1 y가 인자 a의 수준 3에 있을 때
0 기타
y가 인자 a의 수준 2에 있을 때
-1 y가 인자 a의 수준 3에 있을 때
0 기타
y가 인자 a의 수준 1에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
y가 인자 a의 수준 2에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
y가 인자 a의 수준 3에 있을 때
-1 y가 인자 a의 수준 4에 있을 때
0 기타
data ex5_3;
input x1 x2 z1 z2 z3 y;
x1z1=x1*z1;
x1z2=x1*z2;
x1z3=x1*z3;
x2z1=x2*z1;
x2z2=x2*z2;
x2z3=x2*z3;
cards;
1 0 1 0 0 12
1 0 1 0 0 18
0 1 1 0 0 4
0 1 1 0 0 6
-1 -1 1 0 0 20
-1 -1 1 0 0 14
1 0 0 1 0 14
1 0 0 1 0 10
0 1 0 1 0 13
0 1 0 1 0 13
-1 -1 0 1 0 22
-1 -1 0 1 0 22
1 0 0 0 1 12
1 0 0 0 1 4
0 1 0 0 1 6
0 1 0 0 1 12
-1 -1 0 0 1 15
-1 -1 0 0 1 9
1 0 -1 -1 -1 20
1 0 -1 -1 -1 14
0 1 -1 -1 -1 17
0 1 -1 -1 -1 9
-1 -1 -1 -1 -1 17
-1 -1 -1 -1 -1 25
;
run;
proc reg data=ex5_3;
model y=x1 x2 z1 z2 z3 x1z1 x1z2 x1z3 x2z1 x2z2 x2z3;
run;
소개글