Lab 4

ques4.dat=c(7,8,9,10,11,5,6,7,8,9,3,4,5,6,7,11,12,13,14,15)
A=factor(c(rep('a1',10),rep('a2',10)))
B=factor(rep(c(rep('b1',5),rep('b2',5)),2))
sub=1:20
ques4.df=data.frame(sub,A,B,ques4.dat)

ques4.aov=aov(ques4.dat~A*B,data=ques4.df)

ques4.aov

## Call:
##    aov(formula = ques4.dat ~ A * B, data = ques4.df)
## 
## Terms:
##                   A   B A:B Residuals
## Sum of Squares    5  45 125        40
## Deg. of Freedom   1   1   1        16
## 
## Residual standard error: 1.581139
## Estimated effects may be unbalanced

ques4.sum=summary(ques4.aov)

ques4sum=unlist(ques4.sum)

df1=as.numeric(ques4sum['Df1'])
df2=as.numeric(ques4sum['Df2'])
df3=as.numeric(ques4sum['Df3'])
df4=as.numeric(ques4sum['Df4'])

SS1=as.numeric(ques4sum['Sum Sq1'])
SS2=as.numeric(ques4sum['Sum Sq2'])
SS3=as.numeric(ques4sum['Sum Sq3'])
SS4=as.numeric(ques4sum['Sum Sq4'])

ms1=as.numeric(ques4sum['Mean Sq1'])
ms2=as.numeric(ques4sum['Mean Sq2'])
ms3=as.numeric(ques4sum['Mean Sq3'])
ms4=as.numeric(ques4sum['Mean Sq4'])

F1=as.numeric(ques4sum['F value1'])
F2=as.numeric(ques4sum['F value2'])
F3=as.numeric(ques4sum['F value3'])

p1=as.numeric(ques4sum['Pr(>F)1'])
p2=as.numeric(ques4sum['Pr(>F)2'])
p3=as.numeric(ques4sum['Pr(>F)3'])

There was no significant effect of A, F(1, 16)=2, MSe=5, p=0.1764632.

There was a significant wffect of B, F(1, 16)=18, MSe=45, p=6.20655410^{-4}.

There was a sigificant effect of the interaction of A on B, F(1, 16)=50, MSe=125, p=2.645792310^{-6}.