主要内容

帕累托病

带Pareto尾的分段分布

描述

A.帕累托病对象是尾部带有广义帕累托分布(GPD)的分段分布。

A.帕累托病物体由一个或两个gpd在尾部和另一个分布在中心组成。属性可以指定中心的分布类型cdffun的观点帕累托病当你创建一个对象时。有效的值“ecdf”,“内核”,和函数句柄。

帕累托病适合类型的分布cdffun根据观察(x),并找到对应于上尾和下尾累积概率(pl聚氨基甲酸酯分别)。然后,帕累托病将两个GPD安装到下部100*pl观察值的百分比100 * (1-pu)分别占观察值的百分比。如果x在一条尾巴上没有至少两个不同的观察结果吗帕累托病不创建相应的尾段。

使用对象函数边界,,upperparams,lowerparams寻找分布特征。lowerparamsupperparams在尾部返回gpd的参数。边界返回分段分布段之间的边界点,返回包含输入值的分段分布的段片段返回对象中的段数。

使用对象函数cdf,icdf,pdf,随机的评估分布情况。这些函数非常适合copula和其他montecarlo模拟。pdf返回在尾部的GPD密度和在中心的累积分布函数(cdf)的斜率。这些中心的概率密度函数(pdf)值通常不是原始数据潜在密度的良好估计。

创造

使用。创建分段分布对象帕累托病

描述

实例

pd=帕累托(x,pl,聚氨基甲酸酯)返回分段分布对象pd由中心的经验分布和尾部的广义帕累托分布组成。使用较低和较高的尾部累积概率指定尾部的边界pl聚氨基甲酸酯分别地

实例

pd=帕累托(x,pl,聚氨基甲酸酯,cdffun)使用指定中心分布段的类型cdffun

输入参数

全部展开

输入数据,指定为数字向量。

数据类型:

较低的尾部累积概率,指定为范围内的数值标量[0, 1].分位数的pl是下尾观测值的边界。

如果pl0x在下尾翼没有至少两个不同的观测,那么帕累托病将输入数据分为x分成两组,中间和上尾。在这种情况下,拟合的分段分布对象pd由两个部分组成:中心的经验分布和上尾的GPD。

例子:0.1

数据类型:仅有一个的|

上尾累积概率,指定为范围内的数值标量[0, 1].分位数的聚氨基甲酸酯是上尾翼观测的边界。

如果聚氨基甲酸酯1.x在上尾没有至少两个明显的观察结果吗帕累托病将输入数据分为x分为两组,中心和下尾。在这种情况下,拟合的分段分布对象pd由两部分组成:中心的经验分布和下尾的GPD。

例子:0.9

数据类型:仅有一个的|

中心配送段的类型,指定为“ecdf”,“内核”,或函数句柄。

价值 描述
“ecdf”

插值经验cdf。

帕累托病使用值x作为经验cdf垂直步长的中点,并计算x通过线性插值。关于如何找到插值的经验cdf的详细信息,请参见一种分段线性非参数CDF估计

“内核”

插值核平滑估计的cdf。

帕累托病使用ksdensity函数查找范围内的100个点的CDF估计x,并使用线性插值来计算100个点之间的点的估计。

“内核”是否等同于指定一个函数句柄有趣= @ (x) ksdensity (x,“功能”,“提供”);

函数句柄

使用指定函数的插值估计。

帕累托病使用窗体函数的句柄[p, xi] =乐趣(x)它接受输入数据向量x并返回一个向量P的CDF值和一个向量西评估点。值西必须是有序和不同的,但不必须等于值在x.的帕累托病函数为中的值之间的点计算cdf估计值西通过线性插值。

帕累托病使用cdffun计算对应于pl聚氨基甲酸酯

例子:“内核”

性质

全部展开

此属性是只读的。

段的数目,包括中心段和尾部段paretotail对象,指定为标量。数字段如果对象中的尾段数分别为2、1或0,则为3,2或1。

数据类型:

此属性是只读的。

下尾GPD参数,适用于中的下极值观测值x,指定为数字向量。第一个值是图形参数,第二个值是GPD的比例参数。

下尾GPD的位置参数等于的分位数pl. 使用边界函数返回位置参数。例如,运行[p,q]=边界(pd),在那里pd是一个帕累托病对象。问题(1)是位置参数。

数据类型:仅有一个的|

此属性是只读的。

上尾GPD参数,符合上极值观测值x,指定为数字向量。第一个值是图形参数,第二个值是GPD的比例参数。

上尾GPD的位置参数等于的分位数聚氨基甲酸酯. 使用边界函数返回位置参数。例如,运行[p,q]=边界(pd),在那里pd是一个帕累托病对象。问题(2)是位置参数。

数据类型:仅有一个的|

对象的功能

边界 分段分布边界
cdf 累积分布函数
icdf 逆累积分布函数
lowerparams 下帕累托尾参数
片段 分段分布的段数
pdf 概率密度函数
随机的 随机数
包含输入值的分段分布段
upperparams 上帕累托尾参数

例子

全部折叠

生成一个样本数据集,并将带有帕累托尾的分段分布拟合到数据中。通过使用指定中心的经验分布帕累托病使用默认设置。

生成一个样本数据集,其中包含来自T3自由度分布。

rng (“默认”);%的再现性t = trnd (3100 1);

创建一个帕累托病通过将分段分布拟合到T.使用下尾和上尾累积概率指定尾部的边界,以便拟合对象由数据集中80%的经验分布和数据集中10%的gpd组成。

pd = paretotails (t, 0.1, 0.9)
pd =分段分布与3段负无穷到< x < -1.84875 (0 < p < 0.1):低尾巴,加仑日(0.183032,1.00347)-1.84875 < x < 2.07662 (0.1 < p < 0.9):插值经验提供2.07662 < x <正(0.9 < p < 1):上尾巴,加仑日(0.333239,1.19705)

对象显示的每一行显示了每个分段的汇总,包括GPD参数(形状和比例参数)以及分位数和累积概率中的边界值。使用对象函数边界,lowerparams,upperparams返回这些值。

你可以使用片段函数返回段数和函数返回包含输入值的段。

您还可以使用分布函数cdf,icdf,pdf,随机的评估分布并生成随机样本。

画出cdfT产品的分布和cdf帕累托病物体在同一图形上。

x = linspace (5,5);情节(x, tcdf (x, 3),“r——”)举行情节(x, cdf (pd, x)“b -”)

求线段之间的边界点帕累托病对象的使用边界,并在图上标出点。

[p,q]=边界(pd);图(q,p,“波”)传奇(“t分布”,“帕累托尾对象”,边界点的,“位置”,“最佳”)举行

图中包含一个轴对象。axes对象包含3个line类型的对象。这些对象表示t分布、帕累托尾对象、边界点。

生成一个样本数据集,并对数据拟合一个带帕累托尾的分段分布。使用合适的中心部分帕累托病使用函数句柄。

生成包含20%离群值的样本数据集。

rng (“默认”);%的再现性left_tail = -exprnd (1100 1);right_tail = exprnd (5100 1);中心= randn (800 1);x = [left_tail;中心;right_tail];

定义一个函数句柄使用ksdensity指定带宽的非默认值。

myfun1=@(x)ksdensity(x,“带宽”、1。“功能”,“提供”);

创建一个帕累托病通过拟合一个分段分布与指定的核平滑估计x.使用上尾和下尾累积概率指定尾部的边界,以便拟合对象由数据集中80%的核估计量和数据集中10%的gpd组成。

pd1 = paretotails (x, 0.1, 0.9, myfun1)
pd1=3段分段分布-Inf
                   

您还可以为中心段使用参数化分布。定义一个适合数据正态分布的函数并返回cdf值,并在创建中心段时传递函数句柄帕累托病对象。

pd2 = paretotails (x, 0.1, 0.9, @myfun2)
pd2 =分段分布与3段负无穷到< x < -2.70875 (0 < p < 0.1):低尾巴,加仑日(-0.358104,0.831855)-2.70875 < x < 3.52195 (0.1 < p < 0.9):功能:myfun2 < x <正无穷(0.9 < p < 3.52195 1):上尾巴,加仑日(-0.0661815,5.04694)
函数[p,xi] = myfun2(x) pd = fitdist(x)“正常”);ξ= linspace (min (x)马克斯(x),长度(x) * 2);p = cdf (pd, xi);终止
R2007a中引入