Exakte Poisson-Regression mit PROC GENMOD

Das Skript beginnt mit der Erstellung eines Datensatzes namens 'ingots', der die Anzahl der defekten Barren ('Notready') für verschiedene Kombinationen von Heizzeit ('Heat') und Einweichzeit ('Soak') sowie die Gesamtzahl der getesteten Barren ('Total') enthält. Anschließend wird der Logarithmus des Gesamtwerts ('lnTotal') berechnet, um ihn als Offset-Variable im Modell zu verwenden. Die Prozedur PROC GENMOD wird dann aufgerufen, um ein Poisson-Regressionsmodell anzupassen, wobei 'Notready' die abhängige Variable ist. Eine exakte bedingte Analyse wird für die Parameter 'Heat' und 'Soak' angefordert, um präzisere Schätzungen und Tests zu erhalten, was besonders bei kleinen Stichproben nützlich ist.

Datenanalyse

Type : INTERNE_ERSTELLUNG

Die Daten werden direkt im Skript mit einem DATA-Schritt und der Anweisung 'datalines' erstellt.

1 Codeblock

DATA STEP Data

Erklärung :
Dieser Datenblock erstellt die Tabelle 'ingots'. Er liest die Variablen 'Heat', 'Soak', 'Notready' und 'Total' aus den Inline-Daten (datalines) und berechnet dann eine neue Variable 'lnTotal', die der Logarithmus von 'Total' ist. Diese Tabelle ist die Quelle für die statistische Analyse.

Kopiert!

1	DATA ingots;
2	INPUT Heat Soak Notready Total @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
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 Codeblock

PROC GENMOD

Erklärung :
Dieser Block verwendet PROC GENMOD, um eine Poisson-Regression durchzuführen. 'Heat' und 'Soak' werden als Klassenvariablen definiert. Das Modell prognostiziert 'Notready' basierend auf 'Heat' und 'Soak', wobei 'lnTotal' als Offset-Variable zur Normalisierung der Exposition dient. Die Anweisung 'exact' fordert eine exakte bedingte Analyse für die Variablen 'Heat' und 'Soak' an, und 'exactoptions' definiert eine Zeitbegrenzung für den Status dieser Analyse.

Kopiert!

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;

Dieses Material wird von We Are Cas "wie besehen" zur Verfügung gestellt. Es gibt keine ausdrücklichen oder stillschweigenden Garantien hinsichtlich der Marktgängigkeit oder Eignung für einen bestimmten Zweck in Bezug auf die hierin enthaltenen Materialien oder Codes. We Are Cas ist nicht verantwortlich für Fehler in diesem Material, wie es jetzt existiert oder existieren wird, noch bietet We Are Cas technischen Support dafür an.

Urheberrechtsinformationen : SAS SAMPLE LIBRARY, NAME: GENMEX11, TITLE: Example 11 for PROC GENMOD, PRODUCT: STAT

Retour à la liste