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