pouco provavel q o glmRob nao aceite dicotomicos... o problema e' muito provavelmente (pra nao dizer a unica causa possivel) que uma coluna seja combinacao linear de outra.<div><br></div><div>se vc nao estiver disposta a fazer analise descritiva, etc.... pelo menos va' adicionando uma variavel por vez... vai ter um momento q o ajuste vai falhar e vc sabera' qual foi a variavel q causou o problema...</div>
<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">Bom, no meu banco, tenho: <div><br></div><div>Variáveis dicotômicas (0 ou 1) :(DM , INFLAMACAO , HIPERT , SEXOM ,RACACOR) e variáveis quantitativas:(IDADE, IMC, INSULINA,TRIGT, RCQ)</div>
<div>Eu quero ajustar minha resposta (DM) por todos os preditores.</div>
<div>O que acontece é que quando entro com qualquer um dos preditores dicotômicos no modelo, (ex: DM ~ IDADE+IMC+INSULINA+TRIGT+ RCQ+<u> INFLAMACAO</u>) , ele gera o erro citado.</div><div>Está parecendo que o glmRob não aceita preditores dicotômicos. Será que pode ser isso? </div>
<div><br></div><div><br></div><div><br></div><div><br></div><div><div><br></div><div><br></div><div><br><div class="gmail_quote">Em 14 de dezembro de 2012 13:25, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com" target="_blank">beniltoncarvalho@gmail.com</a>></span> escreveu:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">o q vc precisa ver e' o que sao os preditores... e entende-los a fundo... por exemplo, se vc tem uma coluna dizendo:<div>
<br></div><div>sexo masculino? (1 - sim / 0 - nao)</div><div><br></div><div>e outra coluna:</div>
<div><br></div><div>sexo feminino? (1 - sim / 0 - nao)</div><div><br></div><div>apesar de textos diferentes, elas contem a mesma informacao... e vc acaba por obter os problemas reportados....</div><div><br></div><div>definitivamente vc precisa efetuar uma analise descritiva para entender o que esta' acontecendo.</div>
<div><div>
<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">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:<div>
<div><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><div><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" 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><div><div><br><br clear="all">
<div><br></div>-- <br>
Att.,<div>Natália B. Barbieri</div><br>
</div></div></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" 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><div><div class="h5"><br><br clear="all">
<div><br></div>-- <br>
Att.,<div>Natália B. Barbieri</div><br>
</div></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></div>