kronauer: Solve Kronauer's 1990 model

In PabRod/sleepR: Sleep-wake Cycle Models

Description

Solves the Kronauer's 1990 model (equations 2 and 3 in the reference) for the given times, initial condition and parameters

Usage

kronauer(ts, y0, parms = kronauer_default_parms(), ...)

Arguments

ts	Vector of times (in h)
y0	Initial condition
parms	Model parameters (optional, see kronauer_default_parms)
...	Additional arguments passed to the ode integrator

Value

Results of the simulation, including times and states

Author(s)

Pablo Rodríguez-Sánchez (https://pabrod.github.io)

References

Jewett ME, Kronauer RE. Refinement of Limit Cycle Oscillator Model of the Effects of Light on the Human Circadian Pacemaker. J Theor Biol. 1998 Jun 21;192(4):455–65. https://www.sciencedirect.com/science/article/pii/S0022519398906671

See Also

dKronauer, kronauer_default_parms, ode

Examples

y0 <- c(x = 1, xc = 0)
nDays <- 5
ts <- seq(0, nDays*24, length.out = nDays*24*20)
I <- function(t) { 9500 + 0*t }
sol <- kronauer(ts, y0, I)

PabRod/sleepR documentation built on Dec. 21, 2020, 3:27 p.m.