4.5: Análisis Residual
- Page ID
- 149824
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
Para verificar la validez de los supuestos utilizados para desarrollar nuestro modelo, podemos aplicar nuevamente las técnicas de análisis residual que empleamos para examinar el modelo de un factor en la Sección 3.4.
Esta función llama a:
> plot(fitted(int00.lm),resid(int00.lm))produce la gráfica mostrada en la Figura 4.2. Vemos que los residuos parecen estar algo uniformemente dispersos alrededor de cero. Al menos, no vemos ningún patrón obvio que nos lleve a pensar que los residuales no se portan bien. En consecuencia, esta trama no nos da ninguna razón para creer que hemos producido un modelo pobre.
La gráfica Q-Q en la Figura 4.3 se genera usando estos comandos:
> qqnorm(resid(int00.lm))
> qqline(resid(int00.lm))Vemos el que los residuos siguen aproximadamente la línea indicada. En esta trama, podemos ver un poco más de un patrón y algunas no linealidades obvias, llevándonos a ser un poco más cautelosos al concluir que los residuos son
Figura 4.2: Los valores ajustados versus residuales para el modelo multifactor desarrollado a partir de los datos Int2000.
normalmente distribuidos. No debemos necesariamente rechazar el modelo basado en esta prueba, pero los resultados deben servir como recordatorio de que todos los modelos son imperfectos.