#3
chickwts
View(chickwts)
str(chickwts)
#badamy normalność
CA=subset(chickwts, feed=="casein")
head(CA)
shapiro.test(CA$weight)
#pvalue>>0.05, dane z rozkladu normalnego
qqnorm(CA$weight)
qqline(CA$weight)
bartlett.test(weight~feed, data=chickwts)
boxplot(weight~feed, data=chickwts,las=1, col=5)
plot.design(chickwts)
library(doBy)
summaryBy(weight~feed,data=chickwts,FUN=c(mean,var,sum))
#analiza wariancji
H0: uc=uh=ul=um=us=usu
H1~H0
analiza=aov(weight~feed,data=chickwts)
summary(analiza)
pairwise.t.test(chickwts$weight,chickwts$feed)
TukeyHSD(analiza)