#####Read the data "Plant Data"##########

plant <- read.table("Plant experiment.txt",header=T, dec=,)


plant

attach(plant)


summary(weight)

boxplot(weight~group)


fit1 <- lm(weight~group)

summary(fit1)

residuals(fit1)
predict(fit1)

anova(fit1)


#################Two way analysis of Variance Example#######


y <- c(6.8,6.6,5.3,6.1,7.5,7.4,7.2,6.5,7.8,9.1,8.8,9.1)
a <- gl(3,4)
a
b <- gl(2,2, length(a))
b
anova(z <- lm(y~a*b))