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

Обзор

Бета распределение описывает семейство кривых, которые уникальны в этом, они являются ненулевыми только на интервале (0 1). Более общая версия функции присваивает параметры конечным точкам интервала.

Statistics and Machine Learning Toolbox™ обеспечивает несколько способов работать с бета распределением. Можно использовать следующие подходы, чтобы оценить параметры от выборочных данных, вычислить PDF, cdf, и icdf, сгенерировать случайные числа и т.д.

Соответствуйте объекту распределения вероятностей к выборочным данным или создайте объект распределения вероятностей с заданными значениями параметров. Смотрите Using BetaDistribution Objects для получения дополнительной информации.
Работа с вводом данных из матриц, таблиц и массивов набора данных с помощью функций распределения вероятностей. Смотрите Поддерживаемые Дистрибутивы для списка бета функций распределения.
В интерактивном режиме соответствуйте, исследуйте и сгенерируйте случайные числа от распределения с помощью приложения или пользовательского интерфейса.

Для получения дополнительной информации о каждой из этих опций смотрите Работу с Распределениями вероятностей.

Параметры

Бета распределение использует следующие параметры.

Параметр	Описание	Поддержка
`a`	Сначала сформируйте параметр	$a > 0$
`b`	Второй параметр формы	$b > 0$

Функция плотности вероятности

Определение

Функция плотности вероятности (PDF) бета распределения

$y = f (x | a, b) = \frac{1}{B (a, b)} x^{a - 1} {(1 - x)}^{b - 1} I_{[0, 1]} (x)$

где B (·) Бета-функция. Функция индикатора I _(0,1) (x) гарантирует, что только значения x в области значений (0,1) имеют ненулевую вероятность.

График

Этот график показывает, как изменение значения параметров изменяет форму PDF. Постоянный PDF (плоская строка) показывает, что стандартное равномерное распределение является особым случаем бета распределения, которое происходит когда a = b = 1.

X = 0:.01:1;
y1 = betapdf(X,0.75,0.75);
y2 = betapdf(X,1,1);
y3 = betapdf(X,4,4);

figure
plot(X,y1,'Color','r','LineWidth',2)
hold on
plot(X,y2,'LineStyle','-.','Color','b','LineWidth',2)
plot(X,y3,'LineStyle',':','Color','g','LineWidth',2)
legend({'a = b = 0.75','a = b = 1','a = b = 4'},'Location','NorthEast');
hold off

Отношение к другим дистрибутивам

Бета распределение имеет функциональное отношение с распределением t. Если Y является наблюдением от распределения t Студента со степенями свободы ν, то следующее преобразование генерирует X, который является распределенной бетой.

$X = \frac{}{12} + \frac{}{12} \frac{Y}{\sqrt{ν + Y^{2}}}$

Если Y ~t (v), то $X \sim β (\frac{ν}{2}, \frac{ν}{2})$

Это отношение используется, чтобы вычислить значения t cdf и обратной функции, а также генерации, t распределил случайные числа.

Кумулятивная функция распределения

Бета cdf совпадает с неполной бета-функцией.

Пример

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

Один популярный критерий совершенства должен максимизировать функцию правдоподобия. Вероятность имеет ту же форму как бета PDF. Но для PDF, параметры являются известными константами, и переменная является x. Функция правдоподобия инвертирует роли переменных. Здесь, демонстрационные значения (x's) уже наблюдаются. Таким образом, они - фиксированные постоянные. Переменные являются неизвестными параметрами. Оценка наибольшего правдоподобия (MLE) включает вычисление значений параметров, которые дают самую высокую вероятность, учитывая определенный набор данных.

Функциональный betafit возвращает MLEs и доверительные интервалы для параметров бета распределения. Вот пример с помощью случайных чисел от бета распределения с a = 5 и b = 0.2.

rng default  % For reproducibility
r = betarnd(5,0.2,100,1);
[phat, pci] = betafit(r)

phat = 1×2

    7.4911    0.2135

pci = 2×2

    5.0861    0.1744
   11.0334    0.2614

MLE для параметра a 7.4911, по сравнению с истинным значением 5. 95%-й доверительный интервал для a идет от 5,0861 до 11,0334, который не включает истинное значение. В то время как это - маловероятный результат, это действительно иногда происходит при оценке параметров распределения.

Так же MLE для параметра b 0.2135, по сравнению с истинным значением 0,2. 95%-й доверительный интервал для b идет от 0,1744 до 0,2614, который действительно включает истинное значение. В этом искусственном примере вы знаете “истинное значение”. В экспериментировании вы не делаете.

Смотрите также

BetaDistribution

Документация