Thank you, PGStats! The only warning message in the log is WARNING: ODS graphics with more than 5000 points have been suppressed. Use the PLOTS(MAXPOINTS= ) option in the PROC MIXED statement to change or override the cutoff.   I cannot match it with "numerical instability". Is there any way to check for that? Or, you may have other suggestions...   Thank you!  ... View more

Good point, PGStats. I did not think of that. This is the code to create the long data. I always use PROC TRANSPOSE but for this project, I received the code as follows and I am not sure how that works: DATA long; SET wide; measure = 1; *this is the first observation from the repeated measures; Var_1 = Var1; *This is variable that is different in the repeated measures; Var_Y_1 = VarY11; *This is the response variable at first measurement; VAr_Y_2 = VarY12; *I have another response var that is used in a separate regression analysis; output; measure = 2; *this is the second observation from the repeated measures; Var_1 = Var2; *This is variable that is different in the repeated measures; Var_Y_1 = VarY21; VAr_Y_2 = VarY22; output; *the following makes the 3rd record pertaining to the product seen third; measure = 3; Var_1 = Var3; Var_Y_1 = VarY31; VAr_Y_2 = VarY32; output; run; Thank you! ... View more

Rick_SAS, my log does not have rows for 2-60! Is there something wrong with it? I did not check row numbers before... ... View more

Thanks Paige Miller! Do I copy with Ctrl+C? Where is the </> icon? ... View more

1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 61 62 proc mixed data=long; 63 where this_var=0; 64 class Var_1 (ref='a') Var_2 (ref='a') Var_3 (ref='a') 64 ! Var_4 (ref='a') age_group ethnicity education 64 ! income Var_5 Var_6  64 !Var_7 /ref=first; 65 model Var_Y=age_group ethnicity education income  65 !  Var_1 Var_2 Var_3 Var_4 Var_5 Var_6 Var_7/ solution cl; 66 repeated VAr_1/subject=ID type=VC; 67 lsmeans age_group ethnicity education income   67 ! Var_1 Var_2 Var_3 Var_4 Var_5 Var_6 Var_7/diff cl; 68 run;   WARNING: ODS graphics with more than 5000 points have been suppressed. Use the PLOTS(MAXPOINTS= ) option in the PROC MIXED statement to change or override the cutoff. NOTE: Convergence criteria met. NOTE: PROCEDURE MIXED used (Total process time): real time 30.76 seconds cpu time 30.70 seconds     70 71 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 84   ... View more

Thank you Rick_SAS for your idea! I will post my SAS log... ... View more

Thank you, SteveDenham! My concern is why even with the same covariance structure my results are different. I get different results today than my results on Friday. I double-checked everything and I did not change anything. I restarted my computer, I got different results again! ... View more

Thank you, Paige Miller and Rick_SAS for your thoughts on my question.  Rick_SAS, I did not include a random intercept in my model yet (if your question is about that). I am planning to do so though, once I have a better sense of what is going on in my model.  Do you think it has anything with local minimized errors? In some optimization methods there are local and overall minimum errors, I am not sure if this applies to my model.   ... View more

Rick_SAS, may you please provide more details on how I do that? I mean is it a search on the PROC MIXED overall, or in my model. Also, how the LONG data set may be different each time I run program? Thank you! ... View more

Yes! For example when I use a different TYPE for covariance structure, sometimes my model converges, and sometimes it does not! I asked another friend to try the same data and same model, and our results are again different! ... View more

Yes! ... View more