Hi, I'm using the following code. ods graphics on; ods output lsmeans=estimate diffs=diff; proc mixed data = mydata (where=(treatn in (1,2,3))) plots=all; by estremn estrem atypen atype paramcd param avisitn avisit ; class treatn; model aval = treatn base/residual outp=resid; lsmeans treatn / pdiff=control('3'); run; ods output close; ods graphics off; The aim is to get the LS means, SE and the difference of LS means and SEs from plcebo(treatn=3) to each of the treatments (treatn=1,2) as 'placebo - treatment. I believe currently the code is generating the diffrence other way round as (treatment - placebo). Need to know how to reverse the direction of difference to get the difference. Thanks ... View more

Hi, I'm using proc mixed as follows: proc mixed data = mydata (where=(treatn in (1,2,3))) plots=all; by estremn estrem atypen atype paramcd param avisitn avisit ; class treatn; model aval = treatn base/residual outp=resid; lsmeans treatn / pdiff=control('3'); run; The difference of lsmeans and SE for 2 treatments (treatn=1,2) is required to be estimated from placebo (treatn=3) as 'placebo - treatment'. I believe currently pdiff is giving the difference as 'treatment - placebo'. Help is required to reverse the direction of the difference so that I can get the difference of each treatment from placebo. Thanks ... View more

Hi, I'm using the following code.   ods graphics on; ods output lsmeans=estimate_1 diffs=diff_1; proc mixed data = mydata (where=(treatn in (1,2,3))) plots=all; by atype paramcd avisitn avisit; class treatn; model aval = treatn avisitn base; lsmeans treatn / pdiff=control('3'); run; ods output close; ods graphics off; Need to know how to get residual SD/variance and plot studentized residuals using proc mixed ... View more

Hi, I'm using proc mix to fit ANCOVA model to get the following: 1) LS means of each treatment 2) SE of each treatment 3) LS means difference (placebo - each treatment) with placebo as treatn=3 4) SE difference (placebo - each treatment) 5) SD (standard deviation of residuals from the model) 6) Plot studentized residuals versus predicted values, and also check the distribution of (unstandardized) residuals using a plot.   The code is: ods output lsmeans=estimate diffs=diff; proc mixed data = mydata (where=(treatn in (1,2,3))); by treatn treat avisitn avisit; class treatn; model aval = treatn avisitn base; lsmeans treatn / diff=control('3'); run; ods output close;   I'm getting error: "WARNING: Output 'Diff' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used."   Help is needed to resolve the error and guidance to use the proc mix appropriately to achieve above mentioned requirements for ANCOVA. Thanks ... View more