require("spdep")
# directional spatial ARCH process (W is triangular, 1:1 origin)
rho <- 0.5
alpha <- 1
d <- 5
n <- d^2
nblist <- cell2nb(d, d, type = "queen")
W <- nb2mat(nblist)
W[lower.tri(W)] <- 0
y <- sim.spARCH(n = n, rho = rho, alpha = alpha, W = W, type = "spARCH")
out <- qml.spARCH(y ~ 0, W = W, type = "spARCH")
summary(out)
Run the code above in your browser using DataLab