N正态id的乘积的近似分布?特例μ≈0


12

给定 iid和,寻找:X ÑÑμ Xσ 2 Xμ X0N30XnN(μX,σX2)μX0

  1. 精确封闭形式分布近似值 YN=1NXn
  2. 相同乘积的渐近(指数?)逼近

这是一个特殊情况,是一个更一般的问题μX0


1.您是否有有关和?(例如,如果所有都很好。)(2)渐近法线逼近将是可怕的,因为渐近线看起来不会是法线。σ X μ X / σ X » 0μXσXμX/σX0Y
whuber

我刚刚玩了一下。如果您有兴趣,则可以为iid的随机变量的乘积获得精确的封闭形式解。非零的情况会使事情复杂得多。Ñ 0 σ 2μnN(0,σ2)μ
wolfies'2015-4-3

@whuber (1)用不同的和做一些蒙特卡洛之后,我发现分布在和表现得很好;现在,我想为和找到一个不错的表达式,类似于有一些不错的近似值。我通过泰勒展开建立了一些近似值,但是它们的行为不当。 (2)好吧,肯定会“看起来”像是具有卡方平方的法线的总和,因此,如果近似值“证明”了这一点,则可以简化为法线。σ ˚F Ñ > 30 | μ XμσFN>30|μX|10σXμFσFχ2FF
Andrei Pozolotin

3
当,将通过对数正态分布很好地近似(如Barry-Esseen定理在所示)。μX10σXYlog(X)
whuber

@whuber对Barry-Esseen的直接应用给出,的确不错,但是它一些结构:应该为负,应该取决于等,也许有更好的应用方法?μ˚Fσ˚FFN0+1NZμFσFα
Andrei Pozolotin

Answers:


10

在零均值情况下(B部分)可以获得精确的解。

问题

令表示 iid变量,每个变量具有共同的pdf:(X1,,Xn)nN(0,σ2)f(x)

在此处输入图片说明

我们寻找的pdf ,因为i=1nXin=2,3,

两个这样的法线的乘积的pdf就是:

在此处输入图片说明

...我正在使用 mathStatica软件包中的MathematicaTransformProduct函数。支持范围是:

在此处输入图片说明

3、4、5和6法线的乘积是通过迭代应用相同的函数(此处为四次)获得的:

在此处输入图片说明

...其中,MeijerG表示Meijer G函数

通过归纳, iid随机变量的乘积的pdf 为:nN(0,σ2)


1(2π)n2σnMeijerG[{{},{}},{{01,,0n},{}},x22nσ2n] for xR

快速蒙特卡洛检查

这是比较的快速检查:

  • 刚获得的理论pdf(当且):红色虚线σ = 3n=6σ=3
  • 经验蒙特卡洛pdf:弯曲的蓝色曲线

在此处输入图片说明

看起来不错![弯曲的蓝色Monte曲线遮盖了确切的红色虚线]


太好了,谢谢你,科林。现在,我明白了为什么必须购买您的书了:-)也让我想知道看起来更简单。是时候利用我的Wolfram技能了。log(...MeijerG(...))
Andrei Pozolotin
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.