我试图使用lcmm()包从Rplus中的Mplus复制一些生长混合物建模结果。我在将模型规范从Mplus框架转换为R中的线性混合效果模型框架时遇到了一些麻烦。
在Mplus中,一个2级增长混合模型如下所示:
VARIABLE: NAMES ARE y1-y4;
CLASSES = c(2);
ANALYSIS: TYPE = MIXTURE;
MODEL:
%OVERALL%
i s | y1@0 y2@1 y3@2 y4@3;
OUTPUT:
TECH1;
在R中,使用lcmm()包,我认为模型规范看起来应该是这样的:
m1 <- hlme(y ~ time,
mixture = ~ time,
random = ~ time,
subject = 'id',
ng = 2,
idiag = F,
data = dat.long)
但是,我不确定我是否正确指定混合和随机参数。任何建议,将不胜感激!