ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 08.02.2025

Просмотров: 1280

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
ni = 1: Площадь
Рис. 14.2.

^

Если все выборочные значения x1;:::;xn различны, то функцию F (x;~) можно записать в виде

0;

x x(1);

¡

> n

Fn(x) =

8

i

;

x(6i) < x 6 x(i+1); i =

1;n

1;

>

:

1;

x > x(n);

<

т.е. в каждой точке x(i) функция Fn(x) имеет скачок величиной 1=n. График функции Fn(x) изображен на рис. 14.1.

При больших объемах выборки n (n > 50)обычно

производят группирование исходных данных следу-

ющим образом. Промежуток J = [x(1);x(n)], содержа-

щий все выборочные значения, разбивают на m полу-

интервалов J1;:::;Jm, как правило, одинаковой дли-

ны ¢ и таких, что каждый из них, кроме послед-

него, содержит левую границу, а последний содер-

жит обе границы, и подсчитывают число ni элемен-

тов выборки, попавших в i-ый промежуток Ji, i = 1;m,

n = n1 + ¢¢¢ + nm, а результаты представляют в виде

Рис. 14.1.

следующей таблицы 14.2, которую называют интер-

вальным статистическим рядом.

Иногда в верхней строке таблицы 14.3 указывают не интервал, а

J1

J2

. . .

Jm

его середину xi, а в нижней строке вместо частоты ni записывают от-

. . .

m

носительную частоту ni=n.

n1

n2

nm

ni

i=1

промежутков m, на которые разбивают промежуток J =

Число

e

Таблица 14.3.

P

[x(1);x(n)], содержащий все выборочные значения, выбирают в зави-

симости от объема выборки n. Для ориентировочной оценки величины m можно пользоваться следующей приближенной формулой.

m ¼ log2 n + 1;

Определение 14.8 График функции

8 n¢i

x 2 Ji;

(14.3)

pn(x) =

;

n

:

0;

2

<

x = J;

представляющий собой кусочно постоянную функцию называют гистограммой (см. рис. 14.2).

Часто гистограммой называют диаграмму, составленную из прямоугольников с основанием ¢ и высотами ni=(n¢), i = 1;m. Нетрудно увидеть, что суммарная площадь всех прямоугольников, образующих такую диа-

грамму, равна 1, так как Pm nn¢i ¢ = n1 Pm

i=1 i=1

каждого прямоугольника ni=n есть частота попадания элементов выборки в соответствующий интервал Ji статистического ряда.

~

Рассмотрим случайную величину ni(Xn)=n, которая для

~

каждой реализации ~xn случайной выборки Xn ~ равна частоте ni=n. В соответствии с законом больших чисел в форме Бернулли ni(Xn)=n при n ! 1 будет сходиться по вероятности к вероятности попадания случайной величины X в промежуток Ji, i = 1;m, т.е.

~

ni(Xn)

n

©

ª

Ji

P

= Z p(x)dx;

n¡!!1 P X 2 Ji

где p(x) плотность распределения генеральной совокупности X. Если длина ¢ промежутков достаточно мала и объем выборки n велик, то с вероятностью, близкой к 1, можно утверждать, что

nni ¼ p(xei)¢;

или

ni

¼ p(xei);

67


где xei середина промежутка Ji, i = 1;m. Таким образом, при большом объеме выборки n и достаточно малом ¢ с вероятностью, близкой к 1, можно считать, что pn(x) ¼ p(x). Иными словами, функция pn(x) является статистическим аналогом плотности распределения p(x), наблюдаемой в эксперименте случайной величины X, а гистограмма выглядит приблизительно как график плотности случайной величины X.

~

Пусть Xn случайная выборка из генеральной совокупности X с функцией распределения F (x) (и плотностью распределения p(x) в случае непрерывной статистической модели).

Случайную величину

1

n

~

Xi

k

¹k(Xn) =

n

Xi

(14.4)

=1

называют выборочным начальным моментомb

k-го порядка. В частности, выборочный на-

~

чальный момент первого порядка X = ¹ (X ) называют выборочным средним.

Случайную величину

b1

n

1 n

¡

¢

k

~

Xi

ºk(Xn) =

n

=1

Xi ¡ X

(14.5)

называют выборочным центральнымb

моментом k-го порядка. В частности, выборочный

центральный момент 2-го порядка ¾

2

~

~

(Xn) = º2(Xn) называют выборочной дисперсией.

~

) =

2

~

X

¾b(Xn) называют выборочным средним квадратич-

Выборочную характеристику ¾(b n

ным отклонением.

b

q

Случайную величину

b

1

n

~

~

b

Xi

¡

¢¡

¢

(14.6)

K(Xn

;Yn) = n =1

Xi ¡ X

Yi ¡ Y

называют выборочным корреляционным моментом.

Выборочную характеристику

~

~

K(X ;Y )

½(X~n;Y~n) =

¾x(bX~n)n¾y(nY~n)

;

(14.7)

где

nb

¢

b

b

n

¡

¢

1

1

Xi

2

~

2

2

~

2

¾x

(Xn) =

n

i=1

Xi ¡ X

;

¾y

(Yn) =

n

=1

Yi ¡ Y

;

называют выборочнымb

коэффициентом корреляцииb .

Значения выборочных моментов (неслучайные числа) будем для краткости называть теми же терминами, что и сами моменты (случайные величины).

Основное свойство выборочных моментов состоит в том, что при увеличении объема выборки n

они

сходятся по вероятности к соответствующим теоретическим моментам. В частности, при n

! 1

P

P

2

~

имеем X n¡!!1 MX, а ¾b

(Xn) n¡!!1 DX.

68


Лекция 15

Точечные оценки

Одной из задач математической статистики является оценка неизвестных параметров выбранной

параметрической модели.

Предположим, что закон распределения генеральной совокупности принадлежит множеству

~ ~ ~

fF (x;µ) : µ 2 £g, где вид функции распределения задан, а вектор параметров µ = (µ1;:::;µr) неизве-

~

стен. Требуется найти оценку для µ или некоторой функции от него (например, математического ожидания, дисперсии) по случайной выборке (X1;:::;Xn) из генеральной совокупности X.

Например, предположим, что масса X детали имеет нормальный закон распределения, но его параметры µ1 = MX и µ2 = DX неизвестны. Нужно найти приближенное значение параметров по результатам наблюдений x1, . . . , xn, полученным в эксперименте (по реализации случайной выборки).

~

Пусть Xn = (X1;:::;Xn) случайная выборка из генеральной совокупности X, функция распре-

деления F

x;µ которой известна, а µ неизвестный параметр, т.е. рассматривается параметриче-

~ ~

£

(для простоты изложения будем считать пока, что µ скаляр).

ская

модель F (x;µ) : µ

¡ f ¢

2 g

~

Требуется построить статистику µ(X ), которую можно было бы принять в качестве точечной

оценки параметра µ.

b

n

Определение 15.1 Точечной оценкой параметра µ 2 £ назовем любую функцию от наблюдений

^ ~

(т.е. любую статистику) µ(Xn).

b ~

Определение 15.2 Статистику µ(Xn) называют состоятельной оценкой параметра µ 2 £, если с ростом объема выборки n она сходится по вероятности к оцениваемому параметру µ, т.е.

~

)

P

µ:

µ(X

n¡!!1

b

n

b ~

Определение 15.3 Статистику µ(Xn) называют несмещенной оценкой параметра µ, если ее

b ~

математическое ожидание совпадает с µ, т.е. Mµ(Xn) = µ для любого фиксированного n.

Замечание 15.1 Можно показать, что статистика

n

¡

Xi

¡

¢

2

S2

=

1

Xi ¡ X

n

1

=1

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

Замечание 15.2 Можно доказать, что выборочные начальные и центральные моменты являются состоятельными оценками соответствующих теоретических моментов, если только они существуют. Однако эти оценки, кроме X, являются смещенными.

Пример 15.1 Пусть X1, . . . , Xn случайная выборка из генеральной совокупности X, имеющей нормальное распределение с неизвестным средним значением µ и известной дисперсией ¾2.

b b

Оценка µ = µ(X1;:::;Xn) = X1 является несмещенной для µ, ибо MX1 = MX = µ, но не является состоятельной, так как, во-первых, X1 не зависит от объема выборки и, следовательно, ее распределение не меняется с ростом n, а во-вторых,

P jX1

¡ µj < " = ¾p2¼

Z

e¡2¾2

dt 6= 1:

©

ª

2

"

t2

0

69


Метод моментов

Метод моментов был предложен английским статистиком К. Пирсоном и является одним из

первых общих методов оценивания. Он состоит в следующем.

~

;:::;Xn) из генеральной совокупности X, распреде-

Пусть имеется случайная выборка Xn = (X1

~

~

;:::;µr). Требуется найти

ление которой p x

известно с точностью до вектора параметров µ = (µ1

~

~

.

оценку

параметра µ по случайной выборке X

n

¡ ¢

Обозначим

~

k

~

k

¹k(µ) = M(X) ; ºk(µ) = M(X ¡ MX)

начальный и центральный моменты порядка k, k = 1;2;:::. В методе моментов в качестве точечной

~

~

~

~

~

r

µ(X ) = (µ

(X

n

);:::;µ (X ))

вектора~параметров

µ

берут решение системы

уравнений

оценки b

n

1

r

n ~

b

b

¹i® (Xn) =¹i® (µ); ® =

1;k;

(

º (X~

) =º

(µ~); ¯ =

1;l;

;

k + l = r;

(15.1)

b

n

bj¯

относительно неизвестных µ1;:::;µr. Индексы i® и j¯ выбирают так,чтобы система уравнений решалась как можно проще. Можно показать,что при условии непрерывной зависимости решения этой системы от ¹bi® и ºbj¯ , оценка, полученная методом моментов, является состоятельной и имеет асимптотически нормальное распределение, т.е. ее распределение при n ! 1 стремится к нормальному. При этом уравнения (15.1) во многих случаях просты и их решение не вызывает больших вычислительных сложностей.

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

Пример 15.2 Пусть дана случайная выборка (X1;:::;Xn) объема n из генеральной совокупности X, имеющей равномерный закон распределения

p(x) =

1

; x 2 (a; b);

b

a

(

¡0; x = (a; b);

2

с неизвестными параметрами a и b. Найдем методом моментов точечные оценки этих параметров. Известно, что для равномерно распределенной случайной величины X

MX =

a + b

;

DX =

(b ¡ a)2

:

2

12

выборочная дисперсия ¾

2

~

X

Выборочное среднее

и

1 n

2

b~

n

1

n

2

X

Xi

¡ X) :

X = n i=1 Xi;

¾ (Xn) = n =1(Xi

b

Составляем систему двух уравнений

8

2

2

>

a + b

= x;

(b ¡ a)

= ¾2(~x ):

<

12

n

>

b

:

Решая систему, получаем

p

Окончательно имеем

b= x +

3¾b(~xn);

~

p

~

b(Xn) = X + 3¾b(Xn);

p

ba = x ¡ 3¾b(~xn):

p

~ ~

ba(Xn) = X ¡ 3¾b(Xn):

Метод максимального правдоподобия

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

~

из генеральной совокупности X,

Рассмотрим функцию правдоподобия случайной выборки Xn

~

~

£:

p x;µ

точностью до параметра µ

2

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

¡ ¢ которой известна с

n

¢

~

iY ¡

L(X1;:::;Xn;µ) =

p Xi;µ :

=1

70