Взаємозв'язок між сумою гауссових RV та гауссової суміші


13

Я знаю, що сума гауссів - це гаусси. Отже, чим відрізняється суміш гауссів?

Я маю на увазі, що суміш гауссів - це лише сума гауссів (де кожен гаусс множиться на відповідний коефіцієнт змішування) так?


7
Суміш гауссів - це зважена сума гауссових густин , а не зважена сума гауссових випадкових величин.
ймовірністьлогічний

Відповіді:


7

Зважена сума гауссових випадкових величин - випадкова величина Гаусса : if то p i = 1 β i 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)
введіть тут опис зображення

[Джерело: Марін та Роберт, Bayesian Core , 2007]

Випадкова величина з такою щільністю, може бути представлена в вигляді , де і є багаточленом з :X = p i = 1 I ( Z = i ) X i = X Z X iN p ( μ i , σ i ) Z P ( Z = i ) = ω i Z М ( 1 ; ω 1 , , ω pXf(;θ)

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

3

Ось декілька код R, який слід доповнити @ Xi'an:

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


Дякую enthdegree. Я знаю, що наступний приклад за своєю суттю невірний, але це може бути цікаво все одно: скажімо, у нас є особливий вид "суміші" (якщо ми її ще можемо назвати "сумішшю") 2-х гауссових густин, де коефіцієнти змішування обидва відповідають 1, чи буде це однаковою сумою гауссових RV?
njk

Ні, хоча ваша суміш rv буде гауссовою в цьому випадку, якби ви додали два RV з розподілом компонента, сума RV мала б більше дисперсії, ніж суміш RV.
enthdegree

@enthdegree Як суміш rv гауссова? Це все ще може бути бімодальним, якщо засоби не збігаються, правда?
навчання

@learning, так, ти маєш рацію. Коли я написав прев. коментар чомусь я припускав, що вони однакові.
enthdegree
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.