/* schizophrenia example for GEE, lesson 9*/

options nocenter nodate nonumber linesize=72; 
data schiz; 
infile "D:\Statistics\STAT 504\datasets\schiz.dat"; 
input id drug week y; 
sqrtweek = sqrt(week); 
run; 

proc genmod data=schiz; 
class id; 
model y = drug sqrtweek drug*sqrtweek; 
repeated subject=id / type=ind modelse; 
run;