가우스 RV의 합과 가우스 혼합물의 관계


13

가우시안이 가우시안이라는 것을 알고 있습니다. 가우시안 혼합은 어떻게 다른가요?

가우시안 혼합은 가우시안의 합일뿐입니다 (각 가우시안에 각 혼합 계수를 곱한 값)?


7
가우스 혼합은 가우스 랜덤 변수의 가중 합계가 아닌 가우스 밀도 의 가중 합계입니다.
probabilityislogic

답변:


7

가우스 랜덤 변수 의 가중 합 는 가우스 랜덤 변수입니다 . if 다음 P Σ= 1 β 내가 X I ( X 1 , ... , X의 P ) ~ N의 P ( μ , Σ ) β T ( X 1 , ... , X의 P ) ~ N 1 ( β T μ , β T Σ β )X1,,Xp

i=1pβiXi
(X1,,Xp)Np(μ,Σ)
βT(X1,,Xp)N1(βTμ,βTΣβ)

가우스의 혼합물 의 밀도는 가우스의 가중 합으로서 주어진 밀도 갖는다 밀도 : 어느 가우스 밀도와 거의 같지 않습니다. 예를 들어 아래의 청색 추정 혼합물 밀도를 참조하십시오 (노란색 밴드는 추정 혼합물의 변동성 측정 값임).

f(;θ)=i=1pωiφ(;μi,σi)
여기에 이미지 설명을 입력하십시오

[출처 : 마린과 로버트 , 2007 년 베이지안 코어 ]

랜덤 변수 이 밀도, 로서 표현 될 수있는 여기서 및 는 와 함께 :X = p i = 1 I ( Z = i ) X i = X Z X iN p ( μ i , σ i ) Z P ( Z = i ) = ω i Z M ( 1 ; ω 1 , , ω pXf(;θ)

X=i=1pI(Z=i)Xi=XZ
XiNp(μi,σi)ZP(Z=i)=ωi
ZM(1;ω1,,ωp)

3

그리고 @ Xi'an 답변을 보완하는 R 코드가 있습니다.

par(mfrow=c(2,1))
nsamples <- 100000

# Sum of two Gaussians
x1 <- rnorm(nsamples, mean=-10, sd=1)
x2 <- rnorm(nsamples, mean=10, sd=1)
hist(x1+x2, breaks=100)

# Mixture of two Gaussians
z <- runif(nsamples)<0.5 # assume mixture coefficients are (0.5,0.5)
x1_x2 <- rnorm(nsamples,mean=ifelse(z,-10,10),sd=1)
hist(x1_x2,breaks=100)

여기에 이미지 설명을 입력하십시오


1

독립 랜덤 변수의 합의 분포는 분포의 컨벌루션 입니다. 앞서 언급했듯이 두 가우시안의 컨볼 루션은 가우시안입니다.

X,YZXYZ=XZ=Y


감사합니다. 다음 예제는 본질적으로 잘못되었다는 것을 알고 있지만 어쨌든 흥미로울 수 있습니다. 혼합 계수가 2 가우스 밀도 인 특수한 "혼합물"(여전히 "혼합물"이라고 할 수있는 경우)이라고 가정 해 봅시다. 둘 다 1에 해당하며 가우시안 RV의 합과 같습니까?
njk

아니요,이 경우 혼합물 rv는 가우시안이지만 성분 분포를 사용하여 두 개의 RV를 추가하는 경우 합 RV는 혼합물 RV보다 분산이 더 큽니다.
enthdegree

@enthdegree 혼합물 rv 가우시안은 어떻습니까? 수단이 일치하지 않으면 여전히 양손 일 수 있습니다.
학습

@learning, 그렇습니다. 내가 이전을 썼을 때. 어떤 이유로 나는 그들이 동일한 평균을 가지고 있다고 가정했습니다.
enthdegree
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.