coag.times _ scan()
 62 63 68 56
 60 67 66 62
 63 71 71 60
 59 64 67 61
 65 68 63
 66 68 64
 63
 59
 
 diet _ factor(c(rep(LETTERS[1:4],4),rep(LETTERS[1:3],2),c("A","A")))
 split(coag.times,diet)  # check that the factor labels are right