主要内容

このページの翻訳は最新ではありませ。

预测

ガウス过程回帰モデルの予测応答

构文

ypred =预测(gprmdl,xnew)
[ypred,YSD] =预测(gprMdl,Xnew)
[ypred,ysd,yint] =预测(gprmdl,xnew)
[ypred,YSD,yint] =预测(gprMdl,Xnew,名称,值)

说明

ypred=预测(gprMdlXnew的)は,完全またはコンパクトなガウス过程回帰(GPR)モデルgprMdlXnew内の予测子の値について予测した応答ypredを返します。

[ypredYSD] =预测(gprMdlXnew的)は,学习済みGPRモデルから,Xnewに含まれている予测子の値における新しい応答の推定标准偏差も返します。

[ypredYSDyint] =预测(gprMdlXnew的)Xnewの各行に対応する真の応答について95%の予测区间yintも返します。

[ypredYSDyint] =预测(gprMdlXnew名称,价值的)は,1つ以上の名称,价值ペア引数で指定された追加オプションを使用して,予测区间も返します。たとえば,予测区间の信頼度を指定できます。

入力数

开する

ガウス过程回帰モデル。回归(完全)またはCompactRegressionGP(コンパクト)オブジェクトとして指定します。

Fitrgp.でGPRモデルを学习させるために使用する予测子の新しい値。桌子または米行d列の行列として指定します的.mは観测値の数,dは学习データに含まれている予测子変数の数です。

桌子に対してgprMdlを学习させた场合,XnewgprMdlを学习させるために使用した予测子変数がすべて含まれている桌子でなければなりません。

行列に対してgprMdlを学习させた场合,Xnewはd列の数値行列でなければなりません。

データ型:单身的|双倍的|桌子

名前と値のペアの引数

オプションの名称,价值引数のコンマ区切りペアを指定します。名称は引数名で,价值は対応する値です。名称は引用符で囲まなければなりません。name1,value1,...,namen,valuenのように,复数の名前と値のペアの引数を,任意の顺番で指定できます。

予测区间の有意水准。'Α'と0から1の范囲にあるスカラー値から构成されるコンマ区切りのペアとして指定します。

例:'阿尔法',0.01は,99%の予测空间を指定ます。

データ型:单身的|双倍的

出力数

开する

予测した応答の値.N行1列のベクトルとして返されます。

新闻応答値の标准偏差.N行1列のベクトル返されれれ。YSD(ⅰ)には,学习済みGPRモデルから推定した,Xnewの我番目の行に予测予测子の値に対応する新闻応答の推定偏差が格式されます(i = 1,2,...,n)。

Xnewの各行に対応する真の応答値の予测区间.N行2列の行列として返されます。yintの1列目には予测空间のの,2列目には上限が格式され。

开する

标本データを生成します。

N = 10000;RNG(1)再现性的百分比X = linspace(0.5,2.5,N)';Y = SIN(10 * PI * X)./(2. * X)+(X-1)^ 4 + 1.5 *兰特(N,1)。

特性长スケールが各予测子について异なるMatern 3/2カーネル关数を使用し,アクティブセットのサイズを100にして,GPRモデルを近似させます。パラメーター推定には回帰変数サブセット近似法を,予测には完全独立条件法を使用します。

gprMdl = fitrgp(X,Y,'KernelFunction''ardmatern32'......'ActiveSetSize',100,'使用fitmethod''SR''PredictMethod'“同人”);

予测を计算します。

[ypred,〜,YCI] =预测(gprMdl,X);

データを予测および予测区间とともにプロットします。

积(X,Y,'R'。);抓住积(X,ypred);积(X,YCI(:,1),'K--');绘图(x,YCI(:,2),'K--');Xlabel('X');ylabel('是');

图包含一个轴。轴包含型线的4个对象。

标本データを読み込んで桌子に格式します。

加载fisheririsTBL =表(MEAS(:,1),MEAS(:,2),MEAS(:3),MEAS(:,4),物种,......'VariableNames'{'MEAS1''MEAS2''meas3''meas4''物种'});

1番目の测定値を応答て,他のの数をを子として使て,gprモデルを近似さます。

MDL = fitrgp(TBL,'MEAS1');

予测と99%の信息空间を计算ます。

[ypred,〜,YCI] =预测(MDL,TBL,'Α',0.01);

真の応答と予测を予测区间とともにプロットします。

数字();图(mdl.Y,'R'。);抓住;积(ypred);情节(YCI(:,1),'k:');情节(YCI(:,2),'k:');传奇(“真实反应”“GPR预测”......'较低的预测极限'“上预测限制”......'地点''最好的');

图包含一个轴。轴包含型线的4个对象。这些对象表示真响应,GPR的预测,预测下限值,上限值的预测。

标本データを読み込みます。

加载('gprdata.mat');

データには,学习データと検定データが含まれています。学习データには500个,検定データには100个の観测値があります。データには8つの予测子変数があります。このデータは,シミュレーションされたものです。

特性长スケールが各予测子について异なる二乘指数カーネル关数を使用して,GPRモデルを近似させます。学习データの予测子を标准化します。厳密な近似法および予测法を使用します。

gprMdl = fitrgp(Xtrain,ytrain,'基础''持续的'......'使用fitmethod''精确的''PredictMethod''精确的'......'KernelFunction''ardsquaredexponential''标准化',1);

検定データの応答を予测します。

[ytestpred,〜,ytestci] =预测(gprMdl,XTEST);

検定の応答と予测をプロットします。

数字;积(ytest,'r');抓住;积(ytestpred,'B');情节(ytestci(:,1),'k:');情节(ytestci(:,2),'k:');传奇('实际反应'“GPR预测”......'95%“低'95%上”'地点''最好的');抓住离开

图包含一个轴。轴包含型线的4个对象。这些对象表示的实际响应,GPR的预测,95%降低,95%上。

ヒント

  • Fitrgp.で名前と値のペアの引数PredictMethodをををと,GPRモデルを学习方法を选択できます。既定既定予测ます。既定既定の方法は,n≤10000のの合'精确的',それ以外の场合は'bcd'(ブロック座标降下)です.Nは,学习データに含まれている観测値の数です。

  • PredictMethod'bcd'の场合,标准偏差YSDと予测空间yintの计算はサポートされません。

  • gprMdlCompactRegressionGPオブジェクトの场合,PredictMethod'SR'または“同人”に等しいと,标准偏差YSDと予测空间yintを计算できません。PredictMethod'SR'または“同人”に等しい场合にYSDyintを计算するには,完全な回帰オブジェクト(回归)を使用します。

代替方法

resubPredictをを使する,学校データデータ含まれいる観测値学习済み済みモデルモデル予测応答を计算できできでき

拡张机械

R2015Bで导入