Ajuste Regressão Robusta
Boa tarde. Preciso fazer diagnósticos de resíduos para uma regressão robusta. Como por exemplo, plotar a distância de cook para cada observação amostral. No ajuste com "lm", basta rodar: plot(cooks.distance(lm(y~x))) Já na Regressão Robusta, usando o comando "rlm" da biblioteca MASS, o comando gráfico não roda: plot(cooks.distance(rlm(y~x))) Tem algum comando no R? -- ______________________________________________ *MANOEL VITOR DE SOUZA VELOSO----------------------------------------* *ProfessorInstituto de Ciências Sociais AplicadasUniversidade Federal de Alfenas - MGCampus Avançado de Varginha-MGhttp://www.unifal-mg.edu.br/icsa/ <http://www.unifal-mg.edu.br/icsa/>-------Doutor* * em Estatística e ExperimentaçãoUniversidade Federal de Lavras - MG----------------------------------------Tel: 35 - 8865.6116 35 - 3219.8705* P Antes de imprimir pense em sua responsabilidade e compromisso com o* MEIO AMBIENTE* !!!
Caro Manoel, Se existe o método cooks.distance() definido para objetos de classe "rlm", então você facilmente poderia extrair. A minha suspeita é de que não existe, conforme se pode verificar pelo código abaixo. require(MASS) help(rlm, h="html") m0 <- rlm(stack.loss ~ ., stackloss) class(m0) methods(class="rlm") Agora surgem as questões: 1) distância de cook para essa classe de modelos é algo definido e o autor não incluiu na implementação? ou 2) é algo ainda não definido? Tente fazer uma busca com RSiteSearch("robust cooks distance") E termos afins. À disposição. Walmes.
participantes (2)
-
Manoel Vítor de Souza Veloso -
walmes .