library(plot3D)## specify grid of spatial and temporal lagsd <-31st <- expand.grid(0:(d -1),0:(d -1))## get the STCSg14 <- stcfgneiting14(t = st[,1], s = st[,2], a =1/50, c =1/10, alpha =1, beta =1, gamma =0.5, tau =1)## visualize the STCSg14.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))
References
Gneiting, T. (2002). Nonseparable, Stationary Covariance Functions for Space-Time Data, Journal of the American Statistical Association, 97:458, 590-600, tools:::Rd_expr_doi("10.1198/016214502760047113")