对截尾数据拟合分布

版本1.0.0.0 (30.1 KB 列奥尼达斯斑替
当存在左和/或右和/或间隔审查时,适合分布到数据

1.5 k下载

更新2012年9月18日

查看许可协议

根据最大似然来拟合数据x的分布。数据可以向左和/或向右和/或间隔删除。

支持以下发行版:万博1manbetx

1.正常(normfitc)
2.对数正态分布(lognfitc)
3.物流(logistfitc)
4.Log-logistic (loglogistfitc)
5.极值(evfitc)
6.威布尔(wblfitc)
7.指数(expfitc)
8.γ(gamfitc)
9.瑞利(raylfitc)

输入参数:
X:数据的两列矩阵。例如,(对于实数线支持的分布):万博1manbetx

-Inf 5.0000——>被删除
7万Inf- >右删
5.0000 9.0000—>间隔被删除
8.0000 8.0000—>完全符合

也就是说,如果所有的数据都被精确观察到,那么这两列必须是相同的。参见“examples”文件中的说明。

如果要选择一个支持度为正的分布,那么上面的例子是:万博1manbetx

0 5万->被删除
7万Inf- >右删
5.0000 9.0000—>间隔被删除
8.0000 8.0000—>完全符合

可选输入参数:(这些参数可以是根本无法到达的,
或设置为[]以继续进行下一个可选输入参数)

最小化:可以设置为1或2或3。即1:fminsearchbnd例程用于最小化,2:fmincon与Largescale设置为off并使用' spq '算法。3.fmincon的大尺度设置为off,采用“内点”算法。如果设置为[]或根本没有达到,则默认选择最小值1(没有特殊原因)。如果出现收敛问题,请务必首先探索所有三个最小化器。然后您可以在下一个输入参数中研究其他初始值。

init:由用户定义的初始值的行向量。如果您不想定义它并继续查看下一个可选输入参数,请将其设置为[]。如果用户没有提供初始值,则根据精确观测到的数据或基于内置的MATLAB函数进行右截尾。如果只有区间数据可用,那么中点被认为是精确观测到的,并使用内建函数来获得初始值。

options:将在fminsearchbnd或fmincon的优化集中使用的选项,这取决于前面的输入参数选择了哪个最小化器。也就是说,如果最小化设置为1,那么选项将引用fminsearchbnd。如果最小化器设置为2或3,那么选项将引用fmincon。此参数中设置的选项将替换最小化输入参数中提到的默认选项。如果这个输入参数被设置为[]或根本没有达到,则会发生默认值。参见MATLAB文档,了解fmincon或fminsearch的“选项”。

输出参数:

Pars:估计的参数

协方差:估计参数的方差协方差矩阵。这是基于John D’erricos的“自适应鲁棒数值微分”工具在所有发行版中进行评估的。

SE:估计参数的标准误差

Gval:对数似然值

Exitflag:取决于所选的最小化器。在只有精确观测数据和/或正确截尾数据可用的情况下,exitflag设置为NaN,因为估计是基于内置的拟合例程完成的。

这个例程的目的是为MATLAB用户提供一个有用的工具来拟合一些分布(其中许多主要用于生存分析),当左和/或右和/或区间截尾可用时。

这个例程使用John D'Erricos的fminsearchbnd或内置的“fmincon”来实现最小化。

计算协方差矩阵
在数值上使用John D'Erricos的“自适应鲁棒数值微分”工具。可以在这里找到:

//www.tianjin-qmedu.com/matlabcentral/fileexchange/8277

//www.tianjin-qmedu.com/matlabcentral/fileexchange/13490-adaptive-robust-numerical-differentiation

参见“示例”文件中的说明

引用作为

列奥尼达斯·班蒂斯(2022年)。对截尾数据拟合分布(//www.tianjin-qmedu.com/matlabcentral/fileexchange/38226-fit-distributions-to-censored-data), MATLAB中央文件交换。检索

MATLAB版本兼容性
使用R2010a创建
与任何版本兼容
平台的兼容性
窗户 macOS Linux

社区寻宝

在MATLAB中央找到宝藏,并发现社区可以如何帮助你!

开始狩猎!