lognlike

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

Описание

пример

nlogL = lognlike(params,x) возвращает логарифмически нормальную отрицательную логарифмическую правдоподобность параметров распределения (params) учитывая выборочные данные (x). params(1) и params(2) среднее и стандартное отклонение логарифмических значений, соответственно.

nlogL = lognlike(params,x,censoring) задает ли каждое значение в x подвергается цензуре правом или нет. Используйте логический векторный censoring в котором 1 указывает на наблюдения, которые подвергаются цензуре правом, и 0 указывает на наблюдения, которые полностью наблюдаются.

nlogL = lognlike(params,x,censoring,freq) задает частоту или веса наблюдений. Задавать freq не задавая censoring, можно передать [] для censoring.

пример

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

Примеры

свернуть все

Найдите MLEs набора данных с цензурированием при помощи mle, и затем найдите отрицательную логарифмическую правдоподобность MLEs при помощи lognlike.

Сгенерируйте 1 000 случайных чисел от логарифмически нормального распределения параметрами 5 и 2.

rng('default') % For reproducibility
n = 1000; % Number of samples
x = lognrnd(5,2,[n,1]);

Найдите MLEs для параметров распределения (среднее и стандартное отклонение логарифмических значений) при помощи mle.

phat = mle(x,'distribution','LogNormal')
phat = 1×2

    4.9347    1.9969

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

nlogL = lognlike(phat,x)
nlogL = 7.0453e+03

Найдите оценки наибольшего правдоподобия (MLEs) логарифмически нормальных параметров распределения, и затем найдите доверительный интервал соответствующего cdf значения.

Сгенерируйте 1 000 случайных чисел от логарифмически нормального распределения параметрами 5 и 2.

rng('default') % For reproducibility
n = 1000; % Number of samples
x = lognrnd(5,2,n,1);

Найдите MLEs для параметров распределения (среднее и стандартное отклонение логарифмических значений) при помощи mle.

phat = mle(x,'distribution','LogNormal')
phat = 1×2

    4.9347    1.9969

muHat = phat(1);
sigmaHat = phat(2);

Оцените ковариацию параметров распределения при помощи lognlike. Функциональный lognlike возвращает приближение в асимптотическую ковариационную матрицу, если вы передаете MLEs, и выборки раньше оценивали MLEs.

[~,pCov] = lognlike(phat,x)
pCov = 2×2

    0.0040   -0.0000
   -0.0000    0.0020

Найдите cdf значение в 0,5 и его 95%-й доверительный интервал.

[p,pLo,pUp] = logncdf(0.5,muHat,sigmaHat,pCov)
p = 0.0024
pLo = 0.0016
pUp = 0.0037

p cdf значение логарифмически нормального распределения параметрами muHat и sigmaHat. Интервал [pLo,pUp] 95%-й доверительный интервал cdf, оцененного в 0,5, рассматривая неопределенность в muHat и sigmaHat использование pCov. 95% доверительного интервала означают вероятность что [pLo,pUp] содержит истинное cdf значение, 0.95.

Входные параметры

свернуть все

Логарифмически нормальные параметры распределения, заданные как вектор двух числовых значений. params(1) и params(2) среднее и стандартное отклонение логарифмических значений, соответственно. params(2) mustBePositive.

Пример: [0,1]

Типы данных: single | double

Выборочные данные, заданные как вектор.

Типы данных: single | double

Индикатор для цензурирования каждого значения в x, заданный как логический вектор одного размера с x. Используйте 1 в наблюдениях, которые подвергаются цензуре правом и 0 для наблюдений, которые полностью наблюдаются.

Значением по умолчанию является массив 0s, означая, что все наблюдения полностью наблюдаются.

Типы данных: логический

Частота или веса наблюдений, заданных как неотрицательный вектор, который одного размера с x. freq входной параметр обычно содержит неотрицательное целое число, значит соответствующие элементы в x, но может содержать любые неотрицательные значения.

Чтобы получить взвешенную отрицательную логарифмическую правдоподобность для набора данных с цензурированием, задайте веса наблюдений, нормированных к количеству наблюдений в x.

Значением по умолчанию является массив 1 с, означая одно наблюдение на элемент x.

Типы данных: single | double

Выходные аргументы

свернуть все

Отрицательное значение логарифмической правдоподобности параметров распределения (params) учитывая выборочные данные (x), возвратился в виде числа.

Инверсия матрицы информации о Фишере, возвращенной как числовая матрица 2 на 2. aVar основан на наблюдаемой информации о Фишере, учитывая наблюдаемые данные (x), не ожидаемая информация.

Если значения в params MLEs параметров, aVar приближение к асимптотической ковариационной матрице отклонения (также известный как асимптотическую ковариационную матрицу). Чтобы найти MLEs, используйте mle.

Альтернативная функциональность

lognlike функционально-специализированное к логарифмически нормальному распределению. Statistics and Machine Learning Toolbox™ также предлагает родовым функциям mlecov, fitdist, negloglik, и proflik и приложение Distribution Fitter, которые поддерживают различные вероятностные распределения.

  • mlecov возвращает асимптотическую ковариационную матрицу MLEs параметров для распределения, заданного пользовательской функцией плотности вероятности. Например, mlecov(params,x,'pdf',@lognpdf) возвращает асимптотическую ковариационную матрицу MLEs для логарифмически нормального распределения.

  • Создайте LognormalDistribution объект вероятностного распределения путем подбора кривой распределению к данным с помощью fitdist функционируйте или приложение Distribution Fitter. Свойство объекта ParameterCovariance хранит ковариационную матрицу оценок параметра. Чтобы получить отрицательную логарифмическую правдоподобность оценок параметра и профиль функции правдоподобия, передайте объект negloglik и proflik, соответственно.

Ссылки

[1] Эванс, M., Н. Гастингс и Б. Пикок. Статистические Распределения. 2-й редактор Хобокен, NJ: John Wiley & Sons, Inc., 1993.

[2] Беззаконный, J. F. Статистические модели и методы для пожизненных данных. Хобокен, NJ: Wiley-межнаука, 1982.

[3] Более кроткий, W. Q. и Л. А. Эскобар. Статистические методы для данных о надежности. Хобокен, NJ: John Wiley & Sons, Inc., 1998.

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

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

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