RFoptions(seed=0)
#plot shape function of an M3 process of type 'strokorb'
model <- RMbrownresnick(RMfbm(alpha=1.5, s=0.2))
plot(RMstrokorbMono(model))
x <- seq(0, 10, if (interactive()) 0.005 else 1)
z <- RFsimulate(RPsmith(RMstrokorbMono(model), xi=0), x)
plot(z, type="p", pch=20)
RFoptions(seed=NA)
Run the code above in your browser using DataLab