A spatial phase transition between neutral and niche communities: R code and model output

R source code for the determination of critical points in the phase transition, the analysis of richness and Shannon diversity of the model, early warnings and the application to Barro Colorado Island fores plot (BCI) . <br><br>The neutral-niche model C++ program is available at: <br>https://dx.doi.org/10.6084/m9.figshare.969692<br><br><br><b>R Markdown files:</b><br>I usually run Rmd files chunk by chunk because there are long simulations, and save ".RData" to keep the results. Neutral_steadyState.Rmd: Simulations to determine the time to reach steady state<br>Neutral_simulCT.Rmd: Simulations to find critical values, and plots H vs Rho vs Migration<br>Neutral_EarlyWarnings_Smax.Rmd: Early warnings of the critical transition based on Maximun Patch <br>BCIpatchAnalysis.Rmd: Barro Colorado Island analysis you have to obtain the BCI data from http://ctfs.si.edu/webatlas/datasets/bci/Terms%20and%20conditions.htm<br><b>R Files</b><br>Neutral_fun.r: functions used for the analysis and simulations<br>BCI_fun.r: functions used for the analysis of BCI data<br><i>powerlaw.zip: code to fit powerlaws from http://tuvalu.santafe.edu/~aaronc/powerlaws/ (by Cosma Shalizi)</i><i><br></i><i><b>Data files</b></i><i><br>.RData</i> are the results of simulations in data.frame format<br><i><br></i><i>Simul.zip</i> Output of all simulations.<br><br> <br>