Prepares a list of data required for using the SAR model; this is for working with (large) raster data files. For non-raster analysis, see prep_sar_data .
Source
Griffith, Daniel A. (2000). Eigenfunction properties and approximations of selected incidence matrices employed in spatial analyses. Linear Algebra and its Applications 321 (1-3): 95-112. tools:::Rd_expr_doi("10.1016/S0024-3795(00)00031-8") .
prep_sar_data2(row, col, quiet =FALSE)
Arguments
row: Number of rows in the raster
col: Number of columns in the raster
quiet: Controls printing behavior. By default, quiet = FALSE and the range of permissible values for the spatial dependence parameter is printed to the console.
Returns
A list containing all of the data elements required by the SAR model in stan_sar.
Details
Prepare data for the SAR model when your raw dataset consists of observations on a regular tessellation, such as a raster layer or remotely sensed imagery. The rook criteria is used to determine adjacency. This function uses Equation 5 from Griffith (2000) to calculate the eigenvalues for a row-standardized spatial weights matrix of a P-by-Q dimension regular tessellation.
This function can accommodate very large numbers of observations for use with stan_sar; for large N data, it is also recommended to use slim = TRUE or the drop argument. For details, see: vignette("raster-regression", package = "geostan").
Examples
row =100col =120sar_dl <- prep_sar_data2(row = row, col = col)