LogLikelihood for Gaussian Process regression (function: `fitgpr`) for given set of hyperparameter

6 visualizaciones (últimos 30 días)
I am interested in calculating LogLikelihood using Gaussian Process for given hyperparameters and noise parameter i.e. without optimizing for parameters.
In the following example; [3.5, 6.2, 0.2] are provided as parameters and since 'FitMethod' is 'none' fitgpr will not optimize for parameters
load(fullfile(matlabroot,'examples','stats','gprdata2.mat'))
sigma0 = 0.2;
kparams0 = [3.5, 6.2];
gprMdl2 = fitrgp(x,y,'KernelFunction','squaredexponential',...
'FitMethod','none', 'KernelParameters',kparams0,'Sigma',sigma0);
ypred2 = resubPredict(gprMdl2);
but variable gprMdl2.LofLikelihood = [ ], I am interested in LogLikelihood precisely for parameters [3.5, 6.2, 0.2] not for optimized ones.
Thanks
This question is in continuation to this post.

Respuesta aceptada

Gautam Pendse
Gautam Pendse el 6 de Feb. de 2018
Hi Pankaj,
Loglikelihood is not calculated for 'FitMethod','none'. As a temporary workaround, there is an undocumented internal feature that does this calculation:
gprMdl2.Impl.computeLogLikelihoodExact()
Hope this helps,
Gautam

Más respuestas (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by