bri.nonstat: Non-stationary smoothing for Gaussian Process Regression in...

In julianfaraway/brinla: Bayesian Regression with INLA

Non-stationary smoothing for Gaussian Process Regression in 1D

Description

Non-stationary smoothing for Gaussian Process Regression in 1D

Usage

bri.nonstat(x, y, nbasis = 25, sbasis = 5, degree = 2, alpha = 2,
  xout = x, sigma0 = sd(y), rho0 = 0.25 * (max(x) - min(x)))

Arguments

x	the predictor vector
y	the response vector
nbasis	- number of basis functions for the spline (default is 25)
sbasis	- number of basis functions for the smoothing of sigma and rho
degree	- degree for splines (default is 2) - allowable possibilities are 0, 1 or 2.
alpha	- controls shape of the GP kernel (default is 2) - 0 < alpha <=2 is possible
xout	- grid on which posterior will be calculated (default is x)

Value

list consisting of xout, the posterior mean, the lower 95% credibility band, the upper 95% credibility band and the INLA object containing the fit

Author(s)

Julian Faraway, [e-mail]

julianfaraway/brinla documentation built on April 6, 2023, 2:33 p.m.