Гамма-отрицательное логарифмическое правдоподобие
nlogL = gamlike(params,data)
[nlogL,AVAR] = gamlike(params,data)
nlogL = gamlike(params,data) возвращает отрицательное значение гамма-логарифмической правдоподобности параметров, params, учитывая data. params(1)=A, параметры формы и params(2)=B, параметры масштабирования. Параметры в params все должны быть положительными
[nlogL,AVAR] = gamlike(params,data) также возвращает AVAR, которая является асимптотической матрицей дисперсии-ковариации оценок параметров, когда значения в params являются оценками максимального правдоподобия. AVAR является обратной информационной матрицей Фишера. Диагональные элементы AVAR являются асимптотическими дисперсиями их соответствующих параметров.
[...] = gamlike(params,data,censoring) принимает логический вектор того же размера, что и data это 1 для наблюдений, которые подвергаются правой цензуре, и 0 для наблюдений, которые наблюдаются точно.
[...] = gamfit(params,data,censoring,freq) принимает частотный вектор того же размера, что и data. freq обычно содержит целочисленные частоты для соответствующих элементов в data, но может содержать любые неотрицательные значения.
gamlike является полезной функцией для оценки максимального правдоподобия гамма-распределения. С тех пор gamlike возвращает отрицательную функцию гамма-логарифмического правдоподобия, минимизируя gamlike использование fminsearch это то же самое, что максимизировать вероятность.
Вычислить отрицательную логарифмическую вероятность оценок параметров, вычисленных gamfit функция:
a = 2; b = 3; r = gamrnd(a,b,100,1); [nlogL,AVAR] = gamlike(gamfit(r),r) nlogL = 267.5648 AVAR = 0.0788 -0.1104 -0.1104 0.1955