Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2003_-_Gmurman__TV_i_MS

.pdf
Скачиваний:
19
Добавлен:
27.03.2015
Размер:
16.8 Mб
Скачать

Ре w е н и е.

Составим функцию

правдоподсбия, учитывая, что

е =А.:

 

 

 

 

 

 

L=p 1; Лt Р (Х2;

А) '" р (Хn;

А.) =

АХ' е-1.. АХ. е-1..

A.Xn .e-1.. _

A~XI.

x11

xll

 

хnl

 

x1lxl ' .

Найдем логарифмическую функцию правдоподобня:

lп L = (~X,) 10 А-nA-10 (Хl'XI!

••• Хn!).

Найдем первую производную по А.:

 

 

 

d 10 L

~ Х/

n.

 

 

dX -

А.

 

 

Напншем уравненне правдоподобия, для чего приравияем пер­

вую производную иулю:

(~ХIIА)-n=О.

Найдем критнческую точку, для чего решим полученное уравне­

ние относнтельно А:

A=~Xifn=XB'

Найдем вторую производную по А:

da lnL

= -

 

dAa

Легко видеть, что при А = хв вторая производная отрицательна;

следовательно, А=хв-точка макснмума и, значит, в качестве оценка наибольwого праВДОП04обия параметра А распределения Пуассоиа

надо принять выборочкую среднюю л*=хв.

Пример 2. Найти методом наибольшего правдоподоБИЯ оценку

параметра р биномиального

распределения

 

Рn (k) =c~· (l_p)n-k,

 

если D nl независимых испытаниях событие А появилось

Xl = тl раз

и в Па независимых испытаниях событие

А появилось XI = та раз.

Р е w е н и е. Составим

функцию

правдоподобия,

учитывая,

что 6=р:

L=P

n.

(m)Р

n.

)=CmaCmlpml+ml(l_p)[{n.+nll-(m.+m.lJ.

 

1

I

n. n.

Найдем логарифмическую функцию правдоподобия:

10 L= lп (С':.'С':.") +(ml +т.) 10 р + [(nl +nl)- (тl +mJJ lп (I-p).

НаЙдем первую ПРОНЗВОДIJУЮ по р:

d lп L

!'Il +ml

(nl +nl )-(ml +та)

dp

р

I-р'

Напишем уравнение правдопод.обия, д,IIя чего приравняем первую

пронзводную нулю'

(nl +nl) - (тl +та) -О.

I-р

231

Найдем крнтнческую точку, д,IIЯ чеro решим пo.nучеиное урав­

ненне отиосительно р:

р = (ml +m2)/(nl +nl)'

Найдем вторую производную по р:

d? 10 L

ml +m l + (n1+n2)-(ml +ml)

dpl

рl

(l-р)1

Легко убедиться,

что при

р = (ml + m l )/(nl + nl) вторая произ­

водная отрицательна;

следовательно, р = (ml +m2)/(nl +n.) -точка

максимума и, значит, ее надо

принять в качестве оценки наибo.nь­

шего правдоподобия

нензвестной вероятности р биномиальиого рас­

пределения:

р*= (ml +m2)/(nl +nl)'

Б. Непрерывные случайные веJlИЧИНЫ. Пусть Х-не­

п ре рыв н а я с л у чай н а я в е л и ч и н а, которая в ре-

u

зультате n испытании приняла значения Х1X1 , ••• , Хn'

Допустим, что вид плотности распределения f (Х) задан, но не известен параметр В, которым определяется эта

функция.

4>уnкцией nравдоnодо6ия непрерывной случайной вели­ чины Х называют функцию аргумента В:

L 1X1 , ••• , Хn; В) = f 1; В) f (x1 ; В) ..• f n; В),

где Х1' X1 , •••• хn-фиксированные числа.

Оценку наибольшего правдоподобия неизвестного па­

раметра распределения непрерывной случайной величины ищут так же, как в случае дискретной величины.

Пример 3. Найти методом наибо.пьшеro правдоподобия оценку параметра л показательноro распределения

f(х)=Ле-

М

< Х < CIO),

 

 

если в результате n испытаний случайная величина Х, распределен­

ная по показательному закону. приняла значения Xl,

X1 , ••• , Хn

Реш е н и е. Составим функцию правдоподобия,

учитывая,

что 8=Л:

 

L=f(Xl; Л)f(х.; л) ... {(Хn; Л)=(Ле-Ма)(Ле-М1)

(А.е-Ма).

Orсюда

 

Найдем логарифмическую функцию правдоподобия:

 

lп L=n10 л-л ~ Х/.

 

Найдем первую производную по Л:

d 10L

n ~

-1: - ~ xI·

Напишем уравнение правдопод()бия, д,IIя чего прнравняем первую

производную нулю:

(nр..)-}:; Xj=O.

Найдем критическую точку, для чего решнм полученное уравне­ ние относительно А:

A=n!~ Xj= I/(~ Xi/n) = l!xB •

Найдем вторую пронзводную по А:

tJ2 lп L

n

d')..2

А2

Легко вндеть, что при А = '/Хв вторая производная отрицательна;

следовательно, л= l/хв-точка максимума и, значнт, в качестве

оценки нанбольшего правдоподобия параметра л показательного рас­

пределения надо принять величину, обратную выборочной средней:

А*= l/xB

3 а м е ч а н и е. Если плотность распределения f (Х) непрерывной

случайной величины Х опредf'ЛЯется двумя нензвестнымн парамет­

рамн 61 н 61, то

функция правдоподобия явля~ся функцией двух

независимых аргументов 61 и 62:

 

 

 

 

 

 

 

 

L = f (Хl;

61,

(2) f (x1; 61,

 

(2) .•• f n; 61' (2)'

где Хl,

Х2, ••.• Хn-

наблюдавшиеся

 

значения Х.

Далее находят ло­

гарифмическую функцию

правдоподобия и д,IIя отыскаиия ее макси­

мума

составляют н решают систему

 

 

 

 

 

 

 

 

 

 

 

{

д 'П L

 

 

О

 

 

 

 

 

 

 

 

де1

 

 

'

 

 

 

 

 

 

 

 

д lп L_

o

 

 

 

 

 

 

 

 

де

 

 

 

 

 

 

 

 

 

 

 

2

- .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Найти методом наибольшего правдоподобия оценки

пара метров а н о нормального

распределения

 

 

 

 

 

 

()

__

1-==-

 

-(x-a)"/Ial

 

 

 

 

f Х

=

..r - е

 

 

 

 

 

 

 

 

 

О ,. 2п

 

 

 

 

 

 

если

в

результаге

n

испытаний

величина

Х

приняла значения

Хl, Х.,

••• , хn

 

 

 

 

 

 

 

 

 

 

Реш е н н е. Составим функцию правдоподобия, учитывая, что

61 =а и 62=0:

L= 1

о У2п

Orсюда

2ЗЗ

Найдем nогарифмическую функцию праВДОПОДОбия:

I

~ (х;-а)"

IпL=-nlоо+lо (У2п)n

202'

НаАдем частные производные по а и по о:

д 10 L

~ х{-nа

д 'П L

n

~ (x{-а)1

да

(Ji'

до

а+

03 •

Приравняв частные производные нулю и решив полученную си­

стему двух уравненнй относительно а и 02, получим:

а=~ xi!n=Xa; 01 =(~ (Xi-Xa)2)/n =Da .

Итак, искомые оценкн наибольшего правдооодобия: а8 Ха;

08= YDa Заметнм, что первая оценка несмещениая, а вторая сме­

щенная.

§ 23. Другие характеристики вариационного ряда

Кроме выборочной средней и выборочной диспер­

сии применяются и другие характеристики вариационного

ряда. Укажем главные из них.

Модой М. называют варианту, которая имеет наиболь­ шую частоту. Например, для ряда

варианта

... J

4

7

9

частота

.... 5

1

20

6

мода равна 7.

Медианой те называют варианту, которая делит ва­

риационный ряд на две части, равные по числу вариант.

Если число вариант нечетно, т. е. n=2k+ 1, тоте =Х"+l;

при четном n = 2k медиана

те = (х" +х"+ 1)/2.

Например, для ряда 2 3

5 6

7 медиана

равна 5; для

ряда 2 3 5 6 7 9 медиана

равна

(5 +6)/2 =

5,5.

Раэ.4faXо.м варьирования R называют разность между

наибольшей и наименьшей

вариантами:

 

R = xmax-Xmln'

Например, ДЛЯ ряда 1 3 4 5 6 10 размах равен 10-1 = 9.

Размах является простейшей характеристикой рассея­

ния вариационного ряда.

Средни.м абсолютны.м отклонение.м е называют среднее

арифметнческое абсолютных отклонений:

 

=

"

Xj - X-

,...

е

 

<.~ n; I

 

a I)/~ n,.

Например, для ряда

 

 

Х,

1

3

6

16

 

 

nl 4

10

5

1

имеем:

 

 

 

 

 

 

 

Ха

= 4·1 + 10·3+5·6+ 1·16 = 80_ 4.

 

 

4+10+5+1

 

20- ,

-

4· \1-41+10·13-41+5·16-41+1·116-41_22

6 -

 

 

 

20

 

- , .

Среднее

абсолютное

отклонение служит для характерис­

тики рассеяния вариационного ряда.

Коэффициенто.м вариации V наз~вают выраженное в

процентах отношение выборочного среднего квадратичес­

кого отклонения к выборочной средней:

V = ав/х•. 1ОО%.

Коэффициент вариации служит для сравнения величин

рассеяния по отношению к выборочной средней двух

вариационных рядов: тот из рядов имеет большее рас­

сеяние по отношению к выборочной средней, у которого

коэффициент вариации больше. Коэффициент вариации­

безразмерная величина, поэтому он пригоден для срав­

нения рассеяний вариационных рядов, варианты которых

имеют различную размерность, например если варианты

одного ряда выражены в сантиметрах, а другого-в грам­

мах.

3 а м е ч а н и е. Выше предполагалось, что вариационный ряд

составлен по данным выборки, ПО9Тому все описанные характерис·

тики называют 8Ь1борочftbUlU; если вариационный ряд составлен по

данным генеральной совокупности, то характернстикн называют ине·

PaAbftblMU.

Задачи

1. Найти

групповые средние

совокупности, состоящей из

двух групп:

 

 

 

 

 

первая группа • • •

xl

0,1

0,4

0,6

 

 

nj

3

2

5

вторая группа

••

Xj

0,1

0,3

0,4

 

 

n'

10

4

6

Omв. Хl =0,41; xl=O,23.

2. Найти общуlO среднюю по данным задачн 1 двумя способами:

8) объедннить обе группы в одну совокупность; б) нспользовать наЙ· денные в задаче 1 групповые средние.

Omв. х= 0,29.

235

3. Дано распределение статистической совокупности:

Xj

n{

1

4

5

6

11

3

Убедиться, что сумма произведеннй отклонений на соответствуюшне

частоты равна нулю.

4. Дано распределение статистической совокупности:

7

10

15

15

20

5

Найти дисперсию совокупности: а) исходя из определения диспер­

сии; б) пользуясь формулой D=x2-[X]2.

Отв. D =9,84.

5. Найти внутригрупповую, межгрупповую и общую дисперсии

совокупности, состоящеQ из трех групп:

первая группа

.

Xj

1

2

8

 

 

 

 

ni

30

15

5

вторая группа

Xj

1

6

 

 

 

 

 

nj

10

15

 

третья группа

 

••

Х;

3

8

 

 

 

 

 

n,

20

5

 

Omв. Dвигр = 4.6;

Dкежгр = 1;

Dобщ = 5.6.

8. Найти внутригрупповую, межгрупповую и общую дисперсии

совокупности, состоящей из двух групп:

первая группа

••

Х

2

7

 

 

n,

6

4

вторая группа

• • •

xi

2

7

ni 2 8

Отв. Dвигр=5;~Dкежгр=l; Dобщ=6.

7. Найти выборочную и исправленную дисперсни варнационного

ряда, составленного по даиным выборкам:

варианта • • . I

2

5

8

9

частота

. •• 3

4

6

4

3

Отв: о: = 8,4;

S2 = 8,84.

 

 

 

 

В задачах 8-9 даны

среднее

квадратическое отклонение, выбо­

рочная средняя и

объем

выборкн

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

нака. Найти доверительные интервалы для оценки неизвестного мате·

матического ожидания с заданной надежностью.

8. 0=2, хв=5,40, n=10, 1'=0,95.

Оmв. 4,16 < а < 6,64.

9. 0=3, xB=20,1~. n=25, 1'=0,99.

Omв. 18,57 < а < 21,67.

10. Найти миннмальный объем выборки, при котором с надеж­

ностью 0,95 точность оценки математического ожидания нормально

распределенного признака по выборочной средней будет равна 0,2, если среднее квадратическое отклоненне равно 2.

у к а з а н и е. СМ. замечание 2, § 15.

Отв. n=385.

Взадачах 11-12 даны «исправленное» среАнее квадратическое

откnонение, выборочная средняя и объем малой выборки нормально распределенного признака. Найти, пользуясь распределением Стью-

236

дента, доверительные ннтервалы д,IIЯ оценки неизвестного математи­

ческого ожидания с заданной надежностью.

11. s= 1,5, хв = 16,8, n= 12,

1'=0,95.

 

Оmв. 15,85 < а < 17,75.

 

 

 

12. s=2,4, ХВ= 14,2, n=9,

1'=0.99.

 

Ont8. 11.512 < а < 16,888.

 

 

 

13. По данным 16 независимых

равноточных

измерений физичес-

кой величнны найдеиы хв =23.161

и s=O,4oo.

Требуется оценить

истиниое значение а измеряемой величины и точность измерений а

с надежностью 0,95.

Оmв. 22,948 < а< 23.374; 0,224 < 0<0,576.

14. Найти доверительный интервал для оценки нензвестной ве­

роятности р биномиального распределения с надежностью 0.95,

если в 60 нспытанияХ событие появилось 18 раз.

Omв. 0,200 < р < 0,424.

15. Найти методом моментов точечную оценку эксцесса Ek =

= mf,/a4-3 теоретического распределения.

Оmв. ek=тf,/a~-3.

18. Найти методом моментов точечные оценки параметров а. н ~

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

'(х)= ~a+lг~a+l) xCte-Х/р(а.> -1. Р>О, x~O).

у к а з а н и е.

С!1enать подстановку lI=x/fj

и,

используя гамма­

 

 

 

 

ос

 

 

 

функцию

r(n)=~xn-le-xdx,

найти сначала

м (X)=(a.+I)P,

 

 

 

 

о

 

 

 

D (Х) =

+ 1) ~2. а затем приравнять М (Х) =Хв,

D (Х)=DB

01118.

сх* =

-1

 

-

 

 

(xB/D B) - 1; Р* = DB/XB

 

 

17.

Найти

методом наибольшего правдоподобия по выборке хl.

х.' ... ,

хn точечную оценку неизвестного параметра Р гамма-рас­

пределения, еслн параметр а. известен.

 

 

у к а з а н и е.

Использовать

плотность

гамма-распределения.

приведенную в задаче 16.

 

 

 

01118.

~* =хв/(а+ 1).

 

 

 

Глава семнадцатая

МЕТОДЫ РАСЧЕТА СВОДНЫХ ХАРАКТЕРИСТИК

ВЫБОРКИ

§ 1. Условные варианты

Предположим. что варианты выборки располо­

жены в возрастающем порядке. т. е. в виде вариацион­

ного ряда.

Равнооmсmоящuмu называют варианты. которые обра· вуют арифметическую прогрессию с разностью h.

237

УС.ltООНЬUtи называют варианты, определяемые равен­

ством

иi = (x,-C)/h.

где С-ложный нуль (новое начало отсчета); h-шаг,

т. е. разность между любыми двумя соседними первона­

чальными вариантами (новая единица масштаба).

Упрощенные методы расчета сводных характеристик

выборки основаны на замене первоначальных вариант

условнымр.

Покажем. что если вариационный ряд состоит из равно­

отстоящих вариант с шагом h, то условные

варианты

есть ц е л ы е ч и с л а. Действительно, выберем в

качестве

ложного нуля произвольную варианту, например хm

Тогда

и, =

Х'-Хm

=

xl+<l-l)h-[xl+<m-l)h]

.

т.

h

 

h

 

= t -

Так как i

и т-целые числа. то их

разность

i =

=и/-также целое число.

3 а м е ч а н и е 1. В качестве ложного иуля можно принять лю­

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

примерно в середиие вариационного ряда (часто такая варианта

имеет наибольшую частоту).

3 а м е ч а н и е 2. Варианте, которая принята в качестве лож­

ного нуля, соответствует условная варианта, равная нулю.

Пример. Найти условные варианты статнстического распределения:

варианты • • • 23,6

28,6

33,6

38,6

43,6

частоты '"

5

20

50

15

10

Реше и и е. Выберем в качестве ложного нуля варианту 33,6

(па варианта расположена в середине вариационного ряда).

Найдем шаг:

h= 28,6-23,6=5.

Найдем условную варианту:

"1 = (xI-С)/h = (23,6-33,6);5= -2.

Аналогично получим: "1 = - 1, = О, и. = 1, иъ = 2. Мы видим,

'IТO условные вариантынеБОJlьшие целые числа. Разумеется, опе­

рировать с ними проще, чем с первоначальнымн варнантами.

§ 2. Обычные, наЧ8.JIьные и центраJlьные

8Мпирические моменты

Для вычисления сводных характеристик выборки

удобно ПOJlьзоваться эмпирическими моментами, опре..це­

Jlения которых аналогичны определениям соответствую-

2з8

щих теоретнческих моментов (см. гл. VHI, § 10). В от­

личие от теоретических змпирические моменты вычисляют

по данным наблюдений.

Обычны'м эмnирически'м ,Моменто,М порядка k называют

среднее значение k-x степеней разностей Х, -С:

M k= (~ni (Xj-C)k)/n,

где Хj-наблюдаемая варианта, n,-частота варианты,

n = ~nl-объем выборки, С-произвольное постоянное

число (ложный нуль).

Начальным Э'мnирически'м ,Мо,Менто,М порядка k назы­ вают обычный момент порядка k при С = О

M k = (~nIX:)/n'

в частности,

М1 = (~n,XI)/n = Хв,

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

равен выборочной средней.

Центральны'м Э'мnирическим ,Моментом порядка k на-

зывают обычный момент порядка k при С = ХВ

mk = {Dj (Х; - x B )1f)/n.

в частности,

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

Легко выразить центральные моменты через обычные

(рекомендуем читателю сделать это самостоятельно):

та =

M~- (M~)2,

 

тз=М~-ЗМ~М~+2(маа,

}

т" = M~-4M~M~

+ 6M~ {ма2 -З {ма".

 

§ 3. Условные 8мпирические моменты.

Отыскание центраJlЬНЫХ моментов по условным

Вычисление центральных моментов требует до-

вольно громоздких вычислений. Чтобы упростить рас­

четы, заменяют первоначальные варианты условными.

Условны..м Э'мnирически'м ,Мо,Менто,М порядка k называ­

ют начальный момент порядка k, вычисленный для ус-

2э9

ловных вариант:

 

 

 

~

(Xi-C)k

M = ~niuf = ~~i

h

k

n

В частности,

Отсюда

xB=M~h+C.

Таким образом, для того чтобы найти выборочную сред­

нюю, достаточно вычислить условный момент первого

порядка, умножить его на h и к результату прибавить ложный нуль С.

Выразим обычные моменты через условные:

• I

~ ni (Xi-C)k

мА:

M k = hk

n

= hk •

Отсюда

M k= Mkhk

Таким образом, ,идя того чтобы найти обычный момент

порядка k, достаточно условный момент того же порядка

умножить на hk

Найдя же обычные моменты, легко найти централь­

ные моменты по равенствам (**) и (***) предыдущего

параграфа. В итоге получим удобные для вычислений

формулы, выражающие центральные моменты через ус­

ловные:

mll=[M:-(М~)2]hl,

(**)

тз = [м;-зм;м~ + 2 (м~)з] h 3 ,

}

т" = [M:-4M:M~+6М;(М~)2_З (M~)4]h".

(***)

В частности, в силу (**) и соотношения (*)

предыду­

щего параграфа получим формулу для вычисления выбо­

рочной дисперсии по условным моментам первого и вто­

рого порядков

Техника вычислений центральных моментов по услов­

ным описана далее.

240