| GP.summary {BayesGPfit} | R Documentation |
Summary of posterior inference on the Bayesian Gaussian process regression model
Description
Summary of posterior inference on the Bayesian Gaussian process regression model
Usage
GP.summary(GP_fit)
Arguments
GP_fit |
An output object of function GP.Bayes.fit or GP.fast.Bayes.fit. Please refer to them for details. |
Value
A list object consisting of the following elements:
- mean
A list object for posterior mean of the target function,consisting of two elements (f is a vector for function values; x is a vector or matrix for points evaluated).
- work_x
A matrix of real numbers for the standardized grid points for the model fitting. It has the same dimension as "x".
- uci
A list object for 95% upper bound of the creditible interval (uci) of the taget function, consisting of two elements (f is a vector for function values; x is a vector or matrix for points evaluated).
- lci
A list object for 95% lower bound of the creditibel interval (lci) of the taget function, consisting of two elements (f is a vector for function values; x is a vector or matrix for points evaluated).
- sigma2
A vector of posteror mean, the 95% lcl and ucl for variance of the random error.
- tau2
A vector of posterior mean, the 95% lcl and ucl for variance of the target function (hyperparameters).
Author(s)
Jian Kang <jiankang@umich.edu>
Examples
library(BayesGPfit)
library(lattice)
set.seed(1224)
dat = list()
dat$x = GP.generate.grids(d=2,num_grids = 30)
curve = GP.simulate.curve.fast(dat$x,a=0.01,b=0.5,poly_degree=20L)
dat$f = curve$f + rnorm(length(curve$f),sd=1)
fast_fit = GP.fast.Bayes.fit(dat$f,dat$x,a=0.01,b=0.5,poly_degree=20L,progress_bar = TRUE)
reg_fit = GP.Bayes.fit(dat$f,dat$x,a=0.01,b=0.5,poly_degree=20L,progress_bar = TRUE)
sum_fast_fit = GP.summary(fast_fit)
sum_reg_fit = GP.summary(reg_fit)
curves = list(mean_fast = sum_fast_fit$mean,
mean = sum_reg_fit$mean,
lci_fast = sum_fast_fit$lci,
lci = sum_reg_fit$lci,
uci_fast = sum_fast_fit$uci,
uci = sum_reg_fit$uci)
GP.plot.curves(curves,layout=c(2,3))