Methodology from Smyth, A. R., Piehler, M. F. and Grabowski, J. H. (2015), Habitat context influences nitrogen removal by restored oyster reefs. J Appl Ecol, 52: 716–725. doi:10.1111/1365-2664.12435

Statistical analyses were performed using r 2.13.1 (R Foundation for Statistical Computing 2011). Linear mixed-effects models (lme in R nlme package), where habitat nested in sampling location was included as a random effect for the intercept, were used to investigate the effects of oyster reef presence, habitat context, nitrate concentration (ambient vs. elevated) and the interaction between these factors on response variables. Fluxes of N2, NOx ( [math formula] + [math formula] ) [math formula] , denitrification efficiency and SOD were analysed using all three fixed effects. For sediment organic matter, only habitat context and reef presence were included as fixed effects. The effects of ambient vs. elevated nitrate concentration and habitat context on oyster reef-mediated changes in denitrification were also analysed with a mixed-effects model (fixed effects: nitrate concentration × habitat context; random effects: habitat nested in location). Relationships between oyster density and habitat context were made using a mixed-effects model (fixed effects: habitat context; random effects: habitat nested in location). Comparisons were conducted using linear contrasts and judged against an alpha level of 0·05. Interactions were assessed using Tukey's HSD (lsmeans in R lsmeans package). Assumptions of homogeneity were tested using Levene's tests. Regression analyses were used to investigate the effect of oyster density on denitrification. Models with the lowest Akaike's information criterion corrected for small sample sizes (AICc) were chosen.

BCO-DMO Processing Notes:

- column names reformatted to comply with BCO-DMO naming standards.
- lat and lon columns added to correspond with locations.