SNK test

3.11(a)(b)

```{r}
x<-c(24,28,37,30,31)
y<-c(37,44,31,35,36)
z<-c(42,47,52,38,44)
data1<-data.frame(x,y,z)#x is 20g;y is 30g;z is 40g
data1<-stack(data1)
colnames(data1)<-c("values","Does")
data1
anova_data1<-aov(values~Does,data =data1)
summary(anova_data1)
library(agricolae)
out<-SNK.test(anova_data1,3,12,24.2, alpha = 0.05, group=TRUE, main = NULL,console=FALSE)
print(SNK.test(anova_data1,"Does", group=FALSE))
df<-df.residual(anova_data1)
MSerror<-deviance(anova_data1)/df
out <- with(data1,SNK.test(values,Does,df,MSerror, group=TRUE))
print(out$groups)
```

3.24(a)(b)

```{r}
w<-c(17.6,18.9,16.3,17.4,20.1,21.6)
x<-c(16.9,15.3,18.6,17.1,19.5,20.3)
y<-c(21.4,23.6,19.4,18.5,20.5,22.3)
z<-c(19.3,21.1,16.9,17.5,18.3,19.8)
data0<-data.frame(w,x,y,z)#w is type1;x is type2;y is type3;z is type4
data0<-stack(data0)
colnames(data0)<-c("values","TYPE")
data0
anova_data0<-aov(values~TYPE,data =data0)
summary(anova_data0)
library(agricolae)
out<-SNK.test(anova_data0,3,20,3.3, alpha = 0.05, group=TRUE, main = NULL,console=FALSE)
print(SNK.test(anova_data0,"TYPE", group=FALSE))
df<-df.residual(anova_data0)
MSerror<-deviance(anova_data0)/df
out <- with(data0,SNK.test(values,TYPE,df,MSerror, group=TRUE))
print(out$groups)
```