Come creare un grafico residuo in sas

I grafici dei residui vengono spesso utilizzati per valutare se i residui di un modello di regressione sono distribuiti normalmente e se presentano o meno eteroschedasticità .

È possibile utilizzare la seguente sintassi di base per adattare un modello di regressione e produrre un grafico dei residui per il modello in SAS:

 symbol value = circle;

proc reg data=my_data;
    model y = x;
    plot residual. * predicted.;
run ;

L’esempio seguente mostra come utilizzare questa sintassi nella pratica.

Nota : la dichiarazione del simbolo specifica che vogliamo visualizzare i punti residui del grafico come cerchi. Il modulo predefinito è un segno più.

Esempio: creare un grafico residuo in SAS

Supponiamo di avere il seguente set di dati in SAS:

 /*create dataset*/
data my_data;
    input xy;
    datalines ;
8 41
12 42
12 39
13 37
14 35
16 39
17 45
22 46
24 39
26 49
29 55
30 57
;
run ;

/*view dataset*/
proc print data =my_data;

Possiamo utilizzare la seguente sintassi per adattare un semplice modello di regressione lineare a questo set di dati e creare un grafico dei residui per visualizzare i residui rispetto ai valori previsti:

 /*fit simple linear regression model and create residual plot*/
symbol value = circle;
proc reg data =my_data;
   model y = x;
   plot residual. * predicted.;
run ;

Il grafico residuo verrà visualizzato nella parte inferiore dell’output:

L’asse x mostra i valori previsti e l’asse y mostra i residui.

Poiché i residui sono sparsi in modo casuale attorno al valore zero senza una chiara tendenza alla varianza crescente o decrescente, l’ipotesi di omoschedasticità dei residui è soddisfatta.

Nella parte superiore del grafico possiamo anche vedere l’equazione di regressione adattata:

y = 29,631 + 0,7553x

E sul lato destro del grafico possiamo anche vedere i seguenti parametri per il modello di regressione:

• N : numero totale di osservazioni (12)
• Rsq : modello R quadrato (0,6324)
• AdjRsq : R quadrato corretto dal modello (0,5956)
• RMSE : errore quadratico medio del modello (4.4417)

Risorse addizionali

I seguenti tutorial spiegano come eseguire altre attività comuni in SAS:

Come creare istogrammi in SAS
Come creare nuvole di punti in SAS
Come identificare i valori anomali in SAS