MyNixOS website logo
Description

Bayesian Heteroskedastic Gaussian Processes.

Performs Bayesian posterior inference for heteroskedastic Gaussian processes. Models are trained through MCMC including elliptical slice sampling (ESS) of latent noise processes and Metropolis-Hastings sampling of kernel hyperparameters. Replicates are handled efficientyly through a Woodbury formulation of the joint likelihood for the mean and noise process (Binois, M., Gramacy, R., Ludkovski, M. (2018) <doi:10.1080/10618600.2018.1458625>) For large data, Vecchia-approximation for faster computation is leveraged (Sauer, A., Cooper, A., and Gramacy, R., (2023), <doi:10.1080/10618600.2022.2129662>). Incorporates 'OpenMP' and SNOW parallelization and utilizes 'C'/'C++' under the hood.

This repo consists of an in progress interface for the implementation of a Bayesian HetGP

Run R CMD INSTALL bhGP from outside this folder. If already installed and package is updated: remove package, restart R, and then run this install command.

R files include:

  • functions.R : all functions for log likelihood and Sigma contruction

  • mcmc.R : Does mcmc and ess sampling for theta, nugget vector, etc.

  • gibbs.R : Gibbs sampler for MCMC for seperable, iso-tropic, and combination

  • fit.R : Fits bhGP to the data (X,Y)

  • predict.R: Predicts at new locations Xp

  • mcmc_vecchia.R : Does mcmc and ess sampling for theta, nugget vector, w/ vecchia

  • gibbs_vecchia.R : Gibbs sampler for MCMC for seperable, iso-tropic, and combination w/ vecchia

  • gibbs_vdims_N.R: Gibbs for full N version w/ variance changing in some dimenstions

  • predict_vecchia.R: Predicts at new locations Xp w/ vecchia

  • trim.R: trims mcmc object

  • plots.R: Trace plots and prediction plots

  • vecchia.R: Helper vecchia functions

  • general.R: Helper functions

  • checks.R : checks initilizations

  • inits.R: initialization stratergy

src files include:

  • vecchia.cpp: Vecchia related fucntions

  • cov.cpp: Builds cov matrices

  • invdet.c: inv and det fns

  • general.c: used for invdet

  • RcppExports.cpp: created for package.

Metadata

Version

1.0.1

License

Unknown

Platforms (75)

    Darwin
    FreeBSD
    Genode
    GHCJS
    Linux
    MMIXware
    NetBSD
    none
    OpenBSD
    Redox
    Solaris
    WASI
    Windows
Show all
  • aarch64-darwin
  • aarch64-freebsd
  • aarch64-genode
  • aarch64-linux
  • aarch64-netbsd
  • aarch64-none
  • aarch64-windows
  • aarch64_be-none
  • arm-none
  • armv5tel-linux
  • armv6l-linux
  • armv6l-netbsd
  • armv6l-none
  • armv7a-linux
  • armv7a-netbsd
  • armv7l-linux
  • armv7l-netbsd
  • avr-none
  • i686-cygwin
  • i686-freebsd
  • i686-genode
  • i686-linux
  • i686-netbsd
  • i686-none
  • i686-openbsd
  • i686-windows
  • javascript-ghcjs
  • loongarch64-linux
  • m68k-linux
  • m68k-netbsd
  • m68k-none
  • microblaze-linux
  • microblaze-none
  • microblazeel-linux
  • microblazeel-none
  • mips-linux
  • mips-none
  • mips64-linux
  • mips64-none
  • mips64el-linux
  • mipsel-linux
  • mipsel-netbsd
  • mmix-mmixware
  • msp430-none
  • or1k-none
  • powerpc-netbsd
  • powerpc-none
  • powerpc64-linux
  • powerpc64le-linux
  • powerpcle-none
  • riscv32-linux
  • riscv32-netbsd
  • riscv32-none
  • riscv64-linux
  • riscv64-netbsd
  • riscv64-none
  • rx-none
  • s390-linux
  • s390-none
  • s390x-linux
  • s390x-none
  • vc4-none
  • wasm32-wasi
  • wasm64-wasi
  • x86_64-cygwin
  • x86_64-darwin
  • x86_64-freebsd
  • x86_64-genode
  • x86_64-linux
  • x86_64-netbsd
  • x86_64-none
  • x86_64-openbsd
  • x86_64-redox
  • x86_64-solaris
  • x86_64-windows