install.packages(“rmarkdown”)

# Set seed that determines the random number generator
set.seed(48183130)  #use your student ID instead of 12345678

# creates data for your homework
r<-c("A","A","A","A","A","A","B","B","B","B","B","B")
c<-c("M","M","M","F","F","F","M","M","M","F","F","F")
table<-data.frame(r,c)

table$X <-as.numeric(ifelse(table$r=="A" & table$c=="M",round(rnorm(2,mean=9,sd=3),0),
                 ifelse(table$r=="A" & table$c=="F",round(rnorm(2,mean=7,sd=3),0),
                        ifelse(table$r=="B" & table$c=="M",round(rnorm(2,mean=5,sd=3),0),
                               ifelse(table$r=="B" & table$c=="F", round(rnorm(2,mean=15,sd=3),0),"")))))
Row <- table$r

is.numeric(table$X)

## [1] TRUE

table

##    r c  X
## 1  A M 10
## 2  A M 12
## 3  A M 10
## 4  A F  2
## 5  A F  6
## 6  A F  2
## 7  B M  4
## 8  B M  5
## 9  B M  4
## 10 B F 18
## 11 B F 17
## 12 B F 18

# results for Two way ANOVA test
results <- aov(X ~ r*c, data = table)
summary(results)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## r            1   48.0    48.0   26.18 0.000911 ***
## c            1   27.0    27.0   14.73 0.004964 ** 
## r:c          1  320.3   320.3  174.73 1.02e-06 ***
## Residuals    8   14.7     1.8                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

# creates plot
interaction.plot(table$c, Row, table$X, main = "Two Way ANOVA", xlab = "Column", ylab = "X")