主要内容

贝尔替换

β负面的日志可能性

句法

nlogl = betalike(params,data)
[nlogl,avar] = betalike(params,data)

描述

nlogl = betalike(params,data)返回Beta参数的Beta Log-似然函数的负数一种B.在向量中指定参数和列向量中指定的观察结果数据

元素数据必须位于打开间隔(0,1)中,其中定义了测试版。然而,有时也需要将Beta分布符合包含精确零或那些的数据。对于这样的数据,β似然函数是无界的,并且不可能进行标准的最大似然估计。在这种情况下,贝尔替换计算通过处理它们的修改后的可能性,因为它们是已被留下的值sqrt(realmin)或在1-右审查EPS./ 2分别。

[nlogl,avar] = betalike(params,data)也退货avar.,这是如果值的值估计参数估计的渐近方差 - 协方差矩阵参数是最大可能性估计。avar.是Fisher的信息矩阵的倒数。对角线元素avar.是各自参数的渐近差异。

贝尔替换是一个效用函数,可用于最大似然估计测试版分布。可能性假设数据样本中的所有元素都是相互独立的。自从贝尔替换返回负测试曲目日志似然函数,最小化贝尔替换使用fminsearch.与最大化可能性相同。

例子

这个例子继续了Betafit.示例,其计算某些随机生成的Beta分布式数据的Beta参数的估计。

r = betarnd(4,3,100,1);[NLOGL,AVAR] = BETARIKE(Betafit(R),R)NLOGL = -27.5996 AVAR = 0.2783 0.1316 0.1316 0.0867

扩展能力

C / C ++代码生成
使用MATLAB®Coder™生成C和C ++代码。

版本历史记录

在R2006A之前介绍