15) RStudio - One-way ANOVA (Tesadüf parselleri deneme deseninde tek yönlü varyans analizi)

One-way ANOVA Tesadüf parselleri deneme deseninde tek yönlü varyans analizi ANOVA One-way ANOVA One-way classification Tek yönlü varyans analizi Tesadüf parselleri deneme deseni Varyans analizi Lineer model Duncan Tukey Homojenlik testi Levene Bartlett R Kodu ------------ One-way ANOVA Paketleri yükle ---- duncan için ##install.packages("agricolae") library(agricolae) leveneTest için install.packages("car") library(car) plotmeans için install.packages("gplots") library(gplots) Veriler ---- Verileri oku myData = read.table("D:/Belgeler/orhantez.txt", header=T) Veri yapısı str(myData) Analiz ---- # model = aov(AGIRLIK~BREED,data=myData) summary(model) (out = duncan.test(model,"BREED", main="Irk faktörünün canlı ağırlığa etkisi",)) plot(out,variation="IQR", main="Fig.-1: Irk faktörünün canlı ağırlığa etkisi") duncan.test(model,"BREED",alpha=0.01,console=TRUE) TukeyHSD(model, conf.level = 0.95) plot(TukeyHSD(model, conf.level = 0.95), las=1, col = "blue") leveneTest(AGIRLIK ~ BREED, data = myData) bartlett.test(AGIRLIK ~ BREED, data = myData) duncan.test() eski versiyonu df = df.residual(model) MSerror = deviance(model)/df out = with(myData,duncan.test(AGIRLIK,BREED, df,MSerror, group=TRUE)) plot(out,horiz=TRUE,las=1) print(out$groups) lModel = lm(AGIRLIK~BREED, data=myData) anova(lModel) duzey = nlevels(myData$BREED) boxplot(myData$AGIRLIK~myData$BREED, main="Fig.-2: Irk faktörünün canlı ağırlığa etkisi", col= rainbow(duzey)) plotmeans(myData$AGIRLIK~myData$BREED, main="Fig.-3: Irk faktörünün canlı ağırlığa etkisi", ylab = "Canlı ağırlık (kg)", xlab = "Koyun Irkları")