Carrying Out the ANOVA F-test
setwd("/Users/Greg/Documents/Grad School/LAIS 609C/Assignments")
load("flicker.RData")
boxplot(flicker$cff~flicker$color,xlab="Eye Color",ylab="CFF")

tapply(flicker$cff, flicker$color, mean)
## blue brown green
## 28.16667 25.58750 26.92000
tapply(flicker$cff, flicker$color, sd)
## blue brown green
## 1.527962 1.365323 1.843095
cff.aov=aov(cff~color,flicker)
anova(cff.aov)
## Analysis of Variance Table
##
## Response: cff
## Df Sum Sq Mean Sq F value Pr(>F)
## color 2 22.997 11.4986 4.8023 0.02325 *
## Residuals 16 38.310 2.3944
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1