with one-way analyses of vari-
ance for the effects of temperature and treatment. The residuals
were checked for normal distribution graphically as well as with
the Shapiro-Wilk test. In the case of non-normally distributed
residuals, the data were log-transformed, which was done for the
analysis of the treatment effect on cumulative CO2 emission at
25 C. In the case of non-homogeneity of variances, Welch analyses
of variance (oneway.test, R Development Core Team, 2014) were
used followed by pairwise t-tests with Bonferroni corrections
without pooled standard deviations and without the assumption