Movimento Browniano Fracionário

Plote de um movimento Browniano Fracionário

par(mar = c(4, 4, 1.5, 1.5), mfrow=c(2,2))
#install.packages("fractionalBM")
library(fractionalBM)

H <- c(0.2, 0.5, 0.8, 0.9)
N <- rep(1000,4)

Xt <- mapply(FBM, N, H)

sapply(1:4, function(j) {plot( Xt[,j], type="l", main=paste("Hurts ", H[j]), xlab = "Obs", ylab = "Xt" )} )