Высшая математика. Часть III. Теория вероятностей. Математическая статистика. Учебное пособие
.pdf
2) Найдем параметр tэкс ( Zэкс)
Формулы расчета. Для малых выборок ( n1 30, n2 30) |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
x − y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
tэкс= |
|
|
|
|
|
|
|
|
nxny (nx + ny |
− 2) |
, |
||||||||||||||||||
|
(nx −1)Sx2 + |
(ny −1)Sy2 |
nx + ny |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
где nx , ny – объемы выборок. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
Для больших выборок находится параметр tэкс по формулам |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
x |
− y |
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|||||
|
|
|
|
tэкс= |
|
|
|
|
|
|
или |
tэкс= |
|
|
x |
y |
|
|
. |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
S |
2 |
|
|
Sy2 |
|
Sx2 |
+ Sy2 |
||||||||||||||||||
|
|
|
|
|
|
|
|
x |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
nx |
|
ny |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
В |
данной |
задаче nx <30, ny <30, |
поэтому вычисляемtэкс по |
||||||||||||||||||||||||||
формуле 1: |
|
|
90 − |
80 |
|
|
|
|
|
|
12 |
15(12 +15 −2) = |
||||||||||||||||||
|
tэкс= |
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
(12 −1) 5 +(15 −1) 8 |
|||||||||||||||||||||||||||||
|
|
10 |
12 +15 |
|
|
|
|
|
||||||||||||||||||||||
= |
|
|
166,7 = |
|
5 |
|
|
|
12,9 = 5. |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
55 +112 |
|
|
|
|
|
12,9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
3) Найдем критическую точку распределения Стьюдента tкр по уровню значимости 0,05 и числу степеней свободы k = nx +ny – 2 = 12 +
15 – 2 = 25. По таблице 2 Приложения: tкр = tα,k = t0,05;25 = 2,06.
4) Сравниваем tэкс и tкр. Т.к. tэкс = 5 › tкр = 2,06 = 2,06 , то нулевая гипотеза H0 отвергается, (если tэкс<tкр– H0 принимается).
5) Выводы: Различие средних показателей лечения двумя методами существенно на 5% -ном уровне значимости. С доверительной вероятностью 95% можно утверждать, что новая методика лечения более
эффективна, чем старая.
49
9. ЭЛЕМЕНТЫ КОРРЕЛЯЦИОННОГО АНАЛИЗА
Функциональная, статистическая и корреляционная зависимости. Цель и задачи корреляционного анализа.
Для решения многих задач экономики и медицины требуется установить и оценить зависимость между изучаемыми случайными величинами x и y .
Две случайные величины x и y могут быть связаны функцио-
нальной зависимостью, статистической или быть независимыми. Функциональной зависимостью между величинами x и y называ-
ется зависимость, при которой каждому значениюx из множества её допустимых значений соответствует одно или несколько значений
величины y и обозначается y = f (x). Например, зависимость между
|
|
4 |
|
3 |
|
объемом шара и его радиусом V |
= |
|
πR |
|
является функциональной, |
3 |
|
||||
|
|
R |
|
|
|
т.к. каждому значению радиуса шара |
|
|
соответствует определенное |
||
значение его объема V .
Функциональная зависимость между величинами x и y в
экономике и медицине реализуется редко, т.к. обе величины или одна из них подвержены действию случайных факторов. Например, если x - количество вводимого препарата, то его концентрация в крови ( y ) в
произвольный момент времени определяется не только количеством введенного препарата, но и массой больного, скоростью выведения препарата из организма и т. д. В этом случае между величинами возникает статистическая зависимость.
Статистической зависимостью величины y от величины x
называют такую зависимость, при которой каждому значению величины x из множества ее возможных значений, соответствует некоторое множество возможных значений величины y , характеризуемое определенным
законом распределения. При этом изменение величины x приводит к изменению распределения величины y . Например, зависимость между
возрастом ребенка и его весом является статистической, т.к. если у группы детей в возрасте 2-х лет (x = 2)измерить массу тела y , то получим множество различных значений: y1 =12,7 кг, у2 =13,1кг, у3 =13,3 кг
и т. д.
Частным случаем статистической зависимости является кореляционная зависимость.
Корреляционная зависимость проявляется в том, что при изменении одной величины, например x , изменяется условное среднее значение
50
величины y . Например, зависимость между дозой лекарственного препа-
рата и содержанием его в крови, зависимость между возрастом ребенка и массой его тела и т.д.
Корреляционную зависимость y от x можно описать с помощью
уравнения вида
yx = f (x),
где yx - условное среднее значение величины y , соответствующее
данному значениюx ; f (x) - некоторая функция. Уравнение (1) называется уравнением регрессии y на x .
Обратную корреляционную зависимость (если она существует) можно описать уравнением регрессии x на y :
xy = ϕ( y) ,
где |
xy − условное среднее значение |
величины x , соответствующее |
данному значению y ; ϕ(y)- некоторая функция. |
||
x |
Функции f (x) и ϕ(y) называют соответственно регрессиями y на |
|
и x на y , а их графики – линиями регрессии. |
||
|
Целью корреляционного анализа является установление зависимости |
|
(связи) между случайными величинами |
x и y , и если эта зависимость |
|
имеется, то определить числовые характеристики для математического описания этой зависимости, т.е. определить вид функции f (x) или ϕ(y).
Качественная оценка корреляционной зависимости
Пусть изучается корреляционная зависимость между случайными величинами x и y . В результате n независимых опытов получены n
пар значений величин x и y : (x1; y1 ), (x2 , y2 ),…, (xn , yn ), которые представляют выборку из генеральной совокупности. Экспериментальные значения записывают в таблицу.
xi |
x1 |
x |
… |
x |
… |
x |
|
|
|
|
2 |
|
i |
|
|
|
|
|
|
|
|
|
n |
yi |
y1 |
y |
2 |
… |
y |
… |
yn |
|
|
|
|
i |
|
||
Для качественной оценки корреляционной зависимости между величинами x и y строят корреляционное поле: в прямоугольной системе
координат в выбранном масштабе наносят точки с координатами (x1; y1 ),
51
(х2 ; y2 ) ,…, (хn ; yn ) . В результате построения корреляционное поле может иметь один из видов, изображенных на рис.4.
Рис. 4.а |
Рис. 4.б |
Если корреляционное поле имеет вид, изображенный на рис. 4.а, то |
|
можно предположить, что между величинами |
x и y существует прямая |
(или возрастающая) корреляционная зависимость. В случае прямой зависимости обе величины x и y увеличиваются. Например, корреляционная
зависимость между возрастом ребенка и массой его тела является прямой, т.к. с увеличением возраста увеличивается масса тела.
Если корреляционное поле имеет вид, изображенный на рис. 4.б, то можно предположить, что между величинами x и y существует обратная
(или убывающая) корреляционная зависимость. Обратная корреляционная зависимость означает, что одна величина (например y ) уменьшается по
мере увеличения другой ( x ). Например, корреляционная зависимость между процентом охвата населения прививками против гриппа и заболеваемостью гриппом является обратной, т.к. чем больше привито населения против гриппа, тем меньше заболеваемость.
По характеру расположения экспериментальных пар значений (хi ; yi ) на корреляционном поле можно сделать предположение не только
о наличии или отсутствии корреляционной зависимости между величинами x и y , но и о форме этой зависимости. Если экспериментальные точ-
ки расположены приблизительно вдоль прямой линии (рис. 5.а), то можно высказать предположение о существовании линейной корреляционной зависимости между изучаемыми величинами, если точки расположены вдоль ветви параболы (рис. 5.б), - о квадратичной зависимости и т.д.
52
Рис. 5.а |
Рис. 5.б |
Количественная оценка тесноты корреляционной зависимости. Выборочный коэффициент корреляции.
Для оценки тесноты линейной корреляционной зависимости определяют коэффициент корреляции ( r ), который обладает следующими основными свойствами:
1.Коэффициент корреляции двух независимых случайных вели-чин равен нулю.
2.Коэффициент корреляции двух величин, связанных линейной функциональной зависимостью, равен 1 в случае возрастающей зависимости и
−1 в случае убывающей.
3.Абсолютная величина коэффициента корреляции не превышает 1, т.
е. r ≤1 или −1 ≤ r ≤1.
По величине выборочного коэффициента корреляции, используя данные таблицы 1, можно оценить тесноту связи.
|
|
Таблица 1 |
Вид связи |
Характер связи и величина r |
|
Прямая (r > 0) |
Обратная (r < 0 ) |
|
Функциональная |
r = +1 |
r = - 1 |
Сильная |
r [+ 0,7 …+ 1) |
r [- 0,7 …- 1) |
Средняя |
r [+ 0,3 …+ 0,7) |
r [- 0,3 …- 0,7) |
Слабая |
r ( 0,0 …+ 0,3) |
r ( 0,0 …- 0,3) |
Связь |
r = 0 |
r = 0 |
отсутствует |
|
|
53
Выборочные уравнения регрессии. Выборочные коэффициенты регрессии. Метод наименьших квадратов.
Рассмотрим те случаи, когда есть основания предполагать наличие линейной корреляционной зависимости между величинами x и y . В этих
случаях, для описания корреляционной зависимости между величинами x и y вводят выборочные уравнения линейной регрессии y на x и x на y :
yx = ρyx x +b; xy = ρxy y + d ,
где yx - ожидаемое среднее значение величины y , соответствующее данному значению x ; xy−ожидаемое среднее значение величины x , соот-
ветствующее данному значению y ; ρyx , ρxy - выборочные коэффици-
енты регрессии, b и d отрезки отсекаемые линиями регрессии на осях ОY и ОХ соответственно.
Для вычисления выборочного коэффициента корреляции r применяется метод наименьших квадратов, подробно изложенный в учебниках по статистике. В случае линейной корреляционной зависимости y на x и
x на y выборочные коэффициенты регрессии ρyx и ρxy всегда раз-
личны ( ρyx ≠ ρxy ). Вычислив ρyx , ρxy и b , d , можно записать уравнения выборочной регрессии и построить линии выборочной регрессии.
Линии выборочной регрессии y на x и x на y не совпадают,
но обязательно проходят через точку с координатами . (х; y)
На рис. 6 показаны линия регрессии y на x (прямая I) и
линия регрессии x на y (прямая
II). При этом
tgα = ρyx , а tgβ = ρxy
Рис. 6
54
Пример решения задачи.
Задача. Имеются условные данные по районам охвата населения прививкамиxi , % и заболеваемости брюшным тифом yi , %:
Таблица 1
Районы |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
xi |
14,5 |
13,0 |
11,2 |
10,0 |
8,3 |
7,9 |
6,0 |
5,5 |
3,0 |
yi |
3,0 |
3,2 |
4,0 |
4,5 |
5,0 |
5,5 |
6,0 |
7,0 |
8,5 |
1.Проведите качественную оценку корреляционной связи.
2.Определите коэффициент корреляции, характер и размер корреляционной связи.
3.Проверьте достоверность коэффициента корреляции.
4.Найдите коэффициент регрессии заболеваемости брюшным тифом от охвата населения прививками. Объясните смысл.
5.Составьте уравнение регрессии.
6.На корреляционном поле начертите линию регрессии.
7.Определите ожидаемые заболеваемости брюшным тифом при 2 % и 15 % охвате населения прививками.
8.Сделайте выводы.
Решение.
1. По данным табл. 1, по геометрическому положению точек (корреляциионному полю) на рисунке 7 устанавливаем линейная зависимость или нет. Для этого в прямоугольной системе координат x , y откладываем экс-
периментальные точки. По рисунку 7, т.к. корреляционное поле имеет вид эллипса, по хара-ктеру расположения точек делаем вывод о линейной зависимости между x и y .
Примечание: Если корреляционное поле не показывает линейную зависимость, считать зависимость между параметрами x и y линейной.
Таким методом проводим качественную оценку корреляционной связи: прямая (положительная), обратная (отрицательная) или отсутствует линейная связь.
В данной задаче наблюдается линейная обратная зависимость между охватом населения прививками и заболеваемостью брюшным тифом, причем теснота корреляционной связи сильная, т.к. корреляционное поле – вытянутый (узкий) эллипс, диагональ которого расположена под углом 450
коси ОХ.
2.Проводим количественную оценку корреляционной связи. Установив линейную зависимость междуx и y и, учтя, что по условию задачи
объём выборки n 30, коэффициент линейной корреляции рассчитываем по формуле (МНК (см. учебную литературу)):
55
|
|
n |
− x)( yi − y) |
|||
|
|
∑(x2 |
||||
r = |
i=1 |
|
|
|
(1). |
|
n |
|
|
|
|||
|
|
n |
||||
|
∑(xi − x)2 |
∑( yi − y)2 |
||||
|
i=1 |
|
|
i=1 |
|
|
|
|
|
|
|
|
|
Рис. 7
Для определения коэффициента корреляции r по формуле (1) составляем расчетную таблицу 2. Данные заносим в столбцы 1,2,3:
|
|
|
|
|
|
|
|
|
|
Таблица 2 |
Рас- |
x |
y |
i |
x − x |
y − y |
(x − x) |
2 |
( y − y) |
2 |
(x − x)( y − y) |
чет |
i |
|
i |
i |
|
|
|
|||
|
|
|
|
|
i |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
|
7 |
|
8 |
|
1 |
14,5 |
3,0 |
5,68 |
-2,19 |
32,26 |
|
4,80 |
|
-12,44 |
|
2 |
13,0 |
3,2 |
4,18 |
-1,99 |
17,47 |
|
3,96 |
|
-8,32 |
|
3 |
11,2 |
4,0 |
2,38 |
-1,19 |
5,66 |
|
1,42 |
|
-2,83 |
|
4 |
10,0 |
4,5 |
1,18 |
-0,69 |
1,39 |
|
0,48 |
|
-0,81 |
|
5 |
8,3 |
5,0 |
-0,52 |
-0,19 |
0,27 |
|
0,04 |
|
+0,10 |
|
6 |
7,9 |
5,5 |
-0,92 |
+0,31 |
0,85 |
|
0,10 |
|
-0,29 |
|
|
|
|
|
|
|
56 |
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
7 |
6,0 |
6,0 |
-2,82 |
+0,81 |
7,95 |
0,66 |
-2,28 |
8 |
5,5 |
7,0 |
-3,32 |
+1,81 |
11,02 |
3,28 |
-6,01 |
9 |
3,0 |
8,5 |
-5,82 |
+3,31 |
33,87 |
10,96 |
-19,26 |
n=9 ∑=79,4 ∑=46,7 |
|
|
|
∑= 110,74 |
∑=25,7 |
∑=-52,14 |
||||||||||
Средние |
арифметические |
значения |
|
вариант |
получим как |
|||||||||||
|
n |
|
9 |
|
|
|
|
|
|
n |
|
|
|
|
||
|
∑xi |
|
∑xi |
|
79,4 |
|
|
|
∑ yi |
|
46,70 |
|
|
|||
|
|
|
|
|
|
|
|
|||||||||
x = |
i=1 |
= − |
i=1 |
|
= |
|
=8,82, |
у= |
i=1 |
|
= |
|
= 5,19. |
|||
n |
9 |
9 |
n |
9 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||
3. Заполним столбцы 4, 5, 6, 7, 8, производя соответствующие промежуточные расчеты.
Находим суммы столбцов 6,7,8.
4. Расчет коэффициента линейной корреляции по формуле 1:
r = |
−52,14 |
= − |
52,14 |
|
= −0,98 . |
|
110,74 25,7 |
53,35 |
|||||
|
|
|
||||
Примечание: Расчеты по формуле должны производиться с точностью не менее чем до сотых.
5.Вывод. На основании схемы таблицы оценки тесноты корреляционной связи вычисленный коэффициент корреляции r свидетельствует о том, что имеется сильная обратная связь между заболеваемостью брюшным тифом и охватом населения прививками, т.е. чем больше охват населения прививками, тем меньше заболеваемость населения брюшным тифом.
6. Проверка достоверности корреляционной связи.
Для статистических характеристик устанавливается их достоверность, для чего используются один из ниже приведенных способов.
Первый способ.
Для того чтобы убедиться, что коэффициент корреляции, вычисленный по данным выборочного исследования, будет соответствовать размеру связи в генеральной совокупности, необходимо:
1) |
выдвинуть нулевую гипотезу H0 : rB = rГС = 0, |
т.е. полученное |
||
значение выборки случайно и обусловлено ошибками измерений; |
||||
2) |
вычислить параметр tф: tф = |
r |
, где mr - |
средняя ошибка |
|
||||
|
|
mr |
|
|
коэффициента корреляции, вычисляется по формулам:
для n ≤100 |
m |
= ± |
1 |
− r2 |
; |
|
r |
|
n − 2 |
|
|
|
|
|
|
||
57
для n 100 m |
|
= ± |
1 − r2 |
|
. |
||
r |
|
|
n |
|
|||
Для данной задач |
|
||||||
|
|
|
|
||||
m = ± |
1 − |
0,982 |
= ± |
0,006 = ±0,008, |
|||
r |
|
9 |
− 2 |
|
|
|
|
|
|
|
|
|
|||
Тогда tф = |
−0,98 |
= −12,25; |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|||||||||
|
0,08 |
|
|
|
|
|
|
|
|
|
|
|
||
3) |
найти критерий Стьюдента по таблице 2 Приложения выбороч- |
|||||||||||||
ному |
уровню |
значимости |
α = 0,05 и |
|
числу степеней |
свободы |
||||||||
k = n − 2 = 7 |
|
|
|
|
|
|
|
|
|
|
|
|
||
tкр = tα,k = t0,05;7 |
= 2,36. |
|
|
|
|
|
||||||||
4) |
сравнить |
tкр |
|
tф |
|
|
т.к. tкр = 2,36 |
|
tф |
|
=12,25, то |
нулевая |
||
|
|
: |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
гипотеза отвергается, следовательно, коэффициент, вычисленный по данным выборочного исследования, соответствует размеру связи в генеральной совокупности.
Второй способ.
Существует специальная таблица (см. таблицу 1 Приложения) для малых выборок, по которой полученную величину r можно оценить без
расчета mr и tф,а путем сравнения r с критическим коэффициентом корреляции rα,k , рассчитанным в таблице для различной степени
вероятности (или уровня значимости) и различного числа наблюдений.
Здесь k - число коррелируемых пар.
По этой таблице коэффициент корреляции признаётся статистически значимым (достоверным) с вероятностью ошибки меньше 0,05 ,если
r > r0,05,k , и с вероятностью ошибки меньше 0,01, если r > r0,01,k .
В данной задаче по таблице 1 приложения: r0,05;9 = 0,7.Так как rф = −0,98 > r0,05 = 0,666 коэффициент r = −0,98 соответствует
связи в генеральной совокупности.
Примечание. На практике обязательно подвергается проверке на
достоверность коэффициенты корреляции меньше 0,5 ( |
|
r |
|
< 0,5). |
|
||||
|
|
|
|||||||
7)Находим коэффициент регрессии по формуле |
ρyx |
= r |
σy |
||||||
σ |
x |
||||||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||
58
