Rでの反復測定ANOVAでのError()項の指定


10

Rでの2方向反復測定分散分析の誤差項の定義に問題があります。私のデータは、木から抽出されたコアに沿った3つの半径方向位置(内側、中央、および外側)の木材密度推定から構成されています。合計20の木の種、各種の6個体、および各木からの2つのコアがあります。

半径方向の位置が木材の密度に及ぼす影響をテストするために、個人間のばらつきを説明する誤差項を含む次の2元配置分散分析モデルを使用します。

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

ただし、エラー項の指定が適切かどうかはわかりません。コア内の変動性も考慮する必要がありますか?私にとって、この変動性は、私が興味を持っている主な要因である放射状の位置による変動と同じです。

繰り返し測定ANOVAでの誤差項の指定について読むことに少し時間を費やしましたが、実際には誤差項の指定に問題があります。私はこれでいくつかの助けに感謝します。

回答:


8

それはそのようになります

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

これはのサブジェクト内エラーの原因ですRadialposition。他の被験者内要因がある場合は、それらを(相互作用で)Radialpositionエラー分母に入れ+ Error(Individual/(Radialpostion*wiFactorA))ます。それが私の理解です。GLM欠落データがない場合、これはSPSSの繰り返し測定と一致します。


1
+1。詳細なチートシートは、こちらをご覧ください:dwoll.de/rexrepos/posts/anovaMixed.html。そしてここにも:personality-project.org/r/r.guide/r.anova.html
amoeba
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.