Общий устойчивый метод

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

Поэтому автором обобщенных распределений был разработан общий устойчивый метод оценивания параметров [10-13], который по точности не уступает методу наибольшего правдоподобия, но значительно проще последнего.

Здесь так же, как и в случае универсального метода моментов, вводятся два показателя – асимметрии В и островершинности Н, которые зависят от двух параметров формы k =γ/β, u. По этим показателям устанавливается тип выравнивающей кривой распределения и находятся оценки параметров k, u. Оценки двух других параметров рассчитываются по простым формулам.

Достоинством метода является его устойчивость, т.е. он мало чувствителен к выбросам на концах статистического распределения.

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

Если обобщенное распределение задано плотностью р (x), то показатели В, Н равны

, (7.5.1)
где

. (7.5.2)

Исследования показали, что величина Н задана на интервале , а величина В – на интервале –1/4< B <1/4.

Вычислим для разных типов распределений значения показателей В, Н при различных значениях параметров k, u. Далее построим номограмму (Приложение 3). Она справедлива для трех основных систем непрерывных распределений, заданных первыми плотностями. При этом они должны быть приведены к форме плотности р (х).

На номограмме распределения II, II¢ и IV типов представлены кривыми. Типы I, I¢, III, V занимают определенные области. Симметричные распределения IIIc, Vc типов представлены отрезками на оси ОН: для IIIc типов ; для Vc типа . Распределения IVс типа представлены точкой . Распределения IIс типа также представлены точкой .

На номограмме изображены области распределений с левосторонней асимметрией, для которых . Сюда относится часть распределений III-V типов при 0< k <(1–1/ u)/2, а также распределения I, II типов. При этом распределения приведены к форме плотности р (х).

Распределения I¢, II¢ типов, а также часть распределений III-V типов при (1–1/ u)/2< k <1–1/ u имеют правостороннюю асимметрию. Для них
–1/4< B <0, причем для распределений I, II и I¢, II¢ типов справедливы равенства: .

Таким образом, показатели В, Н однозначно определяют тип распределения, приведенного к форме плотности р (х). Более того, с помощью этих показателей могут быть найдены оценки параметров u, k непосредственно из номограммы.

Для распределений III-V типов при В < 0 из номограммы вначале находятся оценки параметров k¢, u (при В > 0), затем вычисляется величина k =1–1/ u – k ¢.

Оценка параметра β для всех типов равна [11]

. (7.5.3)

Тогда γ = k β.

Оценки параметра a для распределений II, II¢ типов и произведения αu для остальных типов равны [10-12]:

(7.5.4)
где в зависимости от типа распределения величины и рассчитываются по формулам:

Типы I, I¢:

(7.5.5)

Типы II, II¢:

(7.5.6)

Типы III-V:

(7.5.7)

Величина

(7.5.8)
может быть вычислена по приближенным формулам:

- при x > 4

(7.5.9)

- при 0< x <4

, (7.5.10)
где

(7.5.11)

Для облегчения расчетов в Приложении 1 приводятся также значения функции g(x).

Для установления типа выравнивающей кривой распределения и нахождения оценок параметров по общему устойчивому методу достаточно найти значения статистических показателей и приравнять их соответствующим теоретическим. Эти показатели для каждой системы непрерывных распределений вычисляются по-своему. Но номограмма применима ко всем трем системам непрерывных распределений.

Оценки статистических показателей в случае выравнивающих распределений, заданных плотностью р (х), вычисляются по формулам:

(7.5.12)
где р_i=m_i /(Mh_i) – эмпирическая плотность распределения; m_i – наблюденная частота случайной величины Х в i -ом интервале - наблюденная частота во всех n интервалах (объем выборки); h_i – ширина i -го интервала; х_i – значение случайной величины Х в середине i -го интервала.

Формулы (7.5.12) можно выразить через абсолютные частоты m_i:

(7.5.13)

Показатель островершинности Н * при h_i = const примет вид

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

Если выравнивающее распределение задано обобщенной плотностью p(t), статистические показатели рассчитываются по формулам:

(7.5.15)

При h_i = const

. (7.5.16)

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

В заключение отметим, что общий устойчивый метод основан на взаимосвязи между законами распределения случайных величин Х и Z.

Запишем обобщенную плотность р (х)

.

Пусть для определенности параметр u > 0.

Введем случайную величину

. (7.5.17)

Тогда плотность р(z) будет равна

.

Поскольку на основании (7.5.17)

,

то

, (7.5.18)

откуда имеем замечательное равенство

β zp(z)=p(x). (7.5.19)

На его базе строится общий устойчивый метод оценивания параметров.

Поскольку плотность р(z) является функцией двух параметров формы , то последняя формула позволяет ввести критерии, зависящие от этих двух параметров.

Запишем на основании формулы (7.5.19) следующее равенство:

.

Введем обозначения

.

Тогда последнее равенство перепишется в виде

. (7.5.20)

Формула (7.5.20) позволяет найти значение параметра β; (например, при r =1), а также получить критерий островершинности, зависящий от двух параметров k, u. Для этого необходимо взять отношение либо . Последнее оказалось наиболее подходящим.

Таким путем был получен показатель островершинности Н.

Показатель асимметрии В найден из условия, чтобы для симметричных распределений он был равен нулю и в то же время использовал ранее введенные величины. Такой показатель может иметь вид

или

.

Покажем, что он зависит от двух параметров k, u.

Поскольку , то

.

По показателям В, Н строится номограмма, позволяющая устанавливать тип выравнивающей кривой распределения и находить оценки параметров k, u. Оценка параметра β вычисляется по величинам . Оценка параметра α или произведения α u вычисляется по тем же формулам, что и в случае универсального метода моментов.

Если в качестве показателей асимметрии и островершинности использовать величины

где x_с – мода, то можно построить аналогичную номограмму для установления типа выравнивающей кривой распределения и нахождения в первом приближении оценок параметров k, u по координатам одной характерной точки С и среднему значению плотности р (х).