exponenta event banner

wbllike

Вейбулл отрицательное логарифмическое правдоподобие

Синтаксис

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 для использования значения по умолчанию.

Отрицательное логарифмическое правдоподобие Вейбулла для неподцензурных данных

(logL) =−log∏i=1f (a, b 'xi) =−∑i=1nlogf (a, 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.

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