Appendix B. Three tables and a figure showing the full results of the analysis of the Linum experiment and a figure showing the results of the analysis of the Lymnaea data set.
|
| FIG. B1. Results of a traditional analysis of the Linum experiment using the log response ratio (Hedges et al. 1999). Bars in the figure are log response ratios for each of the five fertilizer treatments (gm-2 NPK) with error bars are of standard deviation. There were significant fertilization effects on the competitive response of Linum when the log response ratio data were analyzed using ANOVA (SPSS for Windows v.9) with fertilization as a fixed factor and block as a random factor (F4,27 = 3.345, P = 0.024). A subsequent Tukey test only identified significant differences between the 0 and 8 gm-2 fertilization treatments. The conclusions from this analysis of the response ratios analysis are very different from the power curve analysis. This figure suggests that competitive ability in the zero fertilizer treatment was higher than in the 16 and 32 gm-2 while competitive ability was clearly highest in the high-fertilizer treatments in the power curve analysis (Main text Fig. 1). |
TABLE B1. Summaries of the models generated during the analysis of the Linum data set. Global models have a single curve for the entire data set, while separate models have a curve for each treatment. The number of parameters is the number of coefficient values that were estimated from the data for that model. The residual sums of squares and degrees of freedom values are the sum of the residuals from the individual curves for each treatment.
Model |
No. parameters |
Residual Sums of Squares |
Residual df |
AIC |
Global |
2 |
0.02152 |
37 |
-285.9 |
Separate |
10 |
0.01336 |
29 |
-279.4 |
Global, k3=1 |
1 |
0.02922 |
38 |
-276.3 |
Separate, k3=1 |
5 |
0.02180 |
34 |
-277.4 |
Separate, k2 common |
6 |
0.01415 |
33 |
-291.3 |
Separate, k3 common |
6 |
0.01346 |
33 |
-293.3 |
TABLE B2. Results of extra-sums of squares tests to compare models from the Linum experiment. Global models have a single curve for all treatments while separate models have individual curves for each treatment. The first two tests demonstrate that models with a k3 parameter, indicating nonlinearity, are better than linear models. The third test indicates that a model with separate k2 and k3 parameters for each treatment is not significantly better than single model, but the fourth and fifth tests indicate that models with one parameter in common but the other different between treatments are significantly better than the global model. Since one of the separate models with one common parameter is not a nested subset of the other model, they cannot be compared with the F test, and instead should be compared using AIC (Table B1).
Comparison |
SD difference |
df difference |
F ratio |
P |
Global vs. Global, k3=1 |
0.00770 |
1 |
13.24 |
<0.001 |
Separate, k3 common vs. Separate, k3=1 |
0.00834 |
1 |
20.45 |
<0.001 |
Global vs. Separate |
0.00816 |
8 |
2.21 |
0.056 |
Global vs. Separate, k2 common |
0.00734 |
4 |
4.23 |
0.006 |
Global vs. Separate, k3 common |
0.00806 |
4 |
4.94 |
0.003 |
TABLE B3. Pairwise comparisons of models from the Linum experiment. Each test compared a global curve for two fertilization treatments with the k3 parameter value set to 0.2531 (the k3 parameter value from the individual model with a common k3) to individual models with a separate k2 parameter for each treatment. With Bonferroni correction only P values less than 0.005 should be considered significant.
Comparison |
SS difference |
df difference |
F ratio |
P |
0 vs. 4 |
0.00007 |
1 |
0.422 |
0.527 |
0 vs. 8 |
0.00051 |
1 |
3.576 |
0.081 |
0 vs. 16 |
0.00013 |
1 |
0.475 |
0.503 |
0 vs. 32 |
0.00200 |
1 |
2.767 |
0.120 |
4 vs. 8 |
0.00027 |
1 |
2.208 |
0.159 |
4 vs. 16 |
0.00054 |
1 |
2.259 |
0.155 |
4 vs. 32 |
0.00407 |
1 |
6.169 |
0.026 |
8 vs. 16 |
0.00176 |
1 |
7.980 |
0.013 |
8 vs. 32 |
0.00757 |
1 |
11.80 |
0.004 |
16 vs. 32 |
0.00167 |
1 |
2.205 |
0.160 |
|
| FIG. B2. Power curve for Lymnaea fecundity with raw data points for the two species and drying treatments overlain. The gray line is the 1:1 reference line. Any position above that line indicates higher fecundity at higher density, while any position below the line indicates higher fecundity at lower density. Only one line is shown on this figure as these data were best described by a single curve with k3=1. |
LITERATURE CITED
Hedges, L.V., J. Gurevitch, and P.S. Curtis. 1999. The meta-analysis of response ratios in experimental ecology. Ecology 80:11501156.