Buenas,
La instrucción plot admite un parámetro which, su valor por defecto es c(1, 2,
3, 5), que son los cuatro gráficos que salen por defecto.
Usando which = 2 obtienes únicamente el qqplot.
Las observaciones leverage son observaciones que tienen una gran influencia en
los resultados del análisis, de forma que eliminarlas supone
un cambio importante en las estimaciones. No hay que confundirlas con los
outliers o atípicos, pues una observación puede ser atípica y no
leverage, o ser leverage y no atípica.
Un saludo
--
--------------------------------------------------
http://knuth.uca.es/R
--------------------------------------------------
Proyecto R-UCA
--------------------------------------------------
Nombre: Manuel Muñoz Márquez
Departamento: Departamento de Estadística e Investigación Operativa
Institución: Escuela Superior de Ingeniería
Organización: Universidad de Cádiz
--------------------------------------------------
El vie, 16-06-2023 a las 12:41 +0200, Yesica Pallavicini Fernandez
escribió:> Buenos días y muchas gracias por adelantado
> En cuando a probar gráficamente los supuestos de normalidad y
> homocedasticidad de una ANOVA
>
> Lo estoy haciendo de forma gráfica de la siguiente manera:
> Primero hago el nova con: aov()
> luego hago la comprobación de los supuestos con: plot(modelo)
> y me salen 4 gráficos;
> 1 un scatter plot de los "residuals" vs " Fitted",
> 2"standardizez residuals" vs " Fitted",
> 3QQplot,
> 4 Residuals vs leverage
> Pero ¿Qué es el leverage?
> 1)¿teneis algun script que reemplace a plot(modelo) y que solo contenga el
> QQplot y los valores residuales frente a los fitted?
>
> 2) Shapiro.test () no funciona para observaciones mayores a 500 y mis datos
> tienen más de 30000 observaciones. ¿sabeis de algun otro test para este
> caso?
>
> Muchísimas gracias y que acabéis fenomenal la semana
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
>
https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-help-es__;!!D9dNQwwGXtA!S2L6bI_NoCf_ox_hHI7PINlx5SKlSS1KpBOsabN8o1Z248Rp3TRh5axwIvBUT1g9XZtvU7qL25KUdwc$
>