Learn R Programming

GPvecchia

Fast Gaussian-process inference using general Vecchia approximations

For examples of how to use the package, please see the vignettes folder. Please note that GPvecchia is under active development and not stable at this time.

Reporting problems

If you have an issue with GPvecchia, including unsatisfactory runtime, please open a Github ticket.

References

Katzfuss, M., & Guinness, J. (2017). A general framework for Vecchia approximations of Gaussian processes. arXiv:1708.06302.

Katzfuss, M., Guinness, J., Gong, W., & Zilber, D. (2018). Vecchia approximations of Gaussian-process predictions. arXiv:1805.03309.

Zilber, D., & Katzfuss, M. (2019). Vecchia-Laplace approximations of generalized Gaussian processes for big non-Gaussian spatial data. arXiv:1906.07828.

Installation

Note that Rtools is required for compiling C/C++ with OpenMP on Windows systems. When installing Rtools, the system PATH needs to be set so that the C++ compiler included in Rtools can be found by R. Once Rtools is installed, system('g++ -v') can be used to check if the compiler is accessible from R.

Copy Link

Version

Install

install.packages('GPvecchia')

Monthly Downloads

415

Version

0.1.7

License

GPL (>= 2)

Maintainer

Marcin Jurek

Last Published

March 12th, 2024

Functions in GPvecchia (0.1.7)

vecchia_lincomb

linear combination of predictions compute the distribution of a linear combination Hy
vecchia_pred

make spatial predictions using Vecchia based on estimated parameters
order_maxmin_exact

Maximum minimum distance ordering
vecchia_prediction

Vecchia prediction
vecchia_specify

specify a general vecchia approximation
vecchia_laplace_likelihood

Wrapper for VL version of vecchia_likelihood
vecchia_likelihood

evaluation of the likelihood
vecchia_laplace_prediction

Wrapper for VL version of vecchia_prediction
vecchia_laplace_likelihood_from_posterior

Wrapper for VL version of vecchia_likelihood
SelInv

selected inverse of a sparse matrix
V2covmat

compute covariance matrix from V.ord Do not run this function for large n or n.p!!!
ic0

Incomplete Cholesky decomposition of a sparse matrix passed in the compressed sparse row format
getMatCovFromFactorCpp

Calculate the covariance values required by HV for matrix factors passed as sparse matrices
getMatCov

extract the required elements from the covariance matrix
order_middleout

Middle-out ordering
calculate_posterior_VL

Vecchia Laplace extension of GPVecchia for non-Gaussian data
order_coordinate

Sorted coordinate ordering
order_dist_to_point

Distance to specified point ordering
createU

create the sparse triangular U matrix for specific parameters
order_maxmin_exact_obs_pred

Maximum minimum distance ordering for prediction
createL

create the sparse triangular L matrix for specific parameters
ichol

Wrapper for incomplete Cholesky decomposition
GPvecchia

GPvecchia: fast, scalable Gaussian process approximations
order_outsidein

Outside-in ordering
MaternFun

Calculate Matern covariance function
vecchia_estimate

estimate mean and covariance parameters of a Matern covariance function using Vecchia