[R] binomial errors in split-plot design

[Christoph Scherber]

Dear all,

I have proportion data with binomial errors. The problem is that the 
whole experiment was laid out as a split-plot design.

Ideally, what I would like is having a glm with an Error term such as 
glm(y~x+Error(A/B)) but I fear this is not possible. Would using lme be 
an alternative? How could I state the variance structure, then?

I would very much appreciate any suggestions!
Best regards
Christoph