unbal <- read.table("unbal.dat",header=T)
unbal$A <- factor(unbal$A)
unbal$B <- factor(unbal$B)

# full model for unbalanced data
unbal.aov <- aov(y~A*B,unbal)
# Type I SS
print(summary(unbal.aov))
# Type III SS
print(drop1(unbal.aov,formula(unbal.aov)))

# least squares means
library(pda)
print(lsmean(unbal.aov))