Вейбулл отрицательное логарифмическое правдоподобие
nlogL = wbllike(params,data)
[logL,AVAR] = wbllike(params,data)
[...] = wbllike(params,data,censoring)
[...] = wbllike(params,data,censoring,freq)
nlogL = wbllike(params,data) возвращает логарифмическое правдоподобие Вейбулла. params(1) - параметр масштаба, A, и params(2) - параметр формы, B.
[logL,AVAR] = wbllike(params,data) также возвращает AVAR, которая является асимптотической матрицей дисперсии-ковариации оценок параметров, если значения в params являются оценками максимального правдоподобия. AVAR является обратной информационной матрицей Фишера. Диагональные элементы AVAR являются асимптотическими дисперсиями их соответствующих параметров.
[...] = wbllike(params,data,censoring) принимает логический вектор, censoring, того же размера, что и data, что равно 1 для наблюдений, которые подвергаются правой цензуре, и 0 для наблюдений, которые наблюдаются точно.
[...] = wbllike(params,data,censoring,freq) принимает частотный вектор, freq, того же размера, что и data. freq обычно содержит целочисленные частоты для соответствующих элементов в data, но может содержать любые неотрицательные значения. Пройти внутрь [] для censoring для использования значения по умолчанию.
Отрицательное логарифмическое правдоподобие Вейбулла для неподцензурных данных
, b' xi)
где f - pdf Вейбулла.
wbllike является полезной функцией для оценки максимального правдоподобия.
Этот пример продолжает пример из wblfit.
r = wblrnd(0.5,0.8,100,1); [logL, AVAR] = wbllike(wblfit(r),r) logL = 47.3349 AVAR = 0.0048 0.0014 0.0014 0.0040
[1] Патель, Дж. К., К. Х. Кападия и Д. Б. Оуэн. Справочник по статистическим распределениям. Нью-Йорк: Марсель Деккер, 1976.