S3_Barnacle_code.R is R code to run 43 competing models to find the best-supported model for sources of variation in recruitment as presented in Table A2. Models were hierarchical, that is models lacking a particular main effect would not include interactions involving that effect. Code uses lmer() in the lme4 package and then model.sel() in the MuMIn package to obtain corrected AIC, deltas and AIC weights. Code provided is only for log transformed recruitment for barnacles. To run rockweed and mussel recruitment, replace logB_std with logF_std and logM_std, respectively. Code uses data file S3_Recruit_1997-2008_R_ready.csv.

S3_Size_code.R is R code to run competing models that look for thresholds in recruitment as presented in Table A3. Code uses lmer() in the lme4 package and then model.sel() in the MuMIn package to obtain corrected AIC, deltas and AIC weights. Four of the models tested for breaks or thresholds between different sized clearings. For example, the threshold in the first model was set between control plots and clearings of all sizes while the threshold in the fourth model was set between the eight meter clearings and all clearings that were smaller than eight meters. The fifth model tested for no effect due to clearing size (i.e., size was not included as a factor). For rockweeds, we added a sixth model in which tested for an exponential decline because Dudgeon and Petraitis (2001) reported this pattern for rockweeds. All models included time, bays, sites within bays, and interactions as random effects. Code uses data file S3_Recruit_1997-2008_R_ready.csv.

S3_CAP.R is R code to run the constrained analysis of proximities and create Fig. A1. Analysis was run using capscale() in the vegan package. Correlations of the recruitment with the CAP axes were determined using envfit(). Code uses data file S3_Recruit_1997-2008_R_ready.csv.

S3_Recruit_1997-2008_R_ready.csv contains data used in the analyses presented in the current paper and columns with intermediate steps and information about how the data were standardized. The standardizations are number per 1 sq cm per 56 days for A. nodosum, number per 100 sq cm per 56 days for S. balaniodes and number per 40 sq cm per 84 days for M. edulis. Data are in a format that can be easily loaded into R with variables being corrected identified by default as categorical or numeric. R convention of NA for missing data is used. Duration could not be determined for some entries because of missing information. In these cases, duration was inferred either from data from other plots at the same site in the same year or assumed to be the average duration for that year. Note that entry for Y2004-BO-S4-Barnacle is incorrect in Recruitment_data_97-07.txt and has been changed from 308 to 309.