S4_Recruit_cusp.R is R code to run seven competing models to fit percent cover by rockweeds, mussels, and barnacles as a cusp catastrophe. Models used different combinations of recruitment rates are used as parameters. The package allows the parameters and the ecosystem state to be defined as linear functions of several variables, and finds the coefficients of these functions that provide the best fit.
S4_Cusp_data.csv is the data file needed to run R script. Ecosystem state data are from Petraitis et al. (2008; Ecological Archives E089-032) and include understory cover by barnacles and mussels and canopy cover by Ascophyllum and Fucus. Recruitment data are from data file S3_Recruit_1997-2008_R_ready.csv in Supplement 3. Note file has 59 lines of data; site MD for S2 was dropped because of missing data for cover.
