[R-br] Teste Bonferonni-Detecção Outlier

[Edson Lira]

Outra opção é fazer:

m1<-lm(y~x)
par(mfrow=c(2,2))
plot(m1)

Esta rotina vai te mostrar 4 gráficos:

O primeiro, você poderá avaliar a variabilidade dos dados;
O segundo, a normalidade dos dados;
O terceiro, o resíduo padronizado (acho que o ideal é estar dentro de -3 a +3), caso contrário é um problema;
O quarto, você terá a distância de Cook junto com os leverages, aonde os pontos localizados acima de 1 nas retas são considerados outliers. (Acho que para R é 1, é bom checar)

Por outro lado, faça as análises com os pontos que você acha que são outliers, observe R2, erro padrrão, etc e depois retire os "candidatos" e verifique como se comportou o modelo sem os "candidatos".


Espero ter ajudado.
[  ]'s. 

Edson Lira
Estatístico
Manaus-Amazonas



Em Quarta-feira, 13 de Novembro de 2013 9:34, Cesar Rabak <cesar.rabak em gmail.com> escreveu:
 
Você pulou a informação mais relevante do diagnóstico:

"No Studentized residuals with Bonferonni p < 0.05".

Portanto se você aceitar as premissas do teste (p-valor corrigido pela técnica do Bonferroni <=0,05, etc.), nenhum dado nesse teste seria um outlier.

Ademais, o maior valor em módulo do resíduo "studentizado" gera um p-valor de 0,12. A menos que você tenha razões científicas fora do teste de diagnóstico para declarar esse dado como aberrante, ele não o seria pelo teste feito.

HTH
--
Cesar Rabak




2013/11/9 Helder Gramacho <agrohelder em gmail.com>

A saída do teste Bonferonni para detecção de outlier apareceu com os resultados um pouco confusos para mim:
>
>
>> outlierTest(modelo1)
>
>
>No Studentized residuals with Bonferonni p < 0.05
>Largest |rstudent|:
>   rstudent unadjusted p-value Bonferonni p
>25 3.200023          0.0030963      0.12385
>
>
>A interpretação abaixo está correta?
>25: é o número da observação
>3.200023: é o rstudent unadjusted
>0.0030963: p-value
>
>0.12385: Bonferonni p
>
>
>
>No caso, para  a observação 25 o p-value indica que esta observação é um outlier, está correto este entendimento?
>
>Hélder Gramacho 
>agrohelder em gmail.com
>
>
>
>_______________________________________________
>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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131127/f87ae84b/attachment.html>