Exemple d'introduction à PROC ROBUSTREG

Niveau de difficulté
Débutant
Publié le :
Le script commence par créer le jeu de données 'Stack' contenant 21 observations avec 3 variables explicatives (x1, x2, x3) et une variable réponse (y). Il applique ensuite la procédure PROC ROBUSTREG en trois étapes : une analyse standard avec diagnostics et test de variable, une analyse graphique produisant divers plots de diagnostic (histogrammes, QQ-plots), et une analyse utilisant une méthode d'estimation M spécifique avec la fonction de poids 'bisquare'.
Analyse des données

Type : CREATION_INTERNE


Les données sont intégrées directement dans le script (instream data) via l'instruction DATALINES dans l'étape DATA 'Stack'.

1 Bloc de code
DATA STEP Data
Explication :
Création du dataset 'Stack' par lecture des données incluses (datalines). Chaque ligne contient trois variables prédictives, la variable réponse 'y' et un identifiant d'expérience.
Copié !
1DATA Stack;
2 INPUT x1 x2 x3 y exp $;
3 DATALINES;
480 27 89 42 e1 80 27 88 37 e2
575 25 90 37 e3 62 24 87 28 e4
662 22 87 18 e5 62 23 87 18 e6
762 24 93 19 e7 62 24 93 20 e8
858 23 87 15 e9 58 18 80 14 e10
958 18 89 14 e11 58 17 88 13 e12
1058 18 82 11 e13 58 19 93 12 e14
1150 18 89 8 e15 50 18 86 7 e16
1250 19 72 8 e17 50 19 79 8 e18
1350 20 80 9 e19 56 20 82 15 e20
1470 20 91 15 e21
15;
2 Bloc de code
PROC ROBUSTREG
Explication :
Exécution de la régression robuste (estimation M par défaut). Affiche les diagnostics pour les valeurs aberrantes et les points de levier, et effectue un test de significativité pour la variable 'x3'.
Copié !
1PROC ROBUSTREG DATA=stack;
2 model y = x1 x2 x3 / diagnostics leverage;
3 id exp;
4 test x3;
5RUN;
3 Bloc de code
PROC ROBUSTREG
Explication :
Activation des graphiques ODS pour générer des visualisations diagnostiques (RD-plot, DD-plot, histogramme des résidus, QQ-plot) afin d'évaluer la qualité de l'ajustement.
Copié !
1ods graphics on;
2 
3PROC ROBUSTREG DATA=stack plots=(rdplot ddplot histogram qqplot);
4 model y = x1 x2 x3;
5RUN;
6 
7ods graphics off;
4 Bloc de code
PROC ROBUSTREG
Explication :
Réexécution de la régression robuste en spécifiant explicitement la méthode d'estimation M avec la fonction de poids 'bisquare' et un paramètre de réglage c=3.5.
Copié !
1PROC ROBUSTREG method=m(wf=bisquare(c=3.5)) DATA=stack;
2 model y = x1 x2 x3 / diagnostics leverage;
3 id exp;
4 test x3;
5RUN;
Ce matériel est fourni "tel quel" par We Are Cas. Il n'y a aucune garantie, expresse ou implicite, quant à la qualité marchande ou à l'adéquation à un usage particulier concernant le matériel ou le code contenu dans les présentes. We Are Cas n'est pas responsable des erreurs dans ce matériel tel qu'il existe maintenant ou existera, et We Are Cas ne fournit pas de support technique pour celui-ci.
Informations de Copyright : SAS SAMPLE LIBRARY, NAME: rreggs1, PRODUCT: STAT