Я знаю, що сума гауссів - це гаусси. Отже, чим відрізняється суміш гауссів?
Я маю на увазі, що суміш гауссів - це лише сума гауссів (де кожен гаусс множиться на відповідний коефіцієнт змішування) так?
Я знаю, що сума гауссів - це гаусси. Отже, чим відрізняється суміш гауссів?
Я маю на увазі, що суміш гауссів - це лише сума гауссів (де кожен гаусс множиться на відповідний коефіцієнт змішування) так?
Відповіді:
Зважена сума гауссових випадкових величин - випадкова величина Гаусса : if то p ∑ i = 1 β i X i ( X 1 , … , X p ) ∼ N p ( μ , Σ ) β T ( X 1 , … , X p ) ∼ N 1 ( β T μ , β T Σ β )
Суміш гауссових густин має щільність, задану у вигляді зваженої суми гауссових густин : яка майже незмінно не дорівнює гауссова щільності. Дивіться, наприклад, нижню оцінювану щільність суміші нижче (де жовта смуга є мірою змінності розрахункової суміші):
[Джерело: Марін та Роберт, Bayesian Core , 2007]
Випадкова величина з такою щільністю, може бути представлена в вигляді , де і є багаточленом з :X = p ∑ i = 1 I ( Z = i ) X i = X Z X i ∼ N p ( μ i , σ i ) Z P ( Z = i ) = ω i Z ∼ М ( 1 ; ω 1 , … , ω p
Ось декілька код 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)
Розподіл суми незалежних випадкових величин - це згортка їх розподілу. Як ви вже зазначали, згортання двох гауссів відбувається як гауссова.