Das Skript beginnt mit der Erstellung eines 'ossification'-Datensatzes aus eingebetteten Daten (datalines), der die Ergebnisse eines Experiments an Würfen (litters) mit verschiedenen Behandlungen (PHT, TCPO) darstellt. Ein erstes GEE-Modell wird mit einer unabhängigen Arbeitskorrelationsmatrix angepasst. Anschließend werden die Daten in ein binäres Format (ein Datensatz pro Subjekt) im Datensatz 'ossification_b' transformiert. Ein zweites GEE-Modell wird dann an diese transformierten Daten angepasst, diesmal unter Verwendung einer austauschbaren Arbeitskorrelationsmatrix. Die Ergebnisse beider Analysen werden im HTML-Format generiert.
Datenanalyse
Type : CREATION_INTERNE
Die Daten werden vollständig innerhalb des Skripts über eine 'datalines'-Anweisung im ersten DATA-Schritt generiert.
1 Codeblock
DATA STEP Data
Erklärung : Erstellt den Datensatz 'ossification', indem die Daten direkt aus dem Code über 'datalines' gelesen werden. Jede Datenzeile repräsentiert einen Wurf (litter) und die Anzahl der Erfolge (t) bei einer Anzahl von Versuchen (m) für verschiedene Behandlungen (tx). Indikatorvariablen (PHT, TCPO) werden erstellt, um die Behandlungseffekte zu modellieren.
Erklärung : Passt ein verallgemeinertes Schätzungsgleichungsmodell (GEE) an die Daten an. Das Modell verwendet eine Binomialverteilung mit einer Logit-Verknüpfungsfunktion, um den Erfolganteil (t/m) in Abhängigkeit von den Behandlungen zu modellieren. Eine Analyse mit wiederholten Messungen wird an den Würfen (litters) durchgeführt, wobei eine unabhängige Arbeitskorrelationsstruktur (type=ind) angenommen wird.
Kopiert!
ods html;
title "*** Ossification Data -- GEE using GENMOD with Independent Working Correlation Matrix***";
proc genmod data=ossification;
class litters tcpo pht / param=ref;
model t/m = tcpo pht tcpo*pht / dist=binomial link=logit;
repeated subject=litters / type=ind;
run;
ods html close;
1
ods html;
2
title "*** Ossification Data -- GEE using GENMOD with Independent Working Correlation Matrix***";
3
PROC GENMODDATA=ossification;
4
class litters tcpo pht / param=ref;
5
model t/m = tcpo pht tcpo*pht / dist=binomial link=logit;
6
repeated subject=litters / type=ind;
7
RUN;
8
ods html close;
3 Codeblock
DATA STEP Data
Erklärung : Transformiert den aggregierten Datensatz 'ossification' in ein binäres Format. Für jeden Wurf werden 't' Beobachtungen mit einer Antwortvariablen 'y' gleich 1 (Erfolg) und 'm-t' Beobachtungen mit 'y' gleich 0 (Misserfolg) erstellt. Dieses Format ist für bestimmte Arten der Binomialanalyse erforderlich.
Kopiert!
data ossification_b;
set ossification;
do i=1 to t;
y = 1;
output;
end;
do i=1 to m-t;
y = 0;
output;
end;
run;
1
DATA ossification_b;
2
SET ossification;
3
DO i=1 to t;
4
y = 1;
5
OUTPUT;
6
END;
7
DO i=1 to m-t;
8
y = 0;
9
OUTPUT;
10
END;
11
RUN;
4 Codeblock
PROC GENMOD
Erklärung : Passt ein zweites GEE-Modell unter Verwendung des binären Datensatzes 'ossification_b' an. Das Modell ist dem ersten ähnlich, spezifiziert jedoch eine austauschbare Arbeitskorrelationsstruktur (type=exch), die eine konstante Korrelation zwischen allen Beobachtungen innerhalb desselben Wurfs annimmt.
Kopiert!
ods html;
title "*** Ossification Data -- GEE using GENMOD with Exchangeable Working Correlation Matrix***";
proc genmod data=ossification_b descending;
class litters tcpo pht / param=ref;
model y = tcpo pht tcpo*pht / dist=binomial link=logit;
repeated subject=litters / type=exch;
run;
ods html close;
1
ods html;
2
title "*** Ossification Data -- GEE using GENMOD with Exchangeable Working Correlation Matrix***";
3
PROC GENMODDATA=ossification_b descending;
4
class litters tcpo pht / param=ref;
5
model y = tcpo pht tcpo*pht / dist=binomial link=logit;
6
repeated subject=litters / type=exch;
7
RUN;
8
ods html close;
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.
Zugehörige Dokumentation
Aucune documentation spécifique pour cette catégorie.
SAS und alle anderen Produkt- oder Dienstleistungsnamen von SAS Institute Inc. sind eingetragene Marken oder Marken von SAS Institute Inc. in den USA und anderen Ländern. ® zeigt die Registrierung in den USA an. WeAreCAS ist eine unabhängige Community-Site und nicht mit SAS Institute Inc. verbunden.
Diese Website verwendet technische und analytische Cookies, um Ihre Erfahrung zu verbessern.
Mehr erfahren.
