He estado tratando de crear un gráfico QQ usando Excel que imita los gráficos QQ realizados por SAS JMP para datos normalmente distribuidos. Sin embargo, no puedo imitar con precisión las bandas de confianza realizadas por SAS JMP.
He estado calculando los cuantiles normales utilizando las fórmulas proporcionadas en el siguiente enlace:
Detalles estadísticos para el gráfico de cuantiles normales (jmp.com)
Luego, los residuos estandarizados se calculan mediante '(X - µ)/s', donde µ es la media de mis datos, s es la desviación estándar de la muestra.
Luego, las bandas de confianza se calculan utilizando las fórmulas proporcionadas por John Fox como se explica en este enlace:
modelo lineal - Bandas de confianza para la línea QQ - Validación cruzada (stackexchange.com)
Esto me proporciona puntos para las bandas de confianza en cada punto de datos. Luego aproximo el punto de las bandas de confianza utilizando un polinomio de segundo orden. Las bandas de confianza que obtengo son bastante cercanas a las de SAS JMP, pero comienzan a desviarse en conjuntos de datos más grandes y especialmente en conjuntos de datos con valores atípicos, consulte la imagen adjunta.
Tengo curiosidad por saber si es posible arrojar algo de luz sobre cómo SAS JMP calcula las bandas de confianza para el gráfico QQ. Me parece que hay alguna estimación de asíntotas para las bandas que está incorporada en el cálculo, pero no he podido encontrar ninguna información sobre cómo realizar estos cálculos.