在SAS proc中输出预测值混合：禁止性能问题

Question

我已经注意到SAS proc混合的奇怪行为：具有适度大量行数的模型，只需要几秒钟就可以收敛，但是如果我要求输出预测值，则需要花费半小时才能完成运行。残差。看起来不正常的事情是当我使用nlme :: lme（）运行R中的类似模型时，我得到预测值＆amp;残差作为副作用，模型在几秒钟内完成。这让我觉得这不仅仅是我机器的内存限制。

这是一些示例代码。我无法提供我正在看到这个问题的真实数据，但结构是每个主题1-5行，~1500个独特主题，总共约5,000个结果 - 协变量组。

在SAS：

proc mixed data=testdata noclprint covtest;
    class subjid ed gender;
    model outcome = c_age ed gender / ddfm=kr solution residual outp=testpred;
    random int c_age / type=un sub=subjid;
run;

在R：

lme.test <- lme(outcome ~ c_age + ed + gender, data=testdata, 
               random = ~c_age|factor(subjid), na.action=na.omit)

相关统计数据：Win7，SAS 9.4（64位），R 3.3，nlme 3.1-131。