Презентация на тему Оценка параметров распределения по эмпирическим данным (лекция 4)

СВ

Выборка – это конечный набор значений СВ, полученный в

результате наблюдений

Репрезентативная выборка – это выборка, которая достаточно полно характеризует генеральную совокупность

Задача статистических методов – определить свойства СВ в целом на основании анализа выборки

Статистические оценки (mx*, σx*, Dx* и т.д.) – это числовые характеристики СВ, полученные по эмпирическим данным.

f(x1, x2, x3,….,x) – неизвестного параметра G называется состоятельной,

если по мере роста числа наблюдений n она стремиться к оцениваемому значению G, т.е.

ε – сколь угодно малое число
2. Несмещенность. Оценка G* = f(x1, x2, x3,….,x) – неизвестного параметра G называется несмещенной, если при любом объеме выборки n результат ее осреднения по всем возможным выборкам данного объема приводит к точному (истинному) значению оцениваемого параметра, т.е., т.е. M[G*] = G
Несмещенность означает отсутствие систематической погрешности при оценивании параметра
3. Эффективность. Оценка G* = f(x1, x2, x3,….,x) – называется эффективной, если среди всех оценок параметра G она обладает наименьшей мерой случайного разброса относительно истинного значения оцениваемого параметра, т.е. D[G*] = Dmin
Эффективная оценка имеет минимальную случайную погрешность.

обеспеченности, построенная по эмпирическим данным.

Возможны два способа построения

эмпирической кривой обеспеченности.
 
Первый способ – при наличии большого числа наблюдений

Второй способ – при наличии небольшого числа наблюдений

наблюдений
Определяется амплитуда (размах R) колебаний СВ

2. Разбивается амплитуда

колебаний на k равных интервалов. Величина k примерно рассчитывается по формуле k ≈ 5ln(n)

3. Определяется длина расчетного интервала по формуле l = R/k
В левой границы первого интервала принимается значение большее или равное максимальному значению СВ. Тогда значение правой границы первого интервала будет равно разнице между левой границей и длиной расчетного интервала. Левая граница каждого последующего интервала должна быть меньше на 1 соответствующей правой границы интервала

наблюдений (продолжение)
4. Подсчитывается число случаев попадания СВ в каждый

интервал (по этим данным можно построить график эмпирических частот)

5. Рассчитывается относительная частота попадания СВ в каждый интервал по формуле
p* = (mi/n) •100%
 
где mi – число попаданий СВ в i – й интервал
 
Значения p* последовательно суммируются и умножаются на 100% (по этим данным можно построить график эмпирической кривой обеспеченности)
В гидрологии значения СВ принято откладывать по вертикали, а значение вероятности р – по горизонтали.

наблюдений
1. Ряд наблюдений располагается в убывающем или возрастающем порядке.

В гидрологии – в убывающем порядке.
2. Приближенно вычисляется обеспеченность по формуле:

где m - порядковый номер xm в ранжированном ряду;
pm – обеспеченность (в%) m – ного члена ранжированного ряда

При расчете обеспеченности последнего члена ряда по этой формуле получится, что pm = (n/n)100 = 100%
Чтобы этого избежать используют другие формулы
Формула Хансена pm = ((m-0.5)100)/n; Формула Крицкого – Менкеля (Вейбула) pm = 100m/(n+1); Формула Чегодаева pm = (100(m-0.3))/(n+0.4),
Универсальная формула Грингортена
pm = (100(m-a))/(n+1-2a)
при определенных значениях а – получаются все перечисленные формулы. Сам Грингортен предложил определять а по длине ряда (по таблицам).

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

это МО, СКО, Cv, Cs или Cv/Cs

 
Эти методы расчета условно делят на аналитические, графоаналитические и графические.


К числу аналитических методов относятся метод моментов и метод наибольшего правдоподобия.

расчет эмпирических моментов. При этом вместо N → ∞

берется конечное число значений СВ - n, а теоретическая вероятность pi заменяется на расчетную pi = 1/n.
Эмпирический начальный αs* и центральный μs* моменты S го порядка определяются по формулам

Тогда математическое ожидание МО
можно вычислить по формуле
Эта оценка состоятельная и несмещенная.

Дисперсия – это второй центральный момент, поэтому ее можно вычислить так
Эта оценка состоятельная, но смещенная. Поэтому для расчетов используется формула
где SH2 и SC2 – соответственно несмещенная и смещенная оценки дисперсии;
n/(n-1) – поправка на смещенность

SC2 – соответственно несмещенная и смещенная оценки дисперсии;
n/(n-1) –

поправка на смещенность

Тогда σ* и Cv* определяются по формулам

где ki = xi/xср.. – модульный коэффициент

Несмещенная оценка Cs* определяется как

Преимущество: метод не зависит от закона распределения СВ
Недостаток: при больших значениях Cv* (больше 0.5), достоверность оценок ощутимо снижается.

правдоподобия нужно, прежде всего, построить функцию правдоподобия

Для этого делаются

следующие последовательные шаги:

1. Заданная аналитическая функция распределения вероятности логарифмируется. Берется натуральный логарифм.

2. Создается функция правдоподобия путем интегрирования прологарифмированной функции плотности вероятности

3. Затем для каждого параметра распределения создается своя функция правдоподобия путем дифференцирования полученной функции распределения по требуемому параметру и приравнивается к нулю, чтобы найти ее максимум.

4. Из полученных уравнений находятся оценки, например МО и СКО.

нормального распределения

здесь mx и σx - искомые параметры

распределения.
1. Найдем логарифм этой функции

2 Проинтегрируем ее

3. Далее находим уравнение правдоподобия для оценки параметра mx

получаем, что

правдоподобия для оценки σx
Отсюда получим, что
То есть

для нормального распределения оценки параметров, полученные МНП, совпадают с моментными оценками.

Для других функций распределения плотности вероятности система уравнений правдоподобия получается сложной. Поэтому используются численные методы решений, а на их основе строятся номограммы для практического применения МНП.

и Менкеля параметры Cs и Cv определяются с помощью

специально разработанных номограмм как функций вспомогательных статистик λ2 и λ3.

Недостаток: нужно знать аналитическое выражение функции распределения заранее.

Считается, что оценки трех - параметрического гамма – распределения, полученные таким путем, являются состоятельными, эффективными и несмещенными.
При использовании МНП нужно помнить, что наибольший вес придается средним членам выборки, в отличие от метода моментов, где наибольший вклад вносят крайние члены выборки. Однако, это свойство МНП проявляется в случае выборок с большим рядом.