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

318

Version

0.1.8

License

GPL (>= 2)

Maintainer

Marcin Jurek

Last Published

January 31st, 2026

Functions in GPvecchia (0.1.8)

SelInv

selected inverse of a sparse matrix
ic0

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

Calculate Matern covariance function
createU

create the sparse triangular U matrix for specific parameters
V2covmat

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

Vecchia Laplace extension of GPVecchia for non-Gaussian data
getMatCov

extract the required elements from the covariance matrix
GPvecchia

GPvecchia: fast, scalable Gaussian process approximations
getMatCovFromFactorCpp

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

create the sparse triangular L matrix for specific parameters
order_middleout

Middle-out ordering
order_coordinate

Sorted coordinate ordering
order_maxmin_exact_obs_pred

Maximum minimum distance ordering for prediction
order_outsidein

Outside-in ordering
ichol

Wrapper for incomplete Cholesky decomposition
vecchia_laplace_likelihood

Wrapper for VL version of vecchia_likelihood
vecchia_laplace_likelihood_from_posterior

Wrapper for VL version of vecchia_likelihood
order_dist_to_point

Distance to specified point ordering
order_maxmin_exact

Maximum minimum distance ordering
vecchia_lincomb

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

Vecchia prediction
vecchia_laplace_prediction

Wrapper for VL version of vecchia_prediction
vecchia_likelihood

evaluation of the likelihood
vecchia_pred

make spatial predictions using Vecchia based on estimated parameters
vecchia_estimate

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

specify a general vecchia approximation