tomato <- dotna(read.table("tomato.dat",header=T))
tomato$entry <- factor(tomato$entry)
tomato$tg430 <- ordered(tomato$tg430,c("A","H","B"))
tomato$yr <- ordered(tomato$yr)

tomato1 <- tomato[tomato$yr==1&!is.na(tomato$mfwlog)&!is.na(tomato$tg430),]
tomato.aov <- aov(mfwlog~tg430,tomato1)

# reduce to complete data
tomatos <- tomato[!(is.na(tomato$tg430)|is.na(tomato$mfwlog)),]
dimnames(tomatos)[[1]]<- as.character(1:nrow(tomatos))