f1 <- c(17.6,18.9,16.3,17.4,20.1,21.6)
f2 <- c(16.9,15.3,18.6,17.1,19.5,20.3)
f3 <- c(21.4,23.6,19.4,18.5,20.5,22.3)
f4 <- c(19.3,21.1,16.9,17.5,18.3,19.8)
life <- c(f1,f2,f3,f4)
types <- c(rep(1,6),rep(2,6),rep(3,6),rep(4,6))
types <- as.factor(types)
fluid <- cbind(life,types)
fluid <- as.data.frame(fluid)
model <- aov(life~types)
summary(model)
## Df Sum Sq Mean Sq F value Pr(>F)
## types 3 30.17 10.05 3.047 0.0525 .
## Residuals 20 65.99 3.30
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
a <- mean(f1)
b <- mean(f2)
c <- mean(f3)
d <- mean(f4)
plot(model)
