Список функций Statistics Toolbox

В оглавление \ К следующему разделу \ К предыдущему разделу

Синтаксис

phat = gamfit(x)
[phat,pci] = gamfit(x)
[phat,pci] = gamfit(x,alpha)

Описание

phat = gamfit(x) функция служит для расчета вектора точечных оценок phat параметров Гамма распределения a и b по исходной выборке значений x. Значения оценок параметров определяются методом максимального правдоподобия. Выборка x задается как вектор.

[phat,pci] = gamfit(x) функция служит для расчета точечных phat и интервальных pci оценок параметров Гамма распределения a и b по исходной выборке значений x. Интервальные оценки pci параметров a и b задаются в виде матрицы размерностью 2-2. Первый столбец матрицы pci содержит нижнюю и верхнюю границы доверительного интервала параметра a, второй столбец - нижнюю и верхнюю границы доверительного интервала b. Интервальная оценка параметров Гамма распределения соответствует 95% доверительной вероятности.

[phat,pci] = gamfit(x,alpha) в отличие от второго варианта синтаксиса в качестве второго входного параметра задается уровень значимости alpha. Доверительная вероятность для интервальной оценки параметров a, b определяется как (1- alpha).

Примеры использования функции расчета точечных и интервальных оценок параметров Гамма распределения

Расчет точечных оценок параметров Гамма распределения

>> A=4;
>> B=3;
>> x = gamrnd(A,B,100,1);
>> phat = gamfit(x)
phat =
    3.4797    3.3897
>>a= phat(1)
a =
    3.4797
>>b= phat(2)
b =
    3.3897
        
Расчет точечных и интервальных оценок параметров Гамма распределения для уровня значимости 0,05.

>> A=5;
>> B=2;
>> x = gamrnd(A,B,100,1);
>> [phat,pci] = gamfit(x)
phat =
    5.0306    1.9073
pci =
    3.4819    1.3436
    6.5792    2.4711
>>a= phat (1)
a =
    5.0306
>> b= phat (2)
b =
    1.9073
>> a_low= pci (1,1)
a_low =
    3.4819
>> a_high= pci (2,1)
a_high =
    6.5792
>> b_low= pci (1,2)
b_low =
    1.3436
>> b_high= pci (2,2)
b_high =
    2.4711
        
Расчет точечных и интервальных оценок параметров Гамма распределения для уровня значимости 0,01.

>> A=5;
>> B=2;
>> alfa=0.01;
>> x = gamrnd(A,B,100,1);
>> [phat,pci] = gamfit(x, alfa)
phat =
    5.9167    1.8883
pci =
    3.6012    1.0733
    8.2322    2.7034

В оглавление \ К следующему разделу \ К предыдущему разделу