Hello! hier one question about the limits in the number of random effects specified by the random part in lme-Call. could the Call to lme(...) realy fail only becouse of too large random-structure part in lme(..., random=....) ? I really would need many random parameters in the model, to be specified for example by the following: lme(..., random=list(my.Subject=pdBlocked(list(pdIdent(value=diag(vector_1),form=,nam=), pdIdent(value=diag(vector_2),form=,nam=), pdIdent(value=diag(vector_3),form=,nam=), --------------------------------------------------, pdIdent(value=diag(vector_30),form=,nam=) ))), ...) Has anybody got such a problem with lme(...) or glmmPQL(...) or relative problems? Thanks a lot in advance.