mvtcdf

多元Ť累积分布函数

句法

Y = mvtcdf(X,C,DF)
Y = mvtcdf(XL,戌,C,DF)
[Y,ERR] = mvtcdf(...)
[...] = mvntdf(...,选项)

描述

Y = mvtcdf(X,C,DF)返回多元的累积概率Ť用相关性的参数分配C和自由度DF中,在每行计算X。的行ñ-通过-d矩阵X对应于观察或点,和列对应的变量或坐标。ÿ是一个ñ-通过-1向量。

C是对称,正定的,d-通过-d基质,通常是相关矩阵。如果它的对角线元素不为1,mvtcdfC以关联形式。mvtcdf不重新调整XDF是标量,或者与载体ñ元素。

多元Ť在累积概率X被定义为概率的随机矢量Ť,分布多元Ť,将落入半无限矩形与由定义的上限内XPR {T(1)X(1),T(2)X(2),... T(dX(d)}

Y = mvtcdf(XL,戌,C,DF)返回多元Ť评价在矩形累积概率与由下式定义的上限和下限XL, 分别。

[Y,ERR] = mvtcdf(...)返回错误的估计ÿ。对于双变量和三变量的分布,mvtcdf采用自适应正交上的变换Ť密度,基于由GENZ开发的,如在参考文献中描述的方法。对于这些情况下,默认绝对误差容限1E-8。对于四个或更多的维度,mvtcdf使用基于由GENZ和Bretz发展,在参考文献中所述方法拟蒙特卡罗积分算法。对于这些情况下,默认绝对误差容限1E-4

[...] = mvntdf(...,选项)指定的控制参数用来计算数值积分ÿ。这种说法可以通过将呼叫创建statset。的选择statset参数如下:

  • 'TolFun'- 最大绝对误差容限。默认值是1E-8什么时候d<4,或1E-4什么时候d≥4。

  • 'MaxFunEvals'- 允许积评估的最大数量时,d≥4.默认是1E7'MaxFunEvals'被忽略的时候d<4。

  • '显示'- 显示输出的电平。选择是“关”(缺省值),'ITER''最后''显示'被忽略的时候d<4。

例子

全部收缩

计算多元的CDFŤ用相关性的参数分配C = [1 0.4;0.4 1]和2个自由度。

C = [1 0.4;0.4 1];DF = 2;[X1,X2] = meshgrid(linspace(-2,2,25) 'linspace(-2,2,25)');X = [X1(:) X2(:)];P = mvtcdf(X,C,DF);

画出CDF。

数字;冲浪(X1,X2,重塑(P,25,25));

参考

[1] GENZ,A.“矩形二元和三变量正常和叔概率的数值计算。”统计与计算。卷。14,第3号,2004年,第251-260。

[2] GENZ,A.,和F. Bretz。“多元t概率的数值计算与应用多个对比的功率计算。”杂志统计计算与仿真。卷。63,1999年,第361-378。

[3] GENZ,A.,和F. Bretz。“方法的多元t概率的计算的对比。”杂志计算和图形统计。卷。11,第4号,2002年,第950-971。

也可以看看

|

主题

介绍了在R2006a