Exemple détaillé de PROC RSREG

Le script commence par créer un jeu de données initial (DATA A) à partir de données en ligne. Il étend ensuite ce jeu de données (DATA B) en ajoutant des observations avec des valeurs manquantes pour les variables de réponse sur une grille de valeurs pour les prédicteurs, afin de faciliter les prédictions. La procédure PROC RSREG est utilisée une première fois pour ajuster un modèle aux trois variables de réponse (y1, y2, y3) en fonction de deux prédicteurs (x1, x2), et les prédictions sont sauvegardées dans un nouveau jeu de données (OUT=C). Ensuite, ce jeu de données de prédiction est filtré, trié et les 5 premières observations sont affichées (DATA D, PROC SORT, PROC PRINT). Enfin, une deuxième invocation de PROC RSREG est utilisée avec ODS GRAPHICS pour générer et afficher des graphiques de surface pour les relations entre y1, y2 et x1, x2.

Analyse des données

Type : CREATION_INTERNE

Le jeu de données 'a' est créé directement dans le script via une instruction DATALINES. Les jeux de données 'b', 'c' et 'd' sont dérivés de 'a' ou de résultats de PROC RSREG.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée un jeu de données nommé 'a' en lisant des données brutes en ligne (DATALINES). Il initialise les variables x1, x2, y1, y2 et y3 avec 16 observations.

Copié !

1	DATA a;
2	INPUT x1 x2 y1 y2 y3;
3	DATALINES;
4	-1 -1 1.8 1.940 3.6398
5	-1 1 2.6 1.843 4.9123
6	1 -1 5.4 1.063 6.0128
7	1 1 0.7 1.639 2.3629
8	0 0 8.5 0.134 9.0910
9	0 0 3.0 0.545 3.7349
10	0 0 9.8 0.453 10.4412
11	0 0 4.1 1.117 5.0042
12	0 0 4.8 1.690 6.6245
13	0 0 5.9 1.165 6.9420
14	0 0 7.3 1.013 8.7442
15	0 0 9.3 1.179 10.2762
16	1.4142 0 3.9 0.945 5.0245
17	-1.4142 0 1.7 0.333 2.4041
18	0 1.4142 3.0 1.869 5.2695
19	0 -1.4142 5.7 0.099 5.4346
20	;

2 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée le jeu de données 'b' en copiant d'abord toutes les observations de 'a'. Ensuite, à la fin du traitement de 'a' (lorsque 'eof' est vrai), il génère un grand nombre de nouvelles observations en créant une grille de valeurs pour x1 et x2 (de -2 à 2 avec un pas de 0.1) et en définissant les variables de réponse (y1, y2, y3) comme manquantes. Ces observations supplémentaires serviront de points pour la prédiction par PROC RSREG.

Copié !

1	DATA b;
2	SET a END=eof;
3	OUTPUT;
4	IF eof THEN DO;
5	y1=.;
6	y2=.;
7	y3=.;
8	DO x1=-2 to 2 BY .1;
9	DO x2=-2 to 2 BY .1;
10	OUTPUT;
11	END;
12	END;
13	END;
14	RUN;

3 Bloc de code

PROC RSREG Data

Explication :
Cette PROC RSREG exécute une analyse de régression de surface de réponse sur le jeu de données 'b'. Elle modélise les trois variables de réponse (y1, y2, y3) en fonction des prédicteurs x1 et x2. L'option 'predict' demande au programme de calculer les valeurs prédites pour toutes les observations de 'b', y compris celles avec des réponses manquantes ajoutées précédemment. Le jeu de données 'c' est créé pour stocker les résultats, y compris les prédictions.

Copié !

1
2	PROC RSREG
3	DATA=b out=c;
4	model y1 y2 y3=x1 x2 / predict;
5	RUN;
6

4 Bloc de code

DATA STEP et PROC SORT/PRINT Data

Explication :
Ce bloc combine plusieurs étapes de manipulation et d'affichage de données. Le premier DATA STEP 'd' est créé à partir de 'c', en filtrant les observations où y2 est inférieur à 2 ET y3 est inférieur à la somme de y2 et y1. Ensuite, PROC SORT trie ce jeu de données 'd' par ordre décroissant de y1. Le deuxième DATA STEP sur 'd' conserve uniquement les 5 premières observations après le tri. Enfin, PROC PRINT affiche le contenu final du jeu de données 'd'.

Copié !

1	DATA d;
2	SET c;
3	IF y2<2;
4	IF y3
5
6	PROC SORT DATA=d;
7	BY descending y1;
8	RUN;
9
10	DATA d; SET d;
11	IF (_n_ <= 5);
12	PROC PRINT;
13	RUN;

5 Bloc de code

PROC RSREG

Explication :
Ce bloc utilise PROC RSREG pour générer des visualisations graphiques. 'ods graphics on' active la sortie graphique. La procédure modélise spécifiquement y1 et y2 en fonction de x1 et x2 en utilisant le jeu de données original 'a'. L'option 'plots=surface(overlaypairs)' demande la création de graphiques de surface superposés pour visualiser les surfaces de réponse des deux modèles. 'ods graphics off' désactive la sortie graphique après l'exécution de la procédure.

Copié !

1	ods graphics on;
2	PROC RSREG DATA=a plots=surface(overlaypairs);
3	model y1 y2=x1 x2;
4	RUN;
5	ods graphics off;

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 : /****************************************************************/ /* S A S S A M P L E L I B R A R Y */ /* */ /* NAME: RSRDT */ /* TITLE: Details Example for PROC RSREG */ /* PRODUCT: STAT */ /* SYSTEM: ALL */ /* KEYS: response surface regression */ /* PROCS: RSREG */ /* DATA: */ /* */ /* REF: SAS/STAT User's Guide, PROC RSREG chapter */ /* MISC: */ /* */ /****************************************************************/

Retour à la liste