hello all, I would appreciate some help with the syntax for ESTIMATE in the following model:   proc glimmix data=srst2 ; class id IV; effect spl= spline(progress / degree=2 basis=bspline details knotmethod=percentiles(2)); model DV = IV spl IV*spl / noint s; random int / subject=id; store dv; proc plm source=dv; effectplot slicefit(x=progress sliceby=IV); run; This model results in a significant main effect of spl, but no significant main effect of IV or interaction effect. See the plot below.   I would like to, nevertheless to test (using ESTIMATE):   1. The overall linear effect of spl (which is akin to the main effect that I test in the model) 2. The main effect of spl, separately for the two levels of IV (-1 and 1) 3. The effect of spl between -1.5 and 0, separately for the two levels of IV (-1 and 1) 4. The effect of spl between  0 and 1.5, separately for the two levels of IV (-1 and 1)   (*note that my model has both main effects and the interaction)     ... View more

Hello,   With this syntax: proc sgplot data=plot ; vbox value / category=DV group=genre ; run; SGPLOT produces a VBOX with different colors for the two groups defined in GROUP=   See graphs below.    Instead, I would like the two levels of groups (liteary and popular) to have different patterns (for instance, small dots for one and small x for the other).  I have been looking into various options with ATTRMAP, but I do not seem to be able to get it.   Suggestions?     ... View more

Hello all,   I am comparing plots produced by PROC MIXED (+ PROC PLM), GENMOD (using the EFFECTPLOT) and GLIMMIX with spline (+  PROC PLM).  The  first two (MIXED and GENMOD) give identical results while GLIMMIX with spline gives me a different result. This clearly depends on the spline option. As can be seen in the syntax that I include, the plot for GLIMMIX  is with  effect spl=spline(segment);  I tried  effect spl=spline(segment/naturalcubic); Of course the plot changes, but is still very different than what I obtain with MIXED or GENMOD.    The variable segment is standardized, and it has 10 levels (-1.57    -1.22    -0.87    -0.52    -0.17    0.17     0.52     0.87     1.22     1.57) -- I treat it as a continuous variable, of course.   Any suggestion on how to define the spline, to get the same plots as MIXED and GENMOD?   mixed + plm genmod (effectplot) glimmix / spline + plm     proc mixed data=e ; class id iv; model posemo= iv|segment|segment|segment; random id; store mod; run; proc plm source=mod; effectplot slicefit(x=segment sliceby=iv); run; proc genmod data=e ; class id iv; model posemo= iv|segment|segment|segment; repeated subject=id; effectplot slicefit(x=segment sliceby=iv); run; proc glimmix data=e; class id iv; effect spl=spline(segment); model posemo = iv*spl / noint s; random int / subject=id; store mod3; output out=gmxout pred=p; run; proc plm source=mod3; effectplot slicefit(x=segment sliceby=iv); run;     ... View more

hello.  I am using PROC TRANSPOSE to transpose three variables in a dataset, using ID as the BY variable. Data and syntax is included. However, I would like to keep two other variables: age and gender, so that my final file looks like the table below.  I know how to do it with a different method, but at times it would be much more efficient if I could use PROC TRANSPOSE. Any suggestion?   id _NAME_ COL1 gender age 1 v1 3 M 20 1 v2 4 M 20 1 v3 3 M 20 2 v1 4 F 22 2 v2 3 F 22 2 v3 2 F 22 3 v1 5 M 34 3 v2 1 M 34 3 v3 5 M 34     data have; input id age gender$ v1 v2 v3; cards; 1 20 M 3 4 3 2 22 F 4 3 2 3 34 M 5 1 5 ; proc print; proc transpose data=have out=want; var v1 v2 v3; by id; proc print; ... View more