我正在将LCMM包与功能hlme一起使用。
我想知道每个班级都分类了哪些ID。我的一些案例已从模型中删除。因此，当我使用$pprob时，id变量基本上只是从1开始并不断增加。不是原始研究ID。因此我无法将其与原始数据集合并。
有人知道该怎么做吗？
这是代码：
q4r.fv <-hlme（i_rest〜days + days2，主题='ninclu'，ng = 4，idiag = T，混合物=〜days + days2， data = patch，maxiter = 100，returndata = TRUE，classmb =〜tttantal + sexe + atchroniq + hads_anxiI + hads_depI + pcs_totI） 摘要（q4r.fv）
fmq4r.fv <-q4r.fv $ pprob
write.csv（fmq4r.fv，file =“ fmq4r.fv”，row.names = F）
使用csv文件，可以获得以下内容。 ninclu应该是我的ID，但不再与我原来的minclu变量匹配，后者是一个字符串参与者ID变量
print.data.frame（fmq4r.fv） ninclu类prob1 prob2 prob3 prob4 1 1 2 7.416779e-09 9.635142e-01 3.630078e-02 1.850563e-04 2 2 2 5.479232e-02 8.710804e-01 7.412726e-02 9.118352e-16 3 3 1 9.933911e-01 6.607882e-03 9.830514e-07 1.110920e-23 4 4 2 2.620132e-07 9.991825e-01 8.155809e-04 1.631318e-06 5 5 2 4.382259e-04 9.877001e-01 1.186168e-02 1.166050e-11 6 6 3 4.239271e-09 2.361263e-01 7.634902e-01 3.835313e-04
答案 0 :(得分：1)
您提到您的主题ID是字符串参与者ID变量。 如果查看lcmm的源代码，您会发现lcmm需要数字主题：
if(!is.numeric(data[,subject])) stop("The argument subject must be numeric")
我怀疑字符串参与者ID被强制为因数，并且它是默认因数编码，在模型输出预测中显示为“ subj”。
尝试使用数字参与者ID，看看结果是否更有意义。