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