options ls=80 ps=1000 nodate;

/*------------------------------------------------------*/
/*                                                      */
/* Proc Glimmix to fit subject-specific (random effect) */
/* model  for the indonesian children infection disease */
/* data                                                 */
/*                                                      */
/*------------------------------------------------------*/

data indon; 
  infile "indon.dat";
  input id infect intercep age xero cosv sinv sex height stunted
  visit baseage season visitsq;
  time = age - baseage;
run;

title "Random intercept model for infection disease data";
proc glimmix data=indon method=quad;
  class id;
  model infect = season xero age time sex height / dist=bin link=logit s;
  random int / subject=id type=vc;
run;