time = c(41,43,49,51,52,44,50,46,50,46,48,44)
print(time)
rna = c('m1','m1','m2','m2')
method = factor(rep(rna,3))
print(method)
cna = c(rep('A',4),rep('B',4),rep('C',4))
clas = factor(cna)
print(clas)

data = data.frame(time,method,clas)

print(data)

fit = aov(time~method+clas+method:clas,data=data)
summary(fit)