Dear All, I was wondering if it might be possible in my below model to set "D1:Treat" and "D2:Treat" to be equal? lme(value ~ 0 + D1 + D2 + D1:Treat + D2:Treat, random = ~0 + D1 + D2 | Student, data = data) Thank you, Simon [[alternative HTML version deleted]]