dKronauer: Kronauer's 1990 model flow

In PabRod/sleepR: Sleep-wake Cycle Models

Description

Generates the flow (right hand side of the differential equation) for the Kronauer's 1990 model (equations 2 and 3 in the reference)

Usage

dKronauer(time, y, I, parms = kronauer_default_parms())

Arguments

time	The time (in h)
y	The state
I	Illumination function (in lux)
parms	Model's parameters (optional, see kronauer_default_parms)

Value

The flow (right hand side of the differential equation)

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

kronauer, kronauer_default_parms

Examples

## Not run: 
t <- 0
y <- c(x = 1, xc = 1)
dy <- dKronauer(t, y)

## End(Not run)

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