Generate data related to Section 7.2 of Lu et al. (2012)
(code from author). Define n = c(128,
128, 24)
. Then x[[i]] = (1:n[i])/n[i] -
1/2/n[i]
. These are the observed data locations. For
i
, j
, k
spanning the full length
of each element of x
, mu3d[i, j, k] =
x[[1]][i]^2 + x[[2]][j]^2 + x[[3]][k]^2
. Lastly,
data3d = mu3d + 0.5 * rnorm(n[1] * n[2] * n[3])
.