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


12

给定 iid和,寻找:X Ñ ≈ Ñ(μ X,σ 2 X)μ X ≈ 0N≥30Xn≈N(μX,σX2)μX≈0

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

这是一个特殊情况,是一个更一般的问题。μX≈0


1.您是否有有关和?(例如,如果所有都很好。)(2)渐近法线逼近将是可怕的,因为渐近线看起来不会是法线。σ X μ X / σ X » 0μXσXμX/σX≫0Y
— 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定理在所示)。μX≥10σXYlog⁡(X)
— whuber

@whuber对Barry-Esseen的直接应用给出,的确不错,但是它一些结构:应该为负,应该取决于等,也许有更好的应用方法?μ˚Fσ˚FFN≈0+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 x∈R

快速蒙特卡洛检查

这是比较的快速检查:

  • 刚获得的理论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.