find_s: Find the relative amount of self-regulation needed for...

In franzikoch/Feedbackloops: Feedback loop analysis for competition networks

Find the relative amount of self-regulation needed for stability s*

Description

For unstable matrices, diagonal values are incrementally increased by multiplying it with a factor s, until its dominant eigenvalue becomes negative (=the matrix becomes stable).

Usage

find_s(Jacobian, step_size = 0.01, start_s = 1, max_s = 1000)

Arguments

Jacobian	A Jacobian matrix
step_size	Amount by which s is increases in each step
start_s	Starting value of s
max_s	Maximum possible s value. If the matrix is still not stable, NA is returned

Details

This will not work, if the diagonal contains zeros! Replace missing values before calculating s*!

For stable matrices, diagonal values are incrementally decrease, until the dominant eigenvalue becomes positive (= the matrix is unstable)

franzikoch/Feedbackloops documentation built on July 1, 2023, 12:42 p.m.