多元正态分布的尾概率

机器算法验证 r 可能性 分布 数理统计 多元分析
2022-04-04 15:57:16

我想计算不同维度的标准化多元正态分布的尾部概率。例如,在二元法线的情况下,我需要计算下面粘贴的图形上的灰色区域(在这种特殊情况下,该区域由定义。)x>2

图1

在半径为的球形轮廓之外的概率应计算为:xα

Pr(f(x)f(0)α)=P(xTx2ln(α))=02ln(α)fχd2(xTx)=Fχd2(2ln(α))
其中F是累积分布函数。

我发现自己迷失在这个等式中。我恳请计算示例如何计算所需的概率(使用任何编程语言)。

解决方案

我自己找到了解决方案:

pchisq(q = alpha^2, df = d)

alpha代表半径,d代表尺寸。

1个回答

在评论中回答,复制到这里:

R 中的 mvtnorm 包可以做到这一点。查看https://cran.r-project.org/web/packages/mvtnorm/vignettes/MVT_Rnews.pdf 以获取许多示例。