タイプ III エラーと種と生息地内のペアワイズ比較を使用しながら、ホームレンジのサイズに対する種と生息地の影響を比較したいデータセットがあります。
データのサブセットを次に示します。
species<- c("a","b","c","c","b","c","b","b","a","b","c","c","a","a","b","b","a","a","b","c")
habitat<- c("x","x","x","y","y","y","x","x","y","z","y","y","z","z","x","x","y","y","z","z")
homerange<-c(6,5,7,8,9,4,3,5,6,9,3,6,6,7,8,9,5,6,7,8)
data1<-data.frame(cbind(species, habitat, homerange))
data1$homerange<-as.numeric(as.character(data1$homerange))
現在、私は 3 つの種のデータを分割し、それぞれについて別々の ANOVA を実行していますが、1 つの ANOVA で同時に種と生息地について尋ねる方が理にかなっていると思います。ある種に対して実行した ANOVA の例を次に示します。
data.species.a<-subset(data1, species=="a")
fit<-aov(homerange ~ habitat, data=data.species.a)
summary(fit)
TukeyHSD(fit)
aov() はタイプ I エラーを使用しているようです。. . 適切ではないと思います。さらに、テューキーの検定はペアワイズ比較のアプローチとして保守的すぎる可能性があると思います。ホームレンジに対する種と生息地の両方の影響を考慮した1つのANOVAを実行できるアプローチを誰かが助けてくれますか?