RFoptions(seed=0)
model <- RMnsst(phi=RMgauss(), psi=RMfbm(alpha=1), delta=2)
x <- seq(0, 10, if (interactive()) 0.2 else 1)
plot(model, ylim=c(0,1), dim=2)
plot(RFsimulate(model, x=x, y=x))
RFoptions(seed=NA)
Run the code above in your browser using DataLab