Функция обратного кумулятивного распределения Вейбулла
X = wblinv(P,A,B)
[X,XLO,XUP] = wblinv(P,A,B,PCOV,alpha)
X = wblinv(P,A,B) возвращает обратную кумулятивную функцию распределения (cdf) для распределения Вейбулла с параметром масштаба A и параметр формы B, оценивается по значениям в P. P, A, и B могут быть векторами, матрицами или многомерными массивами одинакового размера. Скалярный вход расширяется до постоянного массива того же размера, что и другие входы. Значения по умолчанию для A и B являются оба 1.
[X,XLO,XUP] = wblinv(P,A,B,PCOV,alpha) возвращает доверительные границы для X когда входные параметры A и B являются оценками. PCOV - матрица 2 на 2, содержащая ковариационную матрицу оцененных параметров. alpha имеет значение по умолчанию 0,05 и указывает 100 ( 1 -alpha)% доверительных границ. XLO и XUP массивы того же размера, что и X содержащий нижнюю и верхнюю доверительные границы.
Функция wblinv вычисляет доверительные границы для X используя нормальное приближение к распределению оценки
logqb ^
где q - значение P-й квантиль из распределения Вейбулла с параметрами масштаба и формы, равными 1. Вычисленные границы дают приблизительно требуемый уровень достоверности при оценке mu, sigma, и PCOV из больших выборок, но в меньших выборках другие способы вычисления доверительных границ могут быть более точными.
Обратная величина cdf Вейбулла равна
(1 − p)] 1/b.
Время жизни (в часах) партии лампочек имеет распределение Вейбулла с параметрами a = 200 и b = 6.
Найдите среднее время жизни луковиц:
life = wblinv(0.5, 200, 6) life = 188.1486
Создать 100 случайных значений из этого распределения и оценить 90-й процентиль (с доверительными границами) из случайной выборки
x = wblrnd(200,6,100,1);
p = wblfit(x)
[nlogl,pcov] = wbllike(p,x)
[q90,q90lo,q90up] = wblinv(0.9,p(1),p(2),pcov)
p =
204.8918 6.3920
nlogl =
496.8915
pcov =
11.3392 0.5233
0.5233 0.2573
q90 =
233.4489
q90lo =
226.0092
q90up =
241.1335