if (FALSE) {
data(cola)
fit_gmwmx_1 = estimate_gmwmx(x = cola,
theta_0 = c(0.1,0.1,0.1,0.1),
n_seasonal = 1,
model_string = "wn+matern")
fit_gmwmx_2 = estimate_gmwmx(x = cola,
theta_0 = c(0.1,0.1,0.1),
n_seasonal = 1,
model_string = "wn+powerlaw")
compare_fits(fit_gmwmx_1, fit_gmwmx_2)
}
Run the code above in your browser using DataLab