O que eu posso fazer?<div>Estou num ponto que não sei qual é o próximo passo.</div><div>Meu n=400, e esses preditores binários é<u> mais ou menos</u> n/2=o e n/2= 1 em cada variável.</div><div>Obrigada desde já, </div><div>
Natália<br><br><div class="gmail_quote">Em 14 de dezembro de 2012 12:01, 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">combinacoes lineares nos preditores binarios....<div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">
2012/12/14 Natália Bordin Barbieri <span dir="ltr"><<a href="mailto:nataliabordinbarbieri@gmail.com" target="_blank">nataliabordinbarbieri@gmail.com</a>></span><br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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>
<div><b>Erro em solve.default(V) : </b></div><div><b>  rotina Lapack dgesv: sistema é exatamente singular</b></div><div><br></div></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:<div><div><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>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>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" 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></div></div><span><font color="#888888"><br>


<br clear="all"><div><br></div>-- <br>
Att.,<div>Natália B. Barbieri</div><br>
</font></span></div>
<br>_______________________________________________<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>
</div></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>