根据最大似然来拟合数据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中央文件交换。检索.