betafit

Оценки параметров

Синтаксис

phat = betafit(data)
[phat,pci] = betafit(data,alpha)

Описание

phat = betafit(data) вычисляет максимальные оценки правдоподобия параметров a и b бета распределения data и возвращает вектор-столбец, содержащую оценки a и b, где beta cdf задается как

F(x|a,b)=1B(a,b)0xta1(1t)b1dt

и B (·) является Бета-функцией. Элементы data должен находиться в открытом интервале (0, 1), где задано бета- распределение. Однако иногда также необходимо подгонять бета- распределение к данным, которые включают точные нули или таковые. Для таких данных функция бета-правдоподобия неограниченна, и стандартная максимальная оценка правдоподобия невозможна. В этом случае betafit максимизирует измененную вероятность, которая включает нули или таковые путем их обработки, как если бы они были значениями, которые были подвергнуты цензуре влево в sqrt(realmin) или с правой цензурой в 1- eps/ 2, соответственно .

[phat,pci] = betafit(data,alpha) возвращает доверительные интервалы параметров a и b в матрице 2 на 2 pci. Первый столбец матрицы содержит нижнюю и верхнюю доверительные границы для параметра a, а второй столбец содержит доверительные границы для параметра b. Необязательный входной параметр alpha - значение в области значений [0, 1], задающее ширину доверительных интервалов. По умолчанию alpha является 0.05, что соответствует 95% доверительным интервалам. Доверительные интервалы основаны на нормальном приближении для распределения журналов оценок параметров.

Примеры

Этот пример генерирует 100 бета распределенных наблюдений. Истинные параметры a и b равны 4 и 3, соответственно. Сравните их со значениями, возвращенными в p по бета- подгонке. Обратите внимание, что столбцы ci обе скобки - истинные параметры.

data = betarnd(4,3,100,1);
[p,ci] = betafit(data,0.01)
p =
    5.5328    3.8097
ci =
    3.6538    2.6197
    8.3781    5.5402

Ссылки

[1] Хан, Джеральд Дж., и С. С. Шапиро. Статистические модели в технике. Hoboken, NJ: John Wiley & Sons, Inc., 1994, p. 95.

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ MATLAB ®

.
Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте