两个正态分布之差的分布


20

我有两个正态分布的概率密度函数:

f1(x1|μ1,σ1)=1σ12πe(xμ1)22σ12

f2(x2|μ2,σ2)=1σ22πe(xμ2)22σ22

我正在寻找和之间分离的概率密度函数。我认为这意味着我正在寻找| x_1-x_2 |的概率密度函数。。那是对的吗?我怎么找到那个?x1x2|x1x2|


如果这是家庭作业,请使用self-study标签。我们接受家庭作业问题,但是在这里我们对它们的处理有些不同。
shadowtalker

另外,我不想成为“那个家伙”,但是您尝试过Google吗?“正态分布之间的差异”很快就为我找到了答案。
shadowtalker

@ssdecontrol不,不是家庭作业,但这是一个业余项目,因此如果我走对了路,我不介意自己找一些东西。我曾经尝试过使用google,但是我对此事的掌握程度非常有限,以至于如果它就在我眼前,我可能不会意识到。用引号引起来,我发现很多东西类似于某些x的“正态分布与x的区别”。
Martijn 2015年

Answers:


25

只能通过假设由这些分布控制的两个随机变量和是独立的来回答这个问题。X 2X1X2 X= X 2 - X 1 μ= μ 2 - μ 1 σ 2 = σ 2 1 + σ 2 2 X 1 X 2 这使得它们的差值正常,均值且方差。(以下解决方案可以轻松地推广为任何二元正态分布。)因此,变量X=X2X1μ=μ2μ1σ2=σ12+σ22(X1,X2)

Z=Xμσ=X2X1(μ2μ1)σ12+σ22

具有标准正态分布(即均值和单位方差为零)并且

X=σ(Z+μσ).

表达方式

|X2X1|=|X|=X2=σ(Z+μσ)2

表现出绝对差异作为具有一个自由度和非中心性参数的非中心卡方分布的平方根的缩放版本。具有这些参数的非中心卡方分布具有概率元素λ=(μ/σ)2

f(y)dy=y2πe12(λy)cosh(λy)dyy, y>0.

为写可以在和它的平方根之间建立一一对应的关系,结果是 x > 0 yy=x2x>0y

f(y)dy=f(x2)d(x2)=x22πe12(λx2)cosh(λx2)dx2x2.

对此进行简化,然后按重新缩放,即可得到所需的密度,σ

f|X|(x)=1σ2πcosh(xμσ2)exp(x2+μ22σ2).

此结果得到仿真的支持,例如,的100,000次独立绘制的直方图 (在代码中称为“ x”),其参数为。在其上绘制了,它与直方图值恰好重合。μ 1 = - 1 μ 2 = 5 σ 1 = 4 σ 2 = 1 ˚F | X ||X|=|X2X1|μ1=1,μ2=5,σ1=4,σ2=1f|X|

数字

R此仿真的代码如下。

#
# Specify parameters
#
mu <- c(-1, 5)
sigma <- c(4, 1)
#
# Simulate data
#
n.sim <- 1e5
set.seed(17)
x.sim <- matrix(rnorm(n.sim*2, mu, sigma), nrow=2)
x <- abs(x.sim[2, ] - x.sim[1, ])
#
# Display the results
#
hist(x, freq=FALSE)
f <- function(x, mu, sigma) {
 sqrt(2 / pi) / sigma * cosh(x * mu / sigma^2) * exp(-(x^2 + mu^2)/(2*sigma^2)) 
}
curve(f(x, abs(diff(mu)), sqrt(sum(sigma^2))), lwd=2, col="Red", add=TRUE)

如果我想求平方差,这会有什么不同?例如,如果我想要吗?(f1(.)f2(.))2
user77005

1
@ user77005答案在我的文章中:这是一个非中心的卡方分布。点击链接获取详细信息。
ub

21

我提供的答案与@whuber的答案相辅相成,可以说是非统计学家(例如,对一个自由度的非中心卡方分布了解不多的人)可能写的内容,而且新手可以相对容易地跟随。

借用独立性的假设以及whuber答案的符号,其中 和。因此,对于, ,当然,对于,。关于区别在于 μ = μ 1 - μ 2 σ 2 = σ 2 1 + σ 2 2 X 0 ˚F | Z | x Z=X1X2N(μ,σ2)μ=μ1μ2σ2=σ12+σ22x0

F|Z|(x)P{|Z|x}=P{xZx}=P{x<Zx}since Z is a continuous random variable=FZ(x)FZ(x),
F|Z|(x)=0x<0x
F|ž|XXF|ž|X=[FžX+Fž-X]1个0X=[经验值-X-μ22σ2σ2π+经验值-X+μ22σ2σ2π]1个0X=经验值-X2+μ22σ2σ2π经验值Xμσ2+经验值Xμσ21个0X=1个σ2π科什Xμσ2经验值-X2+μ22σ21个0X
结果与whuber的答案是完全相同的,但得出的结论更透明。

1
+1我总是喜欢看到可以从最基本的原理和假设出发的解决方案。
ub

1

假设X和Y是独立的,则两个正态分布变量X和Y的差值的分布也是正态分布(感谢Mark的评论)。这是派生词:http : //mathworld.wolfram.com/NormalDifferenceDistribution.html

在这里,您根据whuber的答案询问绝对差,并且如果我们假设X和Y的均值差为零,那么它只是密度的两倍的正态分布的一半(感谢Dilip的评论)。


3
您和Wolfram Mathworld隐式地假设2个正态分布(随机变量)是独立的。所不同的是甚至不需要正态分布,如果2个正常随机变量不是二元正常的,如果它们不是独立的..这可能发生
马克·劳伦斯斯通

4
除了Mark指出的假设之外,您还忽略了手段不同的事实。半正常情况仅在因此差的均值为。 0μ1个=μ20
Dilip Sarwate 2015年

谢谢您的意见。现在,我根据您的评论和胡夫的答案修改了我的答案。
yuqian 2015年
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.