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

Биометрия_пособие2

.pdf
Скачиваний:
42
Добавлен:
18.03.2016
Размер:
2.79 Mб
Скачать

Таблица 2.1

Значение t-критерия Стьюдента

Число

0,10

0,05

0,02

0,01

0,002

0,001

степеней

 

 

 

 

 

 

свободы

 

 

 

 

 

 

1

6,314

12,71

31,82

63,66

318,3

636,6

2

2,920

4,303

6,965

9,925

22,33

31,60

3

2,353

3,182

4,541

5,841

10,21

12,92

4

2,132

2,776

3,747

4,604

7,173

8,610

5

2,015

2,571

3,365

4,032

5,893

6,869

6

1,943

2,447

3,143

3,707

5,208

5,959

7

1,895

2,365

2,998

3,499

4,785

5,408

8

1,860

2,306

2,896

3,355

4,501

5,041

9

1,833

2,262

2,821

3,250

4,297

4,781

10

1,812

2,228

2,764

3,169

4,144

4,587

11

1,796

2,201

2,718

3,106

4,025

4,437

12

1,782

2,179

2,681

3,055

3,930

4,318

13

1,771

2,160

2,650

3,012

3,852

4,221

14

1,761

2,145

2,624

2,977

3,787

4,140

15

1,753

2,131

2,602

2,947

3,733

4,073

16

1,746

2,120

2,583

2,921

3,686

4,015

17

1,740

2,110

2,567

2,898

3,646

3,965

18

1,734

2,101

2,552

2,878

3,610

3,922

19

1,729

2,093

2,539

2,861

3,579

3,883

20

1,725

2,086

2,528

2,845

3,552

3,850

21

1,721

2,080

2,518

2,831

3,527

3,819

22

1,717

2,074

2,508

2,819

3,505

3,792

23

1,714

2,069

2,500

2,807

3,485

3,767

24

1,711

2,064

2,492

2,797

3,467

3,745

25

1,708

2,060

2,485

2,787

3,450

3,725

26

1,706

2,056

2,479

2,779

3,435

3,707

27

1,703

2,052

2,473

2,771

3,421

3,690

28

1,701

2,048

2,467

2,763

3,408

3,674

29

1,699

2,045

2,462

2,756

3,396

3,659

30

1,697

2,042

2,457

2,750

3,385

3,646

 

1,645

1,960

2,326

2,576

3,090

3,291

11

Пример: В ходе подсчета числа семян в корзинке 50-ти одуванчиков (Taraxacum officinale), произрастающих на 2-х различных участках города были получены средние показатели количества числа семян и их стандартные ошибки. Так, на первом участке среднее количество числа семян в корзинке и стандартная ошибка составили 187,4 ± 5,8 шт., а на втором - 171,6 ± 6,3 шт. Необходимо ответить на вопрос, достоверна ли разница и уровень значимости этой разница между этим показателем популяций Taraxacum officinale, произрастающих на разных территориях города.

Для ответа на этот вопрос:

1)находим t-критерий Стьюдента по вышеприведенной формуле: t = (187,4 – 171,6) / √((5,82) + (6,32)) = 15,8 / √(33,64 + 39,69) = 15,8 / √ 73,33 = 15,8 / 8,6 = 1,8;

2)определяем число степеней: 50 + 50 – 2 = 98.

3)по таблице определяем уровень достоверности разницы между количеством семян в корзинке Taraxacum officinale, популяции которого произрастают на 2-х различных участках города. Он определяется примерно, как 0,05. То есть разница достоверна на уровне 95%.

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

терию Z проводится по формуле: Z = хср. уср. / ( 2х/nx) + ( 2y/ny). Затем по таблице 2.2, используя показатель Z при числе степеней свободы, равном nx + ny - 2 нужно найти вероятность Р. Если Р 0,95, то расхождение между двумя средними следует признать неслучайным и существенным, то есть различие двух средних арифметических достоверное.

Таблица 2.2

Определение показателя вероятности P по критерию Z

Z

P

Z

P

Z

P

 

 

 

 

 

 

 

 

 

 

 

0,1

0,080

1,4

0,838

2,7

0,993

0,2

0,159

1,5

0,866

2,8

0,995

0,3

0,236

1,6

0,890

2,9

0,996

0,4

0,311

1,7

0,911

3,0

0,997

 

 

 

12

 

 

0,5

0,383

1,8

0,928

3,1

0,9981

0,6

0,452

1,9

0,943

3,2

0,9986

0,7

0,516

2,0

0,954

3,3

0,9990

0,8

0,576

2,1

0,964

3,4

0,9993

0,9

0,632

2,2

0,972

3,5

0,9995

1,0

0,683

2,3

0,979

3,6

0,9997

1,1

0,729

2,4

0,984

3,7

0,9999

1,2

0,770

2,5

0,988

3,8

0,9999

1,3

0,806

2,6

0,991

 

 

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

Контрольные вопросы и задание для самостоятельной работы к теме 2:

1.Из какого предположения исходят, когда сравнивают 2 выборки по t-критерию Стьюдента?

2.Зачем полученное значение t – критерия Стьюдента сравнивают с табличными значениями?

3.Как определяют достоверное отличие 2-х выборок по Z- критерию?

4.При проведении биотестирования проб воды из реки Белой в районе города Стерлитамака были получены коэффициенты роста (количество листецов) для тест-объекта ряска малая (табл. 2.3).

Таблица 2.3

Количество листецов у ряски малой на 7-ой день эксперимента

 

 

Количество

Место отбора проб воды

Повторности

листецов на

 

 

7-й день

Река Белая (р-н турбазы

 

 

«Дружба»)

 

 

 

1

19

 

2

18

 

3

16

 

4

17

13

 

5

18

 

6

19

 

7

15

 

8

18

 

9

17

Река Белая (р-н «Детской коло-

 

 

нии»)

 

 

 

1

17

 

2

17

 

3

16

 

4

19

 

5

18

 

6

21

 

7

22

 

8

19

 

9

18

Определить достоверность отличия коэффициентов роста ряски малой в пробах воды из реки Белой по t-критерию и по Z- критерию.

Тема 3. Дисперсионный анализ

Дисперсионный анализ - статистический метод выявления влияния отдельных факторов на результат эксперимента. Первоначально дисперсионный анализ был предложен английским статистиком Р. Фишером в 1925 г. для обработки результатов агрономических опытов по выявлению условий, при которых испытываемый сорт сельскохозяйственной культуры даёт максимальный урожай.

Современные приложения дисперсионного анализа охватывают широкий круг задач биологии и трактуются обычно в терминах статистической теории выявления систематических различий между результатами непосредственных измерений, выполненных при тех или иных меняющихся условиях. Если значения неизвестных постоянных a1,...., an могут быть измерены с помощью различных методов или измерительных средств M1,..., Мm и в каждом случае систематическая ошибка может

14

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

xij = ai, + bij + δij,

i = 1, 2,..., n; j = 1, 2,..., m, где bij — систематическая ошибка, возникающая при измерении ai по методу Mj, δij — случайная ошибка. Такая модель называется двухфакторной схемой Д. а. (первый фактор — измеряемая величина, второй — метод измерения).

Пример: На учебно-опытном участке изучали влияние различных способов внесения в почву органических удобрений на урожай зеленой массы кукурузы. Опыт проводили на десятиметровых делянках в трех вариантах, не считая контроля. Каждый вариант опыта имел трехкратную повторность. Результаты опыта занесены в таблицу 3.1.

 

 

 

 

Таблица 3.1

Результаты урожайности зеленой массы кукурузы при

различных способах внесения удобрений

 

 

 

 

 

 

 

Варианты опыта

Урожай по

 

Средний

 

повторностям, кг

урожай,

 

 

 

 

 

xiср

 

1

2

 

3

 

Контроль

21,2

28,0

 

31,2

26,8

Удобрения помещали:

 

 

 

 

 

ниже семян на 4 см

23,6

22,6

 

28,0

24,7

в стороне от семян на 4 см

24,0

30,0

 

29,2

27,7

выше заделки семян на 4 см

29,2

28,0

 

27,0

28,1

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

Обозначим фактор, регулируемый в опыте через А, а его градации (варианты опыта) – соответственно через А1, А2, А3 и А4. Для упрощения расчетов вспомогательных величин:

1)уменьшим каждую варианту комплекса на 20;

2)сгруппируем преобразованные данные так, чтобы гра-

15

дации фактора А помещались в верхней части комбинационной таблицы (табл. 3.2).

Таблица 3.2

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

Урожай

 

Градации фактора А

 

Суммы

по по-

 

 

(варианты опыта)

 

 

вторно-

А1

 

А2

А3

 

А4

 

стям X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

3,6

4,0

 

9,2

 

xi - 20

8,0

 

2,6

10,0

 

8,0

а = 4

 

11,2

 

8,0

9,2

 

7,0

 

n

3

 

3

3

 

3

Σn = N =

 

 

 

 

 

 

 

12

Σ xi

20,4

 

14,2

23,2

 

24,2

Σ xi = 82

(Σ xi)2

416,16

 

201,64

538,24

 

585,64

Σ (Σ xi)2 =

 

 

 

 

 

 

 

1741, 68

Σ xi2

190,88

 

88,72

200,64

 

197,64

Σ xi2 =

 

 

 

 

 

 

 

672,88

Определение девиат:

1)находим Н = 822/ 12 = 560,33;

2)находим общую девиату Dy = Σ xi2 – H = 672,88 – 560,33

=112,55;

3)определяем факториальную девиату DA = ((Σ (Σ xi)2) / n) - H = (1741, 68 / 3) - 560,33 = 580,56 - 560,33 = 20,23;

4)вычисляем остаточную девиату DB = Dy - DA = 112,55 - 20,23 = 92,32.

Определение числа степеней свободы:

1)так как общий комплекс содержит 12 вариант, то число степеней свободы для общей дисперсии ky = N – 1 = 12 – 1 = 11;

2)фактор А содержит 4 градации (3 варианта опыта и контроль) и, поэтому, число степеней свободы для факториальной дисперсии kA = a – 1 = 4 – 1 = 3;

3)для внутригрупповой (остаточной) дисперсии число степеней свободы ke = ky - kA = 11 – 3 = 8;

16

Определение дисперсий:

1)находим факториальную дисперсию SA2 = DA / kA = 20,23 / 3 = 6,74;

2)вычисляем остаточную дисперсию Se2 = DB / ke = 92,32 / 8 = 11,54.

Вывод: В данном случае SA2 < Se2, а это означает, что межгрупповая вариация не превышает внутригруппового случайного уровня и, следовательно, считать достоверным влияние фактора на исследуемый признак нет оснований. Это можно проверить и в результате вычисления F-критерия = Se2 / SA2 = 11,54 / 6,74 = 1,71, что значительно меньше Fst = 8,85 для ke = 8 и kA = 3.

Контрольный вопрос и задание для самостоятельной работы к теме 3:

1.Что позволяет выявить дисперсионный анализ?

2.При проведении биотестирования в течении 4-х суток пробы сточной воды в реку Белая (район города Стерлитамака) с использованием лабораторного тест объекта Daphnia magna (дафния магна) были получены данные, которые занесены в таблицу 3.3.

Таблица 3.3

Результаты биотестирования пробы сточной воды в реку Белая

врайоне города Стерлитамака с использованием тест-объекта

D.magna

 

Варианты опыта

Количество выживших дафний,

 

 

 

шт

 

 

 

1

2

3

Контроль

29

28

30

Разбавления сточной воды:

 

 

 

в 2

раза

13

20

15

в 5

раз

20

23

22

в 10 раз

25

28

27

Определить, является ли достоверным влияния качества сточной воды на лабораторный тест-объект?

17

Тема 4. Исследование степени соответствия эмпирических и теоретических данных по различным критериям

Определение соответствия эмпирических и теоретических данных по критерию хи-квадрат (χ2). В биологических ис-

следованиях часто возникает необходимость определить соответствие эмпирического и теоретического распределения или двух или более эмпирических распределений между собой, то есть оценить соответствие фактических данных предлагаемой гипотезе. Это соответствие обычно измеряется особым показателем, обозначаемым как критерий χ2 (хи-квадрат). По этому критерию удается определить, имеется ли согласие или налицо различие в эмпирическом, т. е. полученном исследователем, и теоретическом распределениях.

Можно проследить, например, изменение каких-либо показателей под влиянием каких-то воздействий и т. д. В этих случаях полученные результаты вводят в таблицу в зависимости от количества исследуемых факторов. В связи с этим подсчеты критерия χ2 будут несколько различаться. После получения критерия χ2 исследователю необходимо определить согласие эмпирического распределения с теоретическим или различие их друг от друга. В данном случае необходимо испытание так называемой нулевой гипотезы. По таблице 4.1 при рассчитанном числе степени свободы (см. далее) находят теоретический показатель χ2 при Р<0,05.

 

 

 

 

Таблица 4.1

 

Значения χ2распределения

 

 

 

 

 

 

 

 

 

 

С

0,20

0,05

0,025

 

0,01

 

1

2,71

3,84

5,02

 

6,63

 

2

4,61

5,99

7,38

 

9,21

 

3

6,25

7,81

9,35

 

11,34

 

4

7,78

9,49

11,14

 

13,28

 

5

9,24

11,07

12,83

 

15,09

 

6

10,64

12,59

14,45

 

16,81

 

7

12,02

14,07

16,01

 

18,48

 

8

13,36

15,51

17,53

 

20,09

 

18

9

14,68

16,92

19,02

21,67

10

15,99

18,31

20,48

23,21

11

17,28

19,68

21,92

24,72

12

18,55

21,03

23,34

26,22

13

19,81

22,36

24,74

27,69

14

21,06

23,68

26,12

29,14

15

22,31

25,00

27,49

30,58

16

23,54

26,30

28,85

32,00

17

24,77

27,59

30,19

33,41

18

25,99

28,87

31,53

34,81

19

27,20

30,14

32,85

36,19

20

28,41

31,41

34,17

37,57

21

29,62

32,67

35,48

38,93

22

30,81

33,92

36,78

40,29

23

32,01

35,17

38,08

41,64

24

33,20

36,42

39,36

42,98

25

34,38

37,65

40,65

44,31

26

35,56

38,89

41,92

45,64

27

36,74

40,11

43,19

46,96

28

37,92

41,34

44,46

48,28

29

39,09

42,56

45,72

45,59

30

40,26

43,77

46,98

50,89

Если полученный в ходе расчетов эмпирический показатель равен или меньше этой величины, «нулевая» гипотеза не отвергается, т. е. в данном случае эмпирическое и теоретическое распределения не различаются, если вычисленный показатель больше полученного из таблицы при Р<0,05, значит «нулевая» гипотеза отвергается и эмпирическое распределение статистически достоверно отличается от теоретического, значит данное воздействие статистически достоверно изменяет результат исследования. Вместе с тем критерий χ2 свидетельствует только о наличии согласия или различия между двумя распределениями, но не может служить их мерой. Для определения силы связи следует использовать методы корреляционного, регрессионного или дисперсионного анализов.

19

Вычисление критерия χ2 по четырехпольной таблице.

Например, необходимо решить, оказывает ли действие удобрение на рост растений. Для этого показатели удобренных и неудобренных растений нужно занести в таблицу (табл. 4.2).

Таблица 4.2

Рост удобренных и неудобренных растений

Группа растений

Не усилили рост

 

Усилили рост

 

Неудобренные

30

 

78

 

Удобренные

 

 

15

 

120

 

Если эти цифры заменить буквами, то таблица будет

иметь вид:

 

 

 

 

 

 

а

 

b

 

В полученной таблице: а — не усили-

c

 

d

 

ли рост неудобренные, b — усилили рост

 

 

 

 

неудобренные, с

— не усилили рост

удобренные, d — усилили рост удобренные. По количеству колонок для цифр, а их четыре, такая таблица носит название четырехпольной. Вычисление критерия в этом случае производят по формуле:

χ2 = , где n = a + b + c +

d.

Эта формула действительна, если ни один из показателей а, Ь, с, d не меньше 4. В тех случаях, когда один или несколько из этих показателей меньше 4, в данную формулу вводят поправку Йейтса. Формула приобретает следующий вид:

χ2 = .

После получения показателя критерия χ2 нужно оценить его. Для четырехпольной таблицы число степеней свободы С = 1. Таким образом, при достоверном различии (Р<0,05) критерий χ2 должен быть больше 3,84, а при Р<0,01, этот показа-

20