[R-br] Pacote lcmm: Estimação de variância do resíduo em modelo de mistura com duas classes

[Felipe Buchbinder]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140622/e4ffa951/attachment.html>