主要内容

gamfit

gamma 参数估计

语法

phat = gamfit(data)
[phat,pci] = gamfit(data)
[phat,pci] = gamfit(data,alpha)
[...] = gamfit(data,alpha,censoring,freq,options)

说明

phat = gamfit(data) 返回到给定向量 data 中数据的 gamma 分布参数的最大似然估计 (MLE)。

[phat,pci] = gamfit(data) 返回 MLE 和 95% 置信区间。pci 的第一行是置信区间的下界;最后一行是上界。

[phat,pci] = gamfit(data,alpha) 返回 100(1 - alpha)% 的置信区间。例如,alpha = 0.01 返回 99% 的置信区间。

[...] = gamfit(data,alpha,censoring) 接受与 data 大小相同的布尔向量,对于右删失观测值,该向量为 1,对于精确观测的观测值,该向量为 0。

[...] = gamfit(data,alpha,censoring,freq) 接受与 data 大小相同的频率向量。freq 通常包含 data 中对应元素的整数频率,但可以包含任何非负值。

[...] = gamfit(data,alpha,censoring,freq,options) 接受结构体 options,该结构体指定函数用于计算最大似然估计值的迭代算法的控制参数。gamma 拟合函数接受 options 结构体,该结构体可以使用函数 statset 创建。输入 statset('gamfit') 以查看 gamfitoptions 结构体中接受的参数的名称和默认值。

示例

对 gamma 分布进行从指定 gamma 分布生成的随机数据拟合:

a = 2; b = 4;
data = gamrnd(a,b,100,1);

[p,ci] = gamfit(data)
p =
  2.1990  3.7426
ci =
  1.6840  2.8298
  2.7141  4.6554

参考资料

[1] Hahn, Gerald J., and S. S. Shapiro. Statistical Models in Engineering. Hoboken, NJ: John Wiley & Sons, Inc., 1994, p. 88.

扩展功能

全部展开

版本历史记录

在 R2006a 之前推出