Skip to main content
SpringerLink
Log in

A Fused Lasso Approach to Nonstationary Spatial Covariance Estimation

  • Published:
Journal of Agricultural, Biological, and Environmental Statistics Aims and scope Submit manuscript

Abstract

Spatial data are increasing in size and complexity due to technological advances. For an analysis of a large and diverse spatial domain, simplifying assumptions such as stationarity are questionable and standard computational algorithms are inadequate. In this paper, we propose a computationally efficient method to estimate a nonstationary covariance function. We partition the spatial domain into a fine grid of subregions and assign each subregion its own set of spatial covariance parameters. This introduces a large number of parameters and to stabilize the procedure we impose a penalty to spatially smooth the estimates. By penalizing the absolute difference between parameters for adjacent subregions, the solution can be identical for adjacent subregions and thus the method identifies stationary subdomains. To apply the method to large datasets, we use a block composite likelihood which is natural in this setting because it also operates on a partition of the spatial domain. The method is applied to tropospheric ozone in the US, and we find that the spatial covariance on the west coast differs from the rest of the country.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4

Similar content being viewed by others

References

  • Appel, K., Pouliot, G., Simon, H., Sarwar, G., Pye, H., Napelenok, S., Akhtar, F., and Roselle, S. (2013), “Evaluation of dust and trace metal estimates from the Community Multiscale Air Quality (CMAQ) model version 5.0,” Geoscientific Model Development, 6, 883–899.

  • Banerjee, S., Gelfand, A. E., Finley, A. O., and Sang, H. (2008), “Gaussian predictive process models for large spatial data sets,” Journal of the Royal Statistical Society: Series B, 70, 825–848.

  • Besag, J., York, J., and Mollié, A. (1991), “Bayesian image restoration, with two applications in spatial statistics,” Annals of the Institute of Statistical Mathematics, 43, 1–20.

  • Bien, J. and Tibshirani, R. J. (2011), “Sparse estimation of a covariance matrix,” Biometrika, 98, 807–820.

  • Bornn, L., Shaddick, G., and Zidek, J. V. (2012), “Modeling nonstationary processes through dimension expansion,” Journal of the American Statistical Association, 107, 281–289.

  • Chang, Y.-M., Hsu, N.-J., and Huang, H.-C. (2010), “Semiparametric estimation and selection for nonstationary spatial covariance functions,” Journal of Computational and Graphical Statistics, 19, 117–139.

  • Cressie, N. and Johannesson, G. (2008), “Fixed rank kriging for very large spatial data sets,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), 70, 209–226.

  • Eidsvik, J., Shaby, B. A., Reich, B. J., Wheeler, M., and Niemi, J. (2014), “Estimation and prediction in spatial models with block composite likelihoods,” Journal of Computational and Graphical Statistics, 23, 295–315.

  • Finley, A. O., Sang, H., Banerjee, S., and Gelfand, A. E. (2009), “Improving the performance of predictive process modeling for large datasets,” Computational Statistics & Data Analysis, 53, 2873–2884.

  • Friedman, J. and Popescu, B. E. (2003), “Gradient directed regularization for linear regression and classification”, Tech. rep., Statistics Department, Stanford University.

  • Fuentes, M. (2002), “Spectral methods for nonstationary spatial processes,” Biometrika, 89, 197–210.

  • —— (2007), “Approximate likelihood for large irregularly spaced spatial data,” Journal of the American Statistical Association, 102, 321–331.

  • Furrer, R., Genton, M. G., and Nychka, D. (2006), “Covariance tapering for interpolation of large spatial datasets,” Journal of Computational and Graphical Statistics, 15.

  • Higdon, D. (1998), “A process-convolution approach to modelling temperatures in the North Atlantic Ocean,” Environmental and Ecological Statistics, 5, 173–190.

  • Higdon, D., Swall, J., and Kern, J. (1999), “Non-stationary spatial modeling”, in Bayesian Statistics 6, pp. 761–768.

  • Hsu, N.-J., Chang, Y.-M., and Huang, H.-C. (2012), “A group lasso approach for non-stationary spatial-temporal covariance estimation,” Environmetrics, 23, 12–23.

  • Kaufman, C. G., Schervish, M. J., and Nychka, D. W. (2008), “Covariance tapering for likelihood-based estimation in large spatial data sets,” Journal of the American Statistical Association, 103, 1545–1555.

  • Land, S. R. and Friedman, J. H. (1997), “Variable fusion: A new adaptive signal regression method,” Tech. rep., Department of Statistics, Carnegie Mellon University, Pittsburgh.

  • Lin, Y. and Zhang, H. H. (2006), “Component selection and smoothing in smoothing spline analysis of variance models,” Annals of Statistics, 34, 2272–2297.

  • Lindgren, F., Rue, H., and Lindström, J. (2011), “An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), 73, 423–498.

  • Neto, J. H. V., Schmidt, A. M., and Guttorp, P. (2014), “Accounting for spatially varying directional effects in spatial covariance structures,” Journal of the Royal Statistical Society: Series C (Applied Statistics), 63, 103–122.

  • Nychka, D. and Saltzman, N. (1998), “Design of air-quality monitoring networks,” in Case studies in environmental statistics, Springer, pp. 51–76.

  • Paciorek, C. J. and Schervish, M. J. (2006), “Spatial modelling using a new class of nonstationary covariance functions,” Environmetrics, 17, 483–506.

  • Reich, B. J., Chang, H. H., and Foley, K. M. (2014), “A spectral method for spatial downscaling,” Biometrics, 70, 932–942.

  • Reich, B. J., Eidsvik, J., Guindani, M., Nail, A. J., and Schmidt, A. M. (2011), “A class of covariate-dependent spatiotemporal covariance functions,” The annals of applied statistics, 5, 2265.

  • Sampson, P. D. (2010), ”Constructions for Nonstationary Spatial Processes”, in Handbook of Spatial Statistics, eds. Gelfand, A. E., Diggle, P. J., Fuentes, M., and Guttorp, P., CRC Press, chap. 9.

  • Sampson, P. D. and Guttorp, P. (1992), “Nonparametric estimation of nonstationary spatial covariance structure,” Journal of the American Statistical Association, 87, 108–119.

  • Stein, M. L., Chi, Z., and Welty, L. J. (2004), “Approximating likelihoods for large spatial data sets,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), 66, 275–296.

  • Storlie, C. B., Bondell, H. D., Reich, B. J., and Zhang, H. H. (2011), “Surface estimation, variable selection, and the nonparametric oracle property,” Statistica Sinica, 21, 679.

  • Tibshirani, R., Saunders, M., Rosset, S., Zhu, J., and Knight, K. (2005), “Sparsity and smoothness via the fused lasso,” Journal of the Royal Statistical Society: Series B, 67, 91–108.

  • Tibshirani, R. J. and Taylor, J. (2011), “The solution path of the generalized lasso,” Annals of Statistics, 39, 1335–1371.

Download references

Author information

Authors and Affiliations

  1. SAS Institute, Cary, NC, USA

    Ryan J. Parker

  2. North Carolina State University, Raleigh, NC, USA

    Brian J. Reich

  3. Norwegian University of Science and Technology, Trondheim, Norway

    Jo Eidsvik

Authors
  1. Ryan J. Parker
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Brian J. Reich
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jo Eidsvik
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Ryan J. Parker.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary material 1 (pdf 89 KB)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Parker, R.J., Reich, B.J. & Eidsvik, J. A Fused Lasso Approach to Nonstationary Spatial Covariance Estimation. JABES 21, 569–587 (2016). https://doi.org/10.1007/s13253-016-0251-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s13253-016-0251-8

Keywords

Search

Navigation