> orange.rsm1 <- lm(y2 ~ x1+x2+x3+x11+x22+x33+x12+x13+x23+block1+block2)
> summary(orange.rsm1)

Call:
lm(formula = y2 ~ x1 + x2 + x3 + x11 + x22 + x33 + x12 + x13 + 
    x23 + block1 + block2)

Residuals:
     Min       1Q   Median       3Q      Max 
-13.4608  -4.7507   0.0618   5.8101   9.8099 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)    80.4006   203.7150   0.395 0.703398    
x1             49.5912    10.7830   4.599 0.001758 ** 
x2          -2106.4948   366.3190  -5.750 0.000429 ***
x3            119.0082    20.5177   5.800 0.000405 ***
x11             0.1287     0.2322   0.554 0.594488    
x22           480.4509   294.0095   1.634 0.140872    
x33            -3.3939     0.7106  -4.776 0.001398 ** 
x12            37.5000    10.4608   3.585 0.007138 ** 
x13            -3.7679     0.5230  -7.204 9.21e-05 ***
x23           150.6250    18.3064   8.228 3.57e-05 ***
block1        -88.3570     5.5932 -15.797 2.58e-07 ***
block2         10.6430     5.5932   1.903 0.093560 .  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 10.36 on 8 degrees of freedom
Multiple R-squared: 0.9969,     Adjusted R-squared: 0.9927 
F-statistic: 235.5 on 11 and 8 DF,  p-value: 8.441e-09 

> 
> orange.rsm2 <- lm(y2 ~ x1+x2+x3+x11+x22+x33+x12+x13+x23)
> summary(orange.rsm2)

Call:
lm(formula = y2 ~ x1 + x2 + x3 + x11 + x22 + x33 + x12 + x13 + 
    x23)

Residuals:
   Min     1Q Median     3Q    Max 
-75.77 -55.74  23.58  30.48  41.23 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept)    50.8028  1204.5638   0.042    0.967
x1             49.6119    63.7699   0.778    0.455
x2          -2074.6266  2166.2560  -0.958    0.361
x3            119.0091   121.3401   0.981    0.350
x11             0.1281     1.3732   0.093    0.928
x22           440.6156  1738.5094   0.253    0.805
x33            -3.3939     4.2025  -0.808    0.438
x12            37.5000    61.8642   0.606    0.558
x13            -3.7679     3.0932  -1.218    0.251
x23           150.6250   108.2623   1.391    0.194

Residual standard error: 61.24 on 10 degrees of freedom
Multiple R-squared: 0.8654,     Adjusted R-squared: 0.7442 
F-statistic: 7.143 on 9 and 10 DF,  p-value: 0.002496 

> 
> anova(orange.rsm1,orange.rsm2)
Analysis of Variance Table

Model 1: y2 ~ x1 + x2 + x3 + x11 + x22 + x33 + x12 + x13 + x23 + block1 + 
    block2
Model 2: y2 ~ x1 + x2 + x3 + x11 + x22 + x33 + x12 + x13 + x23
  Res.Df   RSS Df Sum of Sq      F    Pr(>F)    
1      8   858                                  
2     10 37506 -2    -36648 170.87 2.738e-07 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1