# Gaussian process regression - how to fit a basis function but not other parameters

15 visualizzazioni (ultimi 30 giorni)
Bill Tubbs il 23 Feb 2023
Commentato: Vitek Stepien il 6 Set 2023
I'm doing Gaussian process regression with the fitrgp function and want to manually control some parameters such as the kernel parameters and the noise parameter, but automatically fit a basis function (e.g. 'linear').
However, it seems that to fit a basis function the FitMethod parameter cannot be 'none':
1. Manual parameter settings (no parameter fitting) but no basis function
gpr_model = fitrgp(x_d, y_d, ...
'KernelParameters', [1.5; 1], ...
'Sigma', 0.1, ...
'FitMethod', 'none' ...
);
2. To fit the basis function I can't have FitMethod set to 'none' ('exact' is the default) but then it ignores my parameter values:
gpr_model = fitrgp(x_d, y_d, ...
'KernelParameters', [1.5; 1], ...
'Sigma', 0.1, ...
'BasisFunction', 'linear' ...
);
3. If I set FitMethod to none, it ignores the basis function:
gpr_model = fitrgp(x_d, y_d, ...
'KernelParameters', [1.5; 1], ...
'Sigma', 0.1, ...
'BasisFunction', 'linear', ...
'FitMethod', 'none' ...
);
(same a Fig. 1).
Is there a way to fit a basis function but not the other parameters?
##### 3 CommentiMostra 1 commento meno recenteNascondi 1 commento meno recente
mono il 5 Set 2023
Facing a similar problem to perform optimizations of a gaussian process with one specific paramater fixed...
Bill Tubbs il 5 Set 2023
Let me know if you find a solution!

Accedi per commentare.

### Risposte (2)

Vitek Stepien il 5 Set 2023
If I understand correctly, you're trying to programatically figure out whoch fitting function would work best for the data given? If so, then the short answer is no, there isn't a good way - I tried doing that in a much more complex setting and there is no automatic way of trying different functions. The workaround would be to write a script that fits a number of functions to the data (I presume you have a few specific functions that you wnat to check), and check which fit has the best "goodness", for example using SSE as the quantifier. Then the script would return the fucntion with best fit (lowest SSE) and you can go from there. Let me know if that helped!
##### 3 CommentiMostra 1 commento meno recenteNascondi 1 commento meno recente
mono il 6 Set 2023
If I understood the question correctly, @Bill Tubbs and myself are trying to perform parameters optimizations but with one or several specific paramaeters fixed.
E.g. with initial theta as theta0 = [1, 2], now only want to optimizae theta0[1] but with theta0[2] = 2 fixed without optimization.
Vitek Stepien il 6 Set 2023
I see, it looked similar to a problem I had but I guess I don't know how to help in your case. Sorry and good luck finding the answer!

Accedi per commentare.

Alan Weiss il 5 Set 2023
I might misunderstand what you are trying to do, but perhaps you could use the OptimizeHyperparameters name-value argument to find the parameter you are trying to optimize. I am not sure that this can give you what you want, because I don't really understand what you are doing, but you can check whether the argument might be applicable.
Good luck,
Alan Weiss
MATLAB mathematical toolbox documentation
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Categorie

Scopri di più su Gaussian Process Regression in Help Center e File Exchange

R2021b

### Community Treasure Hunt

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

Start Hunting!

Translated by