Simulate a 4-emotion model

This function simulates a 4-emotion model which is nonlinear, bistable, discrete, and (almost) centered to zero. Adapted from the model described by van de Leemput et al. (2014).

Usage

sim_4_emo(time = 200, init = c(1.36, 1.36, 4.89, 4.89), sd = 1)

Arguments

time

The number of time steps to simulate.

init

A vector of initial values for the four variables. Default is c(1.36, 1.36, 4.89, 4.89), which is one of the stable states of the model.

sd

The standard deviation of the noise.

Value

A matrix with the simulated data.

References

van de Leemput, I. A., Wichers, M., Cramer, A. O., Borsboom, D., Tuerlinckx, F., Kuppens, P., ... & Scheffer, M. (2014). Critical slowing down as early warning for the onset and termination of depression. Proceedings of the National Academy of Sciences, 111(1), 87-92.

See also

true_model_4_emo(), compare_4_emo(), quadVAR()