[R-br] Mallows - estimador robusto

Natália Bordin Barbieri nataliabordinbarbieri em gmail.com
Sexta Dezembro 14 12:28:43 BRST 2012


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ália

Em 14 de dezembro de 2012 12:01, Benilton Carvalho <
beniltoncarvalho em gmail.com> escreveu:

> combinacoes lineares nos preditores binarios....
>
>
> 2012/12/14 Natália Bordin Barbieri <nataliabordinbarbieri em 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 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
>>
>>
>> _______________________________________________
>> 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/3d87aa0d/attachment.html>


Mais detalhes sobre a lista de discussão R-br