[R] using lme on multiple datasets in one shot

[maitra at iastate.edu]

Dear list,

I would like to do a huge number of lme's using the same design matrix
(and fixed and random effects). Is it possible to do this efficiently?
Doing otherwise is not an option for my example.

Basically, I am wanting to do the following which is possible using lm:

X <- matrix(rnorm(50),10,5)
Y <- matrix(rnorm(50),10,5)
lm(Y~X)  

with lme. Any suggestions?

Many thanks and best wishes,
Ranjan