Specificare il termine Error () nelle misure ripetute ANOVA in R


10

Sto riscontrando problemi con la definizione dei termini di errore per una misura ripetuta a due vie ANOVA in R. I miei dati consistono in stime della densità del legno per tre posizioni radiali (interna, media ed esterna) lungo un nucleo estratto da un albero. Ci sono un totale di 20 specie di alberi, 6 individui di ogni specie e due nuclei per ogni albero.

Per testare l'effetto della posizione radiale sulla densità del legno, uso il seguente modello ANOVA a due vie con un termine di errore che tiene conto della variabilità tra gli individui:

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

Tuttavia, non sono sicuro che la mia specifica del termine di errore sia adeguata. Devo anche spiegare la variabilità all'interno di un nucleo? Per me questa variabilità è uguale a quella dovuta alla posizione radiale, che è il fattore principale a cui sono interessato.

Sebbene abbia dedicato del tempo alla lettura della specifica del termine di errore nelle misure ripetute ANOVA, ho ancora problemi con la specifica del termine di errore. Apprezzerò un po 'di aiuto in questo.

Risposte:


8

Sarebbe

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

Ciò spiega l'errore interno del soggetto Radialposition. Se hai altri fattori all'interno del soggetto, inseriscili (in un'interazione) con Radialpositionnel denominatore Errore, ad esempio + Error(Individual/(Radialpostion*wiFactorA)). Questa è la mia comprensione. Ciò corrisponde alle misure ripetute di SPSS GLMse non si dispone di dati mancanti.


1
+1. Vedi il cheat sheet dettagliato qui: dwoll.de/rexrepos/posts/anovaMixed.html . E anche qui: personality-project.org/r/r.guide/r.anova.html
amoeba,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.