Here is the code I am using. If I change the "adjust=tukey" to "adjust=t" the code crashes. Am I doing something wrong? I want to test the lsmeans for student's t comparisons. How do I do this? proc import datafile='/home/u48644977/Priyanka/305_All_Dry_Yield_Data_Adjusted.xlsx' out=work.dryYield dbms=xlsx replace; run; proc mixed plot=all covtest data=work.dryYield; class year harvest block treatment; model yield = block treatment|harvest|year /residual outp=resid s; random block*treatment; repeated / type=un subject=plot r; lsmeans treatment|harvest|year/ adjust=tukey; ods output diffs=ppp lsmeans=mmm; run; proc univariate noprint data=resid; histogram resid/ normal; inset n mean std="Std Dev" skewness; ppplot resid/normal (L=1 mu=est sigma=est); qqplot resid /normal (L=1 mu=est sigma=est) square; inset n mean std="Std Dev" skewness; run; %include '/home/u48644977/Priyanka/pdmix800.sas'; %pdmix800(ppp,mmm,alpha=.1,sort=yes); run;
... View more
