Thank you both VERY MUCHfor your help.  Proc mixed seems to work except that I cannot get a main effect for tx_condition from the mianalyze statement -- I only get specific comparisons to the control/reference group. Is there a way to specify things to get that overall F (and p value)  for my categorical predictor tx_condition (with 3 levels). I think this is similar to the question Karen had previously. Amy ... View more

Hi Karen, Can you provide specific code for what you have done to circumvent this problem? I suspect I am encouring a similar issue.  I want to obtain an overall Fand p value for treatment (from the mianalyze step) for a model that has a 3-level categorical predictor (3 different types of treatment). I get a warning message when I run the glm step. proc glm data = XX; class tx_condition; model continuous_dv = continuous_iv tx_condition/inverse; by _imputation_; ods output parameterestimates = glmparms invxpx = glmpxi; run; ***WARNING MESSAGE: 'ParameterEstimates' was not created. **** Then when I turn to step 3 (mianalyze) I get an error. proc mianalyze; class tx_condition; model effects intercept continuous_iv tx_condition; parms = glmparms; xpxi = glmpxi edf = 218; run; *** ERROR MESSAGE: Parms Statement is not valid or it is used out of proper order*** I have read around and suspect it's an issue with the class statement but I'm not sure how to get around this problem. Thanks in advance for any thoughts/suggestions! Amy ... View more