Exemple de PROC NLIN avec Profilage de Paramètres et Bootstrap

Le script ajuste un modèle exponentiel sur un jeu de données. Il utilise la procédure NLIN pour non seulement estimer les paramètres du modèle, mais aussi pour évaluer leur stabilité et leur sensibilité. La déclaration PROFILE est utilisée pour examiner la vraisemblance du profil de deux des paramètres, tandis que la déclaration BOOTSTRAP effectue une analyse de rééchantillonnage pour fournir des intervalles de confiance robustes et des diagnostics graphiques sur les estimateurs.

Analyse des données

Type : CREATION_INTERNE

Le jeu de données 'clarke1987a' est créé directement dans le script à l'aide d'une étape DATA et d'une instruction 'datalines'.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc de code crée la table SAS 'clarke1987a' dans la librairie WORK. La table contient deux variables numériques, 'x' et 'y', dont les données sont lues directement depuis le flux d'entrée via l'instruction 'datalines'.

Copié !

1	DATA clarke1987a;
2	INPUT x y;
3	DATALINES;
4	1 3.183
5	2 3.059
6	3 2.871
7	4 2.622
8	5 2.541
9	6 2.184
10	7 2.110
11	8 2.075
12	9 2.018
13	10 1.903
14	11 1.770
15	12 1.762
16	13 1.550
17	;
18

2 Bloc de code

PROC NLIN

Explication :
Ce bloc active la génération de graphiques ODS, puis exécute la procédure de régression non linéaire (NLIN). Il définit les valeurs initiales des paramètres (parms), spécifie un modèle exponentiel (model), et demande deux analyses avancées : un profilage de la vraisemblance pour 'theta1' et 'theta3' (profile) et une analyse bootstrap avec 2000 rééchantillonnages pour obtenir des intervalles de confiance et des graphiques de diagnostic (bootstrap). Les graphiques ODS sont ensuite désactivés.

Copié !

1	ods graphics on;
2	PROC NLIN DATA=clarke1987a plots(stats=none)=diagnostics;
3	parms theta1=-0.15
4	theta2=2.0
5	theta3=0.80;
6	profile theta1 theta3 / range = -6 to 2 BY 0.2 all;
7	bootstrap / nsamples = 2000 seed=123 bootplots bootci bootcov;
8	model y = theta3 + theta2exp(theta1x);
9	RUN;
10	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

Retour à la liste