active_forcing: Determine when forcing is active for an object of class...

In RobertGM111/havok: Procedures for the Analysis of Nonlinear Dynamical Systems

Determine when forcing is active for an object of class havok.

Description

This function uses a threshold to determine when forcing is active in a fitted havok object.

Usage

active_forcing(x, thresh = stats::sd(x$Vr[, x$r]))

Arguments

x	An object of class havok.
thresh	Numeric; A cutoff value for determining when forcing is active. Defaults to one standard deviation of the forcing term.

Value

A list of forcing values with their corresponding activity status.

Examples

data(ECG_measurements)

xdat <- ECG_measurements[,"channel1"]
dt <- ECG_measurements[2,"time"] - ECG_measurements[1,"time"]

stackmax <- 25
rmax <- 5
lambda <- .001

hav <- havok(xdat = xdat, dt = dt, stackmax = stackmax, lambda = lambda, rmax = rmax)
active_forcing(x = hav)

RobertGM111/havok documentation built on July 8, 2023, 8:23 p.m.