#R
m3.glmm <- glmer(y ~ timebefore + timeafter + x1 + x2 +...+ x11 +
                     (1+timebefore+timeafter|id),
                 data=data, family=binomial(link="logit"), nAGQ=3)

#SAS
proc glimmix data=data METHOD=QUAD(QPOINTS=3) NOCLPRINT ;
  class id x2 x3 x4 x5;
  model y(event='1')=timebefore timeafter x1 x2 x3 x4 x5
        x6 x7  x8 x9 x10 x11 /solution CL link = logit dist = binary;
  random intercept timebefore timeafter/subject = id GCORR SOLUTION;
run;