Exemple d'initiation à PROC MCMC

Le script effectue une analyse de régression linéaire simple en utilisant une approche bayésienne via PROC MCMC. Il définit des distributions a priori pour les coefficients de régression (beta0, beta1) et la variance de l'erreur (sigma2), puis modélise la variable 'weight' en fonction de 'height'. L'échantillonnage de la chaîne de Markov Monte Carlo est configuré avec 10000 itérations et un pas de sélection de 2. Les résultats de la chaîne MCMC sont sauvegardés dans un ensemble de données appelé 'classout'. La production graphique est activée pendant l'exécution de la procédure.

Analyse des données

Type : SASHELP

Les données source proviennent de l'ensemble de données intégré 'SASHELP.CLASS', qui est un jeu de données SAS standard disponible dans l'environnement SAS.

1 Bloc de code

ODS

Explication :
Active le système de destination de sortie (ODS) pour la production de graphiques, permettant à PROC MCMC de générer des visualisations des résultats de l'analyse.

Copié !

1	ods graphics on;

2 Bloc de code

PROC MCMC

Explication :
Ce bloc exécute la procédure PROC MCMC pour une analyse de régression bayésienne. Il utilise 'sashelp.class' comme données d'entrée et crée un ensemble de données de sortie 'classout' contenant les échantillons MCMC. Les paramètres 'beta0', 'beta1' et 'sigma2' sont déclarés avec des valeurs initiales. Des distributions a priori (normale pour les betas, inverse gamma pour sigma2) sont spécifiées. La moyenne linéaire 'mu' est définie comme une fonction de 'height', et le modèle spécifie que 'weight' suit une distribution normale avec cette moyenne et la variance 'sigma2'. Les options 'nmc', 'thin' et 'seed' contrôlent l'échantillonnage MCMC.

Copié !

1	PROC MCMC DATA=sashelp.class outpost=classout nmc=10000 thin=2 seed=246810;
2	parms beta0 0 beta1 0;
3	parms sigma2 1;
4	prior beta0 beta1 ~ normal(mean = 0, var = 1e6);
5	prior sigma2 ~ igamma(shape = 3/10, scale = 10/3);
6	mu = beta0 + beta1*height;
7	model weight ~ n(mu, var = sigma2);
8	RUN;

3 Bloc de code

ODS

Explication :
Désactive la production graphique via le système de destination de sortie (ODS), mettant fin à la génération de graphiques après l'exécution de PROC MCMC.

Copié !

1	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