ドキュメンテーションヘルプ センター
最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。
確率分布のプロファイル尤度関数
[ll,param] = proflik(pd,pnum)
[ll,param] = proflik(pd,pnum,'Display',display)
[ll,param] = proflik(pd,pnum,setparam)
[ll,param] = proflik(pd,pnum,setparam,'Display',display)
[ll,param,other] = proflik(___)
例
[ll,param] = proflik(pd,pnum) は、対数尤度値のベクトル ll と、pnum で示される位置のパラメーターに対応するパラメーター値のベクトル param を返します。
ll
param
pd
pnum
[ll,param] = proflik(pd,pnum,'Display',display) は、対数尤度値と対応するパラメーター値を返し、対数尤度の近似にプロファイル尤度を重ね合わせてプロットします。
display
[ll,param] = proflik(pd,pnum,setparam) は、対数尤度値と、setparam で指定された対応するパラメーター値を返します。
setparam
[ll,param] = proflik(pd,pnum,setparam,'Display',display) は、対数尤度値と、setparam で指定された対応するパラメーター値を返し、対数尤度の近似にプロファイル尤度を重ね合わせてプロットします。
[ll,param,other] = proflik(___) は、前の構文からのいずれかの入力引数を使って、尤度を最大化するその他のパラメーターの値を含む行列 other も返します。
other
すべて折りたたむ
標本データを読み込みます。ワイブル分布をガロンあたりの走行マイル数 (MPG) データに近似して、確率分布オブジェクトを作成します。
MPG
load carsmall pd = fitdist(MPG,'Weibull')
pd = WeibullDistribution Weibull distribution A = 26.5079 [24.8333, 28.2954] B = 3.27193 [2.79441, 3.83104]
分布のパラメーター名を表示します。
pd.ParameterNames
ans = 1x2 cell {'A'} {'B'}
ワイブル分布の場合、A は位置 1 にあり、B は位置 2 にあります。
A
B
位置 pnum = 2 にある B のプロファイル尤度を計算します。
pnum = 2
[ll,param] = proflik(pd,2);
B の推定値に対する対数尤度の値を表示します。
[ll',param']
ans = 21×2 -329.9688 2.7132 -329.4312 2.7748 -328.9645 2.8365 -328.5661 2.8981 -328.2340 2.9597 -327.9658 3.0213 -327.7596 3.0830 -327.6135 3.1446 -327.5256 3.2062 -327.4943 3.2678 ⋮
これらの結果は、B の推定値 3.2678 および 3.3295 の間でプロファイル対数尤度が最大化されることを示しています。これは、対数尤度値 -327.4943 および -327.5178 に対応します。以前の近似から、B の MLE は 3.27193 で、この MLE は予想されたとおりこの区間に含まれています。
標本データを読み込みます。一般化極値分布をガロンあたりの走行マイル数 (MPG) データに近似して、確率分布オブジェクトを作成します。
load carsmall pd = fitdist(MPG,'GeneralizedExtremeValue')
pd = GeneralizedExtremeValueDistribution Generalized Extreme Value distribution k = -0.207765 [-0.381674, -0.0338564] sigma = 7.49674 [6.31755, 8.89603] mu = 20.6233 [18.8859, 22.3606]
ans = 1x3 cell {'k'} {'sigma'} {'mu'}
一般化極値分布の場合、k は位置 1 にあり、sigma は位置 2 にあり、mu は位置 3 にあります。
k
sigma
mu
位置 pnum = 3 にある mu のプロファイル尤度を計算します。計算をパラメーター値 20 ~22 に制限し、プロットを表示します。
pnum = 3
[ll,param,other] = proflik(pd,3,20:.1:22,'display','on');
プロットには、対数尤度を最大化するパラメーター mu の推定値が表示されます。
mu の推定値の対数尤度値と、対応する対数尤度を最大化する他の分布パラメーターの値を表示します。
[ll',param',other]
ans = 21×4 -327.5706 20.0000 -0.1803 7.4087 -327.4971 20.1000 -0.1846 7.4218 -327.4364 20.2000 -0.1890 7.4354 -327.3887 20.3000 -0.1934 7.4493 -327.3538 20.4000 -0.1978 7.4636 -327.3317 20.5000 -0.2023 7.4783 -327.3223 20.6000 -0.2067 7.4932 -327.3257 20.7000 -0.2112 7.5084 -327.3418 20.8000 -0.2156 7.5240 -327.3706 20.9000 -0.2201 7.5399 ⋮
1 列目には、2 列目の mu の推定値に対応する対数尤度値が含まれています。対数尤度は、パラメーター値 20.6000 および 20.7000 の間で最大化されています。これは、対数尤度値 -327.3223 および -327.3257 に対応します。3 列目には、対応する mu の対数尤度を最大化する k の値が含まれています。4 列目には、対応する mu の対数尤度を最大化する sigma の値が含まれています。
確率分布。次のいずれかを使用して作成した確率分布オブジェクトを指定します。
fitdist
プロファイル尤度を計算するパラメーター番号。パラメーター名ベクトルの必要なパラメーターの位置に対応する正の整数値として指定します。たとえば、ワイブル分布にパラメーター名ベクトル {'A','B'} がある場合は、pnum を 2 として指定して、B のプロファイル尤度を計算します。
{'A','B'}
2
データ型: single | double
single
double
パラメーター値の制限。スカラー値またはスカラー値のベクトルとして指定します。setparam を指定しない場合、proflik は、確率分布 pd の既定の信頼区間メソッドに基づいて出力ベクトル param の値を選択します。パラメーターが制限された値のみ取ることができる場合や、信頼区間がその制限に違反している場合は、setparam を使用して有効な値を指定できます。
proflik
例: [3,3.5,4]
[3,3.5,4]
'off'
'on'
表示切り替え。'on' または 'off' として指定します。'on' を指定すると、対数尤度の Wald 近似に正確な対数尤度のプロファイルが重ね合わされて表示されます。'off' を指定すると、表示を省略できます。Wald 近似は、位置 pnum にあるパラメーターの関数またはその対数として、推定されたパラメーター値に対するテイラー級数展開に基づきます。水平方向に点線で表される曲線が交差する部分は、95% 信頼区間の端点としてマークします。
ベクトルとして返される対数尤度値。対数尤度は、位置 pnum のパラメーターが、残りのパラメーターに対して最大化される param の値に設定される尤度の値です。
ll の対数尤度値に対応するパラメーター値。ベクトルとして返されます。setparam を使用してパラメーター値を指定すると、param は setparam と等しくなります。
尤度を最大化する他のパラメーター値。行列として返されます。other の各行には、位置 pnum 内のパラメーターを除くすべてのパラメーターの値が含まれます。
Distribution Fitter | fitdist | mle | negloglik | paramci
mle
negloglik
paramci
この例の変更されたバージョンがシステム上にあります。代わりにこのバージョンを開きますか?
次の MATLAB コマンドに対応するリンクがクリックされました。
コマンドを MATLAB コマンド ウィンドウに入力して実行してください。Web ブラウザーは MATLAB コマンドをサポートしていません。
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list:
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Contact your local office