Выпрямленное распределение Гаусса

редактировать

В теории вероятностей выпрямленное распределение Гаусса является модификацией Гауссово распределение при сбросе его отрицательных элементов в 0 (аналог электронного выпрямителя ). По сути, это смесь дискретного распределения (константа 0) и непрерывного распределения (усеченное распределение Гаусса с интервалом (0, ∞) {\ displaystyle (0, \ infty)}(0, \ infty) ) в результате цензуры.

Содержание
  • 1 Функция плотности
  • 2 Среднее значение и дисперсия
  • 3 Создание значений
  • 4 Применение
  • 5 Расширение до общих границ
  • 6 См. Также
  • 7 Ссылки
Функция плотности

Функция плотности вероятности выпрямленного гауссовского распределения для которые случайные величины X, имеющие это распределение, полученное из нормального распределения N (μ, σ 2), {\ displaystyle {\ mathcal {N}} (\ mu, \ sigma ^ {2}),}{\ displaystyle {\ mathcal {N}} (\ mu, \ sigma ^ {2}),} отображаются как X ∼ NR (μ, σ 2) {\ displaystyle X \ sim {\ mathcal {N}} ^ {\ textrm {R}} (\ mu, \ sigma ^ {2})}X \ sim {\ mathcal {N}} ^ {{{\ textrm {R}}}} (\ mu, \ sigma ^ {2}) , определяется как

f (x; μ, σ 2) = Φ (- μ σ) δ (x) + 1 2 π σ 2 e - (x - μ) 2 2 σ 2 U (x). {\ Displaystyle е (х; \ му, \ sigma ^ {2}) = \ Phi (- {\ frac {\ mu} {\ sigma}}) \ delta (x) + {\ frac {1} {\ sqrt {2 \ pi \ sigma ^ {2}}}} \; e ^ {- {\ frac {(x- \ mu) ^ {2}} {2 \ sigma ^ {2}}}} {\ textrm {U }} (x).}f (x; \ mu, \ sigma ^ {2}) = \ Phi (- {\ frac {\ mu} {\ sigma}}) \ delta (x) + {\ frac {1} {{\ sqrt {2 \ pi \ sigma ^ {2}}}}} \; e ^ {{- {\ frac {(x- \ mu) ^ { 2}} {2 \ sigma ^ {2}}}}} {\ textrm {U}} (x).
Сравнение распределения Гаусса, выпрямленного распределения Гаусса и усеченного распределения Гаусса.

Здесь Φ (x) {\ displaystyle \ Phi (x)}\ Phi (x) - кумулятивная функция распределения (cdf) стандартного нормального распределения :

Φ (x) = 1 2 π ∫ - ∞ xe - t 2/2 dtx ∈ R, {\ displaystyle \ Phi (x) = {\ frac {1} {\ sqrt {2 \ pi}}} \ int _ {- \ infty} ^ {x} e ^ {- t ^ {2} / 2} \, dt \ quad x \ in \ mathbb {R},}\ Phi (x) = {\ frac {1} {{\ sqrt {2 \ pi}}}} \ int _ {{- \ infty}} ^ {x} e ^ {{ -t ^ {2} / 2}} \, dt \ quad x \ in {\ mathbb {R}},

δ (x) {\ displaystyle \ delta (x)}\ delta (x) - дельта-функция Дирака

δ (x) = {+ ∞, Икс знак равно 0 0, Икс ≠ 0 {\ Displaystyle \ Delta (х) = {\ begin {cases} + \ infty, x = 0 \\ 0, x \ neq 0 \ end {cases}}}\ delta (x) = {\ begin {cases} + \ infty, x = 0 \\ 0, x \ neq 0 \ end {case}}

и U (x) {\ displaystyle {\ textrm {U}} (x)}{\ textrm {U}} (x) - это функция единичного шага :

U (x) = {0, х ≤ 0, 1, х>0. {\ displaystyle {\ textrm {U}} (x) = {\ begin {cases} 0, x \ leq 0, \\ 1, x>0. \ end {cases}}}{\textrm {U}}(x)={\begin{cases}0,x\leq 0,\\1,x>0. \ end {cases}} <181 Среднее значение и дисперсия

Так как неустановленное нормальное распределение имеет mean μ {\ displaystyle \ mu}\ mu и поскольку при преобразовании его в исправленное распределение некоторая вероятностная масса был сдвинут на более высокое значение (с отрицательных значений на 0), среднее значение выпрямленного распределения больше μ. {\ displaystyle \ mu.}\ mu.

Поскольку выпрямленное распределение формируется путем перемещения некоторых из вероятностной массы к остальной вероятностной массе исправление представляет собой сокращение с сохранением среднего значения в сочетании с жестким сдвигом распределения, изменяющим среднее значение, и, таким образом, дисперсия уменьшается; поэтому дисперсия выпрямленного распределения меньше σ 2. {\ displaystyle \ sigma ^ {2}. }{\ displaystyle \ sigma ^ {2}.}

Генерация значений

Для генерации значений вычислительным способом можно использовать

s ∼ N (μ, σ 2), x = max (0, s), {\ displaystyle s \ sim {\ mathcal {N}} (\ mu, \ sigma ^ {2}), \ quad x = {\ textrm {max}} (0, s),}{\ displaystyle s \ sim {\ mathcal {N}} (\ mu, \ sigma ^ {2}), \ quad x = {\ textrm {max}} (0, s),}

и затем

x ∼ NR (μ, σ 2). {\ displaystyle x \ sim {\ mathcal {N}} ^ {\ textrm {R}} (\ mu, \ sigma ^ {2}).}{\ displaystyle x \ sim {\ mathcal {N}} ^ {\ textrm {R}} (\ mu, \ sigma ^ {2}).}
Приложение

Выпрямленное распределение Гаусса является полу -сопряжен с гауссовским правдоподобием, и недавно он был применен к факторному анализу или, в частности, (неотрицательному) анализу исправленных факторов. Харва предложил алгоритм вариационного обучения для модели ректифицированных факторов, в котором факторы соответствуют смеси выпрямленных гауссовских; а позже Meng предложил модель бесконечного выпрямленного фактора в сочетании со своим решением выборки Гиббса, где факторы следуют процессу Дирихле смесь выпрямленного распределения Гаусса, и применил ее в вычислительной биологии для реконструкции сети регуляции генов.

Расширение до общих границ

Расширение выпрямленного гауссовского распределения было предложено Палмером и др., Позволяющее исправить положение между произвольными нижними и верхними границами. Для нижней и верхней границ a {\ displaystyle a}aи b {\ displaystyle b}b соответственно cdf, FR (x | μ, σ 2) {\ displaystyle F_ {R} (x | \ mu, \ sigma ^ {2})}{\ displaystyle F_ {R} (x | \ mu, \ sigma ^ {2})} определяется по формуле

FR (x | μ, σ 2) = {0, x < a, Φ ( x | μ, σ 2), a ≤ x < b, 1, x ≥ b, {\displaystyle F_{R}(x|\mu,\sigma ^{2})={\begin{cases}0,x{\ displaystyle F_ {R} (x | \ mu, \ sigma ^ {2}) = {\ begin {cases} 0, x <a, \\\ Phi (x | \ mu, \ sigma ^ {2}), a \ leq x <b, \ \ 1, x \ geq b, \\\ end {case}}}

где Φ (x | μ, σ 2) {\ displaystyle \ Phi (x | \ mu, \ sigma ^ {2})}{\ displaystyle \ Phi (x | \ mu, \ sigma ^ {2})} - cdf нормального распределения с среднее значение μ {\ displaystyle \ mu}\ mu и дисперсия σ 2 {\ displaystyle \ sigma ^ {2}}\ sigma ^ {2} . Среднее значение и дисперсия выпрямленного распределения вычисляются путем первого преобразования ограничений, действующих на стандартное нормальное распределение:

c = a - μ σ, d = b - μ σ. {\ displaystyle c = {\ frac {a- \ mu} {\ sigma}}, \ qquad d = {\ frac {b- \ mu} {\ sigma}}.}{\ displaystyle c = {\ frac {a- \ mu} {\ sigma}}, \ qquad d = {\ frac {b- \ mu} {\ sigma}}.}

Используя преобразованные ограничения, среднее и дисперсия, μ R {\ displaystyle \ mu _ {R}}{\ displaystyle \ mu _ {R }} и σ R 2 {\ displaystyle \ sigma _ {R} ^ {2}}{\ displaystyle \ sigma _ {R} ^ {2}} соответственно, тогда задаются следующим образом:

μ t = 1 2 π (e (- c 2 2) - e (- d 2 2)) + c 2 (1 + erf (c 2)) + d 2 ( 1 - erf (d 2)), {\ displaystyle \ mu _ {t} = {\ frac {1} {\ sqrt {2 \ pi}}} \ left (e ^ {\ left (- {\ frac {c ^ {2}} {2}} \ right)} - ​​e ^ {\ left (- {\ frac {d ^ {2}} {2}} \ right)} \ right) + {\ frac {c} { 2}} \ left (1 + {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}}} \ right) \ right) + {\ frac {d} {2}} \ left (1 - {\ textrm {erf}} \ left ({\ frac {d} {\ sqrt {2}}} \ right) \ right),}{\ displaystyle \ mu _ {t} = {\ frac {1} {\ sqrt {2 \ pi}}} \ left (e ^ {\ left (- {\ frac {c ^ {2}} {2 }} \ right)} - ​​e ^ {\ left (- {\ frac {d ^ {2}} {2}} \ right)} \ right) + {\ frac {c} {2}} \ le ft (1 + {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}}} \ right) \ right) + {\ frac {d} {2}} \ left (1- {\ textrm {erf}} \ left ({\ frac {d} {\ sqrt {2}}} \ right) \ right),}
σ t 2 = μ t 2 + 1 2 ( erf (d 2) - erf (c 2)) - 1 2 π ((d - 2 μ t) e (- d 2 2) - (c - 2 μ t) e (- c 2 2)) + (c - μ t) 2 2 (1 + erf (c 2)) + (d - μ t) 2 2 (1 - erf (d 2)), {\ displaystyle {\ begin {align} \ sigma _ {t} ^ {2} = {\ frac {\ mu _ {t} ^ {2} +1} {2}} \ left ({\ textrm {erf}} \ left ({\ frac {d} {\ sqrt { 2}}} \ right) - {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}}} \ right) \ right) - {\ frac {1} {\ sqrt {2 \ pi}}} \ left (\ left (d-2 \ mu _ {t} \ right) e ^ {\ left (- {\ frac {d ^ {2}} {2}} \ right)} - ​​\ left (c-2 \ mu _ {t} \ right) e ^ {\ left (- {\ frac {c ^ {2}} {2}} \ right)} \ right) \\ + {\ frac { \ left (c- \ mu _ {t} \ right) ^ {2}} {2}} \ left (1 + {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}) }} \ right) \ right) + {\ frac {\ left (d- \ mu _ {t} \ right) ^ {2}} {2}} \ left (1 - {\ textrm {erf}} \ left ({\ гидроразрыва {d} {\ sqrt {2}}} \ right) \ right), \ end {align}}}{\ displaystyle {\ begin {align} \ sigma _ {t} ^ {2 } = {\ frac {\ mu _ {t} ^ {2} +1} {2}} \ left ({\ textrm {erf}} \ left ({\ frac {d} {\ sqrt {2}}) } \ right) - {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}}} \ right) \ right) - {\ frac {1} {\ sqrt {2 \ pi} }} \ left (\ left (d-2 \ mu _ {t} \ right) e ^ {\ left (- {\ frac {d ^ {2}} {2}} \ right)} - ​​\ left (c -2 \ mu _ {t} \ right) e ^ {\ left (- {\ frac {c ^ {2}} {2}} \ right)} \ right) \\ + {\ frac {\ left ( c- \ mu _ {t} \ right) ^ {2}} {2}} \ left (1+ {\ textrm {erf}} \ left ({\ frac {c} {\ sqrt {2}}} \ right) \ right) + {\ frac {\ left (d- \ mu _ {t} \ right) ^ {2}} {2}} \ left (1 - {\ textrm {erf}} \ left ({\ frac {d} {\ sqrt {2}}} \ right) \ right), \ end {align}} }
μ R = μ + σ μ t, {\ displaystyle \ mu _ {R} знак равно \ му + \ сигма \ му _ {т},}{\ displaystyle \ mu _ {R} = \ mu + \ sigma \ mu _ {t},}
σ R 2 = σ 2 σ T 2, {\ displaystyle \ sigma _ {R} ^ {2} = \ sigma ^ {2} \ sigma _ {t} ^ {2},}{\ displaystyle \ sigma _ {R} ^ {2} = \ sigma ^ {2} \ sigma _ {t} ^ {2},}

где erf - это функция ошибок. Это распределение использовали Palmer et al. для моделирования уровней физических ресурсов, таких как количество жидкости в емкости, которое ограничено как 0, так и вместимостью емкости.

Ссылки
Последняя правка сделана 2021-06-03 10:31:13
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте