Microsoft Word - sbe13_anova.docx

Size: px

Start display at page:

Download "Microsoft Word - sbe13_anova.docx"

용전 차
5 years ago
Views:

1 실험설계개요어떤원인이반응에유의한영향을주고있는가를파악하고그영향이양적으로어느정도큰가를알아내고자실시함 ( 추정과검정 적은영향밖에미치지못하는요인 ( 오차 들은전체적으로어느정도도영향을주고있으며, 측정오차는어느정도인가를알아내고자실시함 ( 오차항추정 유의한영향을미치는원인들이어떠한조건을수있는가를알아내기위해서실시함 ( 최적화 가질때가장바람직한반응을얻을 용어 요인 (Factor: 요인은제어가능한변수중하나로반응에미치는영향을의미하며, 실험공정에서연구된다. 요인 ( 인자 은 X로표현된다다. 인자 ( 요인 는온도나시간처럼정량적일수있고또한다른기계, 다른작업자, 또는환경처럼정성적일수있다. 반응 (Response: 주어진수준에서인자들의수준조합에의해해측정되는 ( 생산되는 특성치이다. 수준 (Level: 인자수준이란실험에서조절될수있는인자값을의미한다. 교호작용 (Interaction: 두인자들사이의조합에서일어나는효과 기본원칙 Randomization 랜덤화 관심요인외에기타원인들의영향이실험결과에미치지않게게함 시간에따라변하는인자의효과나경향을실험을시간대에대해균일하게배치함으로써약화시킬수있다. Blocing 블록화 동일한성질을가진단위들의집합 (Bloc 실험전체를시간적혹은공간적으로분할하여 Bloc 을만들어어주면각 Bloc 내에서는실험환경이균일하게되어좀더좋은결과를얻을수있다.

2 Bloc 은실험계획시또다른독립변수( 요인 로취급해야한다. 실험이이틀에걸쳐서수행되었다고하면실험일을 Bloc 이라한다. Replication 반복화 실험조건을처음부터다시 setup 하여실험하는것. 실험의재현성을알아보기위한방법 : 실험결과의신뢰성을높일수있다. 오차변동을계산할수있다. ( 참고 Repetition ( 반복 : 동일실험조건에서하나의관측치로이용 즉시반복실험 => 평균관측치를 Confounding 교락 두개이상인자의효과가함께나타나고그효과를각각의인자에의한효과로분리해낼수없을경우그인자들은교락 (Confounding 되어어있다고한다. 구할필요가없는 인자교호작용이나고차의교호작용 (interaction effect 을블록또는주효과와교락시켜실험의효율을높일수있다. -요인 인설계 (4 개실험조건 요인 A 주효과 : 클럽효과수준 = (New, Old : 수준별평균의차이 요인 B 주효과 : 경기기각수준 = ( 밤, 낮 (AB 교호효과 : 두요인수준의결합조건의평균차이

One-way ANOVAA 일원분산분석개요 요인 ( 분류형, 질적설명변수 이하나 집단이 3

완전임의설계 데이터 기호 i = 처리변수수준첨자, j = 반복 첨자 n i = y : 처리변수

: 처리변수 i 번째수준의종속변수평균 y : 종속변수반응값의 총평균 모형 model MVUE : ˆ ˆ

3 One-way ANOVAA 일원분산분석개요 요인 ( 분류형, 질적설명변수 이하나 집단이 3 개이상인모집단평균차이검정 실험설계에서는 CRD completelyy Randomized Designn 완전임의설계 데이터 기호 i = 처리변수수준첨자, j = 반복 첨자 n i = y : 처리변수 i수준의반복수, n 처리변수 i 번째수준의 j 번째 = 실험총반복수 반복의종속변수결과값 y i. : 처리변수 i 번째수준의종속변수평균 y : 종속변수반응값의 총평균 모형 model MVUE : ˆ ˆ Y y, i y i e i., i,,.,, j,,., n i 가정 assumption e ~ iidn (0, 독립성 / 정규성 / 등분산성 3

4 집단간평균의차이? 귀무가설 :. ( 모든집단평균은동일하다 모든 0 ( 요인의 주효과는유의하지않음 요인의효과는유의하지않다. i 대립가설 : 적어도하나이상의집단 평균은유의하다. y 3 35 변동분해 Variation Decomposition 총변동자승합 SST (Total Sum of Squares SST ni i j ( Y j Y (33 SST ni i j ( Y ni i j Y ( Y Y i. ni i j ( Y i Y. ( ( ( = SSE (Error SS 집단내변동, 오차변동 ( = SSA (Between SS / Treatment SS 집단간변동, 처리변동, 요인 A 주효과 가정 e ~ iidn(0, 으로부터 y ~ iidn (, 이므로 i SSA ~ (, SSE ~ ( n 이고서로독립이다. SSA /( 그러므로 ~ F(, n SS SE /( n 4

5 분산분석표 ANOVA table 변동 자유도 자승합 요인변동 - SSA 오차변동 n-- SSE 총변동 n- SST 평균자승합 SSAA MSA ( SSEE MSE ( n F MSA TSS ~ MSE F( (, n 비즈니스표작성 집단 평균 (M 표준편차 (SD F-검정통계량유의확률포함 사후검정 (Post-hoc test 혹은다중 비교 (multiple comparison 분산분석의 F- 검정은단지귀무가설 검정하는것이다. 그러므로수준별차이 (pairwise: 예 : H 0 : 수준의선형결합대비 (contrast: 예 : 있는데이를사후검정혹은다중비교 ( 대비포함 라한다. 사후검정이므로비록 F- 검정결과와관계없이 ( 귀무가설을채택하더라도 시행하게된다. 이를조정된실험유의수준 (controlled experimental error rate 이라하고 ( 이다. 된다. 여기서 c 는가설수를의미한다. pairwise 다중비교의경우 c Fisher s Least Significant Difference pairwise ( 두수준별평균비교 검정에사용하나이는다중비교에해당되지는않는다. 두모집단평균차이검정에서통합분산 (pooled variance 대신 MSE 사용한다. Tuey HSD(honestly significant difference procedure Y Y j q 분포에대한표제공 MSE / n H H 0 0 : u : u u. u t u 즉 u 3 즉전체적인인차이를 u u 3 가있는지혹은 의차이가있는지검정할필요가 다중비교에서는여러개의가설을동시에검정하므로유의수준을조정해야한다. t ( t / 가 c 5

차이가가장큰것을비교해가면서유의수준을 ( 검정한다. r 은검정단계순서이다. 귀무가설을기각할확률이매우높아자주사용하지않는다.

( 사회과학분야야 Dunnett s procedure 처리효과의수준하나가 control ( 실험집단 인경우 ( 예 : placebo 집단,

xls 유람선운영 S 사는배의크기에따른고객만족도점수의차이가가있는지알아보기위한조사한것이다.

6 가장보수적인 ( 귀무가설기각하지않음 방법이다. ( 자연과학분야 Duncan Multiple range test Tuey 방법과매우유사하나수준별표본평균을크기순으로로나열하여 차이가가장큰것을비교해가면서유의수준을 ( 검정한다. r 은검정단계순서이다. 귀무가설을기각할확률이매우높아자주사용하지않는다. Scheffe s S method r 으로로조정해가면서 대비 (contrast 까지고려한다중 비교방법이다. ( 사회과학분야야 Dunnett s procedure 처리효과의수준하나가 control ( 실험집단 인경우 ( 예 : placebo 집단, 교육을하지않는집단, 이전약투여집단 이집단과다른집단들을을 pairwise 비교할경우사용된다. 예제데이터 SHIPS.xls 유람선운영 S 사는배의크기에따른고객만족도점수의차이가가있는지알아보기위한조사한것이다. 배크기별로 8개유람선임의추출, 각유람선탑승고객의평가점수를평균한측정치이다. One-way ANOVA ods graphics on; proc glm data=ship; classs size; model score=size; means size / scheffe; run;quit; proc tabulate data=ship; class size; var score; table (size all, a (score*(mean std ; run; 고정된경우대응 변인의유의성검증 Type III 는다른 요인의효과가 6

7 F-검정결과배크기에따른만족도의차이는유의함 다중비교결과 ( 소형 > 중형 > ( 중형 > 대형 : 그러므로소형과대형배의승객의만족도차이만유의하다. Two-Way ANOVA 개요 요인 ( 분류형, 질적설명변수 이두개 요인개별효과를주효과 ( 각요인의수준별평균차이 두요인의수준결합에의해평균의차이가나는효과, 사회과학분석에서는일반적으로교호효과를사용하지않는다 ( 이유는설명이용이하지않음. 주효과를보기위한변동분해로만사용한다. 7

8 모형 model Y i b j ( ab e, i,,., a, j,,.,b 가정 assumption ~ iidn (0, e 독립성 / 정규성 / 등분산성성 통계적 가설 귀무가설 : 요인 A 의수준별평균의차이는없다 요인 A 주효과 유의하지 않음 귀무가설 : 요인 B 의수준별평균의차이는없다 귀무가설 : 모든 ( ab 0 요인 (A, B 의교호작용은 유의하지않음 변동분해 Variation Decomposition 총변동자승합 SST (Total Sum of Squares SST ni i j ( Y j Y SST a b i j ( Y Yi. Y. j Y SSE SSA SSB SSAB a b i j ( Yi. Y a b i j ( Y. j Y a b i j ( Y Y Y. j Y 분산분석표 ANOVA table 변동 자유도 자승합 평균자승합 (= 자승합 / 자유도 F 요인 A 주효과 a- SSA MSAA MSA/MSE~F 요인 A 주효과 b- SSB MSBB MSB/MSE~F (AB 교호효과 (a-(b- SSAB MSAB MSAB/MSE~F 오차변동 차이 SSE MSEE Total 총변동 n- SST 비즈니스표작성 요인 A 요인 B A F-F 검정통계량량 b 평균SD 요인 A 주효과과요인 B 주효과과요인 AB 교호효과 8

9 예제데이터 SHIPS.xls Height: 0 if step at the low (5.75" height, if at the high (.5" height Frequency: the rate off stepping. 0 if slow (4 steps/min, if medium ( steps/min, if high (8 steps/min Rest_HR: the resting heart rate of the subject before a trial, in beats per minute HR: the final heart rate of the subject after a trial, in beats per minute 요인 ( 계단높이, 올라가는 속도 가운동후심장박동 (HR 에영향을을미치나? One-way ANOVA ods graphics on; proc glm data=heart; classs frequency height; model hr=frequency height; means frequency / tuey; run;quit; proc tabulate data=ship; class size; var score; table (size all, a (score*(mean std ; run; 대체 빈도, 높이 모두유의함, 교호작용유의하지않음 9

10 높이가높으면심장박동이높아진다 빈도가많으면심장박동이높아진다. 빈도 = 만다른빈도와유의한차이있음 교호작용은유의하지않음 => ( 높이 =, 빈도 = 결합조건에서심장박동이가장높아진다. 예제다음은가솔린종류 (A, B, C 와첨가제 (,, 3, 4 에따른연비의차이가있는알아보기위하여동일차종 4 대를이용하여 ( 물론차에의한차이도있을것이라생각되지만무시한다 측정한결과이다. 가솔린 A B C 첨가제 공분산분석 ANCOVA (Analysis of Covaraince 공변량 (covariate 종속변수값에대한요인들의유의성검정을제대로하기위해해고려되는변량 일반적으로종속변수의실험전값이다. ( 예 교육효과에서의사전점수 공변량은관심의대상이아니라요인의유의성검정을정확하기위하여고려함 모형 Y b ( ab a i j x One-way ANOVA ods graphics on; proc glm data=heart; classs frequency height; model hr=resthr frequency height ; means frequency / tuey; run;quit; 0

11 공변량은유의함 빈도, 높이요인유의함, 교호효과유의하지않음

Microsoft Word - SBE2012_anova.docx

Microsoft Word - SBE2012_anova.docx 실험설계개요어떤원인이반응에유의한영향을주고있는가를파악하고그영향이양적으로어느정도큰가를알아내고자실시함 ( 추정과검정 ) 적은영향밖에미치지못하는요인 ( 오차 ) 들은전체적으로어느정도영향을주고있으며, 측정오차는어느정도인가를알아내고자실시함 ( 오차항추정 ) 유의한영향을미치는원인들이어떠한조건을가질때가장바람직한반응을얻을수있는가를알아내기위해서실시함 ( 최적화 ) 용어 요인