Diversity measure provided in the vegetarian library call bootstrap to estimate standard error for the parameter outputs. Bootstrap generates simulated populations from a multinormal distribution centered around the community matrix provided and applies the function of interest to the simulated data. Requires stats package.

1 |

`abundances` |
Community data as a matrix where columns are individual species and rows are sites. Matrix elements are abundance data (e.g. counts, percent cover estimates). |

`s.sizes` |
A vector of the total abundances at each site. If no vector is provided, then the scores in the abundance matrix will be treated as counts, and sample size of each site will be computed as the row sum. |

`num.iter` |
An integer giving the number of bootstrapping iterations. |

`func` |
The function on which bootstrapping will be performed. The function output must be a scalar. |

`func.arg` |
A list containing any additional arguments that must be fed to the function specified by func. |

`sim.pop` |
A logical specifying whether or not to output the simulated populations for each iteration. |

`sim.par` |
A logical specifying whether or not to output the individual parameters calculated on each simulated population. |

The procedure is outlined in Chao et al. (2008). For each site (row) of the original data set, the total sample size and the relative proportions of each species are used to construct a multinomial distribution. Scores are drawn from this multinomial distribution to generate a simulated population at that site. This is done for each row until an entire simulated community matrix is generated. The supplied function (func) is then applied to this simulated matrix to calculate the focal parameter for that iteration. This process of sampling from the multinormal distributions and applying the function is repeated for as many times as is specified by num.iter. The standard deviation of the bootstrap parameter values serves as the standard error for the parameter estimate based on the real data.

While this function is called internally by other functions in the vegetarian library, only the standard error is reported by those functions. The user may wish to call bootstrap directly to access the simulated populations and measure uncertainty in other ways.

`StdErr` |
Bootstrap Standard Error for parameter estimate |

`Simulated.Populations` |
(optional) Array containing simulated populations |

`Simulated.Parameters` |
(optional) Vector containing parameter for each simulated population |

Noah Charney, Sydne Record

Chao, A, L. Jost, S. C. Chiang, Y.-H Jiang, R. L. Chazdon. 2008. A two-stage probabilistic approach to multiple-community similarity indices. Biometrics 64: 1178-1186.

1 2 3 4 5 6 7 | ```
data(simesants)
bootstrap(simesants[,-1], func=d)
boot_output<-bootstrap(simesants[,-1], func=d, num.iter=10, sim.pop=TRUE, sim.par=TRUE)
boot_output$Simulated.Parameters
boot_output$Simulated.Populations
simesants_freq<-normalize.rows(simesants[,-1])#not a necessary step, but just to show what you would do if your data were in the form of frequencies, not counts
bootstrap(simesants_freq, func=d, num.iter=1000, s.sizes=c(198,186,102,108,187))
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.