Wenfeng Feng & Kazuhiro Takemoto, Heterogeneity in ecological mutualistic networks dominantly determines community stability, 2014.
The code is made available under the GNU GPL license: http://www.gnu.org/licenses/gpl.html
Please feel free to contact Wenfeng Feng if you have any question
When using the codes in your study, please cite the following paper:
Wenfeng Feng & Kazuhiro Takemoto, Heterogeneity in ecological mutualistic networks dominantly determines community stability, 2014.
All the code and data files are constructed using a statistical software R (http://cran.us.r-project.org/).
The following R scripts require the following packages:
igraph, bipartite, plyr, reshape2, doMC, MASS, ggplot2, akima.
plustwins.functions.R
This is a source code for estimating the dominant eigenvalue of a bipartite network and for evaluating network measures using the null models.
plustwins.test.R
Sample scripts for generating the figures in the paper.
datasets.RData
This file includes 40 empirical mutualistic networks collected from the ecological literatures.
> load('datasets.RData')
> source('plustwins.functions.R')
> source('plustwins.test.R')
See the descriptions in the files for details.
Note: Execution of plustwins.test.R is time-consuming; thus, you may need to run the scripts step by step.