Weibull отрицательная логарифмическая правдоподобность
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 использовать его значение по умолчанию.
Отрицательная логарифмическая правдоподобность Вейбула для данных без цензуры
где f - Weibull 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] Patel, J. K., C. H. Kapadia, and D. B. Owen. Справочник по статистическим распределениям. Нью-Йорк: Марсель Деккер, 1976.