Hata () terimini tekrarlanan ölçümlerde belirtme ANOVA R


10

R'de iki yönlü tekrarlanan ölçümler ANOVA için hata terimlerini tanımlamakta sorun yaşıyorum Verilerim, bir ağaçtan çıkarılan bir çekirdek boyunca üç radyal konum (iç, orta ve dış) için ahşap yoğunluk tahminlerinden oluşuyor. Toplam 20 ağaç türü, her türden 6 kişi ve her ağaçtan iki çekirdek vardır.

Radyal konumun ahşap yoğunluğu üzerindeki etkisini test etmek için, bireyler arasındaki değişkenliği açıklayan bir hata terimiyle aşağıdaki iki yönlü ANOVA modelini kullanıyorum:

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

Ancak, hata terimi spesifikasyonumun yeterli olup olmadığından emin değilim. Bir çekirdek içindeki değişkenliği de hesaba katmam gerekir mi? Bana göre bu değişkenlik, ilgilendiğim ana faktör olan radyal konumdan kaynaklanıyor.

Yinelenen önlemler ANOVA'daki hata terimini belirtmek için biraz zaman ayırmış olmama rağmen, hala pratikte hata terimini belirtmekte sorun yaşıyorum. Bununla ilgili bazı yardımları takdir edeceğim.

Yanıtlar:


8

Olurdu

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

Bu konu içindeki hatadan sorumludur Radialposition. Başka bir özne faktörünüz varsa, bunları RadialpositionHata paydasında (etkileşimli olarak) içine atın + Error(Individual/(Radialpostion*wiFactorA)). Bu benim anlayışım. GLMEksik verileriniz yoksa bu, SPSS'nin tekrarlanan ölçümleriyle eşleşir .


1
+1. Ayrıntılı hile sayfasına bakın: dwoll.de/rexrepos/posts/anovaMixed.html . Ve ayrıca burada: personality-project.org/r/r.guide/r.anova.html
amip
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.