A<-rep(seq(1:2),each=6)
B<-rep(seq(1:3),4)
A<-as.factor(A)
B<-as.factor(B)
obs<-c(4.6,10.15,11.01,4.4,10.2,10.58,3.2,9.38,10.81,3.5,10.02,10.60)
model<-aov(obs~A*B)
summary(model)
## Df Sum Sq Mean Sq F value Pr(>F)
## A 1 0.98 0.98 15.26 0.00793 **
## B 2 111.19 55.59 865.16 4.13e-08 ***
## A:B 2 0.58 0.29 4.48 0.06450 .
## Residuals 6 0.39 0.06
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
data.frame(A,B,obs)
## A B obs
## 1 1 1 4.60
## 2 1 2 10.15
## 3 1 3 11.01
## 4 1 1 4.40
## 5 1 2 10.20
## 6 1 3 10.58
## 7 2 1 3.20
## 8 2 2 9.38
## 9 2 3 10.81
## 10 2 1 3.50
## 11 2 2 10.02
## 12 2 3 10.60