[R-sig-ME] Fitting crossed random effects in lme (not lmer)

[Christoph Scherber]

Dear all,

Let´s assume I have an lme model with the following random effects structure:

random=~1|block/plot

Now, further treatments are applied at the next smallest level, but they
are on the same hierarchy. In lmer notation, this could be parameterized
as

...+(1|block/plot/treatment.type1)+(1|block/plot/treatment.type2)

How could I specify the same structure in lme?

I know that I would need to use one of the pdClasses(), but I don´t
understand how to set this up correctly. I tried of course the Pinheiro &
Bates book, but maybe someone can give me a hint. That would be excellent!

All the best and many thanks for any help,
Christoph