02-14-2013 06:57 AM
I need do a plot for my model in proc Mixed.
It is my final model:
proc mixed data=DatosExpIny2010 covtest;
class nest sex;
model Time= treatment date/ddfm=satterth;
*estimate 'volumen huevo' huevovolok 1;
estimate 'orden_p' opuesta 1;
estimate 'trt' trt 1;
*estimate 'Fechatot' Fechatot 1;
*estimate 'Fase Pollo' F_Pollo 1;
And I need a plot with 95% error.
Could you tell me the correct script?
Thank you very much.
All the best
02-14-2013 11:04 AM
Yes, I need to make a graph that represents the model. It should be a plot for Treatment and Time, becaue they are continuous variables. Also I need it with confidence limits. Thanks!
02-14-2013 11:29 AM
As far as I can see, there is no built in plotting with PROC MIXED that makes a plot of the model.
You might want to do this through PROC GLIMMIX followed by PROC PLM with the EFFECTPLOT statement.
02-14-2013 04:39 PM
It sounds great! but... how i can do it? Thanks!
Here's an example using PROC ORTHOREG, maybe this can get you started
02-15-2013 07:49 AM
As I said, it was an example.
I think you will need to use PROC MIXED or PROC GLIMMIX followed by PROC PLM, if you want to include a RANDOM statement in the model.
02-14-2013 01:17 PM
If you want 95% CL for the mean predicted values, I think you can get that in an output data set. On the MODEL statement, use the OUTP= option to get an output data set. The doc says
"The L95 and U95 variables provide a t-type confidence interval for the predicted values, and they correspond to the L95M and U95M variables from the GLM and REG procedures for fixed-effects models. The residuals are the observed minus the predicted values."
You can then use PROC SGPLOT to plot the predicted and 95% CLM.
02-15-2013 08:46 AM
If you want to see examples of usnig PROC PLM, see the paper by Tobias and Cai (2010) : http://support.sas.com/resources/papers/proceedings10/258-2010.pdf