## R-script to the solution of exercise 1.

## 1. ##

blood <-c(62,60,63,59,63,67,71,64,65,66,68,66,71,67,68,68,56,62,60,61,63,64,63,59) 
tr <- c(1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,4,4,4,4,4,4,4,4)
b.data <- data.frame(cbind(blood,tr))
b.data$tr <- as.factor(b.data$tr)

## a)

plot(b.data$tr,b.data$blood)


## d)

fit.blood <- aov(b.data$blood ~ b.data$tr)
summary(fit.blood)


## 3. ##

N2 <- c(19.4,32.6,27,32.1,33,18.2,24.6,25.5,19.4,21.7,20.8,20.7,
        21,20.5,18.8,18.6,20.1,21.3)   
strain <- c(1,1,1,1,1,5,5,5,5,5,5,7,7,7,7,7,7,7)
r.data <- data.frame(cbind(N2,strain))
r.data$strain <- as.factor(b.data$strain)

## a)

plot(r.data$strain,r.data$N2)

## b)

fit.n2 <- aov(r.data$N2 ~ r.data$strain)
summary(fit.n2)

## c)

par(mfrow=c(2,2))
plot(fit.n2)

rr.data <- r.data[-1,]
fit.n2mod <- aov(rr.data$N2~rr.data$strain)
summary(fit.n2mod)

par(mfrow=c(2,2))
plot(fit.n2mod)