lorenz: Lorenz '96 spatPomp simulator

In kidusasfaw/spatPomp: Inference for Spatiotemporal Partially Observed Markov Processes

Lorenz '96 spatPomp simulator

Description

Generate a spatPomp object representing a U-dimensional stochastic Lorenz '96 process with N measurements made at times t_n = n * delta_obs, simulated using an Euler method with time increment delta_t.

Usage

lorenz(
  U = 5,
  N = 100,
  delta_t = 0.01,
  delta_obs = 0.5,
  regular_params = c(F = 8, sigma = 1, tau = 1)
)

Arguments

U	A length-one numeric signifying the number of spatial units for the process.
N	A length-one numeric signifying the number of observations.
delta_t	A length-one numeric giving the Euler time step for the numerical solution.
delta_obs	A length-one numeric giving the time between observations.
regular_params	A named numeric vector containing the values of the F, sigma and tau parameters. F=8 is a common value that causes chaotic behavior.

Value

An object of class ‘spatPomp’ representing a simulation from a U-dimensional Lorenz 96 model

Author(s)

Edward L. Ionides

References

\lorenz

96

\ionides

2021

See Also

Other spatPomp model generators: bm(), bm2(), gbm(), he10(), measles()

Examples

# Complete examples are provided in the package tests
## Not run: 
l <- lorenz(U=5, N=100, delta_t=0.01, delta_obs=1)
# See all the model specifications of the object
spy(l)

## End(Not run)

kidusasfaw/spatPomp documentation built on April 22, 2024, 4:24 p.m.