Para testar, verifiquei colinearidade e se a matriz era inversível. Ambos estão Ok.<div>O que acontece, é que quando coloco os preditores binário, o modelo apresenta o erro já mostrado. Quando é apenas com os preditores quantitativos, o método funciona. </div>
<div>Testando apenas a variável resposta binário com um preditor binário, o erro que apresenta é o seguinte:</div><div><br></div><div><div><b>Initial scale 0 because more than 'h' (=201) observations are identical.</b></div>
<div><b>Erro em solve.default(V) : </b></div><div><b> rotina Lapack dgesv: sistema é exatamente singular</b></div><div><br></div><div>Alguma ideia?</div><div>Att., Natália</div><div><br></div><div class="gmail_quote">Em 14 de dezembro de 2012 09:07, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com" target="_blank">beniltoncarvalho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">nada reproduzivel...<div><br></div><div>verifique se algumas de suas variaveis nao sao combinacoes lineares de outras.</div>
<div><br></div><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On 14 December 2012 10:52, Natália Bordin Barbieri <span dir="ltr"><<a href="mailto:nataliabordinbarbieri@gmail.com" target="_blank">nataliabordinbarbieri@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">Olá, <div><br></div><div>preciso ajustar o seguinte modelo de estimação robusta, através do método de "Mallows' para o modelo de regressão logística que segue:</div>
<div><div><br></div><div><br></div><div>require(robust)</div>
<div>DM.glmRob3<-glmRob(DM ~ INFLAMACAO + IDADE + SEXOM + RACACOR + HIPERT + IMC + RCQ + TRIGT, family=binomial, data=DM, weights=NULL,</div><div> method ="mallows", model = TRUE, control = glmRob.control)</div>
<div><br></div><div>Porém, ao rodar, o que retorna é o seguinte erro:</div><div><br></div><div><div>'The covariance matrix has become singular during</div><div>the iterations of the MCD algorithm.</div><div>There are 284 observations (in the entire dataset of 400 obs.) lying on the hyperplane with equation a_1*(x_i1 - m_1) + ... + a_p*(x_ip - m_p) = 0 with (m_1,...,m_p) the mean of these observations and coefficients a_i from the vector a <- c(1, 0, 0, 0, 0, 0, 0, 0) Erro em solve.default(V) rotina Lapack dgesv: sistema é exatamente singular'</div>
</div><div><br></div><div>Alguém já trabalhou com isso?Sabe me ajudar?</div><div>Obrigada, Natália</div>
</div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
Att.,<div>Natália B. Barbieri</div><br>
</div>