# \dontshow{
x <- borehole$x
y <- borehole$y
theta <- c(
5.767699e+01, 0.000000e+00, 0.000000e+00, 1.433571e-06,
0.000000e+00, 2.366557e-06, 1.695619e-07, 2.454376e-09
)
alpha <- c(
1.110223e-16, 0.000000e+00, 0.000000e+00, 0.000000e+00,
0.000000e+00, 0.000000e+00, 2.494862e-03, 0.000000e+00
)
cor_par <- data.frame(Theta = theta, Alpha = alpha)
rownames(cor_par) <- colnames(borehole$x)
sp_var <- 38783.7
borehole_fit <- GaSPModel(
x = borehole$x, y = borehole$y,
reg_model = ~1, cor_family = "PowerExponential",
cor_par = cor_par, random_error = FALSE,
sp_var = sp_var
)
# }
borehole_cv <- CrossValidate(borehole_fit)
Run the code above in your browser using DataLab