如何找到多元正态的水平曲线?

机器算法验证 多元分析
2022-03-21 06:24:41

我试图对多元分析有一些直觉。认为可视化多元正态的 0.95 概率等值线会很好。

我正在寻找有关在 R 中执行此操作的有效方法的建议。

中均匀分布的 x生成一堆向量,然后尝试操纵 cdf 以找到常数(当任意时)。v=(cos(x),sin(x))[0,2π]kF(kv)=.95μ=0Σ

无法使其正常工作,并想知道是否有更好的方法。我对我尝试的 R 包没有太多运气 - mvtnorm - 不确定“等坐标”分位数是什么意思?

1个回答

好吧,自己搔痒。我们想要水平曲线1αyN2(μ,Σ)

分布函数的水平曲线是以原点为中心,半径为这是因为如果我们考虑从那个分布中抽取1αN2(0,I)σ=χ2,1α2xP(xTxχ2,1α2)=1α

使用它来获得应用标准线性变换的情况yx

通过在上生成一系列角度在单位圆上取一些点并应用,调用通过堆叠这些行向量形成的矩阵 _ϕ[0,2π)(cos(ϕ),sin(ϕ))R

Y~=RσΣ1/2+μ

Y~是必要水平曲线上的一组点。

编辑:整理符号以回应@whuber 的评论,希望它现在不那么不连贯。