I would like to give you some feedback on my question because you took some time to answer it.   In the meantime I was able to find out that MPlus doesn't give empirical coverage rates as a result, but those which are based on a 5% error level for the assumption of a normal distribution. If I look at the empirical coverage rates in MPlus, I get quite exactly the same rates as in SAS.   Nevertheless, it has cost a lot of detective work to look behind the program syntax, which I personally find very unfortunate.   Many thanks again for your detailed help!   Best regards Julia Mang ... View more

Dear Reinhard, dear lotcarrots,    thank you both for your replies!   Both comments have been very helpful in finding the right ods output procedure for my intended analyses.   For documentation reasons, here my final code for this post:   /*ods trace on;*/ proc glimmix data=dataForSAS_200_25_01_Modell1 method=quadrature empirical=classical; class cluster; model ach = / solution dist=normal ; random int/subject=cluster solution ; /*WHERE rep < "10_200_25_01.dat";*/ by rep; ods output ParameterEstimates=param1 CovParms=param2; run; /*ods trace off;*/ ods rtf file="C:\\storage_path\filename.rtf"; proc means data=param1; proc means data=param2; class CovParm; run; ods rtf close;   Thanks again for your help! It was really appreciated by me!   Kind regards,   Julia ... View more