Fast computation of the LBM robustness statistic for uniform extinctions

Source: R/robustness_bm.R

This function provides much faster computation than robustness_lbm for the case where ext_seq = "uniform", when just the AUC statistic is needed and the robustness function is not needed. This is particularly useful if nr gets large.

auc_robustness_lbm(con, pi, rho, nr, nc)

Arguments

con

A matrix, the connectivity parameter

pi

A vector of length nrow(con), the proportion of row blocks

rho

A vector of length ncol(con), the proportion of column blocks

nr

An integer, the number of row (primary) species

nc

An integer, the number of column (secondary) species

Value

The AUC of the LBM Robustness function for uniform primary extinctions sequences.

Examples

con <- matrix(c(.5,.3,.3,.1), 2, 2)
pi  <- c(.25,.75)
rho <- c(1/3, 2/3)
nr <- 50
nc <- 30
auc_robustness_lbm(con, pi, rho, nr, nc)
#> [1] 0.9120858