Prezados, Estou utilizando a função hlme do pacote lcmm para estimar um modelo RMM (*Regression Mixture Model*) com duas classes latentes. Eu gostaria que a função estimasse as variâncias dos resíduos separadamente para cada classe. Porém, o output contém uma única linha para variância dos resíduos. Exemplo:
rmm <- hlme(y ~ x + x2, ng = 2, mixture = y ~ x + x2, subject = 'ID',data = data) summary(rmm)
Fixed effects in the longitudinal model: coef Se Wald p-value intercept class1 -2.26113 0.78246 -2.890 0.00386 intercept class2 0.93273 0.78247 1.192 0.23324 x class1 1.40524 0.32679 4.300 0.00002 x class2 2.09610 0.32679 6.414 0.00000 x2 class1 0.07210 0.02895 2.490 0.01276 x2 class2 0.09123 0.02895 3.151 0.00163 coef se Residual standard error: 0.9408319 0.105216 Como é possível estimar variâncias dos resíduos específicas para cada classe latente? Obrigado! Felipe Buchbinder