# \donttest{
library(spdep)
library(gamlss)
data(oldcol)
W1 <- spdep::nb2mat(COL.nb, style = "W")
W2 <- W1 # In this case, assume the same spatial weights for both
# Fit a SARARgamlss model
result_sarar <- SARARgamlss(formula = CRIME ~ INC + HOVAL,
sigma.formula = ~ INC + pb(HOVAL),
W1 = W1, W2 = W2, data = COL.OLD,
type="SAR")
summary_SAR(result_sarar)
# Example for GEESAR model
result_geesar <- GEESAR(formula = CRIME ~ INC + HOVAL, data = COL.OLD, W = W1)
summary_SAR(result_geesar)
# }
Run the code above in your browser using DataLab