Problems with an anova mehtod to test for differences between groups in a lme ()

Dear all,

 I have problems using anova mehtod to test for differences between groups in
 a lme () model. As described at page 224 and 225 in Mixed-Effects Models in
 S and S-PLUS, by Pinhiro and Bates. Could it that this script doesn't work
 for R?, see below (in yellow).

 Respons: Speed
 Fixed effects: Fuel, CarMod (1,2&3)
 Random effects: Place

 > lme1 <- lme(Speed ~  Fuel*Car, random=~1|Place,data=test)
 > summary (lme1)

 Linear mixed-effects model fit by REML
 Data: test
      AIC      BIC    logLik
 261.2013 275.6996 -121.6007

 Random effects:
 Formula: ~1 | Place
           (Intercept)         Residual
 StdDev: 0.0003238738   5.013858

 Fixed effects: Speed ~ Fuel + Car + Driver + Fuel * Car
                      Value        Std.Error     DF   t-value       p-value
 (Intercept)        -29.33479  12.743084   30   -2.302017   0.0285
 Fuel                10.04684    1.408789    30    7.131542   0.0000
 CarMod2         46.55593    14.192029   7     3.280428   0.0135
 CarMod3          1.65157     18.247158   7     0.090511   0.9304
 Fuel:CarMod2  -5.53264     1.624159    30   -3.406464   0.0019
 Fuel:CarMod3  -0.18452      2.010470   30   -0.091779   0.9275

 Number of Observations: 44
 Number of Groups: 10

 > anova( lme1, L = c(Fuel:CarMod2 = 1, Fuel:CarMod3 = -1) )

 Error: syntax error in "anova( lme1, L = c(Fuel:CarMod2="

 Regards R.S

        [[alternative HTML version deleted]]