(1) In execute_pipeline.R, previously we are using queen-way neighbors (including diagonal, vertical and horizontal neighbors), and we change that to rook way neighbors (only including vertical and horizontal neighbors).
(2) In execute_pipeline.R, the diag previously calculates the downstream neighbors, and we need to change that to upstream neighbors.
(3) In Stan code, the node1 and node2 in the stan is also changed to be consistent with the upstream neighbor sets.
(4) When do_time_slice = FALSE, In execute_pipeline.R, I found that eta_effect contains NA value and results in NA value for residual, so I added if else to accommodate that.
Four things have been changed:
(1) In execute_pipeline.R, previously we are using queen-way neighbors (including diagonal, vertical and horizontal neighbors), and we change that to rook way neighbors (only including vertical and horizontal neighbors).
(2) In execute_pipeline.R, the diag previously calculates the downstream neighbors, and we need to change that to upstream neighbors.
(3) In Stan code, the node1 and node2 in the stan is also changed to be consistent with the upstream neighbor sets.
(4) When do_time_slice = FALSE, In execute_pipeline.R, I found that eta_effect contains NA value and results in NA value for residual, so I added if else to accommodate that.