auc_robustness_lbm {robber}

Fast computation of the LBM robustness statistic for uniform extinctions

Description

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.

Usage

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)