Биостат - помощь / Учебники / Ивантер Коросов Введение в количественную биологию

210 Задача «Найти зависимость между двумя признаками»

Величина (n3 - n) = (103 -10) = 165.

Статистическая ошибка и критерий достоверности отличия коэффициента корреляции от нуля вычисляются по формулам:

Tr = rS /mr = 0.638/0.272 = 2.34.

Величина критерия равна несколько выше критического значения (2.31) для уровня значимости α = 0.05 и числа степеней свободы df = n–2 = 8 (табл. 6П). Казалось бы, это дает основание отвергнуть нулевую гипотезу rS = 0 и с вероятностью P = 95% констатировать достоверность установленной связи. Однако в связи с тем, что при небольших выборках статистические свойства коэффициента Спирмена еще «хуже», чем коэффициент Пирсона, для оценки значимости корреляции лучше воспользоваться специально подготовленной таблицей 17П, аналогичной рассмотренной выше таблице 16П. Чтобы полученный коэффициент можно было считать достоверно отличным от нуля, он должен превышать табличное значение при данном n. В нашем случае (n = 10, α = 0.05) коэффициент r = 0.638 ниже табличного r = 0.64, следовательно, значимо от нуля не отличается. Зависимость численности лисицы и грызунов по приведенным данным достоверно не прослеживается.

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

Степень сопряженности (сочетаемость) двух возможных состояний двух качественных признаков можно измерить с помощью особого коэффициента корреляции – коэффициента контингенции Шарлье.

альтернативными признаками выясняется вопрос о том, чаще ли оба признака одновременно присутствуют или отсутствуют у варианты, чем это могло бы быть по случайным причинам. Достоверность отличия от нуля оценивается по критерию Стьюдента:

самок песцов получены первичные материалы о численности родивших (+) и неродивших (– ) самок из числа хотя бы дважды перекрытых (+) и неперекрытых (–).

9

ЗАДАЧА «КЛАССИФИЦИРОВАТЬ ОБЪЕКТЫ»

Методы многомерного анализа

Методы многомерной статистики – своеобразный отклик математики на запрос современной науки обеспечить, с одной стороны, более полное (многоплановое) количественное описание биологических объектов и окружающей среды (с помощью большого числа переменных), а с другой стороны – представить огромные массивы информации в более наглядном, интегрированном, обобщенном виде. Поиск максимально полной, но интегрированной характеристики каждого объекта привел к идее рассчитывать небольшое число новых признаков, вбирающих в себя почти всю информацию от исходных характеристик; в результате полученные данные «сворачиваются» до размеров, которые в состоянии охватить мысль. Так решается «задача сокращения размерности».

Теоретической основой для методов многомерной статистики служит понятие гиперпространства, или многомерного пространства. В отличие от привычного физического трехмерного пространства, имеющего три ортогональных (взаимно перпендикулярных) оси, многомерное пространство имеет множество осей координат, в качестве которых выступают признаки (переменные) изучаемых объектов. Отдельный объект, охарактеризованный по нескольким признакам, рассматривается как отдельная точка, а множество объектов

– как облако точек. Если объекты (особи разного возраста, пола, органы, пробы, даты, разные популяции, виды, биотопы, местообитания и т. п.) отличаются друг от друга по разным признакам, то они будут занимать разное положение в многомерном пространстве; объекты оказываются рассеянными в нем.

Главной характеристикой объектов становится расстояние между ними в этом гиперпространстве, а главной особенностью всей выборки – форма облака рассеяния со своими пустотами и сгущениями объектов. Методы многомерной статистики изучают информацию, «закодированную» в порядке расположения объектов друг относительно друга. Например, взрослая особь по множеству

размерных признаков превосходит молодую. Она будет находиться в зоне особей с большими размерами, тогда как молодая – в зоне мелких. Исследование относительного места расположения особи в «облаке» других особей раскрывает, «расшифровывает» ее биологический статус.

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

Основы кластерного анализа

Классификация, кластеризация – методы, широко используемые в современной таксономии, – позволяют наглядно представить сходство или различие биологических объектов, охарактеризованных по многим параметрам. Эти подходы можно применять в самых разных областях биологии, в частности, для сравнения условий среды в сериях местообитаний, выявления различий и сходства между сообществами живых организмов, отдельными их популяциями, группами, особями и т. п. Кластерный анализ, как и многие другие многомерные статистические приемы, не имеет достаточно разработанного математического аппарата для статистического оценивания полученных данных; его основная функция – выявление скрытой структуры биологического материала, что позволяет затем целенаправленно ставить и решать конкретные биометрические задачи с помощью простых статистических методов (регрессионного, корреляционного, дисперсионного и др.).

Суть кластерного анализа состоит в

1)определении «расстояний» (меры различия) между объектами по всей совокупности признаков,

2)группировании сходных объектов в кластеры (кластеризация),

3)графическом изображении сети (или древа) расстояний между всеми объектами.

признака (относительного расстояния) на основании нескольких исходных.

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

Один из вариантов такой меры основан на коэффициенте сходства Съёренсена, который используется для сравнения многовидовых сообществ. «Расстояние» вычисляется по формуле:

2 × A ,

B + C

где В и С – число видов в двух сравниваемых сообществах, А – число общих видов.

Рассмотрим в качестве примера анализ биоценотических группировок мелких млекопитающих Приладожья. Видовой состав изучен в 7 основных биотопах: лишайниковых сосняках (A), сосня- ках-зеленомошниках (Б), ельниках (В), спелых лиственных и смешанных лесах (Г), лиственном мелколесье (Д), молодых зарастающих вырубках (Е) и по границе сеяного луга и ольшаника (Ж). Встречаемость (по принципу отсутствие – присутствие) и относительная численность (число особей на 100 ловушко-суток) 14 видов землероек и грызунов показаны в таблицах 9.1 и 9.2. Дальнейшая процедура сводится к следующему.

По данным таблицы 9.1 рассчитывается матрица расстояний между разными биотопами. Например, в ельниках (В) отмечено 12 видов мелких млекопитающих, а на вырубках (Е) – 5; из них общих для обоих биотопов – 5. Отсюда расстояние:

S = 1–10/17 = 0.41.

Смысл следующей операции – собственно кластеризации (от слова «кластер» – гроздь, группа) – состоит в последовательном объединении объектов в кластеры, в группы, внутри которых сходство между объектами выше, чем с другими объектами или кластерами. Вначале объединяются наиболее сходные объекты (с наименьшим расстоянием между собой), затем приближающиеся к ним

по этому показателю и так далее до момента слияния всех объектов в один общий кластер. При этом на промежуточных этапах могут образовываться несколько отдельных кластеров. Уровень каждого объединения фиксируется и затем отображается на графике.

Таблица 9.1

Таблица 9.2

Существует множество вариантов процедуры кластеризации, из них наиболее простым считается метод «ближайшего соседа», не требующий обязательного использования ЭВМ. Сначала по матрице расстояний (табл. 9.2) отыскиваются ближайшие соседи для всех объектов и заносятся в таблицу наименьших расстояний (табл. 9.3). Так, к лишайниковому сосняку (А) ближе всего соснякзеленомошник (Б): SАБ = 0.27, а к сосняку-зеленомошнику (Б) – лиственное мелколесье (Д): SБД = 0.07, (минимальное расстояние из всех изученных биотопов).

Заполнив четыре первые графы, приступают к построению предварительного дендроида расстояний (рис. 9.1, А). Для этого на график наносят индексы объектов и расстояния между ними, соединяют их прямыми линиями. В нашем случае сначала образовалось два отдельных кластера (АБДЕЖ и ВГ), но их может быть и больше. Теперь вновь возвращаемся к таблицам 9.2 и 9.3. В пятой графе против объектов из меньшего кластера следует отметить индекс ближайших объектов из большего кластера, а в шестой – расстояние между ними. Далее выбираем звено наименьшей протяженности – это спелые лиственные леса (Г) и молодняки (Д): SГД = 0.21. На рис. 9.1 соединим кластеры пунктирной линией, и кластеризация завершена.

Последний этап – построение окончательного варианта дендрограммы. Здесь также есть несколько возможностей. Представленное на рис. 9.1, Б «древо минимальной протяженности» строится с учетом единственного условия – соблюдения пропорций расстояний между биотопами-соседями.

Ж

АБ

ДГ

В

Рис. 9.1. А – схема поэтапной кластеризации; Б – «древо минимальной протяженности»; А– Ж – индексы биотопов

Классический вариант дендрограммы приведен на рис. 9.2. По оси абсцисс размещаются объекты в том порядке, который продиктован логикой их связей и субъективными вкусами исследователя, отдельные ветви «древа» при этом не должны пересекаться. По оси ординат откладывается расстояние между ближайшими соседями (рис. 9.2).

Интерпретация полученных результатов подчеркивает достоинства дендрограммы как емкой иллюстрации обобщающих характеристик. Так, в данных по Приладожью кластерный анализ выделил группы биотопов, наиболее близких по условиям обитания и видовому составу зверьков. Наиболее богаты видами еловые и смешанные леса (В и Г). Обедненными териокомплексами, представленными в основном политопными видами, характеризуются сосня- ки-зеленомошники, вырубки и лиственное мелколесье (Б, Е, Д). Население сосняков (Б и А) в общем сходно (табл. 9.1), но в лишайни-

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

0.4

Ж А Г В Е Б Д

Рис. 9.2. Дендрограмма сходства биотопов по видовому составу мелких млекопитающих

При использовании в кластерном анализе количественных признаков применяют евклидову меру расстояния:

лексу количественных признаков, проводят кластеризацию и построение дендрограмм по описанному выше методу. Рассмотрим эту процедуру на уже знакомом примере многовидовых группировок мелких млекопитающих в трех типах биотопов Приладожья (Б – сосняки-зеленомошники, В – ельники, Д – мелколесье), но по данным количественных учетов канавками (экз. на 10 канавко-суток; табл. 9.4).