подобный бете

Бета отрицательная логарифмическая вероятность

Синтаксис

nlogL = betalike(params,data)
[nlogL,AVAR] = betalike(params,data)

Описание

nlogL = betalike(params,data) возвращает отрицание бета логарифмической функции правдоподобия для бета параметров a и b, заданный в векторном params и наблюдениях, заданных в вектор-столбце data.

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

[nlogL,AVAR] = betalike(params,data) также возвращает AVAR, который является асимптотической ковариационной матрицей отклонения оценок параметра, если значения в params являются оценками наибольшего правдоподобия. AVAR является инверсией информационной матрицы Фишера. Диагональные элементы AVAR являются асимптотическими отклонениями своих соответствующих параметров.

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

Примеры

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

r = betarnd(4,3,100,1);
[nlogl,AVAR] = betalike(betafit(r),r)
nlogl =

  -27.5996


AVAR =

    0.2783    0.1316
    0.1316    0.0867

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

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

Представлено до R2006a

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