library(plot3D)
## specify grid of spatial and temporal lags
d <- 31
st <- expand.grid(0:(d - 1),
0:(d - 1))
## get the STCS
g14 <- stcfgneiting14(t = st[, 1],
s = st[, 2],
a = 1/50,
c = 1/10,
alpha = 1,
beta = 1,
gamma = 0.5,
tau = 1)
## visualize the STCS
g14.m <- matrix(g14,
nrow = d)
persp3D(z = g14.m, x = 1: nrow(g14.m), y = 1:ncol(g14.m),
expand = 1, main = "", scale = TRUE, facets = TRUE,
xlab="Time lag", ylab = "Distance", zlab = "STCF",
colkey = list(side = 4, length = 0.5), phi = 20, theta = 120,
resfac = 5, col= gg2.col(100))
Run the code above in your browser using DataLab