Régression de Poisson exacte avec PROC GENMOD

Le script commence par créer un jeu de données nommé 'ingots' qui contient le nombre de lingots défectueux ('Notready') pour différentes combinaisons de temps de chauffage ('Heat') et de temps de trempage ('Soak'), ainsi que le total testé ('Total'). Il calcule ensuite le logarithme du total ('lnTotal') pour l'utiliser comme variable d'offset dans le modèle. La procédure PROC GENMOD est ensuite appelée pour ajuster un modèle de régression de Poisson, où 'Notready' est la variable dépendante. Une analyse exacte conditionnelle est demandée pour les paramètres 'Heat' et 'Soak' afin d'obtenir des estimations et des tests plus précis, ce qui est particulièrement utile avec de petits échantillons.

Analyse des données

Type : CREATION_INTERNE

Les données sont créées directement dans le script à l'aide d'une étape DATA et de l'instruction 'datalines'.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc de données crée la table 'ingots'. Il lit les variables 'Heat', 'Soak', 'Notready' et 'Total' à partir des données en ligne (datalines), puis calcule une nouvelle variable 'lnTotal' qui est le logarithme de 'Total'. Cette table est la source pour l'analyse statistique.

Copié !

1	DATA ingots;
2	INPUT Heat Soak Notready Total @@;
3	lnTotal= log(Total);
4	DATALINES;
5	7 1.0 0 10 14 1.0 0 31 27 1.0 1 56 51 1.0 3 13
6	7 1.7 0 17 14 1.7 0 43 27 1.7 4 44 51 1.7 0 1
7	7 2.2 0 7 14 2.2 2 33 27 2.2 0 21 51 2.2 0 1
8	7 2.8 0 12 14 2.8 0 31 27 2.8 1 22 51 4.0 0 1
9	7 4.0 0 9 14 4.0 0 19 27 4.0 1 16
10	;
11

2 Bloc de code

PROC GENMOD

Explication :
Ce bloc utilise PROC GENMOD pour effectuer une régression de Poisson. 'Heat' et 'Soak' sont définies comme des variables de classe. Le modèle prédit 'Notready' en fonction de 'Heat' et 'Soak', avec 'lnTotal' comme variable d'offset pour normaliser l'exposition. L'instruction 'exact' demande une analyse exacte conditionnelle pour les variables 'Heat' et 'Soak', et 'exactoptions' définit un temps limite pour le statut de cette analyse.

Copié !

1	PROC GENMOD DATA=ingots;
2	class Heat Soak / param=ref;
3	model Notready=Heat Soak / offset=lnTotal dist=Poisson link=log;
4	exact Heat Soak / joint estimate;
5	exactoptions statustime=10;
6	RUN;

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: GENMEX11, TITLE: Example 11 for PROC GENMOD, PRODUCT: STAT

Retour à la liste