Sempre que fazemos uma regressão linear, precisamos verificar se as hipóteses assumidas para o modelo estão corretas. Uma das melhores maneiras de fazer isto é através de gráficos de diagnóstico. Veja o exemplo abaixo:
ajuste <- lm(Petal.Width ~ Petal.Length, data=iris)
library(ggfortify)
autoplot(ajuste)
São quatro os gráficos de diagnóstico produzidos pela função autoplot
. Alguns dos pontos destes gráficos são identificados como desviantes das hipóteses formuladas. Por exemplo, no QQ plot acima, os pontos 115, 135 e 142 são identificados como fora do esperado para os resíduos caso eles fossem distribuídos de acordo com a normal.
Existe alguma maneira de fazer esta identificação de forma automática no R? Como eu poderia pegar o output do autoplot
(ou da própria função nativa plot
do R) e identificar, para cada gráfico plotado, quais pontos violam as hipóteses do modelo?