Câu hỏi của tôi có liên quan rất chặt chẽ với một bài trước Xác định hạn Lỗi () trong các biện pháp lặp đi lặp lại ANOVA trong R . Tuy nhiên, tôi muốn hiểu rõ hơn về cách xác định thuật ngữ lỗi.
Giả sử tôi có ANOVA lặp lại hai chiều, Yếu tố giữa hiệu ứng nhóm là Điều trị (kiểm soát so với giả dược), trong khi Thời gian là hiệu ứng trong nhóm được đo lặp lại hơn 4 lần (T1 ~ T4). ID bệnh nhân được ghi là Chủ đề. Ở đây tôi đã mượn dữ liệu từ một ví dụ từ hướng dẫn trong http://gjkerns.github.io/R/2012/01/20/power-sample-size.html để dữ liệu trông như thế này
Time Subject Method NDI
0min 1 Treat 51.01078
15min 1 Treat 47.12314
48hrs 1 Treat 26.63542
96hrs 1 Treat 20.78196
0min 2 Treat 42.61345
15min 2 Treat 32.77171
Để áp dụng ANOVA:
aovComp <- aov(NDI ~ Time*Method + Error(Subject/Time), theData)
summary(aovComp)
Error: Subject
Df Sum Sq Mean Sq F value Pr(>F)
Method 1 113 112.7 0.481 0.491
Residuals 58 13579 234.1
Error: Subject:Time
Df Sum Sq Mean Sq F value Pr(>F)
Time 3 13963 4654 103.789 < 2e-16 ***
Time:Method 3 1221 407 9.074 1.3e-05 ***
Residuals 174 7803 45
Tôi cũng đã thử thuật ngữ lỗi khác:
aovComp1 <- aov(NDI ~ Time*Method + Error(Subject), theData)
summary(aovComp1)
Error: Subject
Df Sum Sq Mean Sq F value Pr(>F)
Method 1 113 112.7 0.481 0.491
Residuals 58 13579 234.1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Time 3 13963 4654 103.789 < 2e-16 ***
Time:Method 3 1221 407 9.074 1.3e-05 ***
Residuals 174 7803 45
Ai đó có thể giúp tôi giải thích sự khác biệt giữa hai điều khoản lỗi này không? Nếu thuật ngữ đầu tiên là thuật ngữ chính xác, kết quả từ thuật ngữ lỗi thứ hai có nghĩa là gì?
Cập nhật bởi @amoeba: Hai kết quả đầu ra giống nhau nên dường như trong trường hợp này không có sự khác biệt, nhưng câu hỏi vẫn là sự khác biệt về nguyên tắc là gì . Được Error(subject)
và Error(subject/time)
luôn luôn là điều tương tự?
time
), Error(subject)
và Error(subject/time)
sẽ mang lại các giá trị F và p khác nhau cho time
.