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

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

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

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

Добавлен: 11.04.2024

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

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

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

Характеристики этого распределения оказались: х = 394,0 мг; я=41,6 мг\ С1/ = 10,6% и Лх=+0,56. Видно, что перегруппировка вариант сильнее сказалась на величине коэффициента асиммет­ рии и совсем незначительно на других характеристиках ряда. Из приведенных примеров видно, что асимметрия действительно мо­ жет быть кажущейся, возникающей по чисто техническим причи­ нам как следствие группировки совокупности наблюдений в ва­ риационный ряд. Вместе с тем асимметрия может иметь и не случайный, а систематический характер, являясь следствием гетерогенности выборочного материала. В литературе имеются указания на то, что асимметрия и эксцесс могут возникать не только при смешении генетически неоднородного материала, но и вследствие модифицирующих влияний условий среды, в кото­ рых существует популяция. В частности, асимметрия в распреде­ лении некоторых признаков у сельскохозяйственных растений вызывается такими условиями, как засуха и невыравненность агрофона в отношении питания растений (Иогансен, 1933; Кон­ стантинов, 1952 и др.). В качестве примера приводим данные, по­ казывающие, как при прочих равных условиях изменяются ха­ рактеристики распределения озимой ржи по длине колосьев в за­ висимости от густоты стояния растений в посевах этой культуры |табл. 34).

Т а б л и ц а 34

Количество

Измерено

Средняя

Основное от­

Коэффициент

Коэффициент

растений

колосьев (л)

длина колосьев

клонение (а)

вариации (СѴ)

асимметрии

на 1 м ъ

 

( М М ) (X)

 

 

(A s )

197

132.

64,0

19,0

27,7

+ 0,004

222

115

63,3

18,5

29,2

+ 0,410

364

182

49,3

18,2

36,9

4 0,620

С уплотнением массы растений на единице площади умень­ шается средняя длина колосьев, и параллельно с этим увеличи­ вается значение коэффициента асимметрии.

Итак, отклонения от нормального распределения могут быть вызваны разными причинами. По-видимому, не всегда эти при­ чины носят случайный или временный характер. Многие факты говорят о том, что в отношении ряда изученных признаков асим­ метрия является неотъемлемым атрибутом их варьирования. Та­ кого рода отклонения от кривой нормального распределения мо­ гут служить показателем величины тех изменений, которые признак претерпевает в процессе своего развития. Разумеется, помогая исследователю характеризовать ту или иную форму рас­ пределения, биометрия не претендует на причинный анализ асим­ метрии и других отклонений от нормального закона. Она лишь может указывать, в каком направлении следует вести биологиче­

232


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

ОЦЕНКА ЭМПИРИЧЕСКИХ РЯДОВ ПО НОРМАЛЬНОМУ ЗАКОНУ

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

Для нахождения теоретических частот вариационного ряда по нормальному закону Гаусса — Лапласа служит следующая рабо­ чая формула:

р'

п X і

о

 

где р' — вычисленная или теоретическая частота; п —-сумма всех частот, т. е. объем эмпирического вариационного ряда; і — вели­

чина классового интервала эмпирического ряда;

о — среднее

квадратическое отклонение эмпирического ряда; f(t)

— функция

нормированного отклонения, значения которой (ординаты нор­ мальной кривой) приведены в табл. II приложений.

Таким образом, для нахождения теоретических частот эмпи­ рического вариационного ряда достаточно значения вероятности, соответствующие нормированным отклонениям частот эмпириче­

ского ряда умножить на - ^ 1. Продемонстрируем эту методику

о

на примере распределения 267 датских угрей по числу позвонков. Средняя арифметическая этого ряда 114,74 позвонков, а 0 = 1,35 позвонков. Для получения теоретических частот этого распреде­ ления нормируем его, т. е. находим отклонения вариант от их средней арифметической, выражая каждое отклонение в долях среднего квадратического отклонения (табл. 35).

Затем из табл. II приложений выписываем ординаты нормаль-

X X

ной кривой для каждого значения t = ------- эмпирического рас-

0

пределения. Перемножая значения ординат на величину п Х і»


 

 

 

 

 

Т а б л и ц а 35

 

Эмпири­

 

 

 

Теоретііческие

 

Отклонения

 

Ординаты

частоі ы (p')

Варианты

ческие

І Х~ Х

 

 

нормальной

 

 

(X)

частоты

 

 

 

(л--Т)

<3

фактичес­

 

 

(р )

кривой

округленно

 

 

 

 

 

ки

 

 

 

 

 

 

і и

3

- 3 ,7 4

—2,77

0,0086

1,7

2

112

9

- 2 ,7 4

- 2 ,0 3

0,0508

10,0

10

113

31

- 1 ,7 4

- 1 ,2 9

0,1739

34,3

34

114

71

- 0 ,7 4

- 0 ,5 5

0,3429

67,9

68

115

82

+ 0,26

+0,19

0,3918

77,6

78

116

46

+ 1,26

+0,93

0,2589

51,1

50

117

19

+2,26

+ 1,67

0,0989

19,5

20

118

5

+ 3,26

+2,41

0,0219

4,3

4

119

1

+4,26

+ 3,15

0,0028

0,5

1

Сумма

267

267,2

267

 

 

*

267 V 1

igg, получаем теоретически

равную в данном случае ____ _ _=

 

 

 

1,35

 

 

 

вычисленные частоты данного распределения (см. последнюю гра-

Рис. 17. Эмпирическая и выравненная по нормально­ му закону кривые распределения 267 датских угрей по числу позвонков:

на оси абсцисс — количество позвонков, на оси ординат — частоты

фу табл. 35). Из табл. 35 видно, что вычисленные и эмпирические частоты в общем неплохо согласуются между собой, что более наглядно показано на рис. 17.

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

134


кретных распределений. Покажем это на примере распределе­ ния кальция (мг%) в сыворотке крови павианов-гамадрилов с его характеристиками: ж = 11,9 мг% и о=1,2 мг%. В данном слу-

і Х п

0,7 X 100

 

Ход вычислений показан

в

чае ------ = -------—— =58,3 = 58.

о

 

1,2

 

 

 

 

 

 

табл. 36.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

36

 

Эмпири­

 

 

Ординаты

Теоретические час­

Варианты

і Отклонения

Х —Х

тоты (р ’)

 

ческие

(л--7)

<3

нормаль­

 

 

 

(X)

частоты

ной

кри-

фактичес­

округ­

 

ІР)

 

 

вой

1 /(0 1

 

 

 

 

 

 

ки

ленно

 

8,9

2

- 3 ,0

- 2 ,5 0

0,0175

1,0

1

 

9,6

3

- 2 ,3

—1,91

0,0644

3,7

4

 

10,3

9

- 1 ,6

- 1 ,3 3

0,1647

9,6

10

 

11,0

17

—0,9

-0 ,7 5

0,3011

17,6

18

 

11,7

25

- 0 ,2

- 0 ,1 7

0,3932

22,8

23

 

12,4

23

+ 0,5

+0,42

0,3653

21,2

21

 

13,1

10

+ 1,2

-+1,00

0,2420

14,0

14

 

13,8

7

+ 1,9

+ 1,58

0,1145

6,6

7

 

14,5

4

+ 2 ,6

+2,17

0,0379

2,2

2

 

Сумма

100

 

98,6

100

 

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

ОЦЕНКА ЭМПИРИЧЕСКИХ РЯДОВ ПО ЗАКОНУ ПУАССОНА

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

ат

табл. III), в которой помещены значения — е~а ,входящие в co­ rn!

став формулы, определяющей функцию распределения редких событий. Зная среднюю (ж) эмпирического распределения, по указанной таблице находят вероятности, соответствующие значе­ ниям т : 0, 1, 2, 3 и т. д. Затем вероятности умножаются на общее число наблюдений (п) . В результате получаются вычисленные или теоретические частоты (р') для каждого т.

Продемонстрируем эту методику на примере распределения поражаемости клеток альфа-частицами по данным, приведенным в табл. 22 (пятая глава). Средняя арифметическая этого распре­

135


деления я = 1,54. Расчет теоретических частот показан в следую­ щей табл. 37.

Т а б л и ц а 37

 

 

 

Теоретические

частоты

Число клеток

Число случаев

Вероятности

час­

 

(т )

(Р)

тот по формуле

округленно

 

 

Пуассона

фактически

0

112

0,2231

115,34

115

1

168

0,3347

173,04

173

2

130

0,2510

129,77

130

3

68

0,1255

64,88

65

4

32

0,0471

24,35

24

5

5

0,0141

7,29

7

6

1

0,0035

1,81

2

7

1

0,0008

0,41

1

Сумма . .

517

516,90

517

Видно, что вычисленные частоты довольно хорошо согласуются с распределением эмпирических частот этого ряда. Следователь­ но, есть достаточные основания считать, что в данном случае час­ тота поражаемое™ клеток альфа-частицами распределяется по закону Пуассона.

КРИТЕРИИ СООТВЕТСТВИЯ ЭМПИРИЧЕСКИХ И ТЕОРЕТИЧЕСКИХ РАСПРЕДЕЛЕНИЙ

Критерий хи-квадрат (х2)

Статистическая оценка расхождений, наблюдаемых между эмпирическими и теоретическими или ожидаемыми частотами вариационного ряда, производится с помощью особых критериев соответствия, называемых также критериями согласия. Одним из таких критериев, широко используемых в биометрии, является критерий X2 (хи-квадрат), предложенный Пирсоном в 1901 г. Он представляет сумму квадратов отклонений эмпирических частот от частот теоретачееких, отнесенную к теоретическим частотам вариационного ряда, т. е.

(р — р' ) 2

%2 = 2 U- nf L - t (72)

Р

где р — эмпирическая, а р' — соответствующая теоретическая,

т.е. вычисленная, или ожидаемая, частоты, 2 — знак суммы. Обозначив разность между эмпирическими и теоретическими

частотами через d, т„ е. р — p' — d,

формулу 72 можно выразить

в более простом виде:

 

 

d2

L(72a)

Х2 = S

- .

136