Code and Data accompanying the paper:

Wenfeng Feng & Kazuhiro Takemoto, Heterogeneity in ecological mutualistic networks dominantly determines community stability, 2014.

License:

The code is made available under the GNU GPL license: http://www.gnu.org/licenses/gpl.html

Requests:

Please feel free to contact Wenfeng Feng if you have any question

Citation:

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.

Structure:

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.

Usage (in R):

> 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.