Статистика приходит к нам на помощь при решении многих задач, например: когда нет возможности построить детерминированную модель, когда слишком много факторов или когда нам необходимо оценить правдоподобие построенной модели с учётом имеющихся данных. Отношение к статистике неоднозначное. Есть мнение, что существует три вида лжи: ложь, наглая ложь и статистика. С другой стороны, многие «пользователи» статистики слишком ей верят, не понимая до конца, как она работает: применяя, например, тест Стьюдента к любым данным без проверки их нормальности. Такая небрежность способна порождать серьёзные ошибки и превращать «поклонников» теста Стьюдента в ненавистников статистики. Попробуем поставить токи над i и разобраться, какие модели случайных величин должны использоваться для описания тех или иных явлений и какая между ними существует генетическая связь.

В первую очередь, данный материал будет интересен студентам, изучающим теорию вероятностей и статистику, хотя и «зрелые» специалисты смогут его использовать в качестве справочника. В одной из следующих работ я покажу пример использования статистики для построения теста оценки значимости показателей биржевых торговых стратегий.

В работе будут рассмотрены дискретные распределения:

  1. Бернулли;
  2. биномиальное;
  3. геометрическое;
  4. Паскаля (отрицательное биномиальное);
  5. гипергеометрическое;
  6. Пуассона,

а также непрерывные распределения:

  1. Гаусса (нормальное);
  2. хи-квадрат;
  3. Стьюдента;
  4. Фишера;
  5. Коши;
  6. экспоненциальное (показательное) и Лапласа (двойное экспоненциальное, двойное показательное);
  7. Вейбулла;
  8. гамма (Эрланга);
  9. бета.

В конце статьи будет задан вопрос для размышлений. Свои размышления по этому поводу я изложу в следующей статье.

Некоторые из приведённых непрерывных распределений являются частными случаями распределения Пирсона.

Дискретные распределения


Дискретные распределения используются для описания событий с недифференцируемыми характеристиками, определёнными в изолированных точках. Проще говоря, для событий, исход которых может быть отнесён к некоторой дискретной категории: успех или неудача, целое число (например, игра в рулетку, в кости), орёл или решка и т.д.

Описывается дискретное распределение вероятностью наступления каждого из возможных исходов события. Как и для любого распределения ( в том числе непрерывного) для дискретных событий определены понятия матожидания и дисперсии. Однако, следует понимать, что матожидание для дискретного случайного события — величина в общем случае нереализуемая как исход одиночного случайного события, а скорее как величина, к которой будет стремиться среднее арифметическое исходов событий при увеличении их количества.

В моделировании дискретных случайных событий важную роль играет комбинаторика, так как вероятность исхода события можно определить как отношение количества комбинаций, дающих требуемый исход к общему количеству комбинаций. Например: в корзине лежат 3 белых мяча и 7 чёрных. Когда мы выбираем из корзины 1 мяч, мы можем сделать это 10-ю разными способами (общее количество комбинаций), но только 3 варианта, при которых будет выбран белый мяч (3 комбинации, дающие требуемый исход). Таким образом, вероятность выбрать белый мяч: (распределение Бернулли).

Следует также отличать выборки с возвращением и без возвращения. Например, для описания вероятности выбора двух белых мячей важно определить, будет ли первый мяч возвращён в корзину. Если нет, то мы имеем дело с выборкой без возвращения (гипергеометрическое распределение) и вероятность будет такова: — вероятность выбрать белый мяч из начальной выборки умноженная на вероятность снова выбрать белый мяч из оставшихся в корзине. Если же первый мяч возвращается в корзину, то это выборка с возвращением (Биномиальное распределение). В этом случае вероятность выбора двух белых мячей составит .

наверх

Распределение Бернулли



(взято отсюда)

Если несколько формализовать пример с корзиной следующим образом: пусть исход события может принимать одно из двух значений 0 или 1 с вероятностями и соответственно, тогда распределение вероятности получения каждого из предложенных исходов будет называться распределение Бернулли:

Bin_{p,q}\left(x\right)=\begin{cases}q,x=0\\ p, x=1\end{cases}


По сложившейся традиции, исход со значением 1 называется «успех», а исход со значением 0 — «неудача». Очевидно, что получение исхода «успех или неудача» наступает с вероятностью .

Матожидание и дисперсия распределения Бернулли:

E\{Bin_{p,q}\}=p\ \ \ \ \left(1.1.2\right)

D\{Bin_{p,q}\}=pq=p\left(1-p\right)\ \ \ \ \left(1.1.3\right)


наверх

Биномиальное распределение



(взято отсюда)

Количество успехов в испытаниях, исход которых распределен по Бернулли с вероятностью успеха (пример с возвращением мячей в корзину), описывается биномиальным распределением:

B_{n,p}(k)=C^k_np^kq^{n-k}\ \ (1.2.1)

где C^k_n={n!\over{k!(n-k)!}} — число сочетаний из по .

По другому можно сказать, что биномиальное распределение описывает сумму из независимых случайных величин, умеющих распределение Бернулли с вероятностью успеха .
Матожидание и дисперсия:

E\{B_{n,p}\}=np\ \ (1.2.2)

D\{B_{n,p}\}=npq\ \ (1.2.3)

Биномиальное распределение справедливо только для выборки с возвращением, то есть, когда вероятность успеха остаётся постоянной для всей серии испытаний.

Если величины и имеют биномиальные распределения с параметрами и соответственно, то их сумма также будет распределена биномиально с параметрами .

наверх

Геометрическое распределение



(взято отсюда)

Представим ситуацию, что мы вытягиваем мячи из корзины и возвращаем обратно до тех пор, пока не будет вытянут белый шар. Количество таких операций описывается геометрическим распределением. Иными словами: геометрическое распределение описывает количество испытаний до первого успеха при вероятности наступления успеха в каждом испытании . Если подразумевается номер испытания, в котором наступил успех, то геометрическое распределение будет описываться следующей формулой:

Geom_p(n)=q^{n-1}p\ \  (1.3.1)

Матожидание и дисперсия геометрического распределения:

E\{Geom_p\}={1\over{p}}\ \ \ (1.3.2)

D\{Geom_p\}={q\over{p^2}}\ \ \ (1.3.3)

Геометрическое распределение генетически связано с экспоненциальным распределением, которое описывает непрерывную случайную величину: время до наступления события, при постоянной интенсивности событий. Геометрическое распределение также является частным случаем отрицательного биномиального распределения.

наверх

Распределение Паскаля (отрицательное биномиальное рспределение)



(взято отсюда)

Распределение Паскаля является обобщением геометрического распределения: описывает распределение количества неудач в независимых испытаниях, исход которых распределен по Бернулли с вероятностью успеха до наступления успехов в сумме. При , мы получим геометрическое распределение для величины .

NB_{r,p}(k)=C^k_{k+r-1}p^rq^k\ \ (1.4.1)

где C^k_n={n!\over{k!(n-k)!}} — число сочетаний из по .

Матожидание и дисперсия отрицательного биномиального распределения:

E\{NB_{r,p}\}={rq\over{p}}\ \ \ (1.4.2)

D\{NB_{r,p}\}={rq\over{p^2}}\ \ \ (1.4.3)

Сумма независимых случайных величин, распределённых по Паскалю, также распределена по Паскалю: пусть имеет распределение , а . Пусть также и независимы, тогда их сумма будет иметь распределение

наверх

Гипергеометрическое распределение



(взято отсюда)

До сих пор мы рассматривали примеры выборок с возвращением, то есть, вероятность исхода не менялась от испытания к испытанию.

Теперь рассмотрим ситуацию без возвращения и опишем вероятность количества успешных выборок из совокупности с заранее известным количеством успехов и и неудач (заранее известное количество белых и чёрных мячей в корзине, козырных карт в колоде, бракованных деталей в партии и т.д.).

Пусть общая совокупность содержит объектов, из них помечены как «1», а как «0». Будем считать выбор объекта с меткой «1», как успех, а с меткой «0» как неудачу. Проведём n испытаний, причём выбранные объектв больше не будут участвовать в дальнейших испытаниях. Вероятность наступления успехов будет подчиняться гипергеометрическому распределению:

HG_{N,D,n}(k)={C^k_DC^{n-k}_{N-D}\over{C^n_N}}\ \ (1.5.1)

где C^k_n={n!\over{k!(n-k)!}} — число сочетаний из по .

Матожидание и дисперсия:

E\{HG_{N,D,n}\}={nD\over{N}}\ \ (1.5.2)

D\{HG_{N,D,n}\}=n{D\over{N}}{N-D\over{N}}{N-n\over{N-1}}\ \ (1.5.3)

наверх

Распределение Пуассона



(взято отсюда)

Распределение Пуассона значительно отличается от рассмотренных выше распределений своей «предметной» областью: теперь рассматривается не вероятность наступления того или иного исхода испытания, а интенсивность событий, то есть среднее количество событий в единицу времени.

Распределение Пуассона описывает вероятность наступления независимых событий за время при средней интенсивности событий :

P_{\lambda,t}(k)={\left(\lambda t\right)^k\over{k!}}e^{\lambda t}\ \ \ (1.6.1)

Матожидание и дисперсия распределения Пуассона:

E\{P_{\lambda,t}\}=\lambda t\ \ \ (1.6.2)

D\{P_{\lambda,t}\}=\lambda t\ \ \ (1.6.3)

Дисперсия и матожидание распределения Пуассона тождественно равны.

Распределение Пуассона в сочетании с экспоненциальным распределением, описывающим интервалы времени между наступлениями независимых событий, составляют математическую основу теории надёжности.

наверх

Непрерывные распределения


Непрерывные распределения, в отличие от дискретных, описываются функциями плотности (распределения) вероятности , определёнными, в общем случае, на некоторых интервалах.

Если известна плотность вероятности для величины : и определено преобразование , то плотность вероятности для y может быть получена автоматически:

f_y(y)=f\left(g^{-1}(y)\right)\left|{dg^{-1}\over{dy}}(y)\right|\ \ \ (2.0.1)

при условии однозначности и дифференцируемости .

Плотность вероятности суммы случайных величин и () с распределениями и описывается свёрткой и :

h(z)=\int f(t)g(z-t)dt=(f*g)(z)\ \ \ (2.0.2)

Если распределение суммы случайных величин принадлежит к тому же распределению, что и слагаемые, такое распределение называется бесконечно делимым. Примеры бесконечно делимых распределений: нормальное, хи-квадрат, гамма, распределение Коши.

Плотность вероятности произведения случайных величин x и y () с распределениями и может быть вычислена следующим образом:

h(z)=\int f(t)g(z/t)dt\ \ \ (2.0.3)


Некоторые из приведённых ниже распределений являются частными случаями распределения Пирсона, которое, в свою очередь, является решением уравнения:

{df\over{dx}}(x)={a_0+a_1x\over{b_0+2b_1x+b_2x^2}}f(x)\ \ \ (2.0.4)

где и — параметры распределения. Известны 12 типов распределения Пирсона, в зависимости от значений параметров.

Распределения, которые будут рассмотрены в этом разделе, имеют тесные взаимосвязи друг с другом. Эти связи выражаются в том, что некоторые распределения являются частными случаями других распределений, либо описывают преобразования случайных величин, имеющих другие распределения.

На приведённой ниже схеме отражены взаимосвязи между некоторыми из непрерывных распределений, которые будут рассмотрены в настоящей работе. На схеме сплошными стрелками показано преобразование случайных величин (начало стрелки указывает на изначальное распределение, конец стрелки — на результирующее), а пунктирными — отношение обобщения (начало стрелки указывает на распределение, являющееся частным случаем того, на которое указывает конец стрелки). Для частных случаев распределения Пирсона над пунктирными стрелками указан соответствующий тип распределения Пирсона.

Предложенный ниже обзор распределений охватывает многие случаи, которые встречаются в анализе данных и моделировании процессов, хотя, конечно, и не содержит абсолютно все известные науке распределения.

наверх

Нормальное распределение (распределение Гаусса)



(взято отсюда)

Плотность вероятности нормального распределения с параметрами и описывается функцией Гаусса:
f(x)={1\over{\sigma \sqrt{2 \pi}}}e^{(x-\mu)^2\over{2\sigma^2}}\ \ \ (2.1.1)

Если и , то такое распределение называется стандартным.

Матожидание и дисперсия нормального распределения:

E\{N_{\mu,\sigma}\}=\mu\ \ \ (2.1.2)

D\{N_{\mu,\sigma}\}=\sigma^2\ \ \ (2.1.3)

Область определения нормального распределения — множество дествительных чисел.

Нормальное распределение является распределение Пирсона типа VI.

Сумма квадратов независимых нормальных величин имеет распределение хи-квадрат, а отношение независимых Гауссовых величин распределено по Коши.

Нормальное распределение является бесконечно делимым: сумма нормально распределенных величин и с параметрами и соответственно также имеет нормальное распределение с параметрами , где и .

Нормальное распределение хорошо моделирует величины, описывающие природные явления, шумы термодинамической природы и погрешности измерений.

Кроме того, согласно центральной предельной теореме, сумма большого количества независимых слагаемых одного порядка сходится к нормальному распределению, независимо от распределений слагаемых. Благодаря этому свойству, нормальное распределение популярно в статистическом анализе, многие статистические тесты рассчитаны на нормально распределенные данные.

На бесконечной делимости нормального распределении основан z-тест. Этот тест используется для проверки равенства матожидания выборки нормально распределённых величин некоторому значению. Значение дисперсии должно быть известно. Если значение дисперсии неизвестно и рассчитывается на основании анализируемой выборки, то применяется t-тест, основанный на распределении Стьюдента.

Пусть у нас имеется выборка объёмом n независимых нормально распределенных величин из генеральной совокупности со стандартным отклонением выдвинем гипотезу, что . Тогда величина z={\bar{X}-\mu\over{\sigma \sqrt{n}}} будет иметь стандартное нормальное распределение. Сравнивая полученное значение z с квантилями стандартного распределения можно принимать или отклонять гипотезу с требуемым уровнем значимости.

Благодаря широкой распространённости распределения Гаусса, многие, не очень хорошо знающие статистику исследователи забывают проверять данные на нормальность, либо оценивают график плотности распределения «на глазок», слепо полагая, что имеют дело с Гауссовыми данными. Соответственно, смело применяя тесты, предназначенные для нормального распределения и получая совершенно некорректные результаты. Наверное, отсюда и пошла молва про статистику как самый страшный вид лжи.

Рассмотрим пример: нам надо измерить сопротивления набора резистров некоторого номинала. Сопротивление имеет физическую природу, логично предположить, что распределение отклонений сопротивления от номинала будет нормальным. Меряем, получаем колоколообразную функцию плотности вероятности для измеренных значений с модой в окрестности номинала резистров. Это нормальное распределение? Если да, то будем искать бракованные резистры используя тест Стьюдента, либо z-тест, если нам заранее известна дисперсия распределения. Думаю, что многие именно так и поступят.

Но давайте внимательнее посмотрим на технологию измерения сопротивления: сопротивление определяется как отношение приложенного напряжения к протекающему току. Ток и напряжение мы измеряли приборами, которые, в свою очередь, имеют нормально распределенные погрешности. То есть, измеренные значения тока и напряжения — это нормально распределенные случайные величины с матожиданиями, соответствующими истинным значениям измеряемых величин. А это значит, что полученные значения сопротивления распределены по Коши, а не по Гауссу.

Распределение Коши лишь напоминает внешне нормальное распределение, но имеет более тяжёлые хвосты. А значит предложенные тесты неуместны. Надо строить тест на основании распределения Коши или вычислить квадрат сопротивления, который в данном случае будет иметь распределение Фишера с параметрами (1, 1).

к схеме
наверх

Распределение хи-квадрат



(взято отсюда)

Распределение описывает сумму квадратов случайных величин , каждая из которых распределена по стандартному нормальному закону :

\chi^2_n(x)={{\left(1\over 2 \right)}^{k\over 2}\over{\Gamma\left({k\over{2}}\right)}}x^{{k\over 2}-1}e^{-{x\over 2}}\ \ \ (2.2.1)


где — число степеней свободы, x=\sum\limits_{i=1}^n {X^2_i}.

Матожидание и дисперсия распределения :

E\{\chi^2_n\}=n\ \ \ (2.2.2)

D\{\chi^2_n\}=2n\ \ \ (2.2.3)

Область определения — множество неотрицательных натуральных чисел. является бесконечно делимым распределением. Если и — распределены по и имеют и степеней свободы соответственно, то их сумма также будет распределена по и иметь степеней свободы.

является частным случаем гамма-распределения (а следовательно, распределением Пирсона типа III) и обобщением экспоненциального распределения. Отношение величин, распределенных по распределено по Фишеру.

На распределении основан критерий согласия Пирсона. с помощью этого критерия можно проверять достоверность принадлежности выборки случайной величины некоторому теоретическому распределению.

Предположим, что у нас имеется выборка некоторой случайной величины . На основании этой выборки рассчитаем вероятности попадания значений в интервалов (). Пусть также есть предположение об аналитическом выражении распределения, в соответствие с которым, вероятности попадания в выбранные интервалы должны составлять . Тогда величины будут распределены по нормальному закону.

Приведем к стандартному нормальному распределению: ,
где m={1\over n}\sum \limits_{i=1}^n {D_i} и S=\sqrt {{1\over {n-1}}\sum \limits_{i=1}^n {D_i^2}}.

Полученные величины имеют нормальное распределение с параметрами (0, 1), а следовательно, сумма их квадратов распределена по с степенью свободы. Снижение степени свободы связано с дополнительным ограничением на сумму вероятностей попадания значений в интервалы: она должна быть равна 1.

Сравнивая значение с квантилями распределения можно принять или отклонить гипотезу о теоретическом распределении данных с требуемым уровнем значимости.

к схеме
наверх

Распределение Стьюдента (t-распределение)



(взято отсюда)

Распределение Стьюдента используется для проведения t-теста: теста на равенство матожидания выборки стандартно нормально распределённых случайных величин некоторому значению, либо равенства матожиданий двух нормальных выборок с одинаковой дисперсией (равенство дисперсий необходимо проверять f-тестом). Распределение Стьюдента описывает отношение нормально распределённой случайной величины к величине, распределённой по хи-квадрат.

T-тест является аналогом z-теста для случая, когда дисперсия или стандартное отклонение выборки неизвестно и должно быть оценено на основании самой выборки.

Рассмотрим пример проверки равенства матожидания нормальной выборки некоторому значению: пусть нам дана выборка нормальных величин объёмом n из некоторой генеральной совокупности, выдвинем и проверим гипотезу о том, что матожидание этой совокупности равно .

Рассчитаем величину . Эта величина будет иметь распределение хи-квадрат. Тогда величина будет иметь распределение Стьюдента c степенью свободы, где:

T_{n}(x)={\Gamma \left({n+1 \over 2}\right) \over {\sqrt{n \pi}\Gamma \left({n \over 2}\right)\left(1+{x^2 \over n}\right)^{n+1 \over 2}}}\ \ \ (2.3.1)


где — гамма-функция Эйлера.

Полученное значение можно сравнивать с квантилями распределения Стьюдента и принимать либо отклонять гипотезу о равенстве маотожидания значению с требуемым уровнем значимости.

Матожидание и дисперсия распределения Стьюдента:

E\{T_{n}\}=0\ \ \ (2.3.2)

D\{T_{n}\}={n \over {n-2}}\ \ \ (2.3.3)

при .

к схеме
наверх

Распределение Фишера



(взято отсюда)

Пусть и независимые случайные величины, имеющие распределение хи-квадрат со степенями свободы и соответственно. Тогда величина будет иметь распределение Фишера со степенями свободы , а величина — распределение Фишера со степенями свободы .
Распределение Фишера определено для действительных неотрицательных аргументов и имеет плотность вероятности:

F_{n_1,n_2}(x)={\sqrt{ (n_1x)^{n_1}n_2^{n_2}\over {(n_1x+n_2)^{n_1+n_2}}} \over {xB\left({n_1 \over 2},{n_2 \over 2} \right)}}\ \ \ (2.4.1)

Матожидание и дисперсия распределения Фишера:

E\{F_{n_1,n_2}\}={n_2 \over {n_2-2}}\ \ \ (2.4.2)

D\{F_{n_1,n_2}\}={2n_2^2(n_1+n_2-2) \over {n_1(n_2-2)^2(n_2-4)}}\ \ \ (2.4.3)

Матожидание определено для , а диспересия — для .

На распределении Фишера основан ряд статистических тестов, таких как оценка значимости параметров регрессии, тест на гетероскедастичность и тест на равенство дисперсий нормальных выборок (f-тест, следует отличать от точного теста Фишера).

F-тест: пусть имеются две независимые выборки и нормально распределенных данных объёмами и соответственно. Выдвинем гипотезу о равенстве дисперсий выборок и проверим её статистически.

Рассчитаем величину . Она будет иметь распределение Фишера со степенями свободы .

Сравнивая значение с квантилями соответствующего распределения Фишера, мы можем принять или отклонить гипотезу о равенстве дисперсий выборок с требуемым уровнем значимости.

к схеме
наверх

Распределение Коши



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

C_{x_0,\gamma}(x)={1\over{\pi \gamma \left(1+\left({x-x_0\over {\gamma}} \right)^2 \right)}}\ \ \ (2.5.1)


Распределение Коши является бесконечно делимым: сумма независимых случайных величин, распределённых по Коши, также распределена по Коши.

к схеме
наверх


Экспоненциальное (показательное) распределение и распределение Лапласа (двойное экспоненциальное, двойное показательное)



(взято отсюда)

Экспоненциальное распределение описывает интервалы времени между независимыми событиями, происходящими со средней интенсивностью . Количество наступлений такого события за некоторый отрезок времени описывается дискретным распределением Пуассона. Экспоненциальное распределение вместе с распределением Пуассона составляют математическую основу теории надёжности.

Кроме теории надёжности, экспоненциальное распределение применяется в описании социальных явлений, в экономике, в теории массового обслуживания, в транспортной логистике — везде, где необходимо моделировать поток событий.

Экспоненциальное распределение является частным случаем распределения хи-квадрат (для n=2), а следовательно, и гамма-распределения. Так-как экспоненциально распределённая величина является величиной хи-квадрат с 2-мя степенями свободы, то она может быть интерпретирована как сумма квадратов двух независимых нормально распределенных величин.

Кроме того, экспоненциальное распределение является честным случаем распределения Вейбулла.

Дискретный вариант экспоненциального распределения — это геометрическое распределение.

Плотность вероятности экспоненциально распределения:

E_\lambda(x)=\lambda e^{-\lambda x}\ \ \ (2.6.1)

определена для неотрицательных действительных значений .

Матожидание и дисперсия экспоненциального распределения:

E\{E_\lambda\}={1 \over \lambda} \ \ \ (2.6.2)

E\{E_\lambda\}={1 \over \lambda^2} \ \ \ (2.6.3)

Если функцию плотности вероятностей экспоненциального распределения отразить зеркально в область отрицательных значений, то есть, заменить на , то получится распределение Лапласа, также называемое двойным экспоненциальным или двойным показательным.


(взято отсюда)

Для большего обобщения, вводится параметр сдвига, смещающий центр «соединения» левой и правой частей распределения вдоль оси абсцисс. В отличие от экспоненциального, распределение Лапласа, определено на всей действительной числовой оси.

L_{\alpha,\beta}(x)={\alpha \over 2}e^{-\alpha \left|x-\beta\right|} \ \ \ (2.6.4)

где — параметр масштаба, а — параметр сдвига.

Матожидание и дисперсия:
E\{L_{\alpha, \beta}\}=\beta\ \ \ (2.6.5)

D\{L_{\alpha, \beta}\}={2 \over {\alpha^2}}\ \ \ (2.6.6)

Благодаря более тяжёлым хвостам, чем у нормального распределения, распределение Лапласа используется для моделирования некоторых видов погрешностей измерения в энергетике, а также находит применение в физике, экономике, финансовой статистике, телекоммуникации и т.д.

к схеме
наверх

Распределение Вейбулла



(взято отсюда)

Распределение Вейбулла описывается функцией плотности вероятности следующего вида:

W_{k, \lambda}(x)={k \over {\lambda}}\left({x \over {\lambda}}\right)^{k-1}e^{-\left({x \over {\lambda}}\right)^k}\ \ \ (2.7.1)

где ( > 0)- интенсивность событий (аналогично параметру экспоненциального распределения), а — показатель нестационарности (). При , распределение Вейбулла вырождается в экспоненциальное распределение, а в остальных случаях описывает поток независимых событий с нестационарной интенсивностью. При моделируется поток событий с растущей со временем интенсивностью, а при — со снижающейся. Область определения функции распределения плотности вероятностей: неотрицательные действительные числа.

Таким образом, распределение Вейбулла — обобщение экспоненциального распределения на случай нестационарной интенсивности событий. Используется в теории надёжности, моделировании процессов в технике, в прогнозировании погоды, в описании процесса измельчения и т.д.

Матожидание и дисперсия распределения Вейбулла:

E\{W_{k, \lambda}\}=\lambda \Gamma\left(1 + {1 \over k} \right)\ \ \ (2.7.2)

D\{W_{k, \lambda}\}=\lambda^2 \left(\Gamma\left(1 + {2 \over k} \right ) - \Gamma\left(1 + {1 \over k} \right )^2\right)\ \ \ (2.7.3)

где — гамма-функция Эйлера.

к схеме
наверх

Гамма-распределение (распределение Эрланга)



(взято отсюда)

Гамма-распределение является обобщением рапсределения хи-квадрат и, соответственно, экспоненциального распределения. Суммы квадратов нормально распределённых величин, а также суммы величин распределённых по хи-квадрат и по экспоненциальному распределению будут иметь гамма-распределение.

Гамма-распределение является распределением Пирсона III рода. Область определения гамма-распределения — натуральные неотрицательные числа.

Гамма-распределение определяется двумя неотрицательными параметрами — число степеней свободы (при целом значении степеней свободы, гамма-распределение называется распределением Эрланга) и коэффициент масштаба .

Гамма-распределение является бесконечно делимым: если величины и имеют распределения и соответсвенно, то величина будет иметь распределение

G_{k,\theta}(x)=x^{k-1}{e^{-{x \over \theta}}\over \Gamma(k) \theta^k}\ \ \ (2.8.1)

где — гамма-функция Эйлера.

Матожидание и дисперсия:

E\{G_{k,\theta}\}=k\theta\ \ \ (2.8.2)

D\{G_{k,\theta}\}=k\theta^2\ \ \ (2.8.3)

Гамма распределение широко применяется для моделирования сложных потоков событий, сумм временных интервалов между событиями, в экономике, теории массового обслуживания, в логистике, описывает продолжительность жизни в медицине. Является своеобразным аналогом дискретного отрицательного биномиального распределения.

к схеме
наверх

Бета-распределение



(взято отсюда)

Бета-распределение описывает долю суммы двух слагаемых, приходящуюся на каждое из них, если слагаемые являются случайными величинами, имеющими гамма-распределение. То есть, если величины и имеют гамма-распределение, величины и будут иметь бета-распределение.

Очевидно, что область определения бета-распределения . Бета-распределение является распределение Пирсона I типа.

B_{\alpha,\beta}={x^{\alpha -1}(1-x)^{\beta - 1}\over{B(\alpha, \beta)}}\ \ \ (2.9.1)

где параметры и — положительные натуральные числа, — бета-функция Эйлера.

Матожидание и дисперсия:

E\{B_{\alpha, \beta}\}={\alpha \over{\alpha+\beta}}\ \ \ (2.9.2)

D\{B_{\alpha, \beta}\}={\alpha \beta \over{(\alpha+\beta)^2(\alpha+\beta+1)}}\ \ \ (2.9.3)

к схеме
наверх

Вместо заключения


Мы рассмотрели 15 распределений вероятности, которые, на мой взгляд, охватывают большинство наиболее популярных приложений статистики.

Напоследок, небольшое домашнее задание: для оценки надёжности биржевых торговых систем используется такой показатель как профит-фактор. Профит-фактор рассчитывается как отношение суммарного дохода к суммарному убытку. Очевидно, что для системы, приносящей доход, профит-фактор больше единицы, и чем его значение выше, тем система надёжнее.

Вопрос: какое распределение имеет значение профит-фактора?

Свои размышления по этому поводу я изложу в следующей статье.

P.S. Если Вы захотите сослатья на нумерованные формулы из этой статьи, то можете использовать такую сслыку: ссылка_на_статью#x_y_z, где (x.y.z)- номер формулы, на которую Вы ссылаетесь.
Поделиться с друзьями
-->

Комментарии (27)


  1. GeMir
    30.09.2016 11:03
    +1

    Неплохо бы статью перед публикацией показать корректору, который почистит текст от остатков «(3/10) х (2/9)» и «Г(x)» (раз уж в большинстве случаев используется красивый TeX). Нумерация формул кажется лишней, а сама статья скорее похожа на цитату из скрипта лекции по вероятностному исчислению. Справочников по теме предостаточно и как студенты, так и «зрелые специалисты», думается мне, первым делом обратятся к ним, а не к поиску по Хабру.


    1. JamaGava
      30.09.2016 13:45

      Спасибо за комментарий. По поводу формул — согласен, буду дорабатывать.
      А что касается тематики и нумерации формул: эта статья «база» для дальнейших более IT-шных статей по анализу данных (чтобы ставить ссылки на конкретные формулы, собранные вместе).


      1. kxx
        30.09.2016 14:44
        +1

        И есть небольшая неточность: биномиальный коэффициент записывается обычно либо как , либо как , но не как .


        1. JamaGava
          30.09.2016 23:30

          Это действительно ляп с моей стороны… Спасибо.
          Уже исправлено.


    1. parpalak
      01.10.2016 12:37
      +1

      Картинки, кстати, тоже можно сделать красивыми:



      Вот исходник:

      \begin{tikzpicture}[mark=*,mark size=1,only marks]
      \begin{axis}[
          yticklabel style={/pgf/number format/fixed},
          ymax=0.25,
          width=10cm,
          declare function={binom(\k, \n, \p)=(\p^(\k))*((1-\p)^(\n-\k))*factorial(\n)/(factorial(\n-\k)*factorial(\k));}
      ]
      \addlegendimage{blue}
      \addlegendimage{green}
      \addlegendimage{red}
      \foreach \n/\p/\c in {20/0.5/blue, 20/0.7/green, 40/0.5/red} {
          \foreach \k in {0,...,\n} {
              \edef\temp{\noexpand
                  \addplot[\c] coordinates {(\k,{binom(\k, \n, \p)})};
              }\temp
          }
          \edef\temp{\noexpand
              \addlegendentry{$p=\p, n=\n$}
          }\temp
      }
      \end{axis}
      \end{tikzpicture}
      


      1. tomzarubin
        02.10.2016 14:12
        +2

        Непонятно почему вы считаете свой график лучше представленных выше. Точки тут неуместны— чтобы понять о чём график(а потом и форму распределения), надо напрячься.


        1. parpalak
          02.10.2016 14:38
          +2

          Согласен, конечно. Мы тут обсуждали оформление, а не содержание.

          В тексте была откуда-то скопированная png-картинка. Моя svg-картинка сделана в латехе с пакетом tikz. Способ не без проблем: человеку без опыта тяжело сразу готовить такие картинки. Зато исходник можно править на лету, не перерисовывая картинку. Я вот за 5 минут поменял точки на столбцы, и диаграмма стала понятнее:

          image


          1. tomzarubin
            04.10.2016 17:04
            +1

            Постойте, так про оформление и речь. Пример выше, конечно, чище и понятней.


            1. JamaGava
              04.10.2016 20:05

              Такие диаграммы красивее. Я взял материал со страниц википедии (в тексте приведены ссылки). Если для кого-нибудь действительно совсем несложно переделать графики — Вы можете улучшить Вики.
              Сам бы занялся, но не владею технологией построения настолько красивых диаграмм :)
              (Хотя, уже захотелось освоить)


              1. parpalak
                04.10.2016 21:11
                +1

                Переделать — не сложно, но и не просто :)


                Хотите освоить — посмотрите введение к официальной документации tikz. В нем последовательным усложнением строятся полноценные примеры графиков и диагамм. Это хорошая отправная точка.


                1. JamaGava
                  04.10.2016 23:13

                  Спасибо


              1. tomzarubin
                06.10.2016 09:10
                +1

                Я не в претензии)
                Кстати, в R есть прекрасные библиотеки. Тот же «классический уже» ggplot2 или seaborn + несколько библиотек интерактивных графиков.
                Ну и R, по ощущениям, очень прекрасен для статистики и EDA. Самое главное— побороть непонимание синтаксиса самого R.


  1. temas
    30.09.2016 12:47

    Интересно, но ожидал больше примеров. (В дискретных распределениях примеры были)


    1. JamaGava
      30.09.2016 13:43

      Большой пример скоро будет в отдельной статье


  1. Lelushak
    30.09.2016 18:48
    -1

    Зачем захламлять теги отдельным названием каждого распределения? У них другое предназначение.


  1. ratatosk
    30.09.2016 19:03
    +1

    Вот очень хорошая вводная статья про виды распределений: Common Probability Distributions: The Data Scientist’s Crib Sheet.


  1. JamaGava
    01.10.2016 14:52

    Название действительно провокационное, для привлечения внимания. В самой статье я пишу, что предложенные распределения используются в «наиболее часто встречающихся задачах».

    Что касается «дз» это тизер готовящейся к публикации работы. Идея такова: мы не можем знать всё обо всех стратегиях, но может поступить по аналогии с тем, как строятся другие стат.тесты, а именно, построить распределение профитфактора для системы, торгующей случайно. Тогда значение профитфактора реальной системы должно быть таковым, чтобы случайное достижения такого значения являлось маловероятным.


  1. smxfem
    01.10.2016 14:59
    +1

    Название статьи провокационное. Получается, что сингулярные распределения — бесполезная альтернатива и создание этих распределений не продиктовано непригодностью описанных тут распределений к некоторым задачам — «случаям жизни»?

    Домашнее задание смахивает на задачку с подвохом. Требуется смоделировать множество биржевых систем (тут требуются знания предметной области на уровне бога) и напрямую получить распределения аналитически или откуда-то взять статистику по биржевым системам, подобрать для нескольких дающих надежду распределений оптимальные параметры (оптимизационная задача) и выбрать среди полученных конкретных распределений наиболее точные?


  1. ShashkovS
    03.10.2016 00:37
    +4

    Я просто оставлю это здесь :)
    image

    Ну, и с подробностями: http://www.math.wm.edu/~leemis/chart/UDR/UDR.html


    1. JamaGava
      03.10.2016 01:19
      +1

      Вот это действительно: «Вау!» :)
      Спасибо!


      1. JamaGava
        03.10.2016 01:38

        У меня вообще подозрение, что почти все непрерывные распределения должны укладываться в форму , где и — полиномы. Либо являться подстановкой функции (например, нецелой степени либо логарифма) от вместо аргумента в эту формулу.


        Кто-нибудь встречал подобные обобщение?


        1. smxfem
          08.10.2016 10:33
          +2

          Есть обобщенное нормальное и обобщенное гиперболическое распределения, которые содержат многие распределения. Но обобщенный вид затрудняет оценку параметров (не критично), а полу-гуманитарные критерии применимости распределений к реальным задачам слабо распространены на такие виды. Попробуйте убедить, предположим, медиков, что какое-то обобщенное распределение лучше подходит для описания некоторой случайной величины, когда 100 лет уже её описывают в обширных исследованиях, которые считаются фундаментальными, более простым частным распределением. Хотя, если, после оценки, вклад всех параметров в распределение будет существенным, то никто не отвертится и придется принять обобщение.
          Знакомый психолог публиковала в зарубежном журнале статью, где математических претензий не было, а отфутболивали именно из-за отсутствия обоснования применимости использованных распределений, пока не были найдены необходимые подтверждающие ссылки.


    1. kxx
      04.10.2016 22:23
      +1

      Граф впечатляет. Даже не подозревал о некоторых распределениях.


    1. Leo5700
      05.10.2016 01:27

      Ближайшие полчаса путешествую по графу )


  1. JamaGava
    05.10.2016 09:35

    Студентам буду на зачёт задание давать: воспроизвести граф по памяти :)


  1. nickolaym
    06.10.2016 02:00
    +1

    Для плотности f(x) и y=g(x), — формула плотности fy(y) справедлива, только если g-1(y) однозначна. То есть, g(x) биективна. А дифференцируемость g(x) не только недостаточное, но и необязательное условие: легко придумать кучу примеров как с разрывными отображениями
    — например, пила вдоль линии y = int(x)-frac(x)
    так и дифференцируемые отображения, для которых обратные не дифференцируемы
    — например, с седловой точкой, y = x3 (при том, что оно биективно)


    1. JamaGava
      06.10.2016 11:40

      Согласен с Вашим, замечанием: дифференцируемой должна быть обратная к g(x) (опечатка).