Chapter 16, Section 16.5
Marginal Logistic Regression Model
data ecg;
infile 'ecg.dat';
input id seq period trt y;
title1 Marginal Logistic Regression Model;
title2 Crossover Trial on Cerebrovascular Deficiency;
proc genmod descending;
class id;
model y = trt period / d=bin;
repeated subject=id / logor=fullclust;
run;
<Selected Output>
Mixed Effects Logistic Regression Model (Random Intercept)
title1 Mixed Effects Logistic Regression Model (Random Intercept);
title2 Crossover Trial on Cerebrovascular Deficiency;
proc glimmix method=quad(qpoints=100);
class id;
model y = trt period / dist=binomial link=logit solution;
random intercept / subject=id g;
run;