1

時間の経過に伴うサイズの変化を調査しています。年に対する多変量線形モデルに 5 つのサイズ変数があります。car パッケージの Anova() 関数を使用して、各サイズ メジャーの勾配が等しいかどうかをテストしました。それらは大きく異なることが判明したので、サイズ測定値の組み合わせごとにペアワイズ検定を実行したかったのですが (これらは応答変数であることを思い出してください)、その方法を理解することができませんでした。これまでのところ、私のコードは以下のコードとほぼ同じですが、カテゴリ変数の代わりに、説明変数として数値を使用しています。

names(iris) <- c("SL", "SW", "PL", "PW", "SPP")
mod.iris <- lm(cbind(SL, SW, PL, PW) ~ SPP, data=iris)
summary(mod.iris)

manova.iris <- Anova(mod.iris)
summary(manova.iris)

したがって、この例を使用して、SL から SW、SL から PL、SL から PW などの傾きが有意に異なるかどうかをテストし、組み合わせごとにそれを行い、(オプションで) p 値を調整します。多変量回帰の勾配の事後検定を探しています。

ありがとう、

ガボール

4

1 に答える 1

1

勾配 (および切片) は に含まれていmod.iris$coefficientsます。

SL の勾配を SW の勾配と比較するには、次のようにします。

mod.iris$coefficients[1,'SL'] / mod.iris$coefficients[1,'S']

SL と SW の勾配を比較するには、次のようにします。

mod.iris$coefficients[1,'SL'] / iris$SW
于 2014-11-14T08:48:44.377 に答える