Skip to content

wdwang/gptp

Repository files navigation

gptp_multi_output

This package is used to implement multivariate Gaussian process regression (MV-GPR) and multivariate Student-t process regression (MV-TPR). NOTE: This code is proof-of-concept, not optimized for speed.

This package is based on GPML Matlab Code http://www.gaussianprocess.org/gpml/code/matlab/doc/

The main function is gptp_general.m. There are four useful sub-functions, gp_solve_gpml.m, tp_solve_gpml.m, mvgp_solve_gpml.m, mvtp_solve_gpml.m These four functions are used to solve GPR, TPR, MV-GPR and MV-TPR, respectively.

covSEiso.m, covSEard.m and sq_dist.m are collected from GPML Matlab Code.

nu_init.m, Omega_init.m, SE_init.m are used to initialise parameter nu (in TPR and MV-TPR), Omega(in MV-GPR and MV-TPR), and hyperparameters in the specific kernel. Actually, the initial hyperparameters play an important role in the performance of regression. If possible, you can write your own initialisation function.

covSEiso and covSEard are used as default covariance function. More covariance functions can be selected in the GPML Code toolbox.

Notice: if the covariance function is changed, do not forget to replace SE_init.m by the corresponding kernel initialation function.

MultiGamma.m and vec2mat_diag.m are two small functions, which are used in the mvgp_solve_gpml.m and mvtp_solve_gpml.m.

%------------------- 2018/01/19

File added: SimulatedExample.m

Add a simple example for multi-output prediction (MV-GP and MV-TP) compared with independent prediction using GP and TP

Update optimset setting for Matlab 2017 and later in gptp_general.m

%-------------------- 2017/12/04

Add gptp_sample.m file This file is to generate a sample from GP or TP with specificed row and column covariance and zero mean function.

Add mv_gptp_sample.m file This file is to generate a sample from MV-GP or MV-TP with specificed row and column covariance and zero mean function.

Reference:

[1] Chen, Zexun, and Bo Wang. "How priors of initial hyperparameters affect Gaussian process regression models." Neurocomputing 275 (2018): 1702-1710.

[2] Chen, Zexun, Bo Wang, and Alexander N. Gorban. "Multivariate Gaussian and Student $-t $ Process Regression for Multi-output Prediction." arXiv preprint arXiv:1703.04455 (2017).

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages