Кумарасвами Функция плотности вероятности
Кумулятивная функция распределения
Параметры а > 0 { displaystyle a> 0 ,} (реальный) б > 0 { displaystyle b> 0 ,} (реальный) Поддержка Икс ∈ ( 0 , 1 ) { Displaystyle х в (0,1) ,} PDF а б Икс а − 1 ( 1 − Икс а ) б − 1 { displaystyle abx ^ {a-1} (1-x ^ {a}) ^ {b-1} ,} CDF 1 − ( 1 − Икс а ) б { displaystyle 1- (1-x ^ {a}) ^ {b}} Значить б Γ ( 1 + 1 а ) Γ ( б ) Γ ( 1 + 1 а + б ) { displaystyle { frac {b Gamma (1 + { tfrac {1} {a}}) Gamma (b)} { Gamma (1 + { tfrac {1} {a}} + b)} } ,} Медиана ( 1 − 2 − 1 / б ) 1 / а { displaystyle left (1-2 ^ {- 1 / b} right) ^ {1 / a}} Режим ( а − 1 а б − 1 ) 1 / а { displaystyle left ({ frac {a-1} {ab-1}} right) ^ {1 / a}} за а ≥ 1 , б ≥ 1 , ( а , б ) ≠ ( 1 , 1 ) { Displaystyle а geq 1, б geq 1, (а, b) neq (1,1)} Дисперсия (сложный - см. текст) Асимметрия (сложный - см. текст) Ex. эксцесс (сложный - см. текст) Энтропия ( 1 − 1 б ) + ( 1 − 1 а ) ЧАС б − пер ( а б ) { displaystyle left (1 ! - ! { tfrac {1} {b}} right) + left (1 ! - ! { tfrac {1} {a}} right) H_ { b} - ln (ab)}
В вероятность и статистика , то Двойное ограниченное распределение Кумарасвами это семья непрерывные распределения вероятностей определенная на интервале (0,1). Это похоже на Бета-распространение , но намного проще в использовании, особенно в симуляционных исследованиях, поскольку функция плотности вероятности , кумулятивная функция распределения и квантильные функции могут быть выражены в закрытая форма . Это распределение было первоначально предложено Пунди Кумарасвами [1] для переменных, ограниченных сверху и снизу с нулевой инфляцией. Это было распространено на инфляцию с обоими крайними значениями [0,1] дюйма.[2]
Характеристика
Функция плотности вероятности В функция плотности вероятности распределения Кумарасвами без учета инфляции
ж ( Икс ; а , б ) = а б Икс а − 1 ( 1 − Икс а ) б − 1 , где Икс ∈ ( 0 , 1 ) , { displaystyle f (x; a, b) = abx ^ {a-1} {(1-x ^ {a})} ^ {b-1}, { mbox {where}} x в (0,1),} и где а и б неотрицательны параметры формы .
Кумулятивная функция распределения В кумулятивная функция распределения является
F ( Икс ; а , б ) = ∫ 0 Икс ж ( ξ ; а , б ) d ξ = 1 − ( 1 − Икс а ) б . { Displaystyle F (x; a, b) = int _ {0} ^ {x} f ( xi; a, b) d xi = 1- (1-x ^ {a}) ^ {b} . } Квантильная функция Обратная кумулятивная функция распределения (функция квантиля) есть
F ( у ; а , б ) − 1 = ( 1 − ( 1 − у ) 1 б ) 1 а . { Displaystyle F (y; a, b) ^ {- 1} = (1- (1-y) ^ { frac {1} {b}}) ^ { frac {1} {a}}. } Обобщение на поддержку произвольного интервала В простейшей форме распределение имеет носитель (0,1). В более общем виде нормализованная переменная Икс заменяется несмещенной и немасштабированной переменной z куда:
Икс = z − z мин z Максимум − z мин , z мин ≤ z ≤ z Максимум . { displaystyle x = { frac {z-z _ { text {min}}} {z _ { text {max}} - z _ { text {min}}}}, qquad z _ { text {min} } leq z leq z _ { text {max}}. , !} Характеристики
Сырье моменты распределения Кумарасвами дают:[3] [4]
м п = б Γ ( 1 + п / а ) Γ ( б ) Γ ( 1 + б + п / а ) = б B ( 1 + п / а , б ) { displaystyle m_ {n} = { frac {b Gamma (1 + n / a) Gamma (b)} { Gamma (1 + b + n / a)}} = bB (1 + n / a) , б) ,} где B это Бета-функция а Γ (.) обозначает Гамма-функция . Дисперсия, перекос , и избыточный эксцесс можно рассчитать из этих сырых моментов. Например, разница составляет:
σ 2 = м 2 − м 1 2 . { displaystyle sigma ^ {2} = m_ {2} -m_ {1} ^ {2}.} В Энтропия Шеннона (в нац) распределения составляет:[5]
ЧАС = ( 1 − 1 а ) + ( 1 − 1 б ) ЧАС б − пер ( а б ) { displaystyle H = left (1 ! - ! { tfrac {1} {a}} right) + left (1 ! - ! { tfrac {1} {b}} right) H_ {b} - ln (ab)} где ЧАС я { displaystyle H_ {i}} это номер гармоники функция.
Отношение к бета-распределению
Распределение Кумарасвами тесно связано с бета-распределением.[6] Предположим, что Икс а, б Кумарасвами распространяет случайная переменная с параметрами а и б .Потом Икс а, б это а корень -й степени из подходящим образом определенной Бета-распределенной случайной величины. Более формально, Пусть Y 1, б обозначить Бета-версия распространена случайная величина с параметрами α = 1 { Displaystyle альфа = 1} и β = б { displaystyle beta = b} Между ними существует следующая связь. Икс а, б и Y 1, б .
Икс а , б = Y 1 , б 1 / а , { displaystyle X_ {a, b} = Y_ {1, b} ^ {1 / a},} с равным распределением.
п { Икс а , б ≤ Икс } = ∫ 0 Икс а б т а − 1 ( 1 − т а ) б − 1 d т = ∫ 0 Икс а б ( 1 − т ) б − 1 d т = п { Y 1 , б ≤ Икс а } = п { Y 1 , б 1 / а ≤ Икс } . { displaystyle operatorname {P} {X_ {a, b} leq x } = int _ {0} ^ {x} abt ^ {a-1} (1-t ^ {a}) ^ { b-1} dt = int _ {0} ^ {x ^ {a}} b (1-t) ^ {b-1} dt = operatorname {P} {Y_ {1, b} leq x ^ {a} } = operatorname {P} {Y_ {1, b} ^ {1 / a} leq x }.} Можно ввести обобщенные распределения Кумарасвами, рассматривая случайные величины вида Y α , β 1 / γ { displaystyle Y _ { alpha, beta} ^ {1 / gamma}} , с участием γ > 0 { displaystyle gamma> 0} и где Y α , β { displaystyle Y _ { alpha, beta}} обозначает бета-распределенную случайную величину с параметрами α { displaystyle alpha} и β { displaystyle beta} . Необработанный моменты этого обобщенного распределения Кумарасвами даются:
м п = Γ ( α + β ) Γ ( α + п / γ ) Γ ( α ) Γ ( α + β + п / γ ) . { Displaystyle m_ {N} = { гидроразрыва { Gamma ( alpha + beta) Gamma ( alpha + n / gamma)} { Gamma ( alpha) Gamma ( alpha + beta + n / gamma)}}.} Обратите внимание, что мы можем восстановить исходное значение моментов α = 1 { Displaystyle альфа = 1} , β = б { displaystyle beta = b} и γ = а { Displaystyle gamma = а} Однако, как правило, кумулятивная функция распределения не имеет решения в замкнутой форме.
Связанные дистрибутивы
Если Икс ∼ Кумарасвами ( 1 , 1 ) { displaystyle X sim { textrm {Kumaraswamy}} (1,1) ,} тогда Икс ∼ U ( 0 , 1 ) { Displaystyle X сим U (0,1) ,} Если Икс ∼ U ( 0 , 1 ) { Displaystyle X сим U (0,1) ,} (Равномерное распределение (непрерывное) ) тогда ( 1 − ( 1 − Икс ) 1 б ) 1 а ∼ Кумарасвами ( а , б ) { displaystyle {{ Big (} 1 - { left (1-X right)} ^ { tfrac {1} {b}} { Big)}} ^ { tfrac {1} {a}} sim { textrm {Kumaraswamy}} (a, b) ,} Если Икс ∼ Бета ( 1 , б ) { displaystyle X sim { textrm {Beta}} (1, b) ,} (Бета-распространение ) тогда Икс ∼ Кумарасвами ( 1 , б ) { displaystyle X sim { textrm {Kumaraswamy}} (1, b) ,} Если Икс ∼ Бета ( а , 1 ) { Displaystyle X sim { textrm {бета}} (а, 1) ,} (Бета-распространение ) тогда Икс ∼ Кумарасвами ( а , 1 ) { displaystyle X sim { textrm {Kumaraswamy}} (а, 1) ,} Если Икс ∼ Кумарасвами ( а , 1 ) { displaystyle X sim { textrm {Kumaraswamy}} (а, 1) ,} тогда ( 1 − Икс ) ∼ Кумарасвами ( 1 , а ) { Displaystyle (1-Х) сим { textrm {Кумарасвами}} (1, а) ,} Если Икс ∼ Кумарасвами ( 1 , а ) { Displaystyle X sim { textrm {Kumaraswamy}} (1, а) ,} тогда ( 1 − Икс ) ∼ Кумарасвами ( а , 1 ) { displaystyle (1-X) sim { textrm {Kumaraswamy}} (а, 1) ,} Если Икс ∼ Кумарасвами ( а , 1 ) { displaystyle X sim { textrm {Kumaraswamy}} (а, 1) ,} тогда − бревно ( Икс ) ∼ Экспоненциальный ( а ) { displaystyle - log (X) sim { textrm {Exponential}} (а) ,} Если Икс ∼ Кумарасвами ( 1 , б ) { displaystyle X sim { textrm {Kumaraswamy}} (1, b) ,} тогда − бревно ( 1 − Икс ) ∼ Экспоненциальный ( б ) { displaystyle - log (1-X) sim { textrm {Exponential}} (b) ,} Если Икс ∼ Кумарасвами ( а , б ) { displaystyle X sim { textrm {Kumaraswamy}} (a, b) ,} тогда Икс ∼ GB1 ( а , 1 , 1 , б ) { Displaystyle X sim { textrm {GB1}} (а, 1,1, b) ,} , то обобщенное бета-распределение первого рода . пример
Примером использования распределения Кумарасвами является объем хранения резервуара емкости z чья верхняя граница z Максимум а нижняя граница равна 0, что также является естественным примером наличия двух инфляций, поскольку многие резервуары имеют ненулевые вероятности как для пустого, так и для полного состояний резервуара.[2]
Рекомендации
^ Кумарасвами, П. (1980). «Обобщенная функция плотности вероятности для дважды ограниченных случайных процессов». Журнал гидрологии . 46 (1–2): 79–88. Bibcode :1980JHyd ... 46 ... 79K . Дои :10.1016/0022-1694(80)90036-0 . ISSN 0022-1694 . ^ а б Fletcher, S.G .; Поннамбалам, К. (1996). «Оценка дебита коллектора и распределения запасов с использованием анализа моментов». Журнал гидрологии . 182 (1–4): 259–275. Bibcode :1996JHyd..182..259F . Дои :10.1016 / 0022-1694 (95) 02946-х . ISSN 0022-1694 . ^ Лемонте, Артур Дж. (2011). «Улучшенная оценка баллов для распределения Кумарасвами». Журнал статистических вычислений и моделирования . 81 (12): 1971–1982. Дои :10.1080/00949655.2010.511621 . ISSN 0094-9655 . ^ КРИБАРИ-НЕТО, ФРАНЦИСКО; САНТОС, ХССИКА (2019). «Раздутое распределение Кумарасвами» . Anais da Academia Brasileira de Ciências . 91 (2): e20180955. Дои :10.1590/0001-3765201920180955 . ISSN 1678-2690 . PMID 31141016 . ^ Михалович, Йозеф Виктор; Николс, Джонатан М .; Бухольц, Франк (2013). Справочник по дифференциальной энтропии . Чепмен и Холл / CRC. п. 100. ISBN 9781466583177 . ^ Джонс, М. (2009). «Распределение Кумарасвами: распределение бета-типа с некоторыми преимуществами управляемости». Статистическая методология . 6 (1): 70–81. Дои :10.1016 / j.stamet.2008.04.001 . ISSN 1572-3127 .
Дискретный одномерный с конечной опорой Дискретный одномерный с бесконечной поддержкой Непрерывный одномерный поддерживается на ограниченном интервале Непрерывный одномерный поддерживается на полубесконечном интервале Непрерывный одномерный поддерживается на всей реальной линии Непрерывный одномерный с поддержкой, тип которой варьируется Смешанная непрерывно-дискретная одномерная Многовариантный (совместный) Направленный Вырожденный и единственное число Семьи