Отрицательное полиномиальное распределение

редактировать
ОбозначениеNM (x 0, p) {\ displaystyle {\ textrm {NM}} (x_ {0}, \, p)}{\ displaystyle { \ textrm {NM}} (x_ {0}, \, p)}
Параметрыx0∈ N0- количество неудач до остановки эксперимента,. p ∈ R - m-вектор вероятностей «успеха»,.
p0= 1 - (p 1 +… + p m) - вероятность «отказа».
Поддержка xi ∈ {0, 1, 2,…}, 1 ≤ я ≤ м {\ Displaystyle x_ {i} \ in \ {0,1,2, \ ldots \}, 1 \ leq i \ leq m}{\ displaystyle x_ {i} \ in \ {0,1,2, \ ldots \}, 1 \ leq i \ leq m}
PDF Γ (∑ i = 0 mxi) p 0 x 0 Γ (x 0) ∏ i = 1 mpixixi!, {\ displaystyle \ Gamma \! \ left (\ sum _ {i = 0} ^ {m} {x_ {i}} \ right) {\ frac {p_ {0} ^ {x_ {0}}} {\ Гамма (x_ {0})}} \ prod _ {i = 1} ^ {m} {\ frac {p_ {i} ^ {x_ {i}}} {x_ {i}!}},}{ \ Displaystyle \ Gamma \! \ left (\ sum _ {i = 0} ^ {m} {x_ {i}} \ right) {\ frac {p_ {0} ^ {x_ {0}}} {\ Gamma ( x_ {0})}} \ prod _ {i = 1} ^ {m} {\ frac {p_ {i} ^ {x_ {i}}} {x_ {i}!}},} . где Γ (x) - Гамма-функция.
Среднее x 0 p 0 p {\ displaystyle {\ tfrac {x_ {0}} {p_ {0}}} \, p}{\ displaystyle {\ \ tfrac {x_ {0}} {p_ {0}}} \, p}
Дисперсия x 0 p 0 2 pp ′ + x 0 p 0 diag ⁡ (p) {\ displaystyle {\ tfrac {x_ {0}} {p_ {0} ^ {2}}} \, pp '+ { \ tfrac {x_ {0}} {p_ {0}}} \, \ operatorname {diag} (p)}{\displaystyle {\tfrac {x_{0}}{p_{0}^{2}}}\,pp'+{\tfrac {x_{0}}{p_{0}}}\,\operatorname {diag} (p)}
CF (p 0 1 - p ′ eit) x 0 {\ displaystyle {\ bigg (} {\ frac {p_ {0}} {1-p'e ^ {it}}} {\ bigg)} ^ {\! x_ {0}}}{\displaystyle {\bigg (}{\frac {p_{0}}{1-p'e^{it}}}{\bigg)}^{\!x_{0}}}

В теории вероятностей и статистики, отрицательное полиномиальное распределение является обобщением отрицательного биномиального распределения (NB (r, p)) для более чем двух исходов.

Предположим у нас есть эксперимент, который генерирует m + 1≥2 возможных результатов, {X 0,..., X m }, каждый из которых имеет неотрицательные вероятности {p 0,..., p m } соответственно. Если бы выборка продолжалась до тех пор, пока не было сделано n наблюдений, то {X 0,..., X m } было бы полиномиально распределенным. Однако, если эксперимент прекращается, когда X 0 достигает заданного значения x 0, тогда распределение m-кортежа {X 1,..., X m } - отрицательный полином. Эти переменные не распределены полиномиально, потому что их сумма X 1 +... + X m не является фиксированной и является результатом отрицательного биномиального распределения.

Содержание
  • 1 Свойства
    • 1.1 Предельные распределения
    • 1.2 Независимые суммы
    • 1.3 Агрегирование
    • 1.4 Корреляционная матрица
  • 2 Оценка параметров
    • 2.1 Метод моментов
  • 3 Связанные распределения
  • 4 Ссылки
  • 5 Дополнительная литература
Свойства

Предельные распределения

Если m-мерное x разбивается следующим образом

X = [X (1) X (2)] с размерами [n × 1 (m - n) × 1] {\ displaystyle \ mathbf {X} = {\ begin {bmatrix} \ mathbf {X} ^ {(1)} \\\ mathbf { X} ^ {(2)} \ end {bmatrix}} {\ text {with sizes}} {\ begin {bmatrix} n \ times 1 \\ (mn) \ times 1 \ end {bmatrix}}}{\ displaystyle \ mathbf {X} = {\ begin {bmatrix} \ mathbf {X} ^ {(1)} \\ \ mathbf {X} ^ {(2)} \ end {bmatrix}} {\ text {с размерами}} {\ begin {bmatrix} n \ times 1 \\ (mn) \ times 1 \ end {bmatrix}}}

и соответственно p {\ displaystyle {\ boldsymbol {p}}}{\ boldsymbol {p}}

p = [p (1) p (2)] с размерами [n × 1 (m - n) × 1] {\ displaystyle { \ boldsymbol {p}} = {\ begin {bmatrix} {\ boldsymbol {p}} ^ {(1)} \\ {\ boldsymbol {p}} ^ {(2)} \ end {bmatrix }} {\ text {with sizes}} {\ begin {bmatrix} n \ times 1 \\ (mn) \ times 1 \ end {bmatrix}}}{\ displaystyle {\ boldsymbol {p}} = {\ begin {bmatrix} {\ boldsymbol {p}} ^ {(1)} \\ {\ boldsymbol {p}} ^ {(2) } \ end {bmatrix}} {\ text {с размерами}} {\ begin {bmatrix} n \ times 1 \\ (mn) \ times 1 \ end {bmatrix}}}

и пусть

q = 1 - ∑ ipi ( 2) знак равно п 0 + ∑ ipi (1) {\ displaystyle q = 1- \ sum _ {i} p_ {i} ^ {(2)} = p_ {0} + \ sum _ {i} p_ {i} ^ {(1)}}{\ displaystyle q = 1- \ sum _ {i} p_ {i} ^ {(2)} = p_ {0} + \ sum _ {i} p_ {i} ^ {(1)}}

Предельное распределение X (1) {\ displaystyle {\ boldsymbol {X}} ^ {(1)}}{ \ displaystyle {\ boldsymbol {X}} ^ {(1)}} равно NM ( Икс 0, п 0 / q, п (1) / q) {\ displaystyle \ mathrm {NM} (x_ {0}, p_ {0} / q, {\ boldsymbol {p}} ^ {(1)} / q)}{\ displaystyle \ mathrm {NM} (x_ {0}, p_ {0} / q, { \ boldsymbol {p}} ^ {(1)} / q)} . То есть маргинальное распределение также является отрицательным полиномиальным с удалением p (2) {\ displaystyle {\ boldsymbol {p}} ^ {(2)}}{\ displaystyle {\ boldsymbol {p}} ^ {(2)}} , а оставшиеся p правильно масштабированы, чтобы добавить к одному.

Одномерное маргинальное m = 1 {\ displaystyle m = 1}m = 1 - это отрицательное биномиальное распределение.

Независимые суммы

Если X 1 ∼ NM (r 1, p) {\ displaystyle \ mathbf {X} _ {1} \ sim \ mathrm {NM} (r_ { 1}, \ mathbf {p})}{\ displaystyle \ mathbf {X} _ {1} \ sim \ mathrm {NM} (r_ {1}, \ mathbf { p})} и If X 2 ∼ NM (r 2, p) {\ displaystyle \ mathbf {X} _ {2} \ sim \ mathrm {NM} (r_ {2}, \ mathbf {p})}{\ displaystyle \ mathbf {X} _ {2} \ sim \ mathrm {NM} (r_ {2}, \ mathbf {p})} независимы, тогда X 1 + X 2 ∼ NM (r 1 + r 2, p) {\ displaystyle \ mathbf {X} _ {1} + \ mathbf {X} _ {2} \ sim \ mathrm {NM} (r_ {1} + r_ {2}, \ mathbf {p})}{\ displaystyle \ mathbf {X} _ {1} + \ mathbf {X} _ {2} \ sim \ mathrm {NM} (r_ {1} + r_ {2}, \ mathbf {p})} . Аналогично и наоборот, из характеристической функции легко увидеть, что отрицательный многочлен бесконечно делим.

Агрегация

Если

X = (X 1,…, X m) ∼ NM ⁡ (Икс 0, (п 1,…, pm)) {\ displaystyle \ mathbf {X} = (X_ {1}, \ ldots, X_ {m}) \ sim \ operatorname {NM} (x_ {0}, (p_ {1}, \ ldots, p_ {m}))}{\ displaystyle \ mathbf {X} = (X_ {1}, \ ldots, X_ {m}) \ sim \ operatorname {NM} (x_ {0}, (p_ {1}, \ ldots, p_ {m}))}

тогда, если случайные величины с индексами i и j исключены из вектора и заменены их суммой,

X ′ = (X 1,…, X i + X j,…, X m) ∼ NM ⁡ (x 0, (p 1,…, pi + pj,…, pm)). {\ displaystyle \ mathbf {X} '= (X_ {1}, \ ldots, X_ {i} + X_ {j}, \ ldots, X_ {m}) \ sim \ operatorname {NM} (x_ {0}, (p_ {1}, \ ldots, p_ {i} + p_ {j}, \ ldots, p_ {m})).}{\displaystyle \mathbf {X} '=(X_{1},\ldots,X_{i}+X_{j},\ldots,X_{m})\sim \operatorname {NM} (x_{0},(p_{1},\ldots,p_{i}+p_{j},\ldots,p_{m})).}

Это свойство агрегирования можно использовать для получения предельного распределения X i {\ displaystyle X_ {i}}X_ {i} упомянутый выше.

Корреляционная матрица

Элементы корреляционной матрицы :

ρ (X i, X i) = 1. {\ displaystyle \ rho (X_ {i }, X_ {i}) = 1.}\ rho (X_i, X_i) = 1.
ρ (X i, X j) = cov ⁡ (X i, X j) var ⁡ (X i) var ⁡ (X j) = pipj (p 0 + pi) (p 0 + pj). {\ displaystyle \ rho (X_ {i}, X_ {j}) = {\ frac {\ operatorname {cov} (X_ {i}, X_ {j})} {\ sqrt {\ operatorname {var} (X_ { i}) \ operatorname {var} (X_ {j})}}} = {\ sqrt {\ frac {p_ {i} p_ {j}} {(p_ {0} + p_ {i}) (p_ {0 } + p_ {j})}}}.}{\ displaystyle \ rho (X_ {i}, X_ {j}) = {\ frac {\ operatorname {cov} (X_ {i}, X_ {j})} {\ sqrt {\ operatorname {var} (X_ {i}) \ operatorname {var} (X_ { j})}}} = {\ sqrt {\ frac {p_ {i} p_ {j}} {(p_ {0} + p_ {i}) (p_ {0} + p_ {j})}}}. }
Оценка параметра

Метод моментов

Если мы позволим среднему вектору отрицательного многочлена быть

μ = x 0 p 0 p {\ displaystyle {\ boldsymbol {\ mu}} = {\ frac {x_ {0}} {p_ {0}}} \ mathbf {p}}{\ displaystyle {\ boldsymbol {\ mu}} = {\ frac {x_ {0}} {p_ {0}}} \ mathbf {p}}

и ковариационная матрица

Σ знак равно Икс 0 п 0 2 пп '+ Икс 0 п 0 диаг (р) {\ Displaystyle {\ boldsymbol {\ Sigma}} = {\ tfrac {x_ {0}} {p_ {0} ^ {2}}} \, \ mathbf {p} \ mathbf {p} '+ {\ tfrac {x_ {0}} {p_ {0}}} \, \ operatorname {diag} (\ mathbf {p})}{\displaystyle {\boldsymbol {\Sigma }}={\tfrac {x_{0}}{p_{0}^{2}}}\,\mathbf {p} \mathbf {p} '+{\tfrac {x_{0}}{p_{0}}}\,\operatorname {diag} (\mathbf {p})},

тогда это через свойства определителей легко показать, что | Σ | Знак равно 1 п 0 ∏ я знак равно 1 м μ я {\ displaystyle | {\ boldsymbol {\ Sigma}} | = {\ frac {1} {p_ {0}}} \ prod _ {i = 1} ^ {m} {\ mu _ {i}}}{\ displaystyle | { \ boldsymbol {\ Sigma}} | = {\ frac {1} {p_ {0}}} \ prod _ {i = 1} ^ {m} {\ mu _ {i}}} . Отсюда можно показать, что

x 0 = ∑ μ i ∏ μ i | Σ | - ∏ μ я {\ displaystyle x_ {0} = {\ frac {\ sum {\ mu _ {i}} \ prod {\ mu _ {i}}} {| {\ boldsymbol {\ Sigma}} | - \ prod {\ mu _ {i}}}}}{\ displaystyle x_ {0} = {\ frac {\ sum {\ mu _ {i}} \ prod {\ mu _ {i}}} {| {\ boldsymbol {\ Sigma}} | - \ prod {\ mu _ {i}}}}}

и

p = | Σ | - ∏ μ i | Σ | ∑ μ i μ. {\ displaystyle \ mathbf {p} = {\ frac {| {\ boldsymbol {\ Sigma}} | - \ prod {\ mu _ {i}}} {| {\ boldsymbol {\ Sigma}} | \ sum {\ mu _ {i}}}} {\ boldsymbol {\ mu}}.}{\ displaystyle \ mathbf {p} = {\ frac {| {\ boldsymbol {\ Sigma}} | - \ prod { \ mu _ {i}}} {| {\ boldsymbol {\ Sigma}} | \ sum {\ mu _ {i}}}} {\ boldsymbol {\ mu}}.}

Подстановка моментов выборки дает метод моментов оценок

x ^ 0 = (∑ i = 1 mxi ¯) ∏ i = 1 mxi ¯ | S | - ∏ я знак равно 1 мкси ¯ {\ Displaystyle {\ шляпа {х}} _ {0} = {\ гидроразрыва {(\ сумма _ {я = 1} ^ {м} {{\ бар {x_ {i}}})} \ prod _ {i = 1} ^ {m} {\ bar {x_ {i}}}} {| \ mathbf {S} | - \ prod _ {i = 1} ^ {m} {\ bar { x_ {i}}}}}}{\ displaystyle {\ hat {x}} _ {0} = {\ frac {(\ sum _ {i = 1} ^ {m} {{\ bar {x_ {i}}})} \ prod _ {i = 1} ^ {m} {\ bar { x_ {i}}}} {| \ mathbf {S} | - \ prod _ {i = 1} ^ {m} {\ bar {x_ {i}}}}}}}

и

p ^ = (| S | - ∏ i = 1 mx ¯ i | S | ∑ i = 1 mx ¯ i) x ¯ {\ displaystyle {\ hat {\ mathbf {p}}} = \ left ({\ frac {| {\ boldsymbol {S}} | - \ prod _ {i = 1} ^ {m} {{\ bar {x}} _ {i}) }} {| {\ boldsymbol {S}} | \ sum _ {i = 1} ^ {m} {{\ bar {x}} _ {i}}}} \ right) {\ boldsymbol {\ bar {x }}}}{\ displaystyle {\ hat {\ mathbf {p}}} = \ left ({\ frac {| {\ boldsymbol {S}} | - \ prod _ {i = 1} ^ {m} {{\ bar {x}} _ {i}}} {| {\ boldsymbol {S}} | \ sum _ {i = 1} ^ {m} {{\ bar {x}} _ {i}}}} \ right) {\ boldsymbol {\ bar {x}}}}
Связанные распределения
Ссылки
  1. ^Ле Галл, Ф. Режимы отрицательного полиномиального распределения, Статистические и вероятностные письма, Том 76, выпуск 6, 15 марта 2006 г., страницы 619-624, ISSN 0167-7152, 10.1016 / j.spl.2005.09.009.

Уоллер Л.А., Зельтерман Д. (1997). Логлинейное моделирование с отрицательным многочленным распределением. Биометрия 53: 971-82.

Дополнительная литература

Johnson, Norman L.; Коц, Самуэль; Балакришнан, Н. (1997). «Глава 36: Отрицательные полиномиальные и другие полиномиальные распределения». Дискретные многомерные распределения. Вайли. ISBN 978-0-471-12844-1.

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