Статистический анализ в физиологии
.pdf41
К ритерий К олм огорова-См ирн ова осн овы вается н а статистической м одели, котораяпредполагаетн епреры вн остьраспределен ия, так чтовероятн остьсовпаден иявы борочн ы х зн ачен ий равн а н улю . О дн акон а практикекритерий частоприм ен яетсяк сгруппирован ны м дан н ы м идан н ы м вы - борок из дискретн ы х распределен ий. В обоих этих случаях учиты вается возм ож н ость появлен ия равн ы х зн ачен ий н аблю ден ий. В дан н ом случае уровен ьзн ачим ости критериян иж ен ом ин альн ого, и вероятн остьош ибки второгорода возрастает.
К ак критерий χ2, так и критерий К олм огорова-См ирн ова предполагаю т, чтораспределение, фигурирую щ еевн улевой гипотезе, долж н обы ть полн остью определен озаран ее(н априм ер, он ом ож етбы тьн орм альн ы м с н улевы м средн им и единичн ой дисперсией). П ри работес критерием согласияχ2 н еобх одим озаран ееопределитьтипраспределен ия, а парам етры оцениваю тсяповы борочн ы м дан н ы м . К ритерий χ2 легком одифицируется при пом ощ и ум ен ьш ен ия числа степен ей свободы , н о н еизвестн о, какие изм ен ен иядолж н ы бы тьвн есен ы впроцедуру прим енен иякритерияК ол- м огорова-См ирн ова. Н есм отря н а это он ин огда прим ен яется вкачестве критерия при проверке гипотез о закон е распределен ия. Следует учиты - вать, что вэтом случае истин н ы й уровен ь значим ости будет н есколько ниж ен ом инальн огоивозрастетвероятн остьош ибок второгорода. П овсей видим ости, этотэффектн ебудетвелик, есличислооцен иваем ы х парам етровм алопосравн ен ию собъем ом вы борки.
Дан н ы е: вы борка x1, x2, x3, … , xn из n н езависим ы х н аблю ден ий, упорядоченн ы х так, чтоx1< x2< x3< … < xn.
Ст ат ист ическ аям одель: н аблю ден иян езависим ы иберутсяиз генеральн ой совокупн ости, распределение которой предполагается н епреры вн ы м .
Г ипот езы : н улеваягипотеза –фун кцияраспределен ияравн а F(x).
Вы числен ие к р ит ер иальн ой ст ат ист ик и:
1.В ы числяю тсякум улятивн ы еразн ости: D1=1-nF(x1), D2=1-nF(x2) и т.д.
2.Н ах одим |Di|max , н аибольш ееабсолю тн оезн ачен иекум улятивн ы х разн остей.
3.В ы числяем зн ачениекритерияD=(|Di|max)/n
|
П р им |
ер . Д ан ы 20 н аблю ден ий (табл.9) из н еизвестн ой генеральн ой |
|||||
совокупн ости. |
|
|
Т аблица 9 |
||||
|
|
|
|
|
|
||
|
0,33 |
|
-0,52 |
-2,41 |
-1,93 |
0,46 |
|
|
-0,44 |
|
-0,97 |
-0,38 |
0,48 |
1,29 |
|
|
-1,82 |
|
-1,23 |
-0,21 |
2,66 |
-1,22 |
|
|
-0,41 |
|
-0,95 |
1,47 |
-0,83 |
-0,43 |
|
Н еобх одим опроверитьн улевую гипотезу, |
предполагаю щ ую , чтоэто |
норм альн ое распределен ие с н улевы м средн им |
зн ачением и един ичн ой |
42
дисперсией.
Д ля этого н адо первон ачальн о упорядочить наблю ден ия по возрастаю щ ей, а затем –вы числитьн аблю даем ы еиож идаем ы езн ачен ия, которы е м ен ьш екаж догонаблю ден ияилиравн ы ем у. Н еобх одим отакж евы числить абсолю тн оезн ачен иекум улятивной разн ицы вовсех н аблю дениях .
П р им ер . Зафиксирован н ое число н аблю ден ий, которы е н е больш е 0,97 равн о 6. О ж идаем оезначен иеравн о: 20F (-0,97) = 3,32. А бсолю тн ое зн ачен иекум улятивн ой разн остиприx=-0,97 равн о6-3,32=2,68.
П о дан н ы м составлен н ой таблицы |
10 н ах одим м аксим альн оезн аче- |
||
ние|Di| –5,96. Зн ачен иекритериальн ой статистикиравн о: |
|||
965, |
|
||
D = |
|
= |
298 0, |
20 |
П ри пом ощ и таблиц, содерж ащ их сведения о табличны х значениях распределения, известн о, чтокритическоезначениеp=0,05 равн о 0,294. Т а- ким образом , необх одим о отклон ить н улевую гипотезу о том , что дан ное распределениеявляетсян орм альн ы м с нулевы м средним и един ичной дисперсией.
|
|
|
Т аблица 10 |
Н аблю ден ие |
Ч ислон аблю ден ий, м еньш их илиравн ы х xi |
Разн ость |
|
xi |
Н аблю даем ое |
О ж идаем ое |
|Di| |
-2,41 |
1 |
0,16 |
0,84 |
-1,93 |
2 |
0,54 |
1,46 |
-1,82 |
3 |
0,69 |
2,31 |
-1,23 |
4 |
2,19 |
1,81 |
-1,22 |
5 |
2,22 |
2,78 |
-0,97 |
6 |
3,32 |
2,68 |
-0,95 |
7 |
3,42 |
3,58 |
-0,83 |
8 |
4,07 |
3,93 |
-0,52 |
9 |
6,03 |
2,97 |
-0,44 |
10 |
6,60 |
3,40 |
-0,43 |
11 |
6,67 |
4,33 |
-0,41 |
12 |
6,82 |
5,18 |
-0,38 |
13 |
7,04 |
5,96 |
-0,21 |
14 |
8,34 |
5,66 |
0,33 |
15 |
12,59 |
2,41 |
0,46 |
16 |
13,54 |
2,46 |
0,48 |
17 |
13,69 |
3,31 |
1,29 |
18 |
18,03 |
3,03 |
1,47 |
19 |
18,58 |
0,42 |
2,66 |
20 |
19,92 |
0,08 |
43
|
П |
Р О В Е Р КА С ТА ТИ С ТИ Ч Е С КИ Х Г И П О ТЕ З В П А КЕ ТЕ STADIA |
|
С |
И С П О |
Л ЬЗО В А Н И Е М М Е ТО Д О В Н Е П |
А Р А М Е ТР И Ч Е С КО Й С ТА ТИ С ТИ КИ |
|
В |
качествеприм ера используем |
дан н ы е о содерж ан ии х олестерин а |
(м М |
оль/л) всы воротке крови здоровы х м уж чин различны х возрастн ы х |
групп: 1 группа – лица до 19 лет, вторая группа – лица старш е 30 лет
(табл.11).
Т аблица 11
Г ру ппа1
5,50 |
5,53 |
5,51 |
5,46 |
5,53 |
|
5,86 |
|
5,73 |
5,55 |
5,49 |
5,48 |
5,44 |
5,49 |
5,71 |
4,88 |
5,74 |
5,44 |
|
5,48 |
|
5,41 |
5,77 |
5,79 |
5,53 |
|
|
|
|
|
|
Г ру ппа2 |
|
|
|
|
|
||
5,84 |
6,20 |
5,93 |
5,99 |
5,97 |
|
6,31 |
|
6,22 |
6,24 |
6,33 |
6,11 |
6,12 |
6,18 |
6,05 |
5,88 |
5,97 |
5,94 |
|
6,13 |
|
6,21 |
6,18 |
6,22 |
6,24 |
|
П риводим ы й н иж ерисун ок 9 иллю стрируетвведен н ы едан н ы евокнередактора дан н ы х пакета STADIA.
Ф а йл: |
|
|
П е р е м е н н ы х =2 |
|
|
И зм е р е н ий=42 |
|
Var/Cases |
1/21 |
2/21 |
3/0 |
4/0 |
5/0 |
6/0 |
7/0 |
Varname |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
1 |
5.5 |
5.84 |
|
|
|
|
|
2 |
5.53 |
6.2 |
|
|
|
|
|
3 |
5.51 |
5.93 |
|
|
|
|
|
4 |
5.46 |
5.99 |
|
|
|
|
|
5 |
5.55 |
5.97 |
|
|
|
|
|
6 |
5.86 |
6.31 |
|
|
|
|
|
7 |
5.75 |
6.77 |
|
|
|
|
|
8 |
5.55 |
6.24 |
|
|
|
|
|
9 |
5.49 |
6.39 |
|
|
|
|
|
10 |
5.48 |
6.11 |
|
|
|
|
|
11 |
5.44 |
6.17 |
|
|
|
|
|
12 |
5.49 |
6.18 |
|
|
|
|
|
13 |
5.71 |
6.05 |
|
|
|
|
|
14 |
4.88 |
5.88 |
|
|
|
|
|
15 |
5.74 |
5.97 |
|
|
|
|
|
16 |
5.44 |
5.94 |
|
|
|
|
|
17 |
5.48 |
6.13 |
|
|
|
|
|
18 |
5.41 |
6.21 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Б ЛО К РЕ ДАКТО РА Д АННЫ Х |
|
|
||
F1 П ом ощь |
F2 П е ч Э кр |
F3 Ч т е ние F4 |
Запись F5 А р х ив F6 |
Рисуно к F7 |
Оч ист F8 |
П р е о бр F9 Ст ат ис F0 |
Вы х од |
Вво дите в м а тр ицучисла + Enter (р а бо та ю т та кже : Enter/ Ins/ Del/ Tab и F-клю чи |
|
||||||
|
|
|
|
|
|
|
|
Рис.9. В веденн ы едан ны епосодерж ан ию х олестерин а всы вороткекровиу двух группм уж чин
А н алогичн оан ализу вы борки с использован ием м етодовописательной статистики для вы х ода вм ен ю осн овн ы х м етодован ализа н аж м ем клавиш у F9 ивы берем пун кт« 5 = х и-квадрат» . П ривы борепункта « 0 эм - пирич. (одн ородн ость)» н а экран епоявляю тсярезультаты ан ализа (рис.10), свидетельствую щ иеоботсутствииразличий враспределен иивы борочн ы х величин .
44
П ри вы борепун кта « 1 = теоретич. (согласия)» производитсяан ализ согласиявы борочн огои предполагаем огостатистическогораспределен ия. Результаты ан ализа приводятсян а рисун ке11.
Ф а йл: |
П е р е м е н н ы |
х =2 |
И зм е р е н ий=42 |
КРИТЕ РИ Й ХИ -КВАД РАТ. |
|
П е р е м е н н ы |
е : х 1 х 2 |
Хи-ква др а т=6.739 Е -2 |
Зн а чим о сть =1 |
сте п.сво б=20 |
|
Гипо те за 0: Не т р а зличий м е ждудвум я р а спр е де ле н ия м и |
|
||
|
СТАТИ СТИ ЧЕ СКИ Е М Е ТО ДЫ |
||
F10 Вы х о д F1 И н ф о р м а ция |
F2 Экр а н н а пе ча ть / в ф а йл |
Esc Вы х о д спр е р ы ва н ие м |
|
На жм ите Enter=пр о до лжить |
или F2=пе ча ть экр а н а >> |
|
Рис.10. Результаты анализаразличий враспределениивы борочны х величин привы борепун кта« 0 эм пирич. (однородн ость)»
Ф а йл: |
|
П е р е м |
е н н ы |
х =2 |
И зм е р е н ий=42 |
КРИТЕ РИ Й ХИ -КВАД РАТ. |
|
П е р е м е н н ы |
е : х 1 х 2 |
||
Хи-ква др а т=1.227 |
Зн а чим о сть =0.9999 |
сте п.сво б=20 |
|
||
Гипо те за 0: Не т р а зличий м е ждудвум я р а спр е де ле н ия м и |
|
||||
|
|
СТАТИ СТИ ЧЕ СКИ Е М Е ТО ДЫ |
|||
F10 Вы х о д |
F1 И н ф |
о р м а ция F2 Экр а н |
н а пе ча ть / в ф а йл |
Esc Вы х о д спр е р ы ва н ие м |
|
На жм ите Enter=пр о до лжить или F2=пе ча ть |
экр а н а >> |
|
Рис.11. Результаты анализаразличий враспределениивы борочны х величин привы борепун кта « 1 = теоретич. (согласия)»
Д алее проан ализируем используем ы е дан н ы е с прим енен ием критерияВ илкоксон а. Д ляэтоговы берем пун кт« 6 = сдвига (полож ен ия)» разделан епарам етрических тестов. Результаты ан ализа приведен ы н а рисун ке12.
Ф а йл: |
|
|
П е р е м |
е н н ы х =2 |
|
И зм е р е н ий=42 |
КРИТЕ РИ Й СДВИ ГА (П О ЛО ЖЕ НИ Я) |
П е р е м е н н ы |
е : х 1 |
х 2 |
|||
Вилко ксо н =232 |
Z=5.522 |
Зн а чим о сть =0, |
сте п.сво б=21, |
21 |
|
|
Гипо те за 1: Е сть |
р а зличия |
м е ждум е диа н а м и вы бо р о к |
|
|
||
|
|
|
СТАТИ СТИ ЧЕ СКИ Е М Е ТО ДЫ |
|
||
F10 Вы х о д F1 И н ф о р м а ция |
F2 Экр а н н а пе ча ть / в ф а йл |
Esc Вы х о д спр е р ы ва н ие м |
||||
На жм ите Enter=пр о до лжить |
>> |
|
|
|
Рис.12. Результаты ан ализа сприм ен ен ием критерияВ илкоксон а
45
П одан н ы м ан ализа сприм ен ен ием критерияВ илкоксон а м ы долж н ы отвергн утьгипотезу оравен ствем едиан вы борок.
Д ля анализа дан н ы х прим ера с пом ощ ью критерия К олм огороваСм ирн ова длядвух вы борок н еобх одим овглавн ом м ен ю м етодован ализа избратьпун кт« 8 = произвольн ы х альтерн атив». Н а рисун ке13 приведен ы результаты прим ен ен иядан н огокритерияк ан ализируем ы м вы боркам .
Ф а йл: |
|
П е р е м е н н ы х =2 |
|
|
И зм е р е н ий=42 |
|
КРИТЕ РИ Й КО Л М О ГО РО ВА -СМ И РНО ВА |
|
П е р е м е н н ы |
е : х 1 |
х 2 |
||
Ко лм о го р о в-См ир н о в=0.9524, |
Зн а чим о сть =1.068Е -8, сте п.сво б=21, |
21 |
||||
Гипо те за 1: Е сть |
ин те гр а ль н ы |
е р а зличия |
м е ждувы |
бо р ка м и |
|
|
|
|
СТАТИ СТИ ЧЕ СКИ Е М Е ТО ДЫ |
|
|
||
F10 Вы х о д F1 |
И н ф о р м а ция |
F2 Экр а н |
н а пе ча ть / в ф а йл |
Esc Вы х о д спр е р ы ва н ие м |
||
На жм ите Enter=пр о до лжить или F2=пе ча ть |
н а экр а н |
>> |
|
|
Рис.13. Результаты анализа спом ощ ью критерияК олм огорова-См ирн ова
Результаты ан ализа свидетельствую т о н аличии ин тегральн ы х различий м еж дудвум яан ализируем ы м ивы боркам и.
П рим енениен епарам етрических тестовдляпарн ы х вы борок –пун кт м ен ю « 9 = дляпарн ы х вы борок» позволяетпроан ализироватьвы борочн ы е зн ачен иядлявы борок, получен н ы х припарн ом вы борочн ом исследован ии
с использован ием |
как критерия В илкоксон а, |
так и критерия зн аков. Ре- |
|||||||
зультаты ан ализа приведен ы н а рисун ке14. |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
Ф а йл: |
|
|
|
|
П |
е р е м |
е н н ы х =2 |
И зм е р е н ий=42 |
|
П АРНЫ Е |
КРИ ТЕ РИ И СДВИ ГА |
|
П е р е м е н н ы |
е : х 1 х 2 |
||||
|
Вилко ксо н =0 |
Z=-4.015, |
Зн а чим о сть =0, |
сте п.сво б=21 |
|
||||
|
Зн а ко в=0 |
Z=-4.583, |
Зн а чим о сть =0, |
сте п.сво б=21 |
|
||||
|
Гипо те за 1: Е сть |
р а зличия |
м е ждум е диа н а м и вы бо р о к |
|
|||||
|
|
|
|
|
|
СТАТИ СТИ ЧЕ СКИ Е М Е ТО ДЫ |
|||
|
F10 Вы х о д |
F1 |
И н ф |
о р м а ция F2 Экр а н |
н а пе ча ть /в ф а йл Esc Вы х о д спр е р ы ва н ие м |
||||
|
На жм ите Enter=пр о до лжить или F2=пе ча ть |
н а экр а н |
|
Рис.14. Результаты ан ализа припарн ом вы борочн ом исследован иисиспользован ием критерияВ илкоксон а икритериязн аков.
Д анн ы екритериитакж еподтверж даю тналичиеразличий м еж ду м е- дианам ивы борок.
В заклю чен иедан н огораздела необх одим оповторитьприводивш еесяран еепредупреж ден иеоприм ен ен иик одн им дан н ы м н ескольких м ето-
46
дован ализа. Е сли к одн им и тем ж е данн ы м прим ен яю т два различн ы х критериядляпроверки одн ой и той ж ен улевой гипотезы (или двух сх одны х гипотез) ивкаж дом случаеприн им аетсяуровен ьзн ачим ости, равны й, наприм ер, 5%, то вероятн остьтого, что х отябы по одн ом у из критериев н улевая гипотеза будет ош ибочн о отклон ен а, превосх одит 5 %. Следует воспользоватьсялиш ьодн им критерием , ж елательн оболеем ощ н ы м .
|
КО Р Р Е Л Я Ц И О Н Н Ы Й А Н А Л И З |
М |
еж ду явлен иям и сущ ествую т две категории причин н о-следствен - |
ны х связей: ф ункциона льные икорре ля ционные . |
|
Ф |
ункциона льна я свя зь им еетотнош ениек каж дом у отдельном у наблю - |
дению , являетсяобязательной и зависитотфизических и х им ических законов. Н априм ер, при повы ш ении давлениянасы щ енн огопара увеличивается еготем пература; нагреван ием еталловприводитк их расш ирению ит. п.
Связьм еж ду явлениям и, отм ечаю щ аясян евкаж дом отдельн ом случае, а прим ассовом сопоставлен иирассм атриваем ы х признаков, назы ваетсякорре ля ционной. В этом случаекакой-либопризн ак (уровен ьзаболевае- м ости остры м и киш ечн ы м и ин фекциям и) м ож етбы тьсвязан главн ы м образом (н он еобязательн о) сизм ен ен иям и другогопризнака (качествореализуем ы х м олочн ы х продуктов), х отян а интен сивн остьтечен ияэпидем и- ческогопроцесса при остры х киш ечн ы х инфекциях влияю ти другиефакторы : м икробн ая обсем ен ен н ость потребляем ой н аселен ием воды , кулинарны х изделий, м играцион н ы епроцессы идр. К орреляцион н аясвязьпроявляется всредн ем для всей совокупн ости н аблю ден ий, а вотн ош ен ии единичн ы х случаевона неполн а ин еточн а.
К орреляцион наязависим остьм еж дуявлен ием м ож етбы тьпрям ой, когда сповы ш ением одн огопоказателяотм ечаетсяростдругого, и обратной, когда с увеличением одного признака другой ум еньш ается (наприм ер, чем больш енаселен ияох ваченопрививкам и, тем ниж еуровеньзаболеваем ости).
П реж дечем изм ерятьвеличин у коэффициен та корреляции, н еобх о- дим о реш итьвопрос о возм ож н ости причин н о-следствен н ой связи м еж ду изучаем ы м и явлен иям и. П араллельноеизм ен ен иестатистических показателей сам о по себеещ ен еговорито н аличии связи, так как м ож етбы ть обусловлен ослучайн ы м совпаден ием м н огих обстоятельств, н есвязанн ы х другсдругом .
Ц ифровы едан н ы е, подвергаю щ иесякорреляцион н ом у ан ализу, дол- ж н ы бы ть сгруппирован ы с учетом особен н остей изучаем ы х явлений. В противн ом случае зн ачен ие получен н ого коэффициен та будет заведом о ош ибочн ы м . Э тоособен н окасаетсяобластиин фекцион н ой патологии, когда отм ом ен та ин фицирован ия, т. е. действияпричин ы , доклин ического проявлен ия процесса, т. е. следствия, прох одит инкубацион н ы й период, длительн остькоторогоразличн а при отдельн ы х заболеван иях . П родолж и- тельн остьин кубацион н ого периода и является оптим альн ы м сроком вре- м ен н ой группировкицифровы х дан н ы х . П рин еин фекцион н ой патологиив
47
осн ову группировки кладут другие врем ен н ы е признаки, ин тересую щ ие исследователя. Н априм ер, при изучен ии частоты развитияин фаркта м иокарда у лиц одн огои тогож евозраста взависим остиотдлительн ости курен иям атериалм ож н огруппироватьсучетом ин тересую щ ей давн остикурен ия: 1, 2, 3 ит. д. лет.
Ме т од ра нговой корре ля ции н аиболеепростой, требую щ ий н езн ачительн ы х затратврем ен и н а проведениевы числен ий. В тож еврем яон м е- нееточен , даетприблизительн оепредставлен иеох арактереитесн отесвязим еж ду явлен иям и.
П р им ер . О пределитьвозм ож н ую связьм еж ду уровнем заболеваем ости остры м икиш ечны м иинфекциям и(О К И ) икачеством реализуем ы х кулинарны х изделий, определяем огопом икробн ой обсем ененн ости(табл.12).
|
|
|
|
|
|
|
|
|
|
|
Т аблица 12 |
|||
|
|
П роцен тн е- |
П орядковы ен ом ера |
|
|
|
|
|
|
|
|
|
|
|
|
Заболева |
стан дартн ы х |
|
|
|
|
|
|
|
d 2 сн а- |
||||
|
(ран ги) |
|
|
|
|
|
|
|
||||||
Н еделя |
ем ость |
пробкулин ар- |
|
|
|
|
|
|
|
растаю - |
||||
|
|
|
|
|
|
|
|
|
||||||
|
|
d |
d 2 |
|
||||||||||
года |
О К И |
н ы х изделий по |
|
|
|
|||||||||
Заболе- |
% н естан - |
|
|
|
|
|
|
|
щ им ито- |
|||||
|
|
м икробн ой об- |
дартн ой |
|
|
|
|
|
|
|
|
гом |
||
|
|
ваем ости |
|
|
|
|
|
|
|
|
||||
|
|
сем ен ен н ости |
продукции |
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
|
7 |
|
|
8 |
||||
1 |
12,8 |
8,3 |
6 |
2 |
4 |
16 |
|
16 |
||||||
2 |
11,6 |
5,1 |
2 |
1 |
1 |
1 |
|
|
17 |
|||||
3 |
15,4 |
18,0 |
9 |
10 |
1 |
1 |
|
|
18 |
|||||
4 |
17,8 |
14,0 |
11 |
7 |
4 |
16 |
|
34 |
||||||
5 |
12,3 |
11,6 |
4 |
4,5 |
0,5 |
0,25 |
|
34,25 |
||||||
6 |
12,5 |
14,2 |
5 |
8 |
3 |
9 |
|
|
43,25 |
|||||
7 |
14,0 |
24,4 |
7 |
11 |
4 |
16 |
|
59,25 |
||||||
8 |
20,3 |
25,4 |
12 |
12 |
0 |
0 |
|
|
59,25 |
|||||
9 |
15,6 |
11,6 |
10 |
4,5 |
5,5 |
30,25 |
|
89,50 |
||||||
10 |
14,8 |
12,9 |
8 |
6 |
2 |
4 |
|
|
93,50 |
|||||
11 |
11,6 |
14,6 |
3 |
9 |
6 |
36 |
|
139,50 |
||||||
12 |
9,3 |
10,0 |
1 |
3 |
2 |
4 |
|
|
133,50 |
1.О пределяю тсяпорядковы ен ом ера (ран ги) показателей заболевае-
мости (графа 4) и процен та н естандартн ой продукции (графа 5). Ран г 1 присваиваетсян аим ен ьш им показателям . П рин аличиин ескольких , равн ы х
по величин е показателей (11,6 вграфе 3) их порядковы е н ом ера (ран ги) сум м ирую тся(4+5=9), а сум м а делитсян а число один аковы х показателей (вн аш ем прим ерен а 2). О ба показателяим ею тодин итотж еран г(4,5), а следую щ ем у повеличин епоказателю (12,9) присваивается6-й ран г.
2.В ы числяетсяразн ость(d) м еж ду ран гам ивотдельн ы ен еделигода. Результаты зан осятсявграфу 6.
3.Разн остьм еж ду ран гам ивозводитсявквадрат, иполучен н ы едан -
ны есум м ирую тся(графа 7). В н аш ем прим ереåd2 = 133,5.
48
4. П олучен н ы езначен ияподставляю тсявформ улу:
6 ×åd2
1-r =nn(2 -1) ,
гдеρ –коэффициен транговой корреляции; n –числопаркоррелируем ы х рядов(внаш ем прим ере12); åd2 –сум м а разностим еж ду рангам идвух корре-
лируем ы х рядов; 6 –постоянны й коэффициен т. Д лянаш егоприм ера: ρ=0,53. П реж де чем делать заклю чен ие о н аличии корреляцион н ой связи, необх одим оубедитьсяврепрезен тативностикоэффициен та, так как он получен н а осн овевы борочн огоисследован ия. П редставительн остькоэффициента ран говой корреляции определяетсяпутем сопоставлен ияеговеличин ы с критическим и зн ачен иям и, приведен н ы м и вприлож ен ии 3. Е сли получен н аявеличин а при задан н ом числеисследован ий окаж етсябольш е критического уровн я, будут осн ован ияговоритьо наличии корреляцион - ной связи. В н аш ем прим ерекоэффициен травен 0,53, а критическаявеличин а приР<0,05 составляет0,58 дляn=12. Следовательн о, получен н оезн а-
чен иекоэффициента н иж етого, при котором м ож н о бы лобы |
говоритьо |
|||||||||||||||||||||
наличиикорреляцион н ой связи. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
М |
етод лин ейн ой корреляции используется при наличии прям оли- |
|||||||||||||||||||||
нейн ой связи м еж ду взаим освязан н ы м и призн акам и. |
П оследовательн ость |
|||||||||||||||||||||
расчета коэффициен та лин ейн ой корреляции (rх у) рассм отрим |
|
на преды - |
||||||||||||||||||||
дущ ем |
прим ере. Д ля этого зан осим |
дан н ы е граф 1, |
2, 3 таблицы 12 в |
|||||||||||||||||||
таблицу 13. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т аблица 13 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
П роцен тн е- |
О тклон ен иеот |
К вадратотклон е- |
|
|
|
|
|
|
|
|||||||||||
Н еде- |
Заболевае- |
стан дартн ы х |
средн ей арифм е- |
н ияотсредн ей |
|
|
|
|
|
|
|
|||||||||||
|
d x2 ×d y2 |
|
||||||||||||||||||||
ляго- |
м ость |
пробкули- |
|
тической |
арифм етической |
|
|
|||||||||||||||
да |
О К И , (x) |
н арн ы х изде- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
лий (y) |
dx |
d y |
d x2 |
|
d y2 |
|
|
|
|
|
|
|
||||||||
1 |
2 |
3 |
4 |
|
5 |
|
6 |
|
|
7 |
|
|
8 |
|
|
|
||||||
1 |
12,8 |
8,3 |
-1,2 |
-5,9 |
1,44 |
|
34,81 |
|
7,08 |
|
|
|||||||||||
2 |
11,6 |
5,1 |
-2,4 |
-9,1 |
5,76 |
|
82,81 |
|
21,84 |
|
|
|||||||||||
3 |
15,4 |
18,0 |
1,4 |
3,8 |
1,96 |
|
14,4 |
|
5,32 |
|
|
|||||||||||
4 |
17,8 |
14,0 |
3,8 |
-0,2 |
14,44 |
|
0,04 |
|
-0,76 |
|
|
|||||||||||
5 |
12,3 |
11,6 |
-1,7 |
-2,6 |
2,89 |
|
6,76 |
|
4,42 |
|
|
|||||||||||
6 |
12,5 |
14,2 |
-1,5 |
0 |
|
2,25 |
|
0 |
|
|
0 |
|
|
|
||||||||
7 |
14,0 |
24,4 |
0 |
|
10,2 |
0 |
|
|
104,04 |
|
0 |
|
|
|
||||||||
8 |
20,3 |
25,4 |
6,3 |
11,2 |
39,69 |
|
125,44 |
|
70,56 |
|
|
|||||||||||
9 |
15,6 |
11,6 |
1,6 |
-2,6 |
2,56 |
|
6,76 |
|
-4,16 |
|
|
|||||||||||
10 |
14,8 |
12,9 |
0,8 |
-1,3 |
0,64 |
|
1,69 |
|
-1,04 |
|
|
|||||||||||
11 |
11,6 |
14,6 |
-2,4 |
0,4 |
5,76 |
|
0,16 |
|
-0,96 |
|
|
|||||||||||
12 |
9,3 |
10,0 |
-4,7 |
-4,2 |
22,09 |
|
17,64 |
|
19,74 |
|
|
49
1.Сум м ирую тсядан ны еряда х (графа 2), и сум м а делитсян а число наблю дений.
2.Сум м ирую тсядан ны еряда у (графа 3), и сум м а делитсян а число наблю дений.
3.О пределяется разность (d) каж дого показателя ряда х отсредней
арифм етической ( x =14,0) и ряда у от своей средн ей арифм етической ( y =14,2). П олучен н ы езн ачен иязаносятсясоответствен н овграфы 4 и5. Результаты будутим етькак полож ительн ы е, так иотрицательны езн ачения.
4. П олучен н ы евеличин ы dx и dy возводятсявквадрат, а результаты зан осятсявграфы 6 и7, послечегосум м ирую тсядан н ы е.
5. П ерем нож аю тсяпоказатели d x и d y сучетом полож ительн ы х и
отрицательн ы х значен ий (графа 8). П олучен ны е произведен ия сум м иру- ю тся.
6. Ц ифровы едан ны еподставляю тсявформ улу:
|
|
å |
|
|
|
× |
|
|
) d (d |
|||
r = ± |
|
x |
y |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
xy |
|
|
|
|
|
|
|
|
|
|
, гдеrх у –коэффициен тлин ейн ой корреляции. |
|
|
|
åd2x × åd2y |
Д лян аш егоприм ера: rxy= ± 0,616
П оскольку получен ны й коэффициен т определен врезультатевы бо-
рочн ого исследован ия, |
н еобх одим о убедиться встепен и его н адеж н ости. |
П редставительн ость rxy |
определяется по прилож ению 4. К оэффициен т |
считается представительн ы м , если получен н ая величин а превы ш ает критическоезначен иеприР<0,05 изадан н ом числестепен ей свободы (n’=n-2). В наш ем прим ере полученн ая величин а 0,616 превы ш ает критический уровен ь0,576 прип'= 10 иР<0,05.
П ри отсутствии таблицы критических зн ачен ий коэффициен та корреляции его представительн ость м ож ет бы ть определен а по величин е
средн ей ош ибки Srxy |
: |
|
|
|
|
|
|||||||||||
|
|
|
1− r2 |
|
|
|
|
|
|
|
|
|
|||||
Sr |
= ± |
|
|
|
xy |
|
|
|
, причислепарн ы х н аблю ден ий больш е100 или |
||||||||
|
|
|
n |
||||||||||||||
xy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
1− r2 |
|
|
|
|
|
|
|
|
|
||||
Sr |
= ± |
|
|
|
xy |
|
|
|
, причислепарн ы х н аблю ден иим ен ьш е100. |
||||||||
|
|
n − 2 |
|
|
|||||||||||||
xy |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
О тн ош ен иевеличин ы коэффициен та корреляциик величин есредн ей |
|||||||||||||||||
ош ибкипозволяетн айтизн ачен иеt: |
|
|
|
|
|
||||||||||||
|
rxy × |
|
|
|
|
|
|
|
|
|
|
|
|
||||
n |
|
|
|
|
|
|
|
|
|
||||||||
t = |
|
приn>100, или |
t = r |
|
n − 2 |
|
приn<100. |
||||||||||
|
|
|
|
|
|||||||||||||
1- r2 |
|
|
|
− r2 |
|
||||||||||||
|
|
|
|
xy 1 |
|
||||||||||||
|
|
|
|
|
xy |
|
|
|
|
|
|
xy |
|
|
|||
В н аш ем |
прим ере: t=2,471. П о таблице зн ачен ий t-критерия при |
||||||||||||||||
Р<0,05 ичислестепен ей свободы 10 t =2,23. |
|
|
|
50
П ослеустан овлен ияпредставительности получен н ого значениякоэффициента определяетсятесн ота связи, вы раж аю щ аясявеличин ой от0 до 1 (табл.14).
|
Т аблица 14 |
Степен ьсвязи |
В еличин а коэффициен та |
М алая(слабая) |
0,01 –0,30 |
Средн яя(ум еренн ая) |
0,31 –0,70 |
Больш ая(сильная) |
0,71 –1,00 |
Степень« связан н ости» ввариациидвух илинескольких изучаем ы х явлений болееточн оизм еряетсяквадратом коэффициента корреляции(rх у). Э то значит, чтоприrх у= 0,5 25% изм енений одн огопризнака объясняетсявариацией другого. П риrх у= 0,3 изм енчивостьодн огоявлениязакон ом ерн освязана сизм енением другогов10%, при rх у = 0,7 –в50%, приrх у = 0,9 –в81% случаев. Т аким образом , х отякоэффициенткорреляциииуказы ваетна общ - ностьэлем ентоввкоррелируем ы х рядах , ноневсяэта общ ностьобъясняется закон ом ерн ой связью ввариации признаков. И з сказанногоясн о, чтоотесной корреляциим ож н оговоритьлиш ьвслучаях , когдаrх у нениж е0,7.
П равильн аятрактовка коэффициен та корреляции предполагаетн ор- м альн оераспределен иесопряж ен ны х величин коррелируем ы х рядовх иу. О днакоприм алом числен аблю ден ий исравн ительн осильн ой корреляции
(rх у=0,5) распределен иекоэффициен та rх у отличаетсязн ачительн ой асим - м етрией отн орм альн огораспределен ия(рис.15).
rxy=0 rxy=+0,4 rxy=+0,8
Рис.15. К ривы ераспределен ияэм пирическогокоэффициен та корреляции приn=12 дляразличн ы х зн ачений генеральн огопарам етра rxy
(поА .К . М итропольском у, 1971)
Следовательн о, эм пирический коэффициен т корреляции н е будет точн ой оцен кой длявсей ген еральн ой совокупн ости, еслион определен н а м алочислен н ой вы боркеи еговеличин а больш е0,5. В этих случаях коэффициен т rх у целесообразн о зам ен ить преобразован н ой величин ой z (z — преобразован ие Ф иш ера), связан н ой с коэффициен том корреляции сле-