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:
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
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