计算矩阵K在测试输入训练后与fitrgp高斯过程
18视图(30天)
显示旧的评论
如果我训练有素的医生使用训练数据D = {X, y} fitrgp和我获得gprMdl:
gprMdl = fitrgp(数据。X data.Y“KernelFunction”,“squaredexponential”,…
“BasisFunction”,“没有”,“详细”, 1“FitMethod”,“准确”)
[ystar, ysd yint] =预测(gprMdl Xstar)
我怎样才能获得矩阵K (Xstar Xstar) ?我找不到RegressionGP的子函数,计算矩阵K。
谢谢
0评论
接受的答案
Gautam Pendse
2015年11月9日
嗨Umberto,
有一个非法的计算方式你想要什么。这是一个例子:
rng (0,“旋风”);
N = 100;
x = linspace (-10 10 N) ';
y = 1 + x * 5依照+ sin (x)。/ x + 0.2 * randn (N, 1);
探地雷达= fitrgp (x, y,“FitMethod”,“准确”,“PredictMethod”,“准确”);
kfcn = gpr.Impl.Kernel.makeKernelAsFunctionOfXNXM (gpr.Impl.ThetaHat)
K = kfcn (x (1:5:) (1:7,:))
K (i, j)核函数评估(我:)和x (j:)。例如,
K (3、6)
:kfcn (x (3) x (6,:))
我有兴趣了解你为什么要计算K。
希望有所帮助,
Gautam
7评论
英镑贝尔德
2021年7月25日
编辑:英镑贝尔德
2021年7月25日
我想我意识到我要求的是predictExactWithCov函数(见
抽样从探地雷达模型的后验分布fitrgp ()
)。我想我误解了这个函数的“精确”部分。我把它“准确”是指不使用稀疏方法输入数据而不是假设没有噪音(原来我以为是后者)。