用lcmm()指定生长混合物模型

时间:2014-12-17 02:51:52

标签: r

我试图使用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)

但是,我不确定我是否正确指定混合和随机参数。任何建议,将不胜感激!

数据在这里:http://www.statmodel.com/usersguide/chap8/ex8.1.dat

0 个答案:

没有答案