[R-sig-ME] glmmPQL inquiry

[kennedy otwombe]

Dear R users,

I have longitudinal data that is all binary and i have run it in SAS using the following code without a problem:

proc glimmix data=navs;
class id;
model y(event='1')=x1 x2 x3/solution distribution=binary;
random intercept/subject=id type=cs;
run;

I have also written a code in R for the same analysis but i am getting the following error message (iteration 1
Error in switch(mode(x), "NULL" = structure(NULL, class = "formula"),  : 
  invalid formula). 

My code in R reads as follows:

>fit<-glmmPQL(y~x1+x2+x3, random=~1|id, family=binomial, data=navs)
>summary(fit)

I am not sure where the problem lies but i realise that glmmPQL does not seem to cater for binary distributions and i aint sure how to model the probability of the event Y=1 which is my interest in the data i am assuming. My data looks as follows:

  t id y x0 x1 x2 x3
  1  1 0  0  0  1  1
  2  1 0  0  1  0  1
  3  1 0  0  1  0  1
  1  2 0  0  0  1  0
  2  2 0  1  1  1  0
  3  2 0  1  1  1  0

I will appreciate any ideas from this network.

 
Kennedy.N.Otwombe
School of Statistics & Actuarial Science 
University of the Witwatersrand 
Private Bag 3 
Wits 2050 
Johannesburg 
South Africa


      ____________________________________________________________________________________
Looking for last minute shopping deals?