how to fit multivariate linear mixed model with `lme4`

Hello All,

I wonder how to fit multivariate linear mixed model with `lme4`. I fitted
univariate linear mixed models with the following code:

     library(lme4)
     lmer.m1<- lmer(Y1~A*B+(1|Block)+(1|Block:A), data=Data)
     summary(lmer.m1)
     anova(lmer.m1)

     lmer.m2<- lmer(Y2~A*B+(1|Block)+(1|Block:A), data=Data)
     summary(lmer.m2)
     anova(lmer.m2)

I'd like to know how to fit multivariate linear mixed model with `lme4`. The
data is below:

     Block A B    Y1    Y2
      1 1 1 135.8 121.6
      1 1 2 149.4 142.5
      1 1 3 155.4 145.0
      1 2 1 105.9 106.6
      1 2 2 112.9 119.2
      1 2 3 121.6 126.7
      2 1 1 121.9 133.5
      2 1 2 136.5 146.1
      2 1 3 145.8 154.0
      2 2 1 102.1 116.0
      2 2 2 112.0 121.3
      2 2 3 114.6 137.3
      3 1 1 133.4 132.4
      3 1 2 139.1 141.8
      3 1 3 157.3 156.1
      3 2 1 101.2  89.0
      3 2 2 109.8 104.6
      3 2 3 111.0 107.7
      4 1 1 124.9 133.4
      4 1 2 140.3 147.7
      4 1 3 147.1 157.7
      4 2 1 110.5  99.1
      4 2 2 117.7 100.9
      4 2 3 129.5 116.2

Thank in advance for your time and cooperation.

Regards,
Iftikhar Ahmad

	[[alternative HTML version deleted]]