set.seed(123)
N <- 3
nobs <- 20
freqs <- (1:nobs) / nobs
spec <- matrix(NA, nrow = nobs, ncol = N)
for (i in 1:N) {
spec[, i] <- exp(2 * cos(2 * pi * freqs) + rnorm(1, sd = 0.1))
}
data_generater(N = N, nobs = nobs, spec = spec)
Run the code above in your browser using DataLab