Файл: Лакин Г.Ф. Биометрия учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 289

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

п /0,5 X 0,5

,---------

т <і = ^ ^ - г Ѵ - =

УО,0018 = 0,043 и т. д.

'

141

 

Для уровня значимости Р = 0,01 и числа степеней свободы, кото­ рое в таких случаях определяется по формуле k = n — 1= >30 по таблице Стьюдента (см. табл. V приложений) находим tst = 2,6. Так как во всех случаях /ф<7г, нулевая гипотеза сохраняется, вопрос о влиянии возраста производителей на пол потомства остается открытым.

Остановимся на одном из классических опытов Грегора Мен­ деля по скрещиванию двух константных сортов гороха, отличав­ шихся друг от друга по окраске семян: у одного были желтые, а у другого сорта зеленые семена. Во втором поколении от этого скрещивания, т. е. от посева гибридных семян, были получены следующие результаты (табл. 48).

Т а б л и ц а 48

М семьи (отдельные

Окраска

семядолей

Отношение желтых

растения)

желтая

|

зеленая

к зеленым

 

 

 

1

25

 

И

1 : 2 ,2 7

2

32

 

7

1 : 4 , 5 7

3

14

 

5

1 : 2 , 8 0

4

70

 

27

1

: 2 , 5 9

5

24

 

13

1

: 1 ,8 5

6

20

 

6

1

: 3 , 3 3

7

32

 

13

1

: 2 , 4 6

8

44

 

9

1

: 4 , 9 0

9

50

 

14

1 : 3 ,5 7

10

44

 

18

1

: 2 , 4 4

И т о г о . . .

355

 

123

 

1 : 2 , 8 9

В с е г о . . .

 

478

 

 

 

Необходимо установить, соответствуют ли полученные в опыте результаты ожидаемому отношению желтых к зеленым как 3:1. Исходя из теоретической предпосылки, ожидаемое отношение

желтых к зеленым должно быть как

3ДХ478 = 358,5 желтых к

74 X 478= 119,5

зеленым. Откуда доля

зеленых в генеральной

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

=

= 0,26. Разность р—Р 0,26—0,25 = 0,01.

Средняя ошибка этой разности

160


md = ] / 0’25^ -0-7- =

У0,0004 =

0,02.

 

478

 

 

 

Критерий достоверности t — 0,01

0,5. Эта

величина не

пре­

(Г02

 

 

 

вышает £Sf = l,96 для Р —0,05. Следовательно,

отвергнуть

нуле­

вую гипотезу нельзя; в данном случае имеется хорошее совпа­ дение опытных данных с ожидаемыми в Р2 по второму закону Менделя.

Оценка разности между эмпирическими и ожидаемыми частотами альтернатив по критерию хи-квадрат

При оценке различий, наблюдаемых между

эмпирическими

и ожидаемыми или вычисленными частотами

альтернативных

признаков, можно использовать .критерий Пирсона % 2 (хи-квад- рат). Покажем это на только что рассмотренном примере. Ожи­ даемые и полученные в опыте с горохом отношения желтых се­ мян к численности зеленых нам уже известны. Разность между

этими отношениями равна: для желтых семян

р—Р = 355—

—358,5 = —3,5 и для зеленых— 123—119,5= +3,5.

Возведем эту

разность в квадрат, полученный результат отнесем к ожидаемым частотам каждой альтернативы, т. е. к численности желтых и зеленых семян и, сложив полученные отношения, найдем вели­ чину критерия хи-квадрат:

Х 2 _

(~ 3.:5)2+

(+ 3І_5)2 =

о,034 +

0,103 = 0,137 0,14.

Л

358,5 ^

119,5

> -г

>

Остается сравнить полученную в опыте величину % 2 с ее кри­ тическим значением для взятого уровня значимости (Р) и числа степеней свободы, которое при табличной группировке частот определяется по следующей формуле:

k= (г -1 ) (с -1 ),

где г обозначает число граф (столбцов), а с — число строк таб­ лицы, в которой помещены сравниваемые величины; при этом итоговые столбцы и строки таблицы в расчет не принимаются. В данном случае сравниваются две группы признаков — желтые и зеленые семена. Следовательно, число степеней свободы к — = 2—1 = 1. Этому числу для уровня значимости Р = 0,05 соответ­

ствует %2st — 3,8

(см. табл.

VIII приложений).

Так как

%2ф =

= 0 ,1 4 < 3c2s< = 3 ,8 ,

нулевая

гипотеза сохраняется; наблюдаемые

расхождения между ожидаемыми и опытными

данными

носят

случайный характер.

 

 

 

6—2802

161


Применим тот же критерий к оценке результатов дигибридного скрещивания двух сортов гороха, с которыми эксперимен­ тировал Мендель, сравним результаты расщепления признаков во втором поколении с ожидаемым отношением 9 : 3 : 3 : 1 (табл. 49).

Т а б л и ц а 49

 

Данные

Желтые

Зеленые

Желтые

Зеленые

Всего

 

гладкие

гладкие

морщи­

морщи­

 

 

 

 

 

нистые

нистые

 

Наблюдаемые (р)

315

101

108

32

556

Ожидаемые

(р')

313

104

104

35

556

Разница

p')=d

2

3

4

3

Квадрат

(d2)

/ tf2 \

4

9

16

9

Отношение

0,01

0,09

0,15

0,26

0,51

^ — ; j

Приведенные в табл. 49 ожидаемые частоты вычислены умноже­ нием 9/іб желых гладких, 3/іб зеленых гладких, 3/і6 желтых мор­ щинистых и '/іб зеленых морщинистых на 556, т. е. на общее число наблюдений. Так как в данном случае число граф равно 4, а число строк — двум, число степеней свободы /г=(4—1)(2—

—1) =3. Этому числу для Р = 0,05 соответствует %2Sf = 7,8. Видим, что %2 0 <%2 st, следовательно, нулевая гипотеза сохраняется; рас­ хождения между ожидаемыми и полученными в опыте данными носят чисто случайный характер.

Критерий хи-квадрат можно использовать и для оценки раз­ личных состояний одного и того же качественного признака, группируемых в таблицы разной конструкции. Например, при изучении эффекта от добавочного опыления подсолнечника были получены следующие результаты (табл. 50).

 

 

 

Т а б л и ц а 50

 

Количество семян в корзинке

Опыление

заполненных

пустых

Всего

 

 

Естественное ....................

113

42

155

Д обавочное......................

131

11

142

В с е г о . . .

244

53

297

При таком расположении выборочных данных критерий хи-квад­ рат определяется по формуле:

Г1

 

n(ad cby

+

(89)

 

Ь) (с + d) (а + с) (Ь + d)

162


где а, b, с и d — численности

альтернатив, помещенные в клет­

ках четырехпольной таблицы;

n = a + b + c+d — общее число на­

блюдений. В данном случае критерий хи-квадрат, рассчитывае­ мый по указанной формуле, оказывается равным:

2

297(113 X 11 — 131 Х 42)2

тг =

----------------------------------- = 18,9.

А2 4 4 X5 3 X 155X 142

Для уровня значимости Р = 0,05 и числа степеней свободы k = = (2—1)(2—1) = 1 критическое значение критерия %2st = 3,8 (см. табл. VIII приложений). Полученная в опыте величина это­ го критерия почти в пять раз превосходит его стандартное (кри­ тическое) значение, что отвергает нулевую гипотезу. Следова­ тельно, нужно признать, что влияние добавочного опыления на урожай подсолнечника — явление не случайное, а вполне зако­ номерное.

Расчет критерия хи-квадрат по формуле 89, как это видно из приведенного примера, часто связан с перемножением много­ значных чисел, что создает известные неудобства. Чтобы избе­ жать последних, можно изменить группировку выборочного ма­ териала и определить критерий хи-квадрат по формуле 72. Про­ демонстрируем это на том же примере. Расчет критерия хи-квадрат показан в табл. 51. Ожидаемые частоты (р') рассчи­ тываются по итоговым числам четырехпольной таблицы. На­

пример,

частота

127,3,

что

стоит первой

во

второй

строке

табл. 51,

получена на

основании данных табл.

50 следующим

образом:

(244Х'155) : 297=127,3. Следующая за ней ожидаемая

частота

116,7

определена

аналогичным

способом:

(244Х

Х І42) : 297=116,7. Ожидаемая частота 27,7 найдена таким же способом: (53X155) : 297 = 27,7 и т. д. Остальные действия вид­ ны из табл. 51. Результат получился, как и следовало ожидать, тот же, т. е. % 2 = 18,85= 18,9.

 

 

 

 

 

Т а б л и ц а 51

 

 

Заполненные семена

Пустые семена ■

 

Группы сравнения

естест.

добав.

естест.

добав.

Всего

 

 

опыление

опыление

опыление

опыление

 

Данные опыта (р)

113

131

42

п

297

Ожидалось

(/?')

127,3

116,7

27,7

25,3

297

Разница р p' = d

14,3

14,3

14,3

14,3

Квадрат d2

d2

204,5

204,5

204,5

204,5

Отношение

1,6

1,8

7,4

8,1

18,85

-

Р

Как уже сообщалось, распределение хи-квадрат является непрерывным. Поэтому при использовании этого критерия для

6:

163


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

/

п у

[ \ ad — bc\— — ) Х п

у 2 — -------------------------------------------------------

------------------------------------------ ( 9 0 1)

+ Ь) (с + d) + с) (b + d)

Эта формула обеспечивает более точные результаты, особенно на малочисленных выборках.

Применим эту формулу к рассмотренному выше примеру по испытанию действия эндотоксина на выживаемость облученных животных. Соответствующие данные приведены в табл. 46. Мы повторим их здесь для того, чтобы вместе с экспериментальными данными привести и их теоретически вычисленные значения (они заключены в скобки) (табл. 52).

 

 

 

 

 

Т а б л и ц а 52

Группы особей

Выжило

Погибло

Всего

Опытные . . . . . I . . .

23

(18,72)

13

(17,28)

36

Контрольные.....................

3

(7,28)

11

(6,72)

14

В с е г о . . .

 

26

 

24

50

Если вычислить критерий хи-квадрат для этих данных без по­ правки Йейтса, получается следующая величина:

 

50(23 X П - З Х 13)2

2 289 800

 

1

~~ 26 X 24 X 36 X 14 —

314496 _

’ '

Внесение поправки Йейтса значительно снижает ее:

 

50 ( 123 X 11 — 3 X 131— 25)2

1 786 050

 

Х “

26 X 24 X 36 Х~Т4

_ 314 496*

 

* Прямые

скобки, в которых заключена разность ad be,

указывают на

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

164