> summary(fit_all)

Call:
lm(formula = ec ~ bc:year:season:site - 1, data = alldata)

Residuals:
    Min      1Q  Median      3Q     Max
-6.3450 -0.2091  0.0510  0.4306  5.8881

Coefficients: (2 not defined because of singularities)
                                     Estimate Std. Error t value Pr(>|t|)
bc:year2009:seasonA:siteMarylebone   0.757541   0.009116  83.098  < 2e-16 ***
bc:year2010:seasonA:siteMarylebone   0.656476   0.007400  88.710  < 2e-16 ***
bc:year2011:seasonA:siteMarylebone   0.689421   0.007350  93.798  < 2e-16 ***
bc:year2009:seasonSp:siteMarylebone  0.864448   0.008997  96.081  < 2e-16 ***
bc:year2010:seasonSp:siteMarylebone  0.760748   0.010795  70.472  < 2e-16 ***
bc:year2011:seasonSp:siteMarylebone  0.695510   0.009142  76.082  < 2e-16 ***
bc:year2009:seasonSu:siteMarylebone  0.710539   0.008897  79.863  < 2e-16 ***
bc:year2010:seasonSu:siteMarylebone  0.661160   0.008737  75.674  < 2e-16 ***
bc:year2011:seasonSu:siteMarylebone  0.575097   0.007643  75.244  < 2e-16 ***
bc:year2009:seasonW:siteMarylebone   0.637619   0.013366  47.705  < 2e-16 ***
bc:year2010:seasonW:siteMarylebone   0.715393   0.010173  70.323  < 2e-16 ***
bc:year2011:seasonW:siteMarylebone   0.703555   0.008607  81.738  < 2e-16 ***
bc:year2009:seasonA:siteNKensington  0.799219   0.040019  19.971  < 2e-16 ***
bc:year2010:seasonA:siteNKensington  0.994461   0.040281  24.688  < 2e-16 ***
bc:year2011:seasonA:siteNKensington  0.974313   0.037304  26.118  < 2e-16 ***
bc:year2009:seasonSp:siteNKensington 0.955909   0.045654  20.938  < 2e-16 ***
bc:year2010:seasonSp:siteNKensington 1.099792   0.068394  16.080  < 2e-16 ***
bc:year2011:seasonSp:siteNKensington 1.174542   0.063348  18.541  < 2e-16 ***
bc:year2009:seasonSu:siteNKensington 0.700316   0.064641  10.834  < 2e-16 ***
bc:year2010:seasonSu:siteNKensington 0.988856   0.091529  10.804  < 2e-16 ***
bc:year2011:seasonSu:siteNKensington 1.071967   0.102375  10.471  < 2e-16 ***
bc:year2009:seasonW:siteNKensington  0.890885   0.035577  25.041  < 2e-16 ***
bc:year2010:seasonW:siteNKensington  0.923421   0.045026  20.508  < 2e-16 ***
bc:year2011:seasonW:siteNKensington  1.008563   0.062412  16.160  < 2e-16 ***
bc:year2009:seasonA:siteHarwell      0.435643   0.702707   0.620  0.53536
bc:year2010:seasonA:siteHarwell      0.823244   0.138353   5.950 3.11e-09 ***
bc:year2011:seasonA:siteHarwell      0.844239   0.090571   9.321  < 2e-16 ***
bc:year2009:seasonSp:siteHarwell           NA         NA      NA       NA
bc:year2010:seasonSp:siteHarwell     0.786092   0.194016   4.052 5.26e-05 ***
bc:year2011:seasonSp:siteHarwell     0.911042   0.166637   5.467 5.09e-08 ***
bc:year2009:seasonSu:siteHarwell           NA         NA      NA       NA
bc:year2010:seasonSu:siteHarwell     0.830268   0.321735   2.581  0.00993 **
bc:year2011:seasonSu:siteHarwell     0.663322   0.243474   2.724  0.00649 **
bc:year2009:seasonW:siteHarwell      1.853379   0.433678   4.274 2.01e-05 ***
bc:year2010:seasonW:siteHarwell      1.248616   0.139352   8.960  < 2e-16 ***
bc:year2011:seasonW:siteHarwell      0.991226   0.157439   6.296 3.68e-10 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.8445 on 2178 degrees of freedom
  (1073 observations deleted due to missingness)
Multiple R-squared:  0.9735,	Adjusted R-squared:  0.9731
F-statistic:  2351 on 34 and 2178 DF,  p-value: < 2.2e-16