korolev_matematicheskie_osnovy_teorii_ri ska (811435), страница 18
Текст из файла (страница 18)
В некоторых случаях проще иметь дело с другой записью (2.3.1), а именноq(z) = exp½X∞¾aj (z j − 1) ,(2.3.2)j=0где aj = µpj .2.3.1Примеры дискретных обобщенных пуассоновских распределенийРассмотрим некоторые конкретные примеры дискретных обобщенныхпуассоновских распределений.Пример 2.3.1. Прежде всего, это само пуассоновское распределение. Этому распределению в (2.3.1) соответствует p(z) ≡ z.Пример 2.3.2 (продолжение Примера 2.2.4). Пусть в Примере 2.2.4(точнее, в соотношении (2.2.9)) k = 2.
В таком случае представление(2.3.2) принимает видq(z) = exp{a1 (z − 1) + a2 (z 2 − 1)},|z| ≤ 1,(2.3.3)(a1 = p1 µ, a2 = p2 µ1, p1 ≥ 0, p2 ≥ 0, p1 + p + 2 = 1, так что a1 + a2 = µ).Несложно убедиться, что вероятности qn удовлетворяют соотношениям[n/2]qn = e−µXj=0an−2jaj21,(n − 2j)!j!n = 0, 1, 2, . .
. .(2.3.4)2.3 Дискретные обобщенные пуассоновские распределенияРаспределение (2.3.4) известно как распределение Эрмита, см., например, (Kemp and Kemp, 1965). В (Белов, Галкин и Уфимцев, 1985)рассмотрены дальнейшие обобщения распределения Эрмита на случаибольшей множественности k > 2.Пример 2.3.3. Производящая функция (инфекционного) распределения Неймана типа А (см. Пример 1.4.1) может быть записана в виде(2.3.2) с aj = µe−λ λj /j!. Для соответствующих вероятностей мы имеем−λqn = e−µ+µenλnλn Xmn (µe−λ ) = q0(µe−λ )j σn(j) ,n!n! j=0n = 0, 1, 2, .
. . .Здесь q0 = exp{µ(e−λ − 1)} and σn(j) – числа Стирлинга второго рода(см., например, (Riordan, 1968)).Пример 2.3.4. Пусть 0 < p < 1, m = −µ/ log p. Для j ≥ 1 в (2.3.2)положим aj = m(1 − p)j /j. Тогда, как показано в (Quenouille, 1949),(Gurland, 1957), полученное распределение (2.3.2) оказывается отрицательным биномиальным (распределением Паскаля) с вероятностямиnqn = Cn+m−1pm (1 − p)n ,n = 0, 1, 2, . . . .Отметим, что мы уже встречались с распределением, называемым фишеровским распределением логарифмического ряда или логарифмическим распределениемpj = −(1 − p)j,j log pj = 1, 2, . .
. ,с производящей функциейp(z) =log(1 − (1 − p)z)log pв Лемме 2.2.2.Как мы уже отмечали, если в (2.2.1) F (x, y) – функция распределения Пуассона, то независимо от распределения Q(y), распределениеH(x) называется смешанным пуассоновским. Такие распределения более подробно будут рассмотрены в главе 6. Здесь же мы лишь ограничимся замечанием, что пример распределения Неймана типа А убеждает нас в том, что некоторые распределения могут одновременно бытьсмешанными пуассоновскими и пуассоновски-смешанными. Более того,Как мы убедимся в главе 6, если в (2.2.1) F (x, y) – функция распределения Пуассона, а Q(y) – функция гамма-распределения, то итоговаясмесь является отрицательным биномиальным распределением.
Вместес Примером 2.3.4 это дает убедительное свидетельство существованиявероятностных распределений, которые одновременно являются обобщенными пуассоновскими и смешанными пуассоновскими.1031042. Свойства случайных сумм2.3.2Рекуррентные соотношения для дискретныхобобщенных пуассоновских распределенийТочно так же, как мы получили рекуррентные соотношения для моментов обобщенных пуассоновских распределений в разделе 1.5, мыможем получить рекуррентные соотношения для вероятностей, определяющих дискретные обобщенные пуассоновские распределения и факториальных моментов дискретных обобщенных пуассоновских распределений. Если X – случайная величина, то ее факториальный момент порядка k ≥ 1, обозначаемый EX [k] , определяется как EX [k] =EX(X − 1) · . .
. · (X − k + 1). В дальнейшем в качестве X мы рассматриваем случайную величину с производящей функцией q(z) и пусть Y– случайная величина с производящей функцией p(z), участвующей в(2.3.1), так что qj = P(X = j), pj = P(Y = j).Теорема 2.3.1. Для любого n ≥ 0qn =XXµ n−11 n−1(k + 1)qn−k−1 pk+1 =(k + 1)qn−k−1 ak+1 .n k=0n k=0(2.3.5)Если EY n существует при некотором n ≥ 1, тоEX [n] = µn−1XkCn−1EX [n−k−1] · EY [k+1] .(2.3.6)k=0Д о к а з а т е л ь с т в о.
Прежде всего заметим, чтоdq(z)dp(z)dp(z)= µ exp{µ(p(z) − 1)}= µq(z),dzdzdzтак что, применив правило дифференцирования Лейбница для n ≥ 1,мы получимn−1Xdn q(z)dn−k−1 q(z) dk+1 q(z)k=µC·.n−1dz ndz n−k−1dz k+1k=0Теперь, чтобы получить(2.3.7) и формулами(2.3.5),воспользуемся¯¯(2.3.7)соотношениемdn q(z) ¯¯dk p(z) ¯¯n!qn =,k!p=,kdz n ¯z=0dz k ¯z=0справедливыми для любых n ≥ 1 и k ≥ 1. Чтобы получить (2.3.6),воспользуемся соотношением (2.3.7) и формулами¯EX [n] =Теорема доказана.dn q(z) ¯¯,dz n ¯z=1¯EX [k] =dk p(z) ¯¯.dz k ¯z=12.3 Дискретные обобщенные пуассоновские распределения2.3.3105Дискретные безгранично делимые законы какобобщенные пуассоновские распределенияСледуя определению, данному в книге (Феллер, 1984), том 1, главаXII, раздел 3, мы можем определить дискретные безгранично делимыезаконы через их производящие функции.
А именно, дискретное распределение вероятностей, сосредоточенное в неотрицательных целыхточках, называется дискретным безгранично делимым, если для любого n ≥ 1 его производящая функция q(z) может быть представленакак n-я степень некоторой производящей функции. Другими словами,дискретное распределение вероятностей с производящей функцией q(z)безгранично делимо, если для любого n ≥ 1 существует производящаяфункция hn (z) такая, что q(z) ≡ (hn (z))n .
Легко видеть, что если производящая функция q(z) записана в виде (2.3.1), то она представима ввиде q(z) = (hn (z))n с hn (z) = exp{ nµ (p(z) − 1)} и, следовательно, любое дискретное обобщенное пуассоновское распределение безграничноделимо.Обратно, в (Феллер, 1984), том 1, глава XII, раздел 3, доказано,что любое дискретное безграничное распределение представимо в виде(2.3.1) (или (2.3.2)). Это означает, что справедливо следующее утверждение.Теорема 2.3.2. Дискретное распределение, сосредоточенное внеотрицательных целых точках, безгранично делимо тогда и только тогда, когда оно является обобщенным пуассоновским.Следствие 2.3.1.
Пусть X1 , X2 , . . . – одинаково распределенныеслучайные величины. Пусть N – случайная величина с дискретнымбезгранично делимым распределением. Предположим, что случайныевеличины N, X1 , X2 , . . . независимы в совокупности. Тогда распределение случайной суммы S = X1 + . . .
+ XN является обобщенным пуассоновским.Д о к а з а т е л ь с т в о. В соответствии с Теоремой 2.3.2 распределение случайной величины N является обобщенным пуассоновскимс производящей функцией q(z) вида (2.3.1). Как мы видели в разделе 1.4, характеристическая функция g(t) случайной величины S имеетвид q(f (t)), где f (t) – характеристическая функция случайной величины X1 . Но тогда, очевидно,g(t) = q(f (t)) = exp{µ(p(f (t)) − 1)} = exp{µ(ge(t) − 1)},(2.3.8)где ge(t) = p(f (t)) – характеристическая функция (точнее, это – характеристическая функция случайной суммы X1 + .
. . + XM , где M– случайная величина с производящей функцией p(z), независимая от1062. Свойства случайных суммX1 , X2 , . . .). Требуемое утверждение вытекает из представления (2.3.8).2.42.4.1Асимптотическая ноpмальностьпуассоновских случайных суммСходимость распределений пуассоновскихслучайных сумм к нормальному законуПусть X1 , X2 , . . .
– независимые одинаково pаспpеделенные случайныевеличины. Пусть Nλ – случайная величина, имеющая pаспpеделениеПуассона с паpаметpом λ.Пpедположим, что случайные величины Nλ , X1 , X2 , . . . независимыпpи каждом λ. ОбозначимSλ = X1 + . . . + XNλP(для опpеделенности мы будем считать, что 0j=1 = 0). Случайная величина Sλ называется пуассоновской случайной суммой.Из классической теоpии пpедельных теоpем для сумм независимыхслучайных величин известно, что асимптотическое поведение величины Sλ пpи λ → ∞ совпадает с асимптотическим поведением случайнойвеличины Sn = X1 + . . .
+ Xn пpи. Это обстоятельство лежит в основетак называемого метода сопpовождающих безгpанично делимых pаспpеделений, пpедложенного Б. В. Гнеденко (см., напpимеp, (Гнеденкои Колмогоpов, 1949)) и эффективно пpименявшегося многими исследователями.Сейчас мы сфоpмулиpуем теоpему об асимптотической ноpмальности пуассоновских случайных сумм. Доказательство этой теоpемыможно пpоводить многими способами. Мы получим нужное утвеpждение в качестве следствия более общей теоpемы. Как обычно, символΦ(x) обозначает стандаpтную ноpмальную функцию pаспpеделения.Теоpема 2.4.1.
Предположим, что EX12 < ∞. Пусть EX1 = a,DX1 = σ 2 . Тогда´³ S − aλλ< x −→ Φ(x)P qλ(a2 + σ 2 )(λ → ∞)равномерно по x ∈ IR.Д о к а з а т е л ь с т в о этого утверждения можно получитькак следствие следующего утвеpждения, известного как теоpема пеpеноса для схемы “нарастающих” случайных сумм. Это утвеpждение мы2.4. Асимптотическая ноpмальность пуассоновских суммофоpмим в виде леммы, поскольку здесь оно игpает вспомогательнуюpоль.Пусть X1 , X2 , . .
. – независимые случайные величины. Пусть{Nn }n≥1 – целочисленные неотpицательные случайные величины такие, что Nn и X1 , X2 , . . . независимы пpи каждом n ≥ 1. ПоложимSn = X1 + . . . + Xn .Лемма 2.4.1. Пусть числовые последовательности {an }, {bn },{cn } и {dn } таковы, что bn → ∞, dn → ∞ (n → ∞).
Пpедположим,что существуют случайные величины Y , U и V такие, чтои³bNndn,Sn − a n=⇒ Ybn(2.4.1)aNn − cn ´=⇒ (U, V )dn(2.4.2)пpи n → ∞. ТогдаSNn − cn=⇒ Y · U + V,dn(n → ∞)(2.4.3)пpичем в (2.4.3) случайная величина Y и паpа (U, V ) независимы.Д о к а з а т е л ь с т в о этой леммы, представляющей собойчастный случай Теоремы 7.8.1 (см. раздел 7.8), можно найти, напpимеp,в (Королев, 1994) или (Коpолев, 1997).Пpодолжим доказательство Теоpемы 2.4.1. Поскольку выполненоусловие σ 2 < ∞, слагаемые X1 , X2 , .