Я использую пакет LCMM с функцией hlme.
Я хочу знать, какие идентификаторы классифицированы в каждом классе.Некоторые из моих дел исключаются из модели.Поэтому, когда я использую $pprob, переменная id в основном просто начинается с 1 и увеличивается.Это не оригинальный идентификатор исследования.Поэтому я не могу объединить его с исходным набором данных.
Кто-нибудь знает, как это сделать?
Вот код:
q4r.fv <- hlme (i_rest ~ days+ days2, субъект = 'ninclu', ng = 4, idiag = T, смесь = ~ days + days2, данные = патч, maxiter = 100, returndata = TRUE, classmb = ~ tttantal + sexe + atchroniq + hads_anxiI + hads_depI + pcs_totI) сводка (q4r.fv) </p>
fmq4r.fv <-q4r.fv $ pprob </p>
write.csv (fmq4r.fv, file = "fmq4r.fv", row.names = F)
С помощью файла csv я получаю следующее.Ninclu должен быть моим идентификатором, но он больше не соответствует моим исходным переменным minclu, которые являются строковой переменной идентификатора участника
print.data.frame (fmq4r.fv) ninclu class prob1 prob2 prob3 prob4 1 1 27.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-071.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