slope: Calculate slope estimate for time series observations

In lkshrsch/gLVInterNetworks: Inference of interaction networks based on generalised Lotka Volterra dynamics

Usage

slope(abundances)

Arguments

abundances

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (abundances) 
{
    "Get log slopes of measured time points. For discrete LV linear algebraic system"
    "Input: observations without time points == Abundances"
    slopes = matrix(0, nrow = nrow(abundances) - 1, ncol = ncol(abundances))
    for (column in 1:ncol(abundances)) {
        for (row in 1:(nrow(abundances) - 1)) {
            try(slopes[row, column] <- log_slope_luk(abundances[row, 
                column], abundances[row + 1, column]))
        }
    }
    return(slopes)
  }

lkshrsch/gLVInterNetworks documentation built on May 21, 2019, 7:33 a.m.