## Test long memory in time series
library(pracma)
set.seed(123)
data("brown72")
x72 <- brown72 # H = 0.72
xgn <- rnorm(1024) # H = 0.50
xlm <- numeric(1024); xlm[1] <- 0.1 # H = 0.43
for (i in 2:1024) xlm[i] <- 4 * xlm[i-1] * (1 - xlm[i-1])
SLmemory.test(x72)
SLmemory.test(xgn)
SLmemory.test(xlm)
Run the code above in your browser using DataLab