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

книги из ГПНТБ / Мания, Г. М. Статистическое оценивание распределения вероятностей

.pdf
Скачиваний:
21
Добавлен:
20.10.2023
Размер:
7.77 Mб
Скачать

 

 

Отсюда,

при t — 02 — 0Х будем иметь .

 

 

Ш 2- Л У =

 

 

 

 

 

 

 

 

 

 

=

ехр 1

(к -

 

к ?

1 [

 

_

л

[:(4v +

2) х -

к г-

к Г

2 (02 -

0J

 

 

ехр

2 (03-

0а)

 

 

 

 

 

 

7= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

[ (4 v +

4) X + ^ — /2]2

 

 

 

 

 

 

ехр

 

 

 

(1.4.34)

 

 

 

 

 

2 (02

0Х)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(®г

 

 

®i) =

 

 

 

 

 

 

 

 

 

 

— ехр

 

(к ~

 

к?

 

2

ехр

[(4 v +

2) X +

tx +

t2]2

 

2 ( 02

 

0Х)

 

2(02 -

ВО

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7=

0

 

 

 

 

 

-

V

 

ехр Г _

[ ( 4 у +

4 ) А - < , -

t . f i

 

 

(1.4.35)

 

Н

 

 

I

 

 

ЧЧ.-Ч)

J

 

 

 

 

 

 

 

 

 

 

 

 

Из леммы 1.3

и формул (1.4.32), (1.4.33), ,(1.4.24) заклю­

чаем,

 

что

 

 

 

 

 

 

 

 

 

 

и

 

 

 

= Д (02 - 0х) + « 0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P m

=

f%( 0 2 ~ ~ ® i ) +

0 ( 1 ) •

 

 

 

 

 

 

Следовательно, в силу (1.4.21),

 

 

 

 

 

1 Рсs

(Iх; ») =

к (02-0х) +

/2 (02 -

0i) 4- ®п,

(1.4.36))

где еп- » 0 , равномерно относительно tx и /2 из любого ко­ нечного интервала.

Таким образом, мы нашли второй сомножитель общего, члена суммы (1.4.3), а для первого сомножителя имеем

Рп(%, «2) =

п !

6М 02 -

0l)S2“ Sl ( l -

е2)n-Sp

 

Si! («2~ si)! (« -

«а) 1

 

 

_L

А ^ А ^2 •:

 

 

" «

V"0a(02— 91)(1 — ©2)

 

 

Х е х р

_ <*» ~ ^)2

_____ J*____

(1 + йл) >

(1.4.37)»

2 01

2 (02 — 02)

2 ( 1 - 0 ,)

 

 

€0

тде

по-прежнему

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sx =

п

- f itl V

n ,

 

 

— X <

<

X ,

 

 

 

 

 

 

 

s2 =

 

n 02+

t2V~n ,

 

 

- X < 1 2< X ,

 

 

 

 

 

 

 

 

Д tx =

Д t2

=

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Y~n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

an-> 0

при

я с о .

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, в силу (1.4.36) и (1.4.37),

 

 

 

 

 

 

lim Рп (0<?>, 0<jft

X) =

Ф (01э 02; X) =

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

X

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ехр

 

 

 

 

 

-------- ——

} х

2 те ]/ 0i(020х)(1 62),

 

'

2 0*

2(02-

0!)

 

 

 

2 ( 1-

02) ' Л

 

 

 

 

 

 

 

—X —X

 

 

 

 

 

 

 

 

 

 

X - [1 -

/х (02 -

9х) -

h (02

-

0х)1 di t d U .

 

 

 

 

 

 

 

 

Учитывая

(1.4.34)

и (1.4.35)

и раскрывая скобки под

ин­

тегралом,

будем

иметь

 

 

 

 

 

 

 

 

 

 

 

 

Ф (01, 021X) =

2 т е ]/ 01(02 -

01)

(1

02)

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

7

 

 

 

п

 

 

%

-

Ч 2

п

 

 

 

 

 

 

x

j

J « P ( -

 

 

 

 

 

dtр

di*2

 

20!

 

 

2 (0а -

01)

2 ( 1 -

 

0.)

 

 

 

 

—X

—X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

оо

X X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ехр

 

 

J ]

X

 

2 те ]/ 01(03 -

 

01)(1 -

02)

 

 

 

 

 

 

 

v=0 — x — x

 

 

 

2 01.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X ехр j —

2(1

*t

02)

 

1(4v +

2)X — tx -

 

 

 

 

 

 

-

 

 

 

 

2 (02- 0x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— ехр

 

n

 

 

 

 

n

 

 

 

[(4 v -f- 4) X +

^

/2]2

]+

 

2 0i

 

 

2 ( 1 - 0 0

 

 

 

2 (0a -

 

0i)

 

 

 

 

 

 

 

 

 

 

 

 

 

+

exp

 

J l _____

 

f2

 

 

 

,[(4 V +

2) X +

;* +

<,!»

 

 

 

 

fc2

 

 

 

 

2 (02-

 

0i)

 

 

 

 

201-

2(1 -

02)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

61

 

exp

Г

t\

 

 

%

[ ( 4 v + 4 ) X - t i + t , ] r

dt1dt2

 

 

20x

2(1 -

02)

 

2 (0*— Эх)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

X

 

 

 

 

 

 

2 « K

 

 

 

 

 

 

 

 

 

 

 

0x <0* — ©i) (1 - 02)

 

 

 

 

 

 

X

 

X

 

 

 

 

 

 

/l)2

 

 

 

 

X

 

 

 

exp

 

t\

 

 

(t2 -

 

 

d / i d 12 —

 

 

 

 

 

 

 

 

 

 

 

 

■X

— X

 

 

2 0 x

 

 

2 ( 0 2 -

0 x )

 

2 ( 1

- 0 , )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

X.

 

 

 

 

 

2те У

0i (02 -

0i)(A

-

 

 

 

 

 

ea)

 

 

 

 

 

 

 

00

X X

 

 

 

t?

 

t\

 

R4v+ 2 ) X - / x+ *2]2

X 2

> ,

 

;exp

 

 

 

 

' 2

 

 

 

 

26x

2(1 6,)

 

2 (02-

0!)

 

v^O—x —x

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

t\

 

[ (4 v + 4) X — tx— (212

d'ti dt2 .

 

 

20x

2 (1 - 02)

 

2 (0, — 0x)

 

 

 

 

 

 

 

(1.4.38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразуем это

выражение. Полагая,

как и раньше,

 

 

 

 

 

 

h zi К 0х<1 0х) .

 

 

 

 

имеем

 

 

 

t2 =

22У 02(1 — ©i ),-

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

2

75 К

 

©х (02 -

© х ) ( 1

-

е 8)

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

X

 

 

 

 

 

ti

 

 

 

 

 

 

 

 

 

 

 

п

__

 

__

{ /

^l)2

 

 

 

 

 

 

 

___________

<d ^x d 12

X

J

 

J

“ P

i -

2 6x

 

2(1 -

08)

 

2 (02— ©x)

 

— X — X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

X

 

 

 

 

 

 

 

 

 

 

/ 01(1- 01)

/ 0в(1- в 1)'

 

 

 

 

 

' 2 теV "1 -

tf2

J

 

 

 

exp

{; —

у

0- (2j,- Zj)

} dZi dz.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 0/ l - e i )

 

/ е 2( . i - 02)

 

 

 

 

62

где

D _

I

 

 

^2)

л

у \_

I

(«! - 2 R zxz2 + z\).

 

 

■j—

к

V е , ( 1 - в а) ’

e (

l ’ 2' ~

 

 

(1.4.39)

 

Далее, нетрудно видеть, что

 

 

 

 

 

 

 

_____ t\___________ tj

 

_

[ ( 4 у +

2 ) X -

tr - f8]«

 

 

2e2

2(1 -

02)

 

2(es. - e 1)

 

=

_ i

- 0

( Zl, z , ) - 2 ( 2 v +

1)*X*,

 

где

переменные tx и t2 связаны c

zx и г2 следующими соотно­

шениями:

 

 

 

 

 

 

 

 

 

 

к = (4 v + 2)Х вх + гх Т/ 91 (1 - 0г) г

 

 

к =

(4v +

2) X (1 -

02) +

z21/ 02(1 - 02) •

 

Полагая теперь

 

 

 

 

 

 

 

= (4 v + 4) X 0j + zx У (1 — 0г) г

(1.4.40)

 

к =

 

 

 

 

 

z2У02(1

 

 

(4 v +

4) X (1

— 02) +

— 02)

 

будем иметь

 

 

 

 

 

 

 

 

 

_____ к_____ tl

__

[ ( 4 у +

4 )Я, — к — ^ ]2'

 

 

2

0!

2(1 -

0Я)

 

2 (02

0,)

 

=

-1

0 (Zl, z2) -

2 (2 v +

2)*Х*.

 

Преобразовывая интегралы в правой части формулы (1.4.38)

спомощью подстановок (1.4.39) и (L4.40), получим

Ф(0J, 02; X) =

V 0а(10а)

1 0 (zx, z2) \dz1 dz2 —

ехр

2-п У 1 — R2 J

 

X

 

/ M l - Л ) У М 1- У '

 

63

2 л ] Л

- R2

X.— 2 (2 v +

1) X 6,

/ 9 , 0 - 6 , )

X |

exp {— 2 (2 v + l ) 2X2} X

2

V=0

Я— 2 (2 v - f - 1) X (1 — 92)

/02( 1- 02)

j’ expj — -i-0(z1,z2)jd e1dz2 +

— X — 2 (2 v 4 - 1)X0,

— X — 2(2v* + 1) X (1 — 0a)

 

/ 9 , (1- 6,)

 

/ 6 2 ( 1 - 6 2 )

+

 

1

 

2 V exp { — 2 (2 v + 2)2X2} X

2izV

1- Я 2

2

 

 

 

 

 

v-0

 

X — 2 (2 V +

2) X 0,

X

2 (2 V + 2) X (1 0a)

 

/ 0, 0 - 0,)

 

/ 62(1- 62)

X

 

 

 

 

1

 

 

 

 

exp I ------ — 6 (z1,z2)dz1 dz2

— X — 2 (2 V + 2) X 0, — X — 2 (2 v + 2 )X (l-

 

/ 6, ( 1-H ,)

 

/ 62(1- 62)

 

Отсюда

непосредственно следует

Ф (0„ 02; X) =

X

X

 

 

 

 

 

 

 

/ 6 , ( 1 6,) / 0 2 (1 - 02)

 

 

 

X

 

X

 

 

 

/ 9 , 0 - 9 0

/ 62( 1- 62)

 

 

 

 

CD

 

 

г-Х у .’ -

У

г 2

( ~ ‘ )>' 1 е х р 1 ~ 2 *П , , Х

 

 

 

 

fe=l

/

0, 0 - 0,)

X — 2k\(l — (

 

/ 02(1

X

 

 

 

 

exp { -------— 0(2!, zs) } dzxdz2.

— X 2kX 0,

X 2 fe X (1 — 0a)

 

/ 9 , 0 - 9 0

 

/ 62(1- 63)

Тем самым теорема доказана.

64

Теперь рассмотрим частные случаи.

Если

0Х = 0

и

02 =

1,

то

 

 

 

со

 

 

И т Рп (0,

1; X) =

1 -

2 V

( -

l )^1exp { - 2fc2X2} = K(X) .

Л — оо

 

 

 

 

 

Й=1

Если 0Х =

0

и

02 = 0,

то

 

 

 

 

 

 

/

0 ( 1 -

0)

 

 

lim

(0’ , ; 1 ) “

у

к

1 е х р | -

т Н

г

п-+<х>р

 

 

 

 

К 0 ( 1 - 9 )

X— 2feX (16)

 

 

 

 

 

 

 

00

 

 

 

 

/

0(10)

i =

2 ^

( -

 

1)*-* exp { - 2 fe2X2] j

exp { -

]/2

k=\

 

 

 

 

— X — 2 fe X (1 — 0)

 

 

 

 

 

 

 

 

 

 

 

V

0 ( 1 - 0 )

Если,

наконец,

0X =

1

02 = 0,

lim Pn (0,

1

- 0; X)

= Ф (0;

X) =

/!->■ oo

 

X

 

 

X

 

 

 

 

 

 

У 0 ( 1 - 0 )

У 6 ( 1 - 0 )

 

1

Г

 

 

I e x p {

 

 

 

 

2 7C |/ 1 — R2 J

 

 

 

 

X

 

 

X

 

 

/ 0 ( 1 - 0 )

/

0 ( 1 - 0 )

 

1

 

OO

 

 

 

■ 9

(

-

i)* -1 X

2 tzV

1

- - R 2

£ = I

 

 

 

 

 

 

 

 

 

X — 2 * X 0

X — 2 Ы 0

 

 

/ 0 ( 1 - 0 )

У 0 ( 1 - 0 )

to

y 9 ( Z i , г2)

j dz

d

X e x p {— 2& 2 X2} l

 

1

 

d z 1d ?2 ,

j* exp

2 0

(Zi, z2)

 

 

— X — 2fe X0

— X — 2fcX0

/0 ( 1 — 0)

/ 0 ( 1 - 0 )

5. Г. M. Мания

65

где

 

 

 

1 -

0 *

 

 

 

 

 

 

В

приложении к

книге

приведена

таблица

для

функции

Ф (0, X), [25] (см. табл. 2). С

помощью

этой

таблицы при дос­

таточно больших п(п ^

100)

можно указать

верхний и нижний

пределы отклонения Sn(х) — F (%) на интервале

[х: 0 ^

F (х) ^

^ 1

0} ,

гарантируемые

с

наперед

заданной вероятностью.

Пусть,

наприхер,

0

=

0,25,

X =

1, 6;

тогда для доста­

точно больших п с вероятностью

0,,9883

можно

гарантировать

неравенство

 

 

 

 

 

 

 

 

 

 

-- i A < S n( x ) - F ( x ) ^ - b L

У п

У п

для всех х из интервала :0,25 ^ F (х)

0,75) .I

I

J

Г Л А В А II

НЕПАРАМЕТРИЧЕСКАЯ ОЦЕНКА ПЛОТНОСТИ

 

 

РАСПРЕДЕЛЕНИЯ

 

 

 

Первые результаты в области непараметрическсго оценива­

ния плотности распределения принадлежат советским математи­

кам В. И. Г л и в е н к о и Н.

В. С м и р н о в у ,

которые в каче­

стве оценки рассматривали

гистограмму.

В.

И. Г л и в е н к о

установил, что при

некоторых условиях

гистограмма с вероят­

ностью I равномерно сходится к теоретической

плотности

[5], а

Н. В. С м и р н о в

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

мак­

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

граммы от теоретической плотности [58].

 

 

 

 

 

 

 

За последние

годы вопросы

непараметрическсго приближе­

ния плотности получили

ноЕсе

ссЕепение в работах

Р о з е н -

б л а т т а

[96], Н. Н.

 

Ч е н ц о в а

[66],

П а р з е н а

[90],

Э. А. Н а д а р а я

[47], У о т с о н а

и Л и д б е т т е р а

[105]

и др.,

изучивших

свойства „обобщенной гистограммы'1, т.

е.

статистики

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

где Хг, ..., Хп—выборка

из

генеральной

совокупности

X с

не­

прерывной

плотностью

f(x),

ядро

К{х) некоте рая бсрелеЕская

функция,

интегрируемая

 

по мере

Лебега,

=

—последова­

тельность

положительных

чисел,

стремящаяся к нулю,

причем

В § 1

доказывается

теорема Э. А.

Н а д а р а я

[47] об усло­

виях равномерной

сходимости

fn(x)

к

f(x)

с вероятностью

1.

Далее, в

§§ 2 — 3

следуют его же

результаты [48, 49], позво­

ляющие оценить

точность

fn(x)

и при

заданном

коэффициенте

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

67

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

ДЛЯ f(x).

В § 4 исследуется локальная и глобальная точность одной непараметрической оценки плотности двумерного распределения, которая является двумерным аналогом оценки Р о з е н б л а т т а

[96].

 

 

 

§ 1. РАВНОМЕРНАЯ СХОДИМОСТЬ /„ (*) К /(*)

 

 

 

П а р з е н о м

[90] было доказано,

что

при некоторых усло­

виях fn(x)' по вероятности

сходится

к f(x) равномерно на всей

прямой. Условия сходимости с вероятностью 1 устанавливает

 

 

Т е о р е м а

2.1 [47].

Если

 

 

 

 

 

1°

К(х)—функция с ограниченным, изменением и JK ( x ) d x = l ;

(

00

 

 

 

 

 

 

 

2°

ряд

ynfl2

 

 

 

 

 

 

е

сходится при любом у >

 

3° f(x) равномерно непрерывна на всей оси,

 

то с вероятностью 1

 

 

 

 

 

 

 

 

 

 

Vn=

sup

I fn(x)-f(x) |

 

сходится к нулю.

x£Rl

 

 

 

 

 

 

 

 

 

 

 

^

Д о к а з а т е л ь с т в о .

Пусть Vn

sup* I /„ (x) — E fn{x) J.

М /

 

 

 

 

 

 

x^R1

 

Покажем,

что

с вероятностью

1. В самом деле,

применив

формулу интегрирования по частям, находим,, что

 

 

 

sup

I fn(x )-E fn(x) | =

sup

 

 

 

 

 

 

: x£Rl

 

x^R1

 

 

 

 

 

 

 

 

 

 

 

 

 

oo

 

 

-

T . f

 

 

^

sup

Т | | Э Д -

 

 

 

x£Ri

 

 

 

 

 

dK

x —u

<

sup

\ S n( x ) - F ( x ) \

- I j i ,

 

 

 

h

 

 

 

 

 

Vx £ R L

i

'

h

 

 

 

 

 

i

 

 

.

>

(2. 1Л)

68

где [х= Var K(x), Sn(x)- эмпирическая функция распределения, построенная по выборке, a F{x) ~ функция распределения X.

Нетрудно заметить, что

p { d „(0. 4 >

 

- ^

} « : p

{ « (0. i » - ^ }

+

+ P { d - ( 0

,

1 )>

y = r J

= 2 (1 -P * (X )),

(2.1.2)

D~(0,

1) были

определены

в гл. I).

Из

известной

формулы

Н.

В.

С м и р н о в а

[57] (см.

(1.1.3)),

как

показали

Д в о р е ц ­

кий ,

К и ф е р

и В о л ь ф о в и ц [78],

следует неравенство

 

 

 

 

 

1

—aX2

 

(2.1.3)

 

 

 

 

 

 

 

 

где

0 < а ^С 2,

a

ch

1, здесь и ниже означает положительную

константу. Тогда

из (2.1.2) и (2.1.3) имеем

 

 

 

 

 

 

 

 

 

 

 

(2.1.4)

 

В силу (2.1.1) и (2.1.4),

 

 

 

 

 

 

 

 

P { sup |fn(x)~ E fn(x) |>

г] <

 

 

 

 

1

 

x^R1

 

 

 

 

 

 

 

 

D

 

 

 

 

< P

 

sup [ 5

 

 

 

 

 

 

 

{ x^ R 1

 

 

 

 

 

 

 

Неравенство (2.1.5) с применением леммы Бореля-Кантелли

позволяет утверждать, что Vn стремится

к нулю с вероятностью

1 при п—>~со.

 

 

 

 

Для доказательства

теоремы остается

показать,

что

sup ] Е /„(* )—f(x) |->-0.

Пусть S > 0 .

Тогда,

полагая

М —

x^R1

 

 

 

 

= max f(x), имеем

 

 

 

'

x£R'

 

 

 

 

 

 

 

 

Соседние файлы в папке книги из ГПНТБ