[R-br] Mallows - estimador robusto

[Natália Bordin Barbieri]

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 singular*

Alguma ideia?
Att., Natália

Em 14 de dezembro de 2012 09:07, Benilton Carvalho <
beniltoncarvalho em 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 em 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 during
>> the 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 em 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 em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121214/3c256eef/attachment.html>