Загрузите образцы данных. Создайте объект распределения вероятности, установив распределение Вейбулла в милях на галлон (MPG) данные.

load carsmall
pd = fitdist(MPG,'Weibull')

pd = 
  WeibullDistribution

  Weibull distribution
    A = 26.5079   [24.8333, 28.2954]
    B = 3.27193   [2.79441, 3.83104]

Просмотрите имена параметров для распределения.

pd.ParameterNames

ans = 1x2 cell
    {'A'}    {'B'}

Для распределения Вейбулла, A находится в положении 1, и B находится в положении 2.

Вычислить вероятность профиля для B, который находится в положении pnum = 2.

[ll,param] = proflik(pd,2);

Просмотр значений средств к существованию для оценочных значений B.

[ll',param']

ans = 21×2

 -329.9688    2.7132
 -329.4312    2.7748
 -328.9645    2.8365
 -328.5661    2.8981
 -328.2340    2.9597
 -327.9658    3.0213
 -327.7596    3.0830
 -327.6135    3.1446
 -327.5256    3.2062
 -327.4943    3.2678
      ⋮

Эти результаты показывают, что вероятность регистрации профиля максимизируется между оцененными B значения 3.2678 и 3.3295, которые соответствуют значениям логарифма -327.4943 и -327.5178. Из более ранней посадки, MLE B равно 3,27193, что соответствует ожидаемому интервалу.

Загрузите образцы данных. Создайте объект распределения вероятности, установив обобщенное распределение экстремальных значений в милях на галлон (MPG) данные.

load carsmall
pd = fitdist(MPG,'GeneralizedExtremeValue')

pd = 
  GeneralizedExtremeValueDistribution

  Generalized Extreme Value distribution
        k = -0.207765   [-0.381674, -0.0338563]
    sigma =   7.49674   [6.31755, 8.89604]
       mu =   20.6233   [18.8859, 22.3606]

Просмотрите имена параметров для распределения.

pd.ParameterNames

ans = 1x3 cell
    {'k'}    {'sigma'}    {'mu'}

Для обобщенного распределения экстремальных значений k находится в положении 1, sigma находится в положении 2, и mu находится в положении 3.

Вычислить вероятность профиля для mu, который находится в положении pnum = 3. Ограничьте вычисление значениями параметров от 20 до 22 и просмотрите график.

[ll,param,other] = proflik(pd,3,20:.1:22,'display','on');

На графике показано расчетное значение параметра mu это максимизирует средства к существованию.

Просмотр значений средств к существованию для оценочных значений muи значения других параметров распределения, которые максимизируют соответствующие средства к существованию.

[ll',param',other]

ans = 21×4

 -327.5706   20.0000   -0.1803    7.4087
 -327.4971   20.1000   -0.1846    7.4218
 -327.4364   20.2000   -0.1890    7.4354
 -327.3887   20.3000   -0.1934    7.4493
 -327.3538   20.4000   -0.1978    7.4636
 -327.3317   20.5000   -0.2023    7.4783
 -327.3223   20.6000   -0.2067    7.4932
 -327.3257   20.7000   -0.2112    7.5084
 -327.3418   20.8000   -0.2156    7.5240
 -327.3706   20.9000   -0.2201    7.5399
      ⋮

Первый столбец содержит логарифмическое значение правдоподобия, которое соответствует оценке mu во втором столбце. Логарифмическое правдоподобие максимизируется между значениями параметров 20.6000 и 20.7000, соответствующими значениям логарифмического правдоподобия -327.3223 и -327.3257. Третий столбец содержит значение k что максимизирует соответствующее логарифмическое правдоподобие для mu. Четвертый столбец содержит значение sigma что максимизирует соответствующее логарифмическое правдоподобие для mu.