P. S. Petraitis and S. R. Dudgeon. 2015. Variation in recruitment and the establishment of alternative community states. Ecology 96:3186–3196. http://dx.doi.org/10.1890/14-2107.1
R script and data files to run correlations and Mantel tests as shown in Table A7 and Fig. A2.
Ecological Archives E096-274-S5.
File list (downloads)
P. S. Petraitis
Department of Biology
University of Pennsylvania
Philadelphia, PA 19104-6018, USA
S. R. Dudgeon
Department of Biology
California State University
Northridge, CA 91330-8303, USA
S5_Correlation_data.csv (MD5: 0e53bc8dd1f1543dcd1f4f089b86eace)
S5_Mantel_Gower.R (MD5: f4bfc63a0f33630769b6dfdd3308e6aa)
S5_ANpreG.csv (MD5: 8934aa17ffbe40a676ef7da32549f55d)
S5_MEpreG.csv (MD5: 8dcdd84ec700adbb272f54569dd2fc8f)
S5_SBpreG.csv (MD5: 36a277d3087c3665586825f9da2b47ae)
S5_ANpostG.csv (MD5: 91b0c6bf006b68107b0e187fda2ca5e9)
S5_MEpostG.csv (MD5: 449d0c283b64a329fe6c62f6caaf731b)
S5_SBpostG.csv (MD5: d91fc2499bf9b8272fca0f585f82465c)
Jmp Pro 10 was used to calculate the correlations between recruitment in 1997 and the 20112012 averages and between the 19972008 averages and the 20112012 averages. Data can be found in S5_Correlation_data.csv. Data for calculation of averages were extracted from files S3_Recruit_1997-2008_R_ready.csv (see Supplement 3) and S6_Recruit_2011-2012_R_ready.csv (see Supplement 6). Missing data are denoted with a period. The columns of S5_Correlation_data.csv are:
S5_Mantel_Gower.R is the R code for the Mantel tests that tested the concordance in recruitment in individual plots before and after scraping. First, the FD package was used to calculate Gower’s dissimilarities, which handle missing data and standardizes by column. Mantel tests were run using the vegan package. We also checked for concordance between species by comparing mussels versus barnacles, mussels versus rockweeds, and barnacles versus rockweeds.
There are six data files; file names include an identifier for species (AN=Ascophyllum, ME=mussels, SB=barnacles), and if the data are from pre or post scraping. For example, data file for mussel recruitment prior to scraping is named S5_MEpreG.csv. The “pre-scraped” data files were extracted from S3_Recruit_1997-2008_R_ready.csv (see Supplement 3) and the “post-scraped data files were extracted from S6_Recruit_2011-2012_R_ready.csv (see Supplement 6). Data files do not include a column with plot ID; header row identifies Year, e.g., Y1997, Y1998, etc.
There are missing data, which is why Gower dissimilarities were used. For the pre-scraped data, some years are missing entirely, and there are a few missing observations in other years. Post data files have data from 2001 and 2012 for 60 plots (n = 120 possible entries).
S5_ANpreG.csv is missing data for 1999 and 2004, thus has a possible 600 entries (10 years × 60 plots). Total missing entries = 16/600.
S5_MEpreG.csv is missing data for 1999 and has 660 possible entries. Total missing entries = 26/660.
S5_SBpreG.csv is missing data for 1998 and 2002, and has 600 possible entries. Total missing entries = 31/600.
S5_ANpostG.csv is missing 3/120 entries.
S5_MEpostG.csv is missing 0/120 entries.
S5_SBpostG.csv is missing = 4/120 entries.