Chỉ định thuật ngữ Error () trong các số đo lặp lại ANOVA trong R


10

Tôi đang gặp vấn đề với việc xác định các thuật ngữ lỗi cho các biện pháp lặp lại hai chiều ANOVA trong R. Dữ liệu của tôi bao gồm các ước tính mật độ gỗ cho ba vị trí xuyên tâm (bên trong, giữa và bên ngoài) dọc theo lõi được chiết xuất từ ​​cây. Có tổng cộng 20 loài cây, 6 cá thể của mỗi loài và hai lõi từ mỗi cây.

Để kiểm tra ảnh hưởng của vị trí xuyên tâm đến mật độ gỗ, tôi sử dụng mô hình ANOVA hai chiều sau đây với một thuật ngữ lỗi chiếm tỷ lệ thay đổi giữa các cá nhân:

radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual), data=Radpos)

Tuy nhiên, tôi không chắc chắn nếu đặc điểm kỹ thuật của tôi về thuật ngữ lỗi là đủ. Tôi cũng có nghĩa vụ phải tính đến sự thay đổi trong lõi? Đối với tôi sự biến thiên này giống như do vị trí xuyên tâm là yếu tố chính mà tôi quan tâm.

Mặc dù tôi đã dành một chút thời gian để đọc về việc chỉ định thuật ngữ lỗi trong các biện pháp lặp lại ANOVA, tôi vẫn gặp vấn đề với việc chỉ định thực tế thuật ngữ lỗi. Tôi sẽ đánh giá cao một số trợ giúp với điều này.

Câu trả lời:


8

Nó sẽ là

 radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual/(Radialposition)), data=Radpos)
 summary(radpos.aov, type=3)

Đó là tài khoản cho các lỗi trong chủ đề của Radialposition. Nếu bạn có các yếu tố bên trong chủ đề khác, hãy ném chúng vào (trong một tương tác) với Radialpositionmẫu số Lỗi, như + Error(Individual/(Radialpostion*wiFactorA)). Đó là sự hiểu biết của tôi về nó. Điều đó phù hợp với các biện pháp lặp đi lặp lại của SPSS GLMnếu bạn không có dữ liệu bị thiếu.


1
+1. Xem bảng cheat chi tiết tại đây: dwoll.de/rexrepos/posts/anovaMixed.html . Và cũng ở đây: Personal-project.org/r/r.guide/r.anova.html
amip
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.