Estimating Spatial Probit Models in R

In this article we present the Bayesian estimation of spatial probit models in R and provide an implementation in the package spatialprobit. We show that large probit models can be estimated with sparse matrix representations and Gibbs sampling of a truncated multivariate normal distribution with the precision matrix. We present three examples and point to ways to achieve further performance gains through parallelization of the Markov Chain Monte Carlo approach.

Stefan Wilhelm , Miguel Godinho de Matos

CRAN packages used

spBayes, spatial, geoR, sgeostat, spdep, sphet, sna, network, Matrix, sparseM, spatialprobit, McSpatial, LearnBayes, tmvtnorm, mvtnorm, igraph

CRAN Task Views implied by cited packages

Spatial, Bayesian, Distributions, Econometrics, SocialSciences, gR, Multivariate, Optimization, SpatioTemporal, Finance, Graphics, NumericalMathematics, Survival


Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".


For attribution, please cite this work as

Wilhelm & Matos, "The R Journal: Estimating Spatial Probit Models in R", The R Journal, 2013

BibTeX citation

  author = {Wilhelm, Stefan and Matos, Miguel Godinho de},
  title = {The R Journal: Estimating Spatial Probit Models in R},
  journal = {The R Journal},
  year = {2013},
  note = {},
  doi = {10.32614/RJ-2013-013},
  volume = {5},
  issue = {1},
  issn = {2073-4859},
  pages = {130-143}