A nearest‐neighbour Gaussian process spatial factor model for censored, multi‐depth geochemical data

Tilman M. Davies,Sudipto Banerjee,Adam P. Martin,Rose E. Turnbull

DOI: https://doi.org/10.1111/rssc.12565

2022-05-21

Abstract:Journal of the Royal Statistical Society: Series C (Applied Statistics), EarlyView.

What problem does this paper attempt to address?

Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets

Abhirup Datta,Sudipto Banerjee,Andrew O. Finley,Alan E. Gelfand

DOI: https://doi.org/10.48550/arXiv.1406.7343

2016-01-01

Abstract:Spatial process models for analyzing geostatistical data entail computations that become prohibitive as the number of spatial locations become large. This manuscript develops a class of highly scalable Nearest Neighbor Gaussian Process (NNGP) models to provide fully model-based inference for large geostatistical datasets. We establish that the NNGP is a well-defined spatial process providing legitimate finite-dimensional Gaussian densities with sparse precision matrices. We embed the NNGP as a sparsity-inducing prior within a rich hierarchical modeling framework and outline how computationally efficient Markov chain Monte Carlo (MCMC) algorithms can be executed without storing or decomposing large matrices. The floating point operations (flops) per iteration of this algorithm is linear in the number of spatial locations, thereby rendering substantial scalability. We illustrate the computational and inferential benefits of the NNGP over competing methods using simulation studies and also analyze forest biomass from a massive United States Forest Inventory dataset at a scale that precludes alternative dimension-reducing methods.

Methodology
Computationally Scalable Bayesian SPDE Modeling for Censored Spatial Responses

Indranil Sahoo,Suman Majumder,Arnab Hazra,Ana G. Rappold,Dipankar Bandyopadhyay

2024-03-23

Abstract:Observations of groundwater pollutants, such as arsenic or Perfluorooctane sulfonate (PFOS), are riddled with left censoring. These measurements have impact on the health and lifestyle of the populace. Left censoring of these spatially correlated observations are usually addressed by applying Gaussian processes (GPs), which have theoretical advantages. However, this comes with a challenging computational complexity of $\mathcal{O}(n^3)$, which is impractical for large datasets. Additionally, a sizable proportion of the data being left-censored creates further bottlenecks, since the likelihood computation now involves an intractable high-dimensional integral of the multivariate Gaussian density. In this article, we tackle these two problems simultaneously by approximating the GP with a Gaussian Markov random field (GMRF) approach that exploits an explicit link between a GP with Mat\'ern correlation function and a GMRF using stochastic partial differential equations (SPDEs). We introduce a GMRF-based measurement error into the model, which alleviates the likelihood computation for the censored data, drastically improving the speed of the model while maintaining admirable accuracy. Our approach demonstrates robustness and substantial computational scalability, compared to state-of-the-art methods for censored spatial responses across various simulation settings. Finally, the fit of this fully Bayesian model to the concentration of PFOS in groundwater available at 24,959 sites across California, where 46.62\% responses are censored, produces prediction surface and uncertainty quantification in real time, thereby substantiating the applicability and scalability of the proposed method. Code for implementation is made available via GitHub.

Computation,Applications,Methodology
Nearest-Neighbor Mixture Models for Non-Gaussian Spatial Processes

Xiaotian Zheng,Athanasios Kottas,Bruno Sansó

DOI: https://doi.org/10.48550/arXiv.2107.07736

2022-06-27

Abstract:We develop a class of nearest-neighbor mixture models that provide direct, computationally efficient, probabilistic modeling for non-Gaussian geospatial data. The class is defined over a directed acyclic graph, which implies conditional independence in representing a multivariate distribution through factorization into a product of univariate conditionals, and is extended to a full spatial process. We model each conditional as a mixture of spatially varying transition kernels, with locally adaptive weights, for each one of a given number of nearest neighbors. The modeling framework emphasizes the description of non-Gaussian dependence at the data level, in contrast with approaches that introduce a spatial process for transformed data, or for functionals of the data probability distribution. Thus, it facilitates efficient, full simulation-based inference. We study model construction and properties analytically through specification of bivariate distributions that define the local transition kernels, providing a general strategy for modeling general types of non-Gaussian data. Regarding computation, the framework lays out a new approach to handling spatial data sets, leveraging a mixture model structure to avoid computational issues that arise from large matrix operations. We illustrate the methodology using synthetic data examples and an analysis of Mediterranean Sea surface temperature observations.

Methodology
Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

Lu Zhang,Sudipto Banerjee

DOI: https://doi.org/10.1111/biom.13452

2021-02-15

Abstract:Multivariate spatially-oriented data sets are prevalent in the environmental and physical sciences. Scientists seek to jointly model multiple variables, each indexed by a spatial location, to capture any underlying spatial association for each variable and associations among the different dependent variables. Multivariate latent spatial process models have proved effective in driving statistical inference and rendering better predictive inference at arbitrary locations for the spatial process. High-dimensional multivariate spatial data, which is the theme of this article, refers to data sets where the number of spatial locations and the number of spatially dependent variables is very large. The field has witnessed substantial developments in scalable models for univariate spatial processes, but such methods for multivariate spatial processes, especially when the number of outcomes are moderately large, are limited in comparison. Here, we extend scalable modeling strategies for a single process to multivariate processes. We pursue Bayesian inference which is attractive for full uncertainty quantification of the latent spatial process. Our approach exploits distribution theory for the Matrix-Normal distribution, which we use to construct scalable versions of a hierarchical linear model of coregionalization (LMC) and spatial factor models that deliver inference over a high-dimensional parameter space including the latent spatial process. We illustrate the computational and inferential benefits of our algorithms over competing methods using simulation studies and an analysis of a massive vegetation index data set.

Methodology,Computation
An Adaptive Gaussian Process‐based Method for Efficient Bayesian Experimental Design in Groundwater Contaminant Source Identification Problems

Jiangjiang Zhang,Weixuan Li,Lingzao Zeng,Laosheng Wu

DOI: https://doi.org/10.1002/2016wr018598

IF: 5.4

2016-01-01

Water Resources Research

Abstract:Surrogate models are commonly used in Bayesian approaches such as Markov Chain Monte Carlo (MCMC) to avoid repetitive CPU‐demanding model evaluations. However, the approximation error of a surrogate may lead to biased estimation of the posterior distribution. This bias can be corrected by constructing a very accurate surrogate or implementing MCMC in a two‐stage manner. Since the two‐stage MCMC requires extra original model evaluations after surrogate evaluations, the computational cost is still high. If the information of measurement is incorporated, a locally accurate surrogate can be adaptively constructed with low computational cost. Based on this idea, we integrate Gaussian process (GP) and MCMC to adaptively construct locally accurate surrogates for Bayesian experimental design in groundwater contaminant source identification problems. Moreover, the uncertainty estimate of GP approximation error is incorporated in the Bayesian formula to avoid over‐confident estimation of the posterior distribution. The proposed approach is tested with a numerical case study. Without sacrificing the estimation accuracy, the new approach achieves about 200 times of speed‐up compared to our previous work which implemented MCMC in a two‐stage manner.
Statistical Modeling of Right-Censored Spatial Data Using Gaussian Random Fields

Fathima Z. Sainul Abdeen,Akim Adekpedjou,Sophie Dabo Niang

DOI: https://doi.org/10.3390/math12101521

IF: 2.4

2024-05-14

Mathematics

Abstract:Consider a fixed number of clustered areas identified by their geographical coordinates that are monitored for the occurrences of an event such as a pandemic, epidemic, or migration. Data collected on units at all areas include covariates and environmental factors. We apply a probit transformation to the time to event and embed an isotropic spatial correlation function into our models for better modeling as compared to existing methodologies that use frailty or copula. Composite likelihood technique is employed for the construction of a multivariate Gaussian random field that preserves the spatial correlation function. The data are analyzed using counting process and geostatistical formulation that led to a class of weighted pairwise semiparametric estimating functions. The estimators of model parameters are shown to be consistent and asymptotically normally distributed under infill-type asymptotic spatial statistics. Detailed small sample numerical studies that are in agreement with theoretical results are provided. The foregoing procedures are applied to the leukemia survival data in Northeast England. A comparison to existing methodologies provides improvement.

mathematics
Multivariate nearest‐neighbors Gaussian processes with random covariance matrices

Isabelle Grenier,Bruno Sansó,Jessica L. Matthews

DOI: https://doi.org/10.1002/env.2839

2024-01-03

Environmetrics

Abstract:We propose a non‐stationary spatial model based on a normal‐inverse‐Wishart framework, conditioning on a set of nearest‐neighbors. The model, called nearest‐neighbor Gaussian process with random covariance matrices is developed for both univariate and multivariate spatial settings and allows for fully flexible covariance structures that impose no stationarity or isotropic restrictions. In addition, the model can handle duplicate observations and missing data. We consider an approach based on integrating out the spatial random effects that allows fast inference for the model parameters. We also consider a full hierarchical approach that leverages the sparse structures induced by the model to perform fast Monte Carlo computations. Strong computational efficiency is achieved by leveraging the adaptive localized structure of the model that allows for a high level of parallelization. We illustrate the performance of the model with univariate and bivariate simulations, as well as with observations from two stationary satellites consisting of albedo measurements.

environmental sciences,statistics & probability,mathematics, interdisciplinary applications
Bayesian Geostatistical Modeling for Discrete-Valued Processes

Xiaotian Zheng,Athanasios Kottas,Bruno Sansó

DOI: https://doi.org/10.48550/arXiv.2111.01840

2022-03-03

Abstract:We introduce a flexible and scalable class of Bayesian geostatistical models for discrete data, based on the class of nearest neighbor mixture transition distribution processes (NNMP), referred to as discrete NNMP. The proposed class characterizes spatial variability by a weighted combination of first-order conditional probability mass functions (pmfs) for each one of a given number of neighbors. The approach supports flexible modeling for multivariate dependence through specification of general bivariate discrete distributions that define the conditional pmfs. Moreover, the discrete NNMP allows for construction of models given a pre-specified family of marginal distributions that can vary in space, facilitating covariate inclusion. In particular, we develop a modeling and inferential framework for copula-based NNMPs that can attain flexible dependence structures, motivating the use of bivariate copula families for spatial processes. Compared to the traditional class of spatial generalized linear mixed models, where spatial dependence is introduced through a transformation of response means, our process-based modeling approach provides both computational and inferential advantages. We illustrate the benefits with synthetic data examples and an analysis of North American Breeding Bird Survey data.

Methodology
Spatial Hyperspheric Models for Compositional Data

Michael R. Schwob,Mevin B. Hooten,Nicholas M. Calzada

DOI: https://doi.org/10.48550/arXiv.2410.03648

2024-10-09

Abstract:Compositional data are an increasingly prevalent data source in spatial statistics. Analysis of such data is typically done on log-ratio transformations or via Dirichlet regression. However, these approaches often make unnecessarily strong assumptions (e.g., strictly positive components, exclusively negative correlations). An alternative approach uses square-root transformed compositions and directional distributions. Such distributions naturally allow for zero-valued components and positive correlations, yet they may include support outside the non-negative orthant and are not generative for compositional data. To overcome this challenge, we truncate the elliptically symmetric angular Gaussian (ESAG) distribution to the non-negative orthant. Additionally, we propose a spatial hyperspheric regression that contains fixed and random multivariate spatial effects. The proposed method also contains a term that can be used to propagate uncertainty that may arise from precursory stochastic models (i.e., machine learning classification). We demonstrate our method on a simulation study and on classified bioacoustic signals of the Dryobates pubescens (downy woodpecker).

Methodology,Applications,Computation
High-Dimensional Bayesian Geostatistics

Sudipto Banerjee

DOI: https://doi.org/10.48550/arXiv.1705.07265

2017-05-20

Abstract:With the growing capabilities of Geographic Information Systems (GIS) and user-friendly software, statisticians today routinely encounter geographically referenced data containing observations from a large number of spatial locations and time points. Over the last decade, hierarchical spatiotemporal process models have become widely deployed statistical tools for researchers to better understand the complex nature of spatial and temporal variability. However, fitting hierarchical spatiotemporal models often involves expensive matrix computations with complexity increasing in cubic order for the number of spatial locations and temporal points. This renders such models unfeasible for large data sets. This article offers a focused review of two methods for constructing well-defined highly scalable spatiotemporal stochastic processes. Both these processes can be used as "priors" for spatiotemporal random fields. The first approach constructs a low-rank process operating on a lower-dimensional subspace. The second approach constructs a Nearest-Neighbor Gaussian Process (NNGP) that ensures sparse precision matrices for its finite realizations. Both processes can be exploited as a scalable prior embedded within a rich hierarchical modeling framework to deliver full Bayesian inference. These approaches can be described as model-based solutions for big spatiotemporal datasets. The models ensure that the algorithmic complexity has $\sim n$ floating point operations (flops), where $n$ the number of spatial locations (per iteration). We compare these methods and provide some insight into their methodological underpinnings.

Methodology
A Vecchia Approximation for High-Dimensional Gaussian Cumulative Distribution Functions Arising from Spatial Data

Mauricio Nascimento,Benjamin A. Shaby

DOI: https://doi.org/10.48550/arXiv.2007.15195

2020-07-30

Abstract:We introduce an approach to quickly and accurately approximate the cumulative distribution function of multivariate Gaussian distributions arising from spatial Gaussian processes. This approximation is trivially parallelizable and simple to implement using standard software. We demonstrate its accuracy and computational efficiency in a series of simulation experiments and apply it to analyzing the joint tail of a large precipitation dataset using a recently-proposed scale mixture model for spatial extremes. This dataset is many times larger than what was previously considered possible to fit using preferred inferential techniques.

Computation
A multivariate Gaussian random field prior against spatial confounding

Isa Marques,Thomas Kneib,Nadja Klein

DOI: https://doi.org/10.48550/arXiv.2106.03737

2021-06-07

Methodology

Abstract:Spatial models are used in a variety research areas, such as environmental sciences, epidemiology, or physics. A common phenomenon in many spatial regression models is spatial confounding. This phenomenon takes place when spatially indexed covariates modeling the mean of the response are correlated with the spatial random effect. As a result, estimates for regression coefficients of the covariates can be severely biased and interpretation of these is no longer valid. Recent literature has shown that typical solutions for reducing spatial confounding can lead to misleading and counterintuitive results. In this paper, we develop a computationally efficient spatial model in a Bayesian framework integrating novel prior structure that reduces spatial confounding. Starting from the univariate case, we extend our prior structure to case of multiple spatially confounded covariates. In a simulation study, we show that our novel model flexibly detects and reduces spatial confounding in spatial datasets, and it performs better than typically used methods such as restricted spatial regression. These results are promising for any applied researcher who wishes to interpret covariate effects in spatial regression models. As a real data illustration, we study the effect of elevation and temperature on the mean of daily precipitation in Germany.
Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains

Michele Peruzzi,Sudipto Banerjee,Andrew O. Finley

DOI: https://doi.org/10.1080/01621459.2020.1833889

2020-10-05

Abstract:We introduce a class of scalable Bayesian hierarchical models for the analysis of massive geostatistical datasets. The underlying idea combines ideas on high-dimensional geostatistics by partitioning the spatial domain and modeling the regions in the partition using a sparsity-inducing directed acyclic graph (DAG). We extend the model over the DAG to a well-defined spatial process, which we call the Meshed Gaussian Process (MGP). A major contribution is the development of a MGPs on tessellated domains, accompanied by a Gibbs sampler for the efficient recovery of spatial random effects. In particular, the cubic MGP (Q-MGP) can harness high-performance computing resources by executing all large-scale operations in parallel within the Gibbs sampler, improving mixing and computing time compared to sequential updating schemes. Unlike some existing models for large spatial data, a Q-MGP facilitates massive caching of expensive matrix operations, making it particularly apt in dealing with spatiotemporal remote-sensing data. We compare Q-MGPs with large synthetic and real world data against state-of-the-art methods. We also illustrate using Normalized Difference Vegetation Index (NDVI) data from the Serengeti park region to recover latent multivariate spatiotemporal random effects at millions of locations. The source code is available at <a class="link-external link-https" href="https://github.com/mkln/meshgp" rel="external noopener nofollow">this https URL</a>.

Methodology,Computation
Spatiotemporal clustering using Gaussian processes embedded in a mixture model

Jarno Vanhatalo,Scott D. Foster,Geoffrey R. Hosack

DOI: https://doi.org/10.1002/env.2681

2021-05-07

Environmetrics

Abstract:<p>The categorization of multidimensional data into clusters is a common task in statistics. Many applications of clustering, including the majority of tasks in ecology, use data that is inherently spatial and is often also temporal. However, spatiotemporal dependence is typically ignored when clustering multivariate data. We present a finite mixture model for spatial and spatiotemporal clustering that incorporates spatial and spatiotemporal autocorrelation by including appropriate Gaussian processes (GP) into a model for the mixing proportions. We also allow for flexible and semiparametric dependence on environmental covariates, once again using GPs. We propose to use Bayesian inference through three tiers of approximate methods: a Laplace approximation that allows efficient analysis of large datasets, and both partial and full Markov chain Monte Carlo (MCMC) approaches that improve accuracy at the cost of increased computational time. Comparison of the methods shows that the Laplace approximation is a useful alternative to the MCMC methods. A decadal analysis of 253 species of teleost fish from 854 samples collected along the biodiverse northwestern continental shelf of Australia between 1986 and 1997 shows the added clarity provided by accounting for spatial autocorrelation. For these data, the temporal dependence is comparatively small, which is an important finding given the changing human pressures over this time.</p>

environmental sciences,statistics & probability,mathematics, interdisciplinary applications
Gaussian process machine learning and Kriging for groundwater salinity interpolation

Tao Cui,Dan Pagendam,Mat Gilfedder

DOI: https://doi.org/10.1016/j.envsoft.2021.105170

2021-10-01

Abstract:Gaussian processes (GPs) provide statistically optimal predictions in the sense of unbiasedness and maximal precision. Although the modern implementation of GPs as a machine learning technique is more capable and flexible than Kriging, their employment in environmental science is less routine. Their flexibility and capability as a spatial data interpolation technique are demonstrated by applying them to groundwater salinity prediction in a data-sparse region in Australia. By learning from multiple data sources, including AEM and DEM data, GPs have generated groundwater salinity maps with rich local details and quantified uncertainty to support risk-based decision making. The results demonstrate the great worth of nonpoint data with regional spatial coverage to provide more realistic heterogeneity in aquifer properties that are critical for many studies such as contaminant transport. GPs should be further encouraged in groundwater science for data interpolation and prediction, especially when point measurements are sparse and multiple predictors are available.

environmental sciences,engineering, environmental,water resources,computer science, interdisciplinary applications
Sparse Approximate Inference for Spatio-Temporal Point Process Models

Botond Cseke,Andrew Zammit Mangion,Tom Heskes,Guido Sanguinetti

DOI: https://doi.org/10.48550/arXiv.1305.4152

2015-07-06

Abstract:Spatio-temporal point process models play a central role in the analysis of spatially distributed systems in several disciplines. Yet, scalable inference remains computa- tionally challenging both due to the high resolution modelling generally required and the analytically intractable likelihood function. Here, we exploit the sparsity structure typical of (spatially) discretised log-Gaussian Cox process models by using approximate message-passing algorithms. The proposed algorithms scale well with the state dimension and the length of the temporal horizon with moderate loss in distributional accuracy. They hence provide a flexible and faster alternative to both non-linear filtering-smoothing type algorithms and to approaches that implement the Laplace method or expectation propagation on (block) sparse latent Gaussian models. We infer the parameters of the latent Gaussian model using a structured variational Bayes approach. We demonstrate the proposed framework on simulation studies with both Gaussian and point-process observations and use it to reconstruct the conflict intensity and dynamics in Afghanistan from the WikiLeaks Afghan War Diary.

Machine Learning
Efficient Modeling of Spatial Extremes over Large Geographical Domains

Arnab Hazra,Raphaël Huser,David Bolin

DOI: https://doi.org/10.48550/arXiv.2112.10248

2024-04-30

Abstract:Various natural phenomena exhibit spatial extremal dependence at short spatial distances. However, existing models proposed in the spatial extremes literature often assume that extremal dependence persists across the entire domain. This is a strong limitation when modeling extremes over large geographical domains, and yet it has been mostly overlooked in the literature. We here develop a more realistic Bayesian framework based on a novel Gaussian scale mixture model, with the Gaussian process component defined by a stochastic partial differential equation yielding a sparse precision matrix, and the random scale component modeled as a low-rank Pareto-tailed or Weibull-tailed spatial process determined by compactly-supported basis functions. We show that our proposed model is approximately tail-stationary and that it can capture a wide range of extremal dependence structures. Its inherently sparse structure allows fast Bayesian computations in high spatial dimensions based on a customized Markov chain Monte Carlo algorithm prioritizing calibration in the tail. We fit our model to analyze heavy monsoon rainfall data in Bangladesh. Our study shows that our model outperforms natural competitors and that it fits precipitation extremes well. We finally use the fitted model to draw inference on long-term return levels for marginal precipitation and spatial aggregates.

Methodology
Spatial Latent Gaussian Modelling with Change of Support

Erick A. Chacón-Montalván,Peter M. Atkinson,Christopher Nemeth,Benjamin M. Taylor,Paula Moraga

2024-03-13

Abstract:Spatial data are often derived from multiple sources (e.g. satellites, in-situ sensors, survey samples) with different supports, but associated with the same properties of a spatial phenomenon of interest. It is common for predictors to also be measured on different spatial supports than the response variables. Although there is no standard way to work with spatial data with different supports, a prevalent approach used by practitioners has been to use downscaling or interpolation to project all the variables of analysis towards a common support, and then using standard spatial models. The main disadvantage with this approach is that simple interpolation can introduce biases and, more importantly, the uncertainty associated with the change of support is not taken into account in parameter estimation. In this article, we propose a Bayesian spatial latent Gaussian model that can handle data with different rectilinear supports in both the response variable and predictors. Our approach allows to handle changes of support more naturally according to the properties of the spatial stochastic process being used, and to take into account the uncertainty from the change of support in parameter estimation and prediction. We use spatial stochastic processes as linear combinations of basis functions where Gaussian Markov random fields define the weights. Our hierarchical modelling approach can be described by the following steps: (i) define a latent model where response variables and predictors are considered as latent stochastic processes with continuous support, (ii) link the continuous-index set stochastic processes with its projection to the support of the observed data, (iii) link the projected process with the observed data. We show the applicability of our approach by simulation studies and modelling land suitability for improved grassland in Rhondda Cynon Taf, a county borough in Wales.

Methodology
A Low-Rank Bayesian Approach for Geoadditive Modeling

Bryan Sumalinab,Oswaldo Gressani,Niel Hens,Christel Faes

DOI: https://doi.org/10.48550/arXiv.2407.05854

2024-07-08

Abstract:Kriging is an established methodology for predicting spatial data in geostatistics. Current kriging techniques can handle linear dependencies on spatially referenced covariates. Although splines have shown promise in capturing nonlinear dependencies of covariates, their combination with kriging, especially in handling count data, remains underexplored. This paper proposes a novel Bayesian approach to the low-rank representation of geoadditive models, which integrates splines and kriging to account for both spatial correlations and nonlinear dependencies of covariates. The proposed method accommodates Gaussian and count data inherent in many geospatial datasets. Additionally, Laplace approximations to selected posterior distributions enhances computational efficiency, resulting in faster computation times compared to Markov chain Monte Carlo techniques commonly used for Bayesian inference. Method performance is assessed through a simulation study, demonstrating the effectiveness of the proposed approach. The methodology is applied to the analysis of heavy metal concentrations in the Meuse river and vulnerability to the coronavirus disease 2019 (COVID-19) in Belgium. Through this work, we provide a new flexible and computationally efficient framework for analyzing spatial data.

Methodology
Spatial Factor Models for High-Dimensional and Large Spatial Data: An Application in Forest Variable Mapping

Daniel Taylor-Rodriguez,Andrew O. Finley,Abhirup Datta,Chad Babcock,Hans-Erik Andersen,Bruce D. Cook,Douglas C. Morton,Sudipto Banerjee

DOI: https://doi.org/10.48550/arXiv.1801.02078

2018-11-09

Abstract:Gathering information about forest variables is an expensive and arduous activity. As such, directly collecting the data required to produce high-resolution maps over large spatial domains is infeasible. Next generation collection initiatives of remotely sensed Light Detection and Ranging (LiDAR) data are specifically aimed at producing complete-coverage maps over large spatial domains. Given that LiDAR data and forest characteristics are often strongly correlated, it is possible to make use of the former to model, predict, and map forest variables over regions of interest. This entails dealing with the high-dimensional ($\sim$$10^2$) spatially dependent LiDAR outcomes over a large number of locations (~10^5-10^6). With this in mind, we develop the Spatial Factor Nearest Neighbor Gaussian Process (SF-NNGP) model, and embed it in a two-stage approach that connects the spatial structure found in LiDAR signals with forest variables. We provide a simulation experiment that demonstrates inferential and predictive performance of the SF-NNGP, and use the two-stage modeling strategy to generate complete-coverage maps of forest variables with associated uncertainty over a large region of boreal forests in interior Alaska.

Applications

A nearest‐neighbour Gaussian process spatial factor model for censored, multi‐depth geochemical data

Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets

Computationally Scalable Bayesian SPDE Modeling for Censored Spatial Responses

Nearest-Neighbor Mixture Models for Non-Gaussian Spatial Processes

Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

An Adaptive Gaussian Process‐based Method for Efficient Bayesian Experimental Design in Groundwater Contaminant Source Identification Problems

Statistical Modeling of Right-Censored Spatial Data Using Gaussian Random Fields

Multivariate nearest‐neighbors Gaussian processes with random covariance matrices

Bayesian Geostatistical Modeling for Discrete-Valued Processes

Spatial Hyperspheric Models for Compositional Data

High-Dimensional Bayesian Geostatistics

A Vecchia Approximation for High-Dimensional Gaussian Cumulative Distribution Functions Arising from Spatial Data

A multivariate Gaussian random field prior against spatial confounding

Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains

Spatiotemporal clustering using Gaussian processes embedded in a mixture model

Gaussian process machine learning and Kriging for groundwater salinity interpolation

Sparse Approximate Inference for Spatio-Temporal Point Process Models

Efficient Modeling of Spatial Extremes over Large Geographical Domains

Spatial Latent Gaussian Modelling with Change of Support

A Low-Rank Bayesian Approach for Geoadditive Modeling

Spatial Factor Models for High-Dimensional and Large Spatial Data: An Application in Forest Variable Mapping