Help! My lme function had been running for more than 24 hours...

My random effect part is similar to those in Professor Bate's draft book.

However, I have the following heteroskedacity structure:

 MyCovStructure=varIdent(form=~1|MyGroups),

And there are around 500 groups...

Any pointers about how to improve the speed?

Thanks a lot!

	[[alternative HTML version deleted]]

