# NOT RUN {
n <- 1000
d <- 2
times <- seq(n)
delta_time <- 0.01
noise <- cbind(ghyp::rghyp(n, ghyp::ghyp()), ghyp::rghyp(n, ghyp::ghyp()))
data <- construct_path(
diag(d),
noise = noise, y_init = rep(0, d), delta_time = delta_time
)
levy_recovery(adj = diag(d), data = data, times = times)
# }
Run the code above in your browser using DataLab