options ps=54 ls=76; data wheat; infile '../../data/wintwheat.dat'; input id variety yield moist; run; proc mixed; class variety; model yield = moist / solution; random int moist / type=un sub=variety solution; run; quit;