options ls=75 nodate; data one; *infile "recent2.txt" firstobs=2; input Season $ Plant $ Rep Location Results; day=rep; *sy=sqrt(results); *sy2=exp(0.2*log(results)); *sy3=exp(0.1*log(results)); y=results; if results=0 then y=15000+80*normal(123); *if day=1 and location=1 then y=y+100000*abs(normal(123)); *if day=1 then y=y+10000*abs(normal(123)); if day=1 then y=y+5000*abs(normal(12)); y=round(y,100); ly=log(y); sy=sqrt(y); sy2=exp(0.2*log(y)); if season="Winter"; site=location; /* Season Plant Rep Location Results */ cards; Winter Goldkist 1 1 42000 Winter Goldkist 1 1 37000 Winter Goldkist 1 1 45000 Winter Goldkist 1 1 95000 Winter Goldkist 1 1 0 Winter Goldkist 1 1 30000 Winter Goldkist 1 1 36000 Winter Goldkist 1 1 285000 Winter Goldkist 1 1 49000 Winter Goldkist 1 1 27000 Winter Goldkist 2 1 35000 Winter Goldkist 2 1 195000 Winter Goldkist 2 1 38000 Winter Goldkist 2 1 0 Winter Goldkist 2 1 32000 Winter Goldkist 2 1 0 Winter Goldkist 2 1 88000 Winter Goldkist 2 1 29000 Winter Goldkist 2 1 220000 Winter Goldkist 2 1 92000 Winter Goldkist 3 1 68000 Winter Goldkist 3 1 110000 Winter Goldkist 3 1 27000 Winter Goldkist 3 1 345000 Winter Goldkist 3 1 41500 Winter Goldkist 3 1 0 Winter Goldkist 3 1 27000 Winter Goldkist 3 1 195000 Winter Goldkist 3 1 33000 Winter Goldkist 3 1 91000 Winter Goldkist 1 2 45000 Winter Goldkist 1 2 68000 Winter Goldkist 1 2 14000 Winter Goldkist 1 2 18000 Winter Goldkist 1 2 110000 Winter Goldkist 1 2 12000 Winter Goldkist 1 2 0 Winter Goldkist 1 2 12000 Winter Goldkist 1 2 71000 Winter Goldkist 1 2 77000 Winter Goldkist 2 2 54000 Winter Goldkist 2 2 45000 Winter Goldkist 2 2 62000 Winter Goldkist 2 2 0 Winter Goldkist 2 2 82000 Winter Goldkist 2 2 39000 Winter Goldkist 2 2 54000 Winter Goldkist 2 2 67000 Winter Goldkist 2 2 49000 Winter Goldkist 2 2 53000 Winter Goldkist 3 2 11000 Winter Goldkist 3 2 35000 Winter Goldkist 3 2 6200 Winter Goldkist 3 2 18500 Winter Goldkist 3 2 46000 Winter Goldkist 3 2 77000 Winter Goldkist 3 2 55000 Winter Goldkist 3 2 26000 Winter Goldkist 3 2 0 Winter Goldkist 3 2 42000 Winter Goldkist 1 3 4000 Winter Goldkist 1 3 8000 Winter Goldkist 1 3 6000 Winter Goldkist 1 3 7000 Winter Goldkist 1 3 5000 Winter Goldkist 1 3 3000 Winter Goldkist 1 3 4000 Winter Goldkist 1 3 0 Winter Goldkist 1 3 6000 Winter Goldkist 2 3 0 Winter Goldkist 2 3 4000 Winter Goldkist 2 3 7000 Winter Goldkist 2 3 0 Winter Goldkist 2 3 6000 Winter Goldkist 2 3 4000 Winter Goldkist 2 3 4000 Winter Goldkist 2 3 0 Winter Goldkist 2 3 6000 Winter Goldkist 2 3 5000 Winter Goldkist 3 3 6000 Winter Goldkist 3 3 0 Winter Goldkist 3 3 2000 Winter Goldkist 3 3 1000 Winter Goldkist 3 3 4000 Winter Goldkist 3 3 4000 Winter Goldkist 3 3 6000 Winter Goldkist 3 3 0 Winter Goldkist 3 3 7000 Winter Goldkist 3 3 2000 Winter Goldkist 1 3 4000 Winter Goldkist 1 4 2000 Winter Goldkist 1 4 3000 Winter Goldkist 1 4 0 Winter Goldkist 1 4 0 Winter Goldkist 1 4 0 Winter Goldkist 1 4 5000 Winter Goldkist 1 4 0 Winter Goldkist 1 4 3000 Winter Goldkist 1 4 9000 Winter Goldkist 1 4 6000 Winter Goldkist 2 4 3000 Winter Goldkist 2 4 0 Winter Goldkist 2 4 7000 Winter Goldkist 2 4 5000 Winter Goldkist 2 4 3000 Winter Goldkist 2 4 0 Winter Goldkist 2 4 5000 Winter Goldkist 2 4 0 Winter Goldkist 2 4 0 Winter Goldkist 2 4 4000 Winter Goldkist 3 4 3000 Winter Goldkist 3 4 6000 Winter Goldkist 3 4 0 Winter Goldkist 3 4 0 Winter Goldkist 3 4 3000 Winter Goldkist 3 4 0 Winter Goldkist 3 4 0 Winter Goldkist 3 4 4000 Winter Goldkist 3 4 3000 Winter Goldkist 3 4 5000 Fall Perdue 1 1 97000 Fall Perdue 1 1 95000 Fall Perdue 1 1 36000 Fall Perdue 1 1 0 Fall Perdue 1 1 60000 Fall Perdue 1 1 40000 Fall Perdue 1 1 24000 Fall Perdue 1 1 110000 Fall Perdue 1 1 77000 Fall Perdue 1 1 33000 Fall Perdue 2 1 18000 Fall Perdue 2 1 60000 Fall Perdue 2 1 74000 Fall Perdue 2 1 55000 Fall Perdue 2 1 80000 Fall Perdue 2 1 0 Fall Perdue 2 1 0 Fall Perdue 2 1 13500 Fall Perdue 2 1 25300 Fall Perdue 2 1 37500 Fall Perdue 3 1 118000 Fall Perdue 3 1 90000 Fall Perdue 3 1 23000 Fall Perdue 3 1 128000 Fall Perdue 3 1 50000 Fall Perdue 3 1 0 Fall Perdue 3 1 96000 Fall Perdue 3 1 13800 Fall Perdue 3 1 40000 Fall Perdue 3 1 20000 Fall Perdue 1 2 9000 Fall Perdue 1 2 0 Fall Perdue 1 2 7000 Fall Perdue 1 2 7500 Fall Perdue 1 2 36000 Fall Perdue 1 2 41000 Fall Perdue 1 2 36000 Fall Perdue 1 2 20000 Fall Perdue 1 2 0 Fall Perdue 1 2 22000 Fall Perdue 2 2 30000 Fall Perdue 2 2 40000 Fall Perdue 2 2 0 Fall Perdue 2 2 7000 Fall Perdue 2 2 30000 Fall Perdue 2 2 6000 Fall Perdue 2 2 0 Fall Perdue 2 2 30000 Fall Perdue 2 2 15000 Fall Perdue 2 2 20000 Fall Perdue 3 2 0 Fall Perdue 3 2 20000 Fall Perdue 3 2 50000 Fall Perdue 3 2 20000 Fall Perdue 3 2 30000 Fall Perdue 3 2 20000 Fall Perdue 3 2 20000 Fall Perdue 3 2 9000 Fall Perdue 3 2 0 Fall Perdue 3 2 5000 Fall Perdue 1 3 8000 Fall Perdue 1 3 11000 Fall Perdue 1 3 0 Fall Perdue 1 3 14000 Fall Perdue 1 3 21000 Fall Perdue 1 3 0 Fall Perdue 1 3 18000 Fall Perdue 1 3 15000 Fall Perdue 1 3 12000 Fall Perdue 1 3 0 Fall Perdue 2 3 21000 Fall Perdue 2 3 0 Fall Perdue 2 3 24000 Fall Perdue 2 3 0 Fall Perdue 2 3 18000 Fall Perdue 2 3 16000 Fall Perdue 2 3 11000 Fall Perdue 2 3 22000 Fall Perdue 2 3 15000 Fall Perdue 2 3 19000 Fall Perdue 3 3 26000 Fall Perdue 3 3 31000 Fall Perdue 3 3 0 Fall Perdue 3 3 18000 Fall Perdue 3 3 22000 Fall Perdue 3 3 28000 Fall Perdue 3 3 0 Fall Perdue 3 3 21000 Fall Perdue 3 3 0 Fall Perdue 3 3 18000 Spring Perdue 1 1 77000 Spring Perdue 1 1 98000 Spring Perdue 1 1 44000 Spring Perdue 1 1 56000 Spring Perdue 1 1 0 Spring Perdue 1 1 68000 Spring Perdue 1 1 59000 Spring Perdue 1 1 0 Spring Perdue 1 1 66000 Spring Perdue 1 1 71000 Spring Perdue 2 1 75000 Spring Perdue 2 1 66000 Spring Perdue 2 1 0 Spring Perdue 2 1 58000 Spring Perdue 2 1 0 Spring Perdue 2 1 78000 Spring Perdue 2 1 55000 Spring Perdue 2 1 65000 Spring Perdue 2 1 82000 Spring Perdue 2 1 61000 Spring Perdue 3 1 55000 Spring Perdue 3 1 49000 Spring Perdue 3 1 77000 Spring Perdue 3 1 0 Spring Perdue 3 1 62000 Spring Perdue 3 1 59000 Spring Perdue 3 1 46000 Spring Perdue 3 1 52000 Spring Perdue 3 1 54000 Spring Perdue 3 1 82000 Spring Perdue 1 2 12000 Spring Perdue 1 2 11000 Spring Perdue 1 2 0 Spring Perdue 1 2 15000 Spring Perdue 1 2 9000 Spring Perdue 1 2 19000 Spring Perdue 1 2 17000 Spring Perdue 1 2 0 Spring Perdue 1 2 15000 Spring Perdue 1 2 14000 Spring Perdue 2 2 21000 Spring Perdue 2 2 0 Spring Perdue 2 2 26000 Spring Perdue 2 2 13000 Spring Perdue 2 2 9000 Spring Perdue 2 2 16000 Spring Perdue 2 2 31000 Spring Perdue 2 2 14000 Spring Perdue 2 2 0 Spring Perdue 2 2 33000 Spring Perdue 3 2 0 Spring Perdue 3 2 22000 Spring Perdue 3 2 18000 Spring Perdue 3 2 9000 Spring Perdue 3 2 0 Spring Perdue 3 2 14000 Spring Perdue 3 2 13000 Spring Perdue 3 2 11000 Spring Perdue 3 2 14000 Spring Perdue 3 2 11000 Spring Perdue 1 3 9000 Spring Perdue 1 3 0 Spring Perdue 1 3 14000 Spring Perdue 1 3 12000 Spring Perdue 1 3 8000 Spring Perdue 1 3 12000 Spring Perdue 1 3 11000 Spring Perdue 1 3 14000 Spring Perdue 1 3 0 Spring Perdue 1 3 0 Spring Perdue 2 3 9000 Spring Perdue 2 3 0 Spring Perdue 2 3 15000 Spring Perdue 2 3 11000 Spring Perdue 2 3 0 Spring Perdue 2 3 9000 Spring Perdue 2 3 0 Spring Perdue 2 3 8000 Spring Perdue 2 3 7000 Spring Perdue 2 3 11000 Spring Perdue 3 3 9000 Spring Perdue 3 3 6000 Spring Perdue 3 3 0 Spring Perdue 3 3 9000 Spring Perdue 3 3 10000 Spring Perdue 3 3 0 Spring Perdue 3 3 0 Spring Perdue 3 3 11000 Spring Perdue 3 3 8000 Spring Perdue 3 3 11000 ; run; /*proc print; var Season day site y; run; endsas;*/ proc sort; by season; run; proc means mean std n; by season; class location; var results; run; /*proc freq; by season; table ypos*location; run; endsas; */ proc glm; by season; class day location; *model results sy=location|day; model y sy ly =location|day; lsmeans location/pdiff adj=tukey; output out=two p=p ps pl r=r rs rl; run; /*proc gplot; by season; plot ly*day; plot ly*location; plot r*p; plot rs*ps; plot rl*pl; *plot r3*p3; run; endsas;*/ proc univariate normal plot; by season; var r rs rl; histogram r rs rl; run; proc mixed method=type3; class day location; model ly=location/ddfm=satterth; random day day*location; lsmeans location/adj=bon; run; proc mixed method=type3; class day location; parms / nobound; model ly=location/ddfm=satterth; random day day*location; lsmeans location/adj=bon; run; options ps=500; proc print data=two; var y ly day location p pl r rl; run; proc varcomp; by season; class day location; model sy=location|day/fixed=1; run;