After running proc nlmixed, I can make predictions (PREDICT_1) using the entire estimated equation (along with the random error term u) and conduct some post-estimation analyses as in the example below. Is there a way of doing this in PROC GLIMMIX and how would I do it. Thank you for any help you can give me.
Yes, you can obtain predicted values from the GLIMMIX procedure. However, it is probably not quite as easy to write the code in GLIMMIX because you need to write an ESTIMATE statement for every subject. Thus, you would need code like the following:
Note that each estimate statement specifies the same linear combination of the fixed effects but picks up a different column of the random effect q03. We need a separate ESTIMATE statement for each level of q03. In addition, in order to predict the expected probability value (rather than eta), we need the ILINK option on each ESTIMATE statement.
Note, too, the use of the vertical bar (|) before the random effect q03 in ESTIMATE statement. The vertical bar indicates that the terms which follow are found in the random effects structure.
Additionally, it should be observed that you will not get the same results out of the ESTIMATE statement if you specify the RANDOM statement as
random intercept / subject=r03;
The GLIMMIX procedure makes model specification much easier than does the NLMIXED procedure. However, obtaining predicted values conditional on levels of the random effects is considerably easier in NLMIXED.