4

死亡率データに対してテューキー検定を実行しようとしています。そこでは、死亡率が銅の量 (一元配置分散分析) と銅と温度の組み合わせ (二元配置分散分析) によって影響を受けるかどうかをテストしたいと考えています。 . これらは私の式です：

lm2<-lm(Mortality~Cu) 
anova(lm2) 
TukeyHSD(aov(Mortality~Cu)) 

lm2<-lm(Mortality~Cu+Temp+Cu:Temp) 
anova(lm2) 
TukeyHSD(aov(Mortality~Cu+Temp+Cu:Temp))

ANOVA は問題ありませんが、両方の Tukey について、次のエラー メッセージが表示されます。

    Error in TukeyHSD.aov(aov(Mortality ~ Cu + Temp + Cu:Temp)) : 
      no factors in the fitted model 
    In addition: Warning messages: 
    1: In replications(paste("~", xx), data = mf) : non-factors ignored: Cu 
    2: In replications(paste("~", xx), data = mf) : non-factors ignored: Temp 
    3: In replications(paste("~", xx), data = mf) : 
      non-factors ignored: Cu, Temp

どこかに要因があるはずだと他の投稿で読んだことがありますが、私のデータはすべて数値です! 私はかなり困惑しており、次に何をすべきかわかりません。

よろしくお願いします。

ランディル

4

1 に答える 1

6

TukeyHSDカテゴリ変数でのみ機能するため、式で因子を探します。Cu個々のバストの値は、数値としてコード化されているだけですか? もしそうなら、それを使用してください

fCu<-factor(Cu)
TukeyHSD(aov(Mortality~fCu))

それとも観察は継続的ですか？次に、それを間隔に分割して因数に変換する必要があります。あなたができる

fCu<-cut(Cu, breaks=4)
TukeyHSD(aov(Mortality~fCu))

方程式の右辺の変数のサブセットのみが因子である場合は、それらを のwhichパラメーターで明示的に指定する必要がありますTukeyHSD。したがって、カテゴリを使用しfCu、Temp連続数値変数である場合、次のことができます

TukeyHSD(aov(Mortality~fCu+Temp), which="fCu")

それでも他の列に関する警告が表示されるため、結果を解釈する方法がわかりません

于 2014-05-09T18:12:03.417 に答える