ฉันจะได้รับ ANOVA โดยรวมที่สำคัญได้อย่างไร


18

ฉันแสดงด้วย ANOVA R และฉันก็มีความแตกต่างที่สำคัญ อย่างไรก็ตามเมื่อตรวจสอบว่าคู่ไหนมีความแตกต่างอย่างมีนัยสำคัญโดยใช้ขั้นตอนของ Tukey ฉันไม่ได้รับเลย สิ่งนี้จะเป็นไปได้อย่างไร

นี่คือรหัส:

fit5_snow<- lm(Response ~ Stimulus, data=audio_snow)
anova(fit5_snow)

> anova(fit5_snow)
Analysis of Variance Table

Response: Response
          Df Sum Sq Mean Sq F value  Pr(>F)  
Stimulus   5  73.79 14.7578  2.6308 0.02929 *
Residuals 84 471.20  5.6095                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

df<-df.residual(fit5_snow)
MSerror<-deviance(fit5_snow)/df

comparison <-  HSD.test(audio_snow$Response, audio_snow$Stimulus, df, MSerror, group=FALSE)

> comparison <-  HSD.test(audio_snow$Response, audio_snow$Stimulus, df, MSerror, group=FALSE)

Study:

HSD Test for audio_snow$Response 

Mean Square Error:  5.609524 

audio_snow$Stimulus,  means

                audio_snow.Response   std.err replication
snow_dry_leaves            4.933333 0.6208034          15
snow_gravel                6.866667 0.5679258          15
snow_metal                 6.333333 0.5662463          15
snow_sand                  6.733333 0.5114561          15
snow_snow                  7.333333 0.5989409          15
snow_wood                  5.066667 0.7713110          15

alpha: 0.05 ; Df Error: 84 
Critical Value of Studentized Range: 4.124617 

Comparison between treatments means

                              Difference   pvalue sig        LCL      UCL
snow_gravel - snow_dry_leaves  1.9333333 0.232848     -0.5889913 4.455658
snow_metal - snow_dry_leaves   1.4000000 0.588616     -1.1223246 3.922325
snow_sand - snow_dry_leaves    1.8000000 0.307012     -0.7223246 4.322325
snow_snow - snow_dry_leaves    2.4000000 0.071587   . -0.1223246 4.922325
snow_wood - snow_dry_leaves    0.1333333 0.999987     -2.3889913 2.655658
snow_gravel - snow_metal       0.5333333 0.989528     -1.9889913 3.055658
snow_gravel - snow_sand        0.1333333 0.999987     -2.3889913 2.655658
snow_snow - snow_gravel        0.4666667 0.994348     -2.0556579 2.988991
snow_gravel - snow_wood        1.8000000 0.307012     -0.7223246 4.322325
snow_sand - snow_metal         0.4000000 0.997266     -2.1223246 2.922325
snow_snow - snow_metal         1.0000000 0.855987     -1.5223246 3.522325
snow_metal - snow_wood         1.2666667 0.687424     -1.2556579 3.788991
snow_snow - snow_sand          0.6000000 0.982179     -1.9223246 3.122325
snow_sand - snow_wood          1.6666667 0.393171     -0.8556579 4.188991
snow_snow - snow_wood          2.2666667 0.103505     -0.2556579 4.788991

คุณให้ข้อมูลได้หรือไม่
ttnphns

1
ฉันพบการตอบกลับคำถามนี้ในภายหลังstats.stackexchange.com/questions/74174/… (ทำเครื่องหมายว่าทำซ้ำเธรดนี้) เป็นประโยชน์อย่างยิ่ง
อะมีบาพูดว่า Reinstate Monica

คำตอบ:


2

ทำไมมันเป็นไปไม่ได้?

การทดสอบโดยรวมและการทดสอบแบบจับคู่ถามคำถามที่แตกต่างกันเพื่อให้พวกเขาสามารถได้รับคำตอบที่แตกต่างกัน


1
คุณช่วยพูดเพิ่มอีกหน่อยได้ไหม
rolando2

2
ANOVA โดยรวมถามคำถามเกี่ยวกับตัวแปรอิสระทั้งหมดและความสัมพันธ์ของมัน (หรือขาดมัน) กับตัวแปรตาม การเปรียบเทียบค่าตามจำนวนคู่ถามเกี่ยวกับความแตกต่างระหว่างคู่ จากนั้นค่า p จะดูที่ sig เชิงสถิติ ของแต่ละรายการด้วยการปรับค่าจำนวนคู่สำหรับการเปรียบเทียบหลาย ๆ แบบ (ในกรณีนี้ใช้วิธี HSD ของ Tukey)
Peter Flom - Reinstate Monica

1
ขอบคุณปีเตอร์ อาจน้อยกว่าที่พวกเขาถาม "คำถามที่แตกต่าง" และการปรับปรุงเพิ่มเติมสำหรับการเปรียบเทียบหลายอย่างที่บัญชีสำหรับผลลัพธ์ที่แตกต่างกัน
rolando2

17

นี่คือสาเหตุหลักมาจากความไวของ ANOVA (มากกว่าความไวการทดสอบแบบคู่) จากนั้น ANOVA จะตรวจจับความแปรปรวนที่ต่ำกว่ารอบ ๆ ค่าเฉลี่ยเมื่อการทดสอบแบบคู่ตามลำดับนั้นแยกความแตกต่างระหว่างค่าเฉลี่ยของทั้งคู่ไม่ได้ การวิเคราะห์จะต้องมุ่งเน้นไปที่ความแตกต่างและคุณสามารถมีความยืดหยุ่นมากขึ้นในการวิเคราะห์โพสต์เฉพาะกิจโดยทราบว่าคุณเพิ่งพบว่ามีความแตกต่างในค่าเฉลี่ย อย่าลืมตรวจสอบสมมติฐานของ ANOVA

ในทางกลับกันมีบางหัวข้อที่เกี่ยวข้องกับการใช้การทดสอบแบบจับคู่โดยไม่ใช้ ANOVA: เราจำเป็นต้องมีการทดสอบระดับโลกก่อนการทดสอบแบบโพสต์เฉพาะกิจหรือไม่?


ไม่จำเป็นต้องมีการทดสอบทั่วโลกก่อนทำการเปรียบเทียบ Tukey hsd เนื่องจาก hsd จะควบคุมอัตราความผิดพลาด Type I ฉันเกลียดที่จะเรียกพวกเขาว่าโพสต์ hocs เพราะพวกเขาควรได้รับการวางแผนมาก่อน
David Lane
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.