Dear SAS community,   I am running a simulation study with different hierarchical models and different weighting scenarios with proc glimmix. I use e.g. for the nullmodell the following code:   proc glimmix data=data method=quadrature empirical=classical; class cluster; model ach = / dist=normal solution; random int / subject=cluster weight=schwt solution; by rep; ods output ParameterEstimates=param1 CovParms=param2; run;   Data have been simulated in R having for this scenario variable ach ~N(500,100) and ICC=0.1 in each of 1000 data sets.   When I look at the estimated outcomes for each single "by" replication, I recover a very bad coverage rate, i.e. the estimated parameters are badly included in the confidence intervall with given population values from the simulation (e.g. for the intercept of this nullmodell the coverage rate is 71%; this appears by the way for all different weighting scenarios).    I use the same simulated data and procedures also with the programm MPLUS and very high coverage rates (over 95%) appear. From my research, MPLUS is using the same estimation technique as I stated in the SAS syntax with method=quadrature and empirical=classical.   I assume that I am specifing still something wrong, otherwise I could not explain those bad coverage rates. Can you imaging what can couse this problem? Or should I use PROC MIXED instead of PROC GLIMMIX, but I though the latter one is preferred nowadays?   Thank you in advance for your help and kind regards,   Julia Mang ... View more

Dear All,   I am using the following SAS Syntax:   proc glimmix data=dataForSAS_200_25_01_Modell1 method=quadrature empirical=classical; class cluster; model ach = / solution dist=normal ; random int/subject=cluster solution ; by rep; ods output ParameterEstimatesStart=parms; run;   I am analyzing different datasets defined through the "by" statement. The different analyses performed well (and output for each dataset is provided), but when I want to summarize those results (ods output ParameterEstimatesStart=parms;), I get the following error message:   WARNING: Output 'ParameterEstimatesStart' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used.   For my analyses, I am reffering to Appendix E from the paper   Lee, S. (2015). Implementing a Simulation Study Using Multiple Software Packages for Structural Equation Modeling. SAGE Open, 5(3), 215824401559182. https://doi.org/10.1177/2158244015591823.   I assume I am missing a special code for the proc glimmix insted of proc calis (described in the paper). Does anybody know what I am doing wrong? I am using SAS 9.4.   Thank you in advance for your support. I am looking forward to an answer.   Kind regards, Julia ... View more