o q vc precisa ver e' o que sao os preditores... e entende-los a fundo... por exemplo, se vc tem uma coluna dizendo:sexo masculino? (1 - sim / 0 - nao)e outra coluna:sexo feminino? (1 - sim / 0 - nao)apesar de textos diferentes, elas contem a mesma informacao... e vc acaba por obter os problemas reportados....definitivamente vc precisa efetuar uma analise descritiva para entender o que esta' acontecendo.2012/12/14 Natália Bordin Barbieri <nataliabordinbarbieri@gmail.com>
O que eu posso fazer?Estou num ponto que não sei qual é o próximo passo.Meu n=400, e esses preditores binários é mais ou menos n/2=o e n/2= 1 em cada variável.Obrigada desde já,NatáliaEm 14 de dezembro de 2012 12:01, Benilton Carvalho <beniltoncarvalho@gmail.com> escreveu:
combinacoes lineares nos preditores binarios....2012/12/14 Natália Bordin Barbieri <nataliabordinbarbieri@gmail.com>
Para testar, verifiquei colinearidade e se a matriz era inversível. Ambos estão Ok.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.Testando apenas a variável resposta binário com um preditor binário, o erro que apresenta é o seguinte:Initial scale 0 because more than 'h' (=201) observations are identical.Erro em solve.default(V) :rotina Lapack dgesv: sistema é exatamente singularAlguma ideia?Att., NatáliaEm 14 de dezembro de 2012 09:07, Benilton Carvalho <beniltoncarvalho@gmail.com> escreveu:
nada reproduzivel...verifique se algumas de suas variaveis nao sao combinacoes lineares de outras.On 14 December 2012 10:52, Natália Bordin Barbieri <nataliabordinbarbieri@gmail.com> wrote:
_______________________________________________Olá,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:require(robust)DM.glmRob3<-glmRob(DM ~ INFLAMACAO + IDADE + SEXOM + RACACOR + HIPERT + IMC + RCQ + TRIGT, family=binomial, data=DM, weights=NULL,method ="mallows", model = TRUE, control = glmRob.control)Porém, ao rodar, o que retorna é o seguinte erro:'The covariance matrix has become singular duringthe iterations of the MCD algorithm.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'Alguém já trabalhou com isso?Sabe me ajudar?Obrigada, Natália
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
--
Att.,Natália B. Barbieri
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
--
Att.,Natália B. Barbieri
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.