options ps=54 ls=76; goptions reset=all device=pdf gsfname=output gsfmode=replace; data amoeba; retain time_order 0; infile 'entozamoeba.dat'; input trt yield; time_order=time_order+1; run; proc glm; class trt; model yield = trt; means trt / bon tukey hovtest; output out=amoebaout r=amoebares; run; filename output 'amoeba1.pdf'; proc gchart; vbar amoebares / midpoints=-67.5 to 67.5 by 15; run; quit; filename output 'amoeba2.pdf'; symbol v=star i=join; proc gplot; plot amoebares*time_order; run; proc npar1way wilcoxon; class trt; var yield; run; quit;