philrob: Solve Phillip and Robinson's model

In PabRod/sleepR: Sleep-wake Cycle Models

Description

Solves the Phillip and Robinson's model for the given times, initial condition and parameters

Usage

philrob(ts, y0, parms = philrob_default_parms(), tStabil = 0, ...)

Arguments

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

Value

Results of the simulation, including times, states and asleep/awake status

Author(s)

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

References

Phillips AJK, Robinson PA. A Quantitative Model of Sleep-Wake Dynamics Based on the Physiology of the Brainstem Ascending Arousal System. J Biol Rhythms. 2007 Apr 29;22(2):167–79. Available from: http://journals.sagepub.com/doi/10.1177/0748730406297512

See Also

dPhilrob, philrob_default_parms, ode

Examples

y0 <- c(Vv = -13, Vm = 1, H = 10)
nDays <- 30
ts <- seq(0, nDays*24, length.out = nDays*24*20)
sol <- philrob(ts, y0)

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