Perché R sta tramando i residui standardizzati rispetto ai quantili teorici in un diagramma QQ?

In R, perché le impostazioni predefinite qqplot(linear model)utilizzano i residui standardizzati sull'asse y? Perché R non usa i residui "normali"?

Quando si utilizzano i residui standardizzati, il valore atteso dei residui è zero e la varianza è (approssimativamente) di uno. Questo ha due vantaggi:

1. Se ridimensioni una delle tue variabili (ad esempio cambia chilometri in miglia), i grafici residui rimangono invariati.
2. Nel qqplot, i residui dovrebbero trovarsi sulla linea y = x
3. Ti aspetti che il 95% dei tuoi residui sia compreso tra -1,96 e 1,96. Questo rende più facile individuare i valori anomali.

I residui teorici in un modello lineare sono indipendenti identicamente normalmente distribuiti. Tuttavia i residui osservati non sono indipendenti e non presentano varianze uguali. Quindi la standardizzazione dei residui si divide per la deviazione standard stimata associata a quel residuo rendendoli più uguali nelle loro varianze (utilizzando le informazioni dalla matrice del cappello per calcolare questo). Questo è un residuo più significativo da guardare nel qqplot.

Inoltre, stai davvero eseguendo qqplot sul modello montato? o è questo il qqplot dal diagramma di esecuzione sul modello?