model { for (i in 1:N) { Y[i] ~ dnorm(mu[i], error.prec) mu[i] <- beta[1] + beta[2]*X[Site[i]] + V[Time[i]] + W[Site[i]] } error.prec ~ dgamma(0.10,0.10) for (i in 1:p) {beta[i] ~ dnorm(0.0, 0.0001)} tausq <- 1/error.prec for (i in 1:NSITES) {muW[i] <- 0.0} W[1:NSITES] ~ spatial.exp(muW[], x[], y[], spat.prec, phi, 1) phi ~ dunif(0,10) for (i in 1:NTIME) {muV[i] <- 0.0} V[1:NTIME] ~ spatial.exp(muV[],x.t[],y.t[],temporal.prec,phi.t,1) phi.t ~ dunif(0,10) spat.prec ~ dgamma(0.10, 0.10) sigmasq <- 1/spat.prec temporal.prec ~ dgamma(0.10,0.10) sigmasq.temp <- 1/temporal.prec } #Initial Values list(beta=c(0,0),phi=1.0,phi.t=1.0,spat.prec=1.0,error.prec=1.0,temporal.prec=1.0,W=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),V=c(0,0,0,0,0,0,0,0,0,0,0,0)) #Data list(N=600,p=2,NSITES=50,NTIME=12,Y = c(0.07, 0.25, 0.02, -0.02, -0.22, 0.79, 0.18, 0.53, -0.49, 0.48, 0.73, -0.28, 0.13, -0.06, 0.52, 0.11, -0.13, -0.1, 0.65, -0.73, 0.08, 0.45, -0.16, -0.02, -0.44, 0.45, 0.41, 0.26, 0.88, 0.22, 0.49, 0.44, -0.04, 0.05, 0.34, 0.23, 0.1, -0.25, -0.38, 0.14, 0.62, -0.06, 0.56, 0.43, 0.91, -0.11, -0.26, 0.24, 0.06, -0.27, 0.22, 0.51, 0.36, 0.25, 0.07, 0.8, 0.37, 0.74, -0.33, 0.62, 0.88, 0, 0.21, 0.21, 0.82, 0.34, 0.08, -0.09, 0.86, -0.66, 0.36, 0.77, -0.06, 0.21, -0.27, 0.56, 0.78, 0.3, 0.94, 0.27, 1.08, 0.64, 0.02, 0.09, 0.57, 0.41, 0.47, -0.07, -0.17, 0.31, 0.71, 0.06, 0.57, 0.52, 0.99, 0.15, 0.04, 0.6, 0.37, -0.01, 1.08, 1.15, 1.2, 1.1, 0.68, 1.85, 1.23, 1.61, 0.27, 1.55, 1.84, 0.76, 0.95, 0.81, 1.61, 1.16, 0.71, 0.86, 1.69, 0.07, 1.01, 1.57, 0.68, 0.92, 0.33, 1.2, 1.53, 1.3, 1.85, 1.06, 1.95, 1.5, 0.78, 0.69, 1.44, 1.15, 1.32, 0.58, 0.62, 1.12, 1.47, 0.98, 1.32, 1.43, 1.91, 1.13, 0.78, 1.62, 1.17, 0.73, 1.14, 1, 1.01, 1.18, 0.67, 1.63, 1.04, 1.54, 0.18, 1.42, 1.49, 0.64, 0.84, 0.75, 1.58, 1.28, 0.71, 1.07, 1.56, 0.22, 0.83, 1.55, 0.73, 0.58, 0.53, 0.95, 1.39, 1.13, 1.58, 0.88, 1.85, 1.47, 0.75, 0.53, 1.33, 0.96, 1.28, 0.47, 0.54, 1.02, 1.26, 0.97, 1.03, 1.12, 1.59, 1.27, 0.88, 1.4, 1.16, 0.63, 2.08, 1.75, 1.71, 2.14, 1.43, 2.47, 1.85, 2.31, 0.89, 2.25, 2.47, 1.38, 1.55, 1.55, 2.52, 2.16, 1.47, 1.89, 2.4, 1.08, 1.47, 2.53, 1.61, 1.22, 1.44, 1.79, 2.18, 1.86, 2.46, 1.8, 2.76, 2.26, 1.64, 1.29, 2.25, 1.88, 2.1, 1.21, 1.19, 1.76, 2.23, 1.99, 1.88, 2, 2.52, 2.27, 1.68, 2.34, 2.05, 1.34, 2.53, 2.24, 2.27, 2.71, 2.01, 3, 2.21, 2.93, 1.52, 2.85, 2.98, 2.11, 2.17, 2.14, 3.07, 2.67, 2.07, 2.5, 2.92, 1.62, 1.97, 2.98, 2.27, 1.76, 2.06, 2.4, 2.72, 2.41, 2.92, 2.29, 3.24, 2.92, 2.24, 1.95, 2.71, 2.4, 2.69, 1.71, 1.92, 2.39, 2.78, 2.55, 2.46, 2.59, 3.05, 2.86, 2.3, 2.86, 2.58, 2.06, 2.83, 2.55, 2.56, 2.96, 2.29, 3.36, 2.49, 3.34, 1.79, 3.21, 3.48, 2.31, 2.41, 2.36, 3.38, 2.94, 2.34, 2.61, 3.33, 1.88, 2.24, 3.24, 2.56, 1.97, 2.32, 3.06, 3.01, 2.83, 3.37, 2.61, 3.44, 3.35, 2.41, 2.27, 2.93, 2.91, 2.97, 1.96, 2.21, 2.64, 3.08, 2.91, 2.75, 3.01, 3.46, 3.11, 2.7, 3.24, 2.89, 2.23, 2.61, 2.39, 2.37, 2.84, 2.13, 3.13, 2.38, 3.04, 1.61, 2.86, 3.06, 2.11, 2.28, 2.21, 3.16, 2.76, 2.22, 2.44, 2.98, 1.78, 2.01, 3, 2.38, 1.91, 2.16, 2.67, 2.87, 2.54, 3.08, 2.49, 3.29, 3.01, 2.29, 2.02, 2.82, 2.56, 2.76, 1.79, 1.99, 2.47, 2.81, 2.75, 2.54, 2.64, 3.04, 2.98, 2.47, 2.93, 2.78, 2.04, 2.27, 2.18, 2.15, 2.34, 1.73, 2.89, 2.16, 2.81, 1.41, 2.64, 2.81, 1.86, 2.07, 1.83, 2.67, 2.41, 1.79, 2.11, 2.68, 1.31, 1.77, 2.74, 2.02, 1.75, 1.86, 2.19, 2.59, 2.37, 2.83, 2.26, 2.99, 2.77, 2.09, 1.83, 2.4, 2.26, 2.46, 1.44, 1.79, 2.26, 2.5, 2.24, 2.31, 2.44, 2.88, 2.43, 2.23, 2.68, 2.27, 1.78, 1.53, 1.54, 1.37, 1.64, 1.03, 2.05, 1.52, 2.07, 0.62, 1.79, 1.97, 1.07, 1.05, 1.15, 1.92, 1.65, 1.11, 1.57, 1.97, 0.45, 1.07, 1.85, 1.28, 1.07, 1.13, 1.47, 1.9, 1.62, 2.19, 1.41, 2.23, 1.93, 1.33, 1.16, 1.8, 1.62, 1.77, 0.82, 0.99, 1.47, 1.88, 1.52, 1.76, 1.63, 1.84, 1.73, 1.5, 1.87, 1.61, 0.98, 0.78, 0.82, 0.66, 0.68, 0.19, 1.03, 0.77, 1.14, -0.23, 0.87, 1.06, 0.17, 0.3, 0.27, 1.03, 0.89, 0.22, 0.52, 0.98, -0.41, 0.43, 1.16, 0.34, 0.3, 0.22, 0.56, 1.07, 0.73, 1.31, 0.55, 1.46, 0.98, 0.44, 0.29, 0.91, 0.7, 0.78, -0.04, 0.02, 0.55, 0.98, 0.54, 0.88, 0.68, 0.95, 0.72, 0.54, 0.89, 0.79, 0.1, 0.14, 0.38, 0.18, 0.04, -0.04, 0.51, 0.34, 0.53, -0.51, 0.25, 0.29, -0.28, -0.08, -0.09, 0.44, 0.34, -0.04, -0.08, 0.54, -0.71, 0.08, 0.6, -0.02, -0.08, -0.42, 0.25, 0.42, 0.23, 0.63, 0.16, 0.73, 0.38, 0.03, 0.02, 0.27, 0.38, 0.16, -0.33, -0.26, -0.15, 0.59, -0.06, 0.55, 0.19, 0.32, -0.07, -0.06, 0.51, 0.23, -0.26 ), X= c(2.156, 2.464, 2.561, 1.806, 2.652, 1.552, 2.332, 0.903, 3.384, 1.115, 1.214, 2.707, 3.003, 2.561, 1.409, 1.976, 2.485, 2.379, 1.498, 3.024, 3.232, 1.841, 2.471, 3.116, 2.705, 2.158, 2.074, 2.396, 1.326, 2.583, 1.51, 1.028, 2.377, 2.901, 1.884, 1.833, 1.982, 3.152, 3.098, 2.475, 1.611, 1.903, 2.135, 1.695, 1.57, 1.61, 2.318, 1.418, 1.903, 2.866), x=c(-109.33, -107.65, -106.657, -108.08, -106.969, -104.02, -107.58, -101.53, -105.87, -101.7, -102.28, -106.67, -105.341, -106.781, -108.58, -109.03, -107.167, -107.17, -104.65, -107.9667, -105.25, -109.48, -106.27, -105.122, -106.97, -104.87, -106.73, -105.48, -102.97, -105.28, -108.73, -101.38, -107.47, -105.18, -107.58, -104.65, -106.83, -107.294, -106.953, -105.53, -105, -107.65, -105.32, -103.7, -103.32, -108.8, -106.13, -104.7, -107.9, -106.158), y=c(37.87, 37.38, 36.956, 40.52, 40.848, 37.63, 37.37, 40.05, 39.646, 39.37, 37.4, 40.08, 36.699, 40.534, 39.07, 37.82, 41.117, 38.47, 38.18, 41.15, 37, 37.62, 40.73, 37.214, 38.87, 39.1, 36.6, 38.13, 36.9, 36.62, 36.77, 39.78, 40.03, 38.75, 38.85, 41.15, 39.67, 40.108, 38.894, 37.72, 39.75, 41.03, 39.63, 39.18, 37.45, 40.08, 40.05, 40.42, 40.03, 39.718), Time=c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12 ), Site=c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50), x.t=c(1,2,3,4,5,6,7,8,9,10,11,12),y.t=c(0,0,0,0,0,0,0,0,0,0,0,0))