more fixes to loo2-non-factorized.Rmd

Author: jgabry

vignettes/loo2-non-factorized.Rmd

@@ -463,9 +463,9 @@ loop over the observations.
 ```{r fit_dummy, cache = TRUE}
 # see help("mi", "brms") for details on the mi() usage
 fit_dummy <- brm(
-  CRIME | mi() ~ INC + HOVAL, 
+  CRIME | mi() ~ INC + HOVAL + sar(COL.nb, type = "lag"), 
   data = COL.OLD,
-  autocor = cor_lagsar(COL.nb), 
+  data2 = list(COL.nb = COL.nb),
   chains = 0
 )
 ```
@@ -494,7 +494,7 @@ for (i in seq_len(N)) {
     y_miss_i <- y
     y_miss_i[i] <- p$Ymi
     eta <- p$b_Intercept + p$b_INC * fit_i$data$INC + p$b_HOVAL * fit_i$data$HOVAL
-    W_tilde <- diag(N) - p$lagsar * fit_i$autocor$W
+    W_tilde <- diag(N) - p$lagsar * spdep::nb2mat(COL.nb)
     Cinv <- t(W_tilde) %*% W_tilde / p$sigma^2
     g <- Cinv %*% (y_miss_i - solve(W_tilde, eta))
     cbar <- diag(Cinv);