Hello, I'm using SAS studio and I am having a peculiar problem with proc glm. I am using the below command to perform an ANCOVA looking at the difference between 2 SubjectTypes (cases-214 and controls-87) on the dependent variable (AVG_skel) controlling for age,sex,and site. proc glm data= sublist plots=meanplot(CL) ;
class SubjectType;
model AVG_Faskel = SubjectType Age_cal sexd sited;
means SubjectType;
lsmeans SubjectType/ stderr pdiff cov out=adjmeans;
run;
quit; For some reason when I run the command the procedure only uses 9/214 cases, but all of the controls. The log has no errors when running the procedure. I have tried doing proc freq to get tables of the SubjectTypes and it reports the correct number in each group so it is not a problem with the class variable. I have also run proc means for the AVG_skel variable and it reports the correct means and number of people in each group. What's even stranger is that proc glm does work with an earlier version of this file in which the AVG_skel variable has identical values. Any ideas on what could be the problem? I've also included the log output: 133 134 proc glm data= sublist plots=meanplot(CL) ; WARNING: The SAS/STAT product with which GLM is associated will be expiring soon, and is currently in warning mode to indicate this upcoming expiration. Please run PROC SETINIT to obtain more information on your warning period. 135 class group; 136 model AVG_Faskel = group Age_cal sexd sited; 137 means group; 138 lsmeans group/ stderr pdiff cov out=adjmeans; 139 run; NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. 140 quit; NOTE: The data set WORK.ADJMEANS has 2 observations and 7 variables. NOTE: PROCEDURE GLM used (Total process time): real time 9.10 seconds cpu time 0.92 seconds
... View more