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