负の対数关数
统计和机器学习工具箱™の分布サポートれいるの対数尤度关数はすべてExply
のように像
でます关数は,分布パラメーター族を表します。引数はは,,データデータデータがが従う分布分布族のの特定特定ののメンバーをを指定するする一连一连一连の値値値値の対数を返します。
最(mle)ををには,の尤度尤度をを化化问题问题の目的关数关数としてとして,,,,®关数fminsearch
または优化工具箱™や全局优化工具箱の关数し解くことができます。のにより探索アルゴリズムの选択选択とと,,アルゴリズムの実行に対する低低水准水准ののの制御制御mle
と,NORMFIT
やGamfit
などの合身
で终わる关数,あらかじめされいるアルゴリズム使用し,そのオプションは关数Statset
によって设定れるもの限られます。
分布の族は确率密度(PDF)f(x |θ)を使用指定,,,,,,は确率,,,,,θは分布パラメーターをそれぞれますますますますf(x |θ)をのx xに対するθの关数见る见る,关数关数关数关数f(x |θ)は结果结果结果结果パラメーターθのの尤度。一分布に従う独立独立ししセットセットセットセットセットセットセットセットセットセットセットパラメーターパラメーターパラメーターθθθのののの尤度
x xのについて,はなななななにおけるl(θ)ををををを化ししますますます。。。数値数値数値数値数値アルゴリズムアルゴリズムアルゴリズムアルゴリズムアルゴリズムアルゴリズムアルゴリズムアルゴリズムでででででは关数关数关数关数关数关数。は値になると思われるのをを対数対数和和ににに変换変换しししししますますます。でででで,,,计算计算计算计算においてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいてにおいて度关数はこののの“负”の値。これは最适アルゴリズムではに最大値値ではなく最最小値をを求める。
负の尤度を使用したたたの
この例で,关数赌注
およびfminsearch
Mleをを求める方法を示し。。。。
关数Gamrnd
を使用て,特定のガンマ分布から无作为を生成し。。
RNG默认;%可再现性a = [1,2];x = gamrnd(a(1),a(2),1e3,1);
关数赌注
を使用て,与えたたX
に対する一种
の近傍尤度を可视化し。。
网格= 50;delta = 0.5;a1 = linspace(a(1)-delta,a(1)+delta,网格);a2 = linspace(a(2)-delta,a(2)+delta,网格);logl = zeros(mesh);%preallocate内存为了i = 1:网眼为了j = 1:网格logl(i,j)= gamlike([a1(i),a2(j)],x);结尾结尾[a1,a2] = meshgrid(a1,a2);Surfc(A1,A2,Logl)
关数fminsearch
を使用,尤度の最小値求め。。。
ll = @(u)gamlike([u(1),u(2)],x);给定x的可能性%mles = fminsearch(LL,[1,2])
mles =1×20.9980 2.0172
关数Gamfit
によって返される値とMLE
を比较し。
ahat = gamfit(x)
ahat =1×20.9980 2.0172
MLE
と一顶帽子
の间各パラメーターのはは1E-4
より小さくています。
mleをプロットにし。。
抓住在Plot3(MLE(1),MLE(2),LL(MLE),,,'ro',,,,“标记”,5,“ MarkerfaceColor”,,,,'r')查看([-60 40])%旋转以显示最低
参考
nenloglik
|Statset
|fminsearch
|冲浪