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

Дуплякин В.М. Статистический анализ

.pdf
Скачиваний:
43
Добавлен:
16.03.2015
Размер:
1.3 Mб
Скачать

1.3. Обобщённая теорема Чебышева

Теорема П.Л.Чебышева обобщается на более сложный случай независимых опытов в переменных условиях, а именно когда закон распределения случайной величины Х от опыта к опыту изменяется.

В этом случае мы имеем дело со средним арифметическим Y =

1

n

åXi случайных

 

n i=1

величин X1, X2, ... , Xn с различными в каждом опыте математическими ожиданиями mx1, mx2 ,...,mxn и дисперсиями Dx1, Dx2 ,..., Dxn .

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

число L , т.е.

 

Dxi £ L при i =1, 2, ..., n .

(1.21)

В такой постановке обобщённая теорема Чебышева сводится к следующему: "При возрастании числа независимых опытов n в переменных условиях среднее арифметическое наблюдаемых значений величин X1, X2, ... , Xn сходится по вероятности к среднему арифметическому их математических ожиданий".

 

n

 

n

 

 

 

 

å Xi

 

åmxi

 

< ε ) >1-δ , если D £ L (i =1,2,.. < n) ,

 

P (

i=1

-

i=1

 

(1.22)

 

 

 

n

 

n

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

где ε и δ − произвольные малые положительные числа.

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

Рассмотрим величину Y = 1 ån Xi . n i=1

Найдём для этой величины математическое ожидание и дисперсию, пользуясь известными

теоремами о числовых характеристиках

 

1

n

 

 

 

 

 

 

 

 

1

 

n

 

 

 

 

 

my =

åmx i

и

 

 

Dy =

åDxi .

 

 

(1.23)

n

 

 

2

 

 

 

i=1

 

 

 

 

 

 

 

 

n

i=1

 

 

 

Применим к величине Y неравенство Чебышева

P(

 

Y - m

y

 

³ ε ) £

Dy

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε 2

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åXi

 

åmx i

³ ε ) £

å Dx i

 

 

 

После подстановки получим

P (

i=1

-

i=1

 

 

 

 

 

i=1

 

.

 

(1.24)

n

 

n

 

 

 

n2ε 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

При замене каждой дисперсии Dxi

 

на гарантированно превосходящую величину L

неравенство может только усиливаться, поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åXi

 

 

 

 

åmxi

 

 

³ ε ) £

L

 

 

 

 

P (

 

 

i=1

 

-

i=1

 

 

 

.

(1.25)

 

 

 

 

 

n

 

 

n

n ε 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как бы ни было мало ε , можно выбрать число опытов n настолько большим, чтобы

выполнялось неравенство

L

£ δ , где δ − произвольное малое положительное число, тогда

n ε 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åXi

 

 

 

 

 

åmx i

 

 

³ ε ) £ δ .

 

 

 

 

P (

 

 

i=1

 

-

i=1

 

 

 

(1.26)

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходя к противоположному событию с вероятностью P = 1− δ , получим

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åXi

 

 

 

åmxi

 

 

< ε ) > 1-δ ,

 

 

 

P (

i=1

 

 

 

-

i=1

 

 

(1.27)

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что и представляет собой доказательство обобщённой теоремы Чебышева.

1.4. Теорема Маркова

Эта теорема представляет собой наиболее общий случай закона больших чисел в следующей формулировке:

"Если имеются статистически зависимые случайные величины X1, X2, ... , Xn для совокупности которых при неограниченном увеличении числа опытов n → ∞ соблюдается

 

1

 

n

 

 

условие

D êé

åXi úù

® 0 , то среднее арифметическое наблюдаемых значений сходится по

2

 

n

ë i=1

û

 

вероятности к среднему арифметическому их математических ожиданий".

Марков Андрей Андреевич (1856-1922)ближайший ученик П.Л.Чебышева, который внес особенно большой вклад в становление и развитие теории случайных процессов, специфическая разновидность которых вошла в

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

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

зависимость включает в себя как частный случай функциональную зависимость и поэтому является более широким понятием.

11

Формализованное оформление теоремы Маркова выглядит следующим образом

 

n

 

n

 

 

 

 

 

 

 

å Xi

 

åmxi

 

 

1

n

 

 

P (

i=1

-

i=1

 

< ε ) >1-δ , если при n ® ¥ имеем

D êéå Xi úù

® 0 . (1.28)

n

n

2

 

 

 

 

n

ë i=1

û

 

 

 

 

 

 

 

 

 

 

 

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

Рассмотрим среднее арифметическое наблюдаемых значений Y = 1 ån Xi . n i=1

Найдём для этой величины математическое ожидание и дисперсию, пользуясь известными

 

my = 1

n

 

 

 

1

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что

åmx i

и

 

Dy =

åDxi .

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

n

i=1

 

 

 

n i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применим к величине Y неравенство Чебышева P(

 

Y - m

y

 

³ ε ) £

Dy

.

 

 

 

(1.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

 

По условию теоремы при n → ∞ неограниченно убывает величина

 

D êé

åXi úù

® 0 , т.е.

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

ë i=1

û

 

неограниченно

уменьшается

дисперсия Dy

® 0 ,

поэтому

 

неравенство

Чебышева можно

представить в

виде

P(

 

Y - my

 

³ ε ) £ δ ,

где

δ − произвольная

малая

положительная

 

 

величина.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходя к противоположному событию, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åXi

 

 

åmxi

< ε ) > 1-δ ,

 

 

 

 

 

 

 

 

 

 

 

 

P(

i=1

-

i=1

 

 

 

 

(1.30)

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что и требовалось доказать.

1.5. Теорема Я. Бернулли*

Известная теорема Якоба Бернулли устанавливает связь между частотой появления наблюдаемого события и его вероятностью:

"При неограниченном увеличении числа независимых опытов в постоянных условиях,

т.е. при , n → ∞ частота события А, а именно p* сходится по вероятности к его истинной вероятности p ".

* Якоб Бернулли (Jakob Bernoulli, 1654-1705)знаменитый швейцарский математик.

12

Частота случайного события p* здесь, как это принято в теории вероятностей,

понимается в виде отношения числа опытов m благоприятствующих появлению данного события к общему числу опытов n , т.е.

 

 

 

 

p* = m .

(1.31)

 

 

 

 

 

n

 

Формально теорема Я.Бернулли записывается в виде неравенства, которое

обеспечивается при увеличении числа наблюдений ( n → ∞ )

 

 

 

P(

 

p* - p

 

< ε ) > 1- δ

(1.32)

 

 

 

 

для произвольных малых положительных чисел

 

 

ε и δ .

 

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

 

 

 

 

 

 

 

 

Обозначим

случайный

результат i го

опыта как Xi .

Эта дискретная случайная

величина может

принимать

два значения: Xi1 = 0, когда опыт не привёл к желаемому

результату и Xi2 = 1, если i й опыт закончился с интересующим нас результатом.

Как и для всякой дискретной величина, её распределение характеризуется статистическим рядом, который в данном случае имеет вид

Xi k

0

1

pk

1− p

p

 

 

 

Математическое ожидание величины Xi в отдельно взятом опыте вычисляется как

2

 

mxi = åxi k × pk = p .

(1.33)

k =1

 

Частота представляет собой среднее арифметическое величин X1, X2 , ...,

Xi ,..., Xn

n

åXi

p* =

i=1

.

(1.34)

n

 

 

 

13

Согласно теореме П.Л.Чебышева при независимых опытах в постоянных условиях среднее арифметическое сходится по вероятности к математическому ожиданию, что в

нашем случае записывается следующим образом

P(

p* p

< ε ) > 1− δ .

(1.35)

Именно это и требовалось доказать.

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

1.6. Теорема Пуассона*

При изменяющихся условиях опыта так же имеет место устойчивость частоты появления события, что формулируется в тереме Пуассона:

"Если производится n независимых опытов и вероятность появления события А в

i м опыте равна pi , то при увеличении числа опытов частота появления события А

сходится по вероятности к среднему арифметическому вероятностей pi

".

Формализованная запись теоремы Пуассона имеет следующий вид

 

p* 1

n

 

 

 

 

 

 

P(

å pi

 

< ε ) > 1−δ ,

(1.36)

 

n

i=1

 

 

для любых малых положительных величин ε и δ , если n → ∞.

Доказательство этой теоремы аналогично доказательству предыдущей теоремы, но при этом используется обобщённая теорема Чебышева, применяя которую необходимо убедиться в ограниченности дисперсий.

Обобщённая теорема Чебышева декларирует сходимость по вероятности среднего

арифметического наблюдаемых значений величин X1, X2, ... , Xn к среднему

арифметическому их математических ожиданий при независимых опытах в переменных условиях и ограниченности дисперсий

 

n

 

n

 

 

 

å Xi

åmxi

 

< ε ) >1−δ , если D L (i =1,2,.. < n) , (1.37)

P (

i=1

i=1

 

 

 

 

n

 

n

 

xi

 

 

 

 

 

 

 

 

 

 

где ε и δ − произвольные малые положительные числа.

* Simeon Denis Poison (1781-1840)выдающийся французский математик и физик, известное распределение

которого и в настоящее время широко используется в теории массового обслуживания при моделировании потоков событий.

14

Применительно к частоте события на основании обобщённой теоремы Чебышева можно записать

 

 

1

n

 

 

P (

p* -

å pi

< ε ) >1-δ , если Dpi £ L (i =1,2,.. < n) .

(1.38)

 

 

n

i=1

 

 

Таким образом, что бы доказать сходимость частоты события к его вероятности, в

данном случае нужно доказать ограниченность дисперсий статистических оценок

вероятностей, т.е. убедиться в выполнении условия Dp £ L (i = 1,2,.. < n) .

 

 

 

 

 

 

i

 

 

 

 

 

 

 

é

- mp i

)

2

ù

Обратимся к анализу дисперсии оценки вероятности Dpi = M ë( pi

 

û.

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

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

 

 

 

 

0 £ pi £ 1, 0 £ mp £ 1,

 

 

 

 

 

 

 

i

 

 

 

 

é

- mp i )

2

ù

 

 

 

 

то из этого следует ограничение Dpi = M ë( pi

 

û £ 0,25 .

 

 

 

 

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

15

2. ВЫБОРОЧНЫЕ ОЦЕНКИ И ИХ СВОЙСТВА

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

Любая из выборочных оценок случайна и при её использовании возможны ошибки, поэтому желательно выбрать такие оценки искомых характеристик, что бы эти ошибки были по возможности минимальными.

2.1. Требования к выборочным оценкам

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

Выполнено n независимых опытов и получены значения рассматриваемой случайной

величины

X1, X2 , ..., Xi ,..., Xn .

Выборочная оценка a* является функцией наблюдаемых случайных величин и

поэтому сама является случайной величиной

a* = a(X1, X2 , ..., Xi ,..., Xn ) .

Закон распределения величины a* зависит от самого a , от наблюдаемых значений

X1, X2 , ..., Xi ,..., Xn , от числа опытов n и может быть найден известными методами теории вероятностей, что представляет сложную в математическом плане задачу, но для многих практически значимых случаев имеются известные решения.

Исходя из общей постановки задач математической статистики, к выборочным оценкам a* предъявляются требования состоятельности, несмещённости и эффективности.

16

1. Состоятельность сходимость по вероятности к истинному значению при увеличении числа опытов:

P( a* - a < ε ) >1-δ ,

(2.1)

для любых малых положительных величин ε и δ , если n ® ¥.

2. Несмещённость отсутствие систематическое ошибки (постоянное отклонение в одну сторону), что выражается в равенстве математического ожидания выборочной оценки

истинному значению данной характеристики

é

* ù

= a .

(2.2)

M ëa

û

3. Эффективность минимальное рассеивание (минимальная дисперсия)

é * ù

® min .

(2.3)

D ëa û

2.2. Свойства выборочных оценок математического ожидания

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

В качестве статистической оценки математического ожидания используем среднее арифметическое наблюдаемых значений случайной величины Х

 

1

n

 

mх* =

åxi .

(2.4)

 

 

n n−1

 

Тем не менее, в этом следует убедиться, проверив состоятельность, несмещённость и эффективность предлагаемой статистической оценки дисперсии.

2.2.1. Состоятельность поскольку речь об оценке математического ожидания,

декларируется следующим образом

P(

m* - m

< ε ) >1-δ ,

 

 

x x

 

(2.5)

для любых малых положительных величин ε и δ .

 

Это свойство выборочных оценок математического ожидания фактически устанавливается теоремой Чебышева, в которой утверждается, что среднее арифметическое

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

 

n

 

 

P(

1 åXi - mx

< ε ) > 1- δ.

(2.6)

 

n i=1

 

 

2.2.2. Несмещённость – M[m*x ] = mx .

(2.7)

17

Для случайной величины, которая является линейной функцией других случайных

 

1

n

 

 

 

 

 

 

величин, а именно mх* =

åxi

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

 

 

n n−1

 

 

 

 

 

 

вероятностей известно что M êé

n

 

n

 

, но, поскольку рассматриваются опыты в

1 åxi úù =

1 åmxi

 

 

ën n−1

û

n n−1

 

 

 

неизменных условиях, т.е. когда mx i - const,

то mx i = mx , то в итоге имеем доказательство

свойства несмещённости выборочной оценки математического ожидания

 

 

 

 

 

M êé

n

 

= 1 n ×mx = mx .

 

 

 

 

 

1 åxi úù

(2.8)

 

 

 

 

ën n−1

û

n

 

2.2.3. Эффективность D[m*x ] ® min .

 

(2.9)

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

соотношениям

é

1

n

ù

 

1

n

1

 

D

 

D ê

 

åxi ú

=

 

 

åDxi =

 

 

n × Dx =

x .

(2.10)

 

n

2

n

2

ën n−1

û

 

 

n−1

 

 

n

 

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

2.3. Свойства выборочных оценок дисперсии

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

 

1

n

1

n

 

Dx* =

å(xi - mx* )2 , где mх* =

åxi .

(2.11)

 

n i=1

n n−1

 

2.3.1. Состоятельность т.к. речь идёт об оценке дисперсии, то её состоятельность

формулируется следующим образом

P(

D* - D

³ ε ) £ δ ,

(2.12)

 

x x

 

 

для любых малых положительных величин ε и δ при неограниченном числе опытов.

Попробуем убедиться в состоятельности оценки дисперсии, обратившись к теореме

Чебышева для некоторой величины Y

 

1

n

 

 

P(

åYi - my

³ ε ) £ δ .

(2.13)

 

n i=1

 

 

18

В качестве Y возьмём второй начальный момент Y = α2x . Выборочная оценка этого

момента вычисляется как

α2*x =

1

n

 

åxi2 .

(2.14)

 

n

i=1

 

Как видно, выборочный второй начальный момент находится как среднее арифметическое квадратов отдельных значений рассматриваемой случайной величины,

потому данная оценка сходится по вероятности к своему истинному значению по теореме Чебышева, т.е.

P(

α2*x -α2x

³ ε ) £ δ ,

(2.15)

для любых малых положительных величин ε и δ при неограниченном увеличении числа опытов.

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

P(

(m* )2

- (m )2

³ ε ) £ δ ,

(2.16)

 

x

x

 

 

для любых малых положительных величин ε и δ при неограниченном увеличении числа опытов.

Объединяя два предыдущих неравенства в одно и выполнив перегруппировку, и

учитывая соотношения

D* = α*

- (m* )2

и D =

2x

- (m )2

,

(2.17)

x 2x

x

x

x

 

 

получим следующее неравенство показывающее сходимость статистической оценки дисперсии к своему истинному значению

P(

Dx* - Dx

³ ε ) £ δ .

(2.18)

2.3.2. Несмещённость а именно то, что M[Dx*] = Dx .

Выразим статистическую оценку дисперсии через оценки второго начального момента и математического ожидания

 

 

 

*

 

 

 

 

*

 

 

-

*

2

=

1

n

 

 

2

-

 

æ 1

 

n

 

ö2

 

 

 

 

 

 

Dx

 

= α2 x

 

(mx )

 

 

n

åxi

 

ç

 

 

 

åxi ÷ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

è n

 

i=1

 

ø

 

 

 

 

(2.19)

= 1

n

 

1

 

n

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

n -1

 

 

n

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åxi2 -

åxi - 2

åxi xj =

 

åxi2 - 2

åxi xj .

 

2

2

2

2

 

n

i=1

 

n

 

 

 

i=1

 

 

 

 

 

n

 

 

i< j

 

 

 

 

 

 

 

 

n

 

 

 

i=1

 

 

 

 

 

n

i< j

 

Найдём математическое ожидание полученного выражения для оценки дисперсии

 

M éD*

ù

= n -1 n M

éx2 ù

- 2

1

 

 

 

 

 

 

M éx x

ù.

 

(2.20)

 

 

 

 

2

å

 

 

 

ë

 

x

û

 

n

2

å

 

 

 

ë

i

û

 

 

 

 

 

n

 

ë

i

 

j û

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i< j

 

 

 

 

 

 

 

 

 

Перейдём к центрированным значениям x&i

= xi - mx ,

 

 

 

так как дисперсия не зависит от

того, в какой точке выбрать начало координат

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

éD*

ù

 

= n -1 n M éx&2

ù - 2

 

 

1

 

 

å

M

éx&

 

x&

j

ù.

 

(2.21)

 

 

 

 

 

 

2

 

 

 

 

 

ë

x

 

û

 

 

n

2

å

 

 

ë

i

û

 

 

 

n

 

 

 

ë i

 

 

û

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

i< j

 

 

 

 

 

 

 

 

 

 

Очевидно что

n

 

M éx&

 

 

 

 

= n × D

,

 

 

 

M

éx&

x&

 

=K

 

.

 

 

(2.22,

2.23)

å

2

ù

 

å

ù

i j

 

 

 

 

 

 

 

 

ë

i

û

 

 

 

 

x

 

 

 

 

 

 

 

ë

 

 

i

 

 

 

j û

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i< j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19