create_Lmatrix {ATNr}

Make L matrix

Description

Make L matrix

Usage

create_Lmatrix(BM, nb_b, Ropt = 100, gamma = 2, th = 0.01)

Arguments

BM

float vector, body mass of species.

nb_b

integer, number of basal species.

Ropt

numeric, consumer/resource optimal body mass ratio.

gamma

numeric, code for the width of the Ricker function.

th

float, the threshold below which attack rates are considered = 0.

Details

The L matrix contains the probability for an attack event to be successful based on allometric rules and a Ricker function defined by Ropt and gamma. If at least one species has not resource or consumer (i.e. it is an isolated species), another food web is generated, until a maximum of 100 iterations.

Value

A numeric matrix with the probability for an attack event between two species to be successful.

Examples

set.seed(123)
mass <- sort(10 ^ runif(30, 2, 6))
L <- create_Lmatrix(mass, nb_b = 10, Ropt = 100)
image(t(L))