create_Lmatrix: Make L matrix

In ATNr: Run Allometric Trophic Networks Models

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

ATNr documentation built on June 8, 2025, 10:49 a.m.