pdf("ex0605.pdf")

ex0605dat <- read.fwf("http://www.stat.ufl.edu/~winner/data/biostat/ex0605.dat", width=c(8,8,8),
col.names=c("ethnic", "gender", "clrnc"))

# Create qualitative factor variable for ethnic, and assign names to levels

fethnic <- factor(ex0605dat$ethnic, levels=1:2)
levels(fethnic) <- c("hispanic", "anglo")

# Create qualitative factor variable for gender, and assign names to levels

fgender <- factor(ex0605dat$gender, levels=1:2)
levels(fgender) <- c("female", "male")


ex0605 <- data.frame(clrnc=ex0605dat$clrnc, fethnic, fgender)

attach(ex0605)


# Run 2 way ANOVA with main effects for fethnic and fgender, and interaction
#   by typing fethnic*fgender as "independent variable", r includes main effects and interactions

ex0605.aov <- aov(clrnc ~ fethnic*fgender)

# Print out the ANOVA table and F-tests

anova(ex0605.aov)


dev.off()