proflik

概率分布的轮廓似然函数

描述

例子

参数) = proflik (pdpnum返回一个向量对数似然值和向量参数的对应参数值,为所指示位置的参数pnum

参数) = proflik (pdpnum“显示”,显示返回对数似然值和相应的参数值,并绘制覆盖在对数似然近似上的剖面似然值。

参数) = proflik (pdpnumsetparam返回指定的对数似然值和相应的参数值setparam

例子

参数) = proflik (pdpnumsetparam“显示”,显示返回指定的对数似然值和相应的参数值setparam,并在对数似然的近似上绘制轮廓似然重叠图。

例子

参数其他) = proflik (___也返回一个矩阵其他使用前面语法中的任何输入参数,包含使可能性最大化的其他参数的值。

例子

全部折叠

加载示例数据。通过将威布尔分布拟合到每加仑英里数(英里/加仑)数据。

负载carsmallpd = fitdist(英里/加仑,“威布尔”
Weibull分布A = 26.5079 [24.8333, 28.2954] B = 3.27193 [2.79441, 3.83104]

查看分布的参数名称。

pd。ParameterNames
ans =1 x2单元格{A} {B}

对于威布尔分布,一个在1号位置,然后呢B在2号位置。

计算轮廓的可能性B,已经就位pnum = 2

(噢,param) = proflik (pd, 2);

显示的估计值的对数似然值B

(会,参数)
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。由早期的拟合可知,的MLEB是3.27193,与预期的一样在这个区间内。

加载示例数据。通过拟合每加仑英里的广义极值分布来创建概率分布对象(英里/加仑)数据。

负载carsmallpd = fitdist(英里/加仑,“GeneralizedExtremeValue”
pd =广义dextremevaluedistribution广义极值分布k = -0.207765 [-0.381674, -0.0338564] sigma = 7.49674 [6.31755, 8.89603] mu = 20.6233 [18.8859, 22.3606]

查看分布的参数名称。

pd。ParameterNames
ans =1 x3单元格{“k”}{“σ”}{“亩”}

对于广义极值分布,k在1号位置,σ在2号位置,然后呢μ在3号位置。

计算轮廓的可能性μ,已经就位pnum = 3.将计算限制为从20到22的参数值,并显示绘图。

(噢,参数,其他)= proflik (pd 3 20: .1:22,“显示”“上”);

图中显示了参数的估计值μ这使对数似然最大化。

显示的估计值的对数似然值μ,以及使相应的对数似然最大化的其他分布参数的值。

(会,参数,其他)
ans =21日×4-327.5706 20.0000 -0.1803 7.4087 -327.4971 20.1000 -0.1846 7.4218 -327.4364 20.2000 - 1890 7.4354 -327.3887 20.3000 - 1934 7.4493 -327.3538 20.4000 - 1978 7.4636 -327.3317 20.5000 -0.2023 7.4783 -327.3223 20.6000 - 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⋮

第一列包含对应于的估计的对数似然值μ在第二列。参数值20.6000 ~ 20.7000之间的对数似然值最大,对应的对数似然值为-327.3223和-327.3257。第三列包含值k使对应的对数可能性最大化μ.第四列包含值σ使对应的对数可能性最大化μ

输入参数

全部折叠

概率分布,指定为使用下列方法之一创建的概率分布对象。

功能或应用程序 描述
fitdist 拟合一个概率分布对象到样本数据。
分布更健康 使用交互式distribution Fitter应用程序拟合样本数据的概率分布,并将拟合对象导出到工作空间。

用于计算轮廓可能性的参数编号,指定为一个正整数值,对应于所需参数在参数名称向量中的位置。例如,Weibull分布有一个参数名称向量{' A ', ' B '},所以指定pnum作为2来计算轮廓的可能性B

数据类型:|

参数值限制,指定为标量值或此类值的向量。如果没有指定setparamproflik选择输出向量的值参数基于默认置信区间法的概率分布pd.如果参数只能取受限制的值,并且置信区间违反了该限制,则可以使用setparam指定有效值。

例子:(3, 3.5, 4)

显示开关,指定为其中之一“上”“关闭”.指定“上”以显示叠加在对数似然的Wald近似上的准确对数似然的轮廓。指定“关闭”省略显示。Wald近似是基于估计参数值周围的泰勒级数展开,作为位置参数的函数pnum或其对数。曲线与水平虚线的交点是95%置信区间的端点。

输出参数

全部折叠

Loglikelihood值,作为向量返回。对数似然是参数处于位置时的似然值pnum设置为中的值参数,在其余参数上最大化。

中的对数似然值对应的参数值,返回为向量。如果使用setparam,然后参数等于setparam

其他最大化可能性的参数值,以矩阵的形式返回。每一行的其他包含除位置参数外的所有参数的值pnum

介绍了R2013a