## AHW3
## March 3, 2008

###
rm(list = ls())
##

## population paramaters
muA <- 2
muB <- 6
muC <- 7
muD<- 10

sigma <- 2

## generate data
yA <- rnorm(4, muA, sigma)
yB <- rnorm(6, muB, sigma)
yC <- rnorm(6, muC, sigma)
yD <- rnorm(8, muD, sigma)

y <- c(yA, yB, yC, yD)
x <- c(rep("A", 4), rep("B", 6), rep("C", 6), rep("D", 8))

## ANOVA table
anova(lm(y~x))