14

私はアンバランスなデザイン/サンプルを扱っており、もともと学んだaov(). lm()ANOVA テストでは、 を使用するのではなく、 を使用してフィッティングを使用するタイプ III 平方和を使用する必要があることがわかりましたaov()。

問題は、 を使用して事後テスト (特に Tukey の HSD) を取得することlm()です。私が行ったすべての調査simintでは、multcompパッケージでの使用は機能すると言われていましたが、更新されたので、そのコマンドは利用できないようです. aov()また、計算するために通過することに依存しているようです。

基本的に、R について私が見つけたすべての Tukey HSD テストaov()は、lm(). タイプ III の二乗和を取得するには、不均衡な設計に必要な以下を使用する必要があります。

mod<-lm(Snavg~StudentEthnicity*StudentGender)

Anova(mod, type="III")

を使用して mod で Tukey HSD テストを使用するにはどうすればよいlm()ですか? または逆に、Type III を使用して ANOVA を計算しても、Tukey HSD テストを実行できますか?

ありがとう！

4

4 に答える 4

9

試してHSD.testみるagricolae

library(agricolae)
data(sweetpotato)
model<-lm(yield~virus, data=sweetpotato)
comparison <- HSD.test(model,"virus", group=TRUE,
main="Yield of sweetpotato\nDealt with different virus")

出力

Study: Yield of sweetpotato
Dealt with different virus

HSD Test for yield 

Mean Square Error:  22.48917 

virus,  means

      yield  std.err replication
cc 24.40000 2.084067           3
fc 12.86667 1.246774           3
ff 36.33333 4.233727           3
oo 36.90000 2.482606           3

alpha: 0.05 ; Df Error: 8 
Critical Value of Studentized Range: 4.52881 

Honestly Significant Difference: 12.39967 

Means with the same letter are not significantly different.

Groups, Treatments and means
a        oo      36.9 
ab       ff      36.33333 
 bc      cc      24.4 
  c      fc      12.86667
于 2011-10-12T00:51:35.293 に答える
3

また、使用しているまたはモデルHSD.test()の構築方法についても非常に細心の注意を払っていることがわかりました。lm()aov()

HSD.test()次のコーディングのアイデアを使用した場合、データからの出力はありませんでしたlm():

    model<-lm(sweetpotato$yield ~ sweetpotato$virus)  
    out <- HSD.test(model,"virus", group=TRUE, console=TRUE)

出力は次のとおりでした。

    Name:  virus 
    sweetpotato$virus

同じロジックを使用した場合、出力は同様に悪かったaov()

    model<-aov(sweetpotato$yield ~ sweetpotato$virus)

の出力を取得するにHSD.test()はlm() (またはaov()モデルに使用する場合も)、MYaseen208 の回答に示されているロジックを厳密に使用して構築する必要があります。

    model <- lm(yield~virus, data=sweetpotato)

これが、 から適切な出力を得られない人に役立つことを願っていますHSD.test()。

于 2016-10-07T11:09:19.593 に答える
1

私は HSD.test が何も出力しないという同じ問題に悩まされていました。console=TRUE関数内に配置する必要があるため、自動的に出力されます。

例えば：

HSD.test(alturacrit.anova, "fator", console=TRUE).
Hope it helps!
于 2014-10-28T03:33:08.620 に答える