Jak utworzyć wykres resztkowy w sas-ie

Wykresy reszt są często używane do oceny, czy reszty modelu regresji mają rozkład normalny i czy wykazują heteroskedastyczność .

Możesz użyć następującej podstawowej składni, aby dopasować model regresji i utworzyć wykres resztowy dla modelu w SAS-owym stylu życia:

 symbol value = circle;

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

Poniższy przykład pokazuje, jak zastosować tę składnię w praktyce.

Uwaga : Instrukcja symbol określa, że chcemy wyświetlić pozostałe punkty wykresu jako okręgi. Domyślną formą jest znak plus.

Przykład: Utwórz działkę resztkową w SAS-ie

Załóżmy, że mamy następujący zestaw danych w SAS-ie:

 /*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;

Możemy użyć następującej składni, aby dopasować prosty model regresji liniowej do tego zbioru danych i utworzyć wykres reszt w celu wizualizacji reszt w stosunku do przewidywanych wartości:

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

Wykres pozostałości zostanie wyświetlony na dole wyniku:

Oś x wyświetla przewidywane wartości, a oś y wyświetla reszty.

Ponieważ reszty są losowo rozproszone wokół wartości zerowej, bez wyraźnego trendu wzrostu lub spadku wariancji, założenie o homoskedastyczności reszt jest spełnione.

Na górze wykresu możemy również zobaczyć dopasowane równanie regresji:

y = 29,631 + 0,7553x

Po prawej stronie wykresu możemy zobaczyć także następujące metryki modelu regresji:

• N : Całkowita liczba obserwacji (12)
• Rsq : model R-kwadrat (0,6324)
• AdjRsq : skorygowany model R-kwadrat (0,5956)
• RMSE : Pierwiastek błędu średniokwadratowego modelu (4,4417)

Dodatkowe zasoby

Poniższe samouczki wyjaśniają, jak wykonywać inne typowe zadania w SAS-ie:

Jak tworzyć histogramy w SAS-ie
Jak tworzyć chmury punktów w SAS-ie
Jak zidentyfikować wartości odstające w SAS