Hello all,   I am conducting a factor analysis on a matrix of tetrachoric correlations (since my vars are binary). I create the matrix using PROC CORR, and then use that as input dataset for PROC FACTOR. I would like, however, compute factor scores for the two factors that I extract, and have them added to the original dataset (DATA). Proc factor does not do that when the input is not raw data but a matrix of correlations. I can of course extract coefficients, and then do it "manually".  Is there a way to do this automatically? Syntax below. Thank you all in advance. Eman   proc corr data = DATA polychoric out=TETCORR; var x1-x10; proc factor data=TETCORR (type=corr) nfactors=2 method=prin rotate=oblimin ; var x1-x10; ... View more

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