R：一元配置ANOVA後のTUKEY

Question

私はRで一方向ANOVAを行ったが、Tukey post-hocを実行してお互いに異なる治療法を見つけようとするとエラーメッセージが表示され続ける。 （私は、AB、B、BCD ...など（結果がランク付けされることを希望）

データの詳細：

データ=

のx = "abh2" 6回の治療： "治療"

Y =水分読取値 "湿った"（処置当たりn = 63、合計= 378）

＃の結果は、私が事後（P < 0.05）に移動できることを示し

anov <- anova(lm(moist~treatment, data=abh2)) 

：：私はHSDテューキーを選んだ

Analysis of Variance Table 

Response: moist 
      Df Sum Sq Mean Sq F value Pr(>F)  
treatment 5 1706.3 341.27 25.911 < 2.2e-16 *** 

し、それを実行しようとした。私は、一方向ANOVAを実行しました2つの方法ではなく、両方のエラーメッセージを取得：

内蔵R機能：

私はMSerrorが

1だった

HSD.test(anov, "treatment", group=TRUE, console=TRUE) 
    # ERROR : Error in HSD.test(anov, "treatment", group = TRUE, console = TRUE) : 
argument "MSerror" is missing, with no default 

） "＃古いバージョンHSD.test（）"（しかし、私はちょうどagricolaeパッケージを更新しました：

Agricolaeパッケージに）

2）MSerror < -deviance（モデル）/

dfをだから私は試してみました：

HSD.test(anov, "treatment", MSerror=deviance(moist)/5, group=TRUE, console=TRUE) 
*but still* # ERROR: $ operator is invalid for atomic vectors 

ここから先に進んでください。それはかなり単純な問題のようですが、私はこれに数時間を費やしました！

多くの感謝:)

Answer 1

おかげで、それはその右のトラックに私を得ましたRはデータを認識していませんでした。

でも、私はこのコードを使用して問題を解決：

library(agricolae) 

model<-aov(moist~treatment, data=abh2) 

out <- HSD.test(model,"treatment", group=TRUE,console=TRUE) 

を私が最初に（Rの基本パッケージから）を使用したANOVAコマンドは、単に私がしようとしていたTukey検定によって理解されていなかったことが表示されます後に実行する（agricolaeパッケージで）。

私のためのテイクホームメッセージ：同じパッケージで関連する一連の分析を行うようにしてください！

p.s.これはp.valueを与えません。その情報のANOVA自体を実行します。

anova(model) 

Answer 2

次のコードを要因として、あなたの治療を指定してみてください：フィードバックアニー・クロードため

abh2$treatment <- factor(abh2$treatment)