# set up effects for two-factor model
ntukey.eff <- list(geno=c(-17, 14, 3), lat=c(-4, 1, 3))
names(ntukey.eff$geno) <- c("A","H","B")
names(ntukey.eff$lat) <- c("25","35","45")

# generate artificial data
ntukey <- expt(eff = ntukey.eff, grand=100,tukey = 0.2, rep = 1,anova=F)
ntukey$geno <- ordered(ntukey$geno,c("A","H","B"))
ntukey$reps <- NULL