Математико-статистический анализ данных по результатам морфофизиологического исследования ерша Лахтинской губы Онежского озера
Автор: Клюкина Елена Александровна, Дзюбук Ирина Михайловна
Журнал: Ученые записки Петрозаводского государственного университета @uchzap-petrsu
Рубрика: Физико-математические науки
Статья в выпуске: 6 (111), 2010 года.
Бесплатный доступ
Анализ данных, морфофизиологические параметры, ерш, онежское озеро, выборка, сравнение двух независимых выборок, компонентный анализ
Короткий адрес: https://sciup.org/14749772
IDR: 14749772
Текст статьи Математико-статистический анализ данных по результатам морфофизиологического исследования ерша Лахтинской губы Онежского озера
С 1999 года сотрудниками лаборатории экологических проблем Севера Петрозаводского государственного университета на базе Шелтозерской учебно-производственной станции проводятся ихтиологические исследования в рамках работ по мониторингу юго-западного района Онежского озера (Лахтинская губа), целью которых является изучение динамики состояния экосистемы, выявление возможностей использования и условий сохранения ее биоресурсов. В ходе исследований проводятся сбор и обработка материалов по различным видам рыб. В уловах в больших количествах встречается ерш, что дало возможность собрать репрезентативный материал.
Цель настоящего исследования – провести математико-статистический анализ данных по результатам морфофизиологического исследования ерша Лахтинской губы Онежского озера.
МАТЕРИАЛ И МЕТОДЫ ИССЛЕДОВАНИЯ
Отлов рыбы проводился в летний период (июль, 2002 год) с помощью ставных сетей. Были исследованы и проанализированы 94 особи возраста от 2+ до 8+, из них 75 самок и 19 самцов (табл. 1). Обработка ихтиологического материала осуществлялась по общепринятой методике [12].
При математико-статистической обработке результатов исследования применялся дескриптивный (описательный) анализ, методы сравнения двух независимых выборок и компонентный анализ. Оценка достоверности различий между средними арифметическими осуществлялась по критерию Стьюдента (t-тест), а между диспер- сиями - по критерию Фишера (F-тест) (при нормальном распределении признака). Нормальность распределения признака определяли по тесту Шапиро – Уилкса (W-тест). Оценка по t-, F- и W-тестам проводилась для уровня значимости α = 0,05. Все математические операции, графические построения и статистическая обработка данных проведены с использованием стандартных программных пакетов Word и Excel и компьютерной программы StatGraphics Centurion XV.
Таблица 1
Количество материала |
(n = |
94) |
||||||
Пол |
Возраст |
Всего |
||||||
2+ |
3+ |
4+ |
5+ |
6+ |
7+ |
8+ |
||
Самки |
4 |
16 |
22 |
15 |
6 |
6 |
6 |
75 |
Самцы |
3 |
4 |
9 |
3 |
0 |
0 |
0 |
19 |
Всего |
7 |
20 |
31 |
18 |
6 |
6 |
6 |
94 |
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
Выборка состоит из самок и самцов. В данном случае мы имеем дело с альтернативным признаком (самка - самец). Таким образом, из 94 проб 75 содержат самок (значение 1), а 19 -не содержат (значение 0). Таким образом, выборка состоит из 75 единиц и 19 нулей. Доля вариант со значением 1 (доля самок) составляет:
p =
= 0,8
(то есть 80 %), тогда доля самцов составляет q = 0,2 (то есть 20 %). Отметим, что доля самок сильно отличается от доли самцов (в 4 раза), поэтому наиболее точные границы доверительного
интервала для альтернативного признака «пол ерша» даст φ-преобразование Фишера [7]. Итак, найдем границы доли самок ерша p = 0,8 при уровне значимости α = 0,05:
Флев. = 2,214 -1,96 ■ = 2,012, фправ..
= 2,214 + 1,96 ■ = 2,416, тогда 94
Р лев . (2,012) = 71,4%, Р прав . (2,416) = 87,4%.
То есть доля самок в генеральной совокупности (популяции ершей) составляет минимум 71,4 %, а максимум - 87,4 %.
Самки ерша были представлены возрастными группами от 2+ до 8+, самцы – от 2+ до 5+. Наибольшее количество особей как самок, так и самцов было возраста 4+ (пятилетки) (рис. 1). В отдельных возрастных группах соотношение полов было сдвинуто в сторону самок (рис. 2).
Размеры ерша были в пределах 20,4–82,0 г и 13,9–18,7 см (табл. 2), что значительно превышает размеры ерша из других районов Онежского озера. Например, размеры ерша Кондопожской губы (2+…7+) составляют 16,6–30,6 г и 11,7–14,3 см, а ерша района Кузаранды (1+…5+) – 7,9–14,6 г и 9,8–11,9 см [9]. Известно, что в большинстве карельских водоемов обитает мелкий ерш – 5–10 см и до 15 г, но в некоторых озерах он может достигать веса 50 и даже 100 г [5]. Ерш – медленнорастущий вид [6], [11]. Однако для ерша Лахтинской губы показатели роста (прирост длины, массы) значительны. Так, величины абсолютного прироста ерша Лахтинской губы находятся в диапазоне от 1,6 до 24,3 г и от 0,2 до 2,2 см. Максимальный абсолютный прирост массы и длины ерша характерен для возрастного периода 7+…8+ - 24,3 г, 2,2 см. Величина относительного прироста массы ерша наи- большая в возрасте 5+…6+ - 1,5 г, относительный прирост длины для всех возрастных периодов у ерша находится в пределах 1,0–1,1.
Потенциальные возможности роста ерша велики, и их реализация возможна в условиях хорошей обеспеченности пищей. Вероятно, что в условиях Лахтинской губы сложились благоприятные кормовые условия для ерша. Это подтверждают и результаты по упитанности ерша.
Таблица 2
Размерно-весовые показатели ерша Лахтинской губы (июль, 2002 год) (n = 94)
Параметры |
Возраст |
|||||||
2+ |
3+ |
4+ |
5+ |
6+ |
7+ |
8+ |
||
Масса, г |
♀♂ |
20,4 ± 3,8 |
29,1 ± 3,1 |
30,7 ± 1,4 |
35,6 ± 1,5 |
52,2 ± 2,4 |
57,7 ± 4,2 |
82,0 ± 6,0 |
♀ |
23,8 ± 6,2 |
28,7 ± 3,8 |
31,8 ± 1,7 |
36,6 ± 1,6 |
52,2 ± 2,4 |
57,2 ± 4,2 |
82,0 ± 6,0 |
|
♂ |
15,8 ± 2,2 |
30,5 ± 2,1 |
28,1 ± 2,2 |
30,3 ± 2,0 |
– |
– |
– |
|
Длина тела, см |
♀♂ |
13,9 ± 1,2 |
14,1 ± 0,5 |
14,3 ± 0,3 |
15,3 ± 0,3 |
16,0 ± 0,2 |
16,5 ± 1,0 |
18,7 ± 0,4 |
♀ |
13,4 ± 1,4 |
13,9 ± 0,7 |
14,6 ± 0,3 |
15,4 ± 0,3 |
16,0 ± 1,2 |
16,5 ± 1,0 |
18,7 ± 0,4 |
|
♂ |
14,6 ± 2,5 |
14,8 ± 0,2 |
13,5 ± 0,5 |
14,8 ± 0,6 |
– |
– |
– |
|
Упитанность |
♀♂ |
1,1 ± 0,2 |
1,3 ± 0,1 |
1,4 ± 0,1 |
1,2 ± 0,1 |
2,1 ± 0,7 |
1,9 ± 0,7 |
1,5 ± 0,1 |
♀ |
1,3 ± 0,1 |
1,3 ± 0,1 |
1,3 ± 0,1 |
1,3 ± 0,1 |
2,1 ± 0,7 |
1,9 ± 0,7 |
1,5 ± 0,1 |
|
♂ |
0,9 ± 0,4 |
1,3 ± 0,2 |
1,5 ± 0,1 |
1,1 ± 0,1 |
– |
– |
– |
Примечание. ♀♂ – общее, ♀ – самки, ♂ – самцы.

Рис. 1. Возрастной состав ерша Лахтинской губы Онежского озера (июль, 2002 год)

Рис. 2. Соотношение самок и самцов ерша Лахтинской губы Онежского озера (июль, 2002 год)
Упитанность (u) рассчитывали по следующей формуле:
100m u= , d3
где m - масса тела в г, d - длина тела без хвостового плавника в см. Величина u была в пределах 1,1–2,1, наиболее упитанными были рыбы (самки) возраста 6+…7+, наименее упитанными – возраста 2+ (табл. 2). В возрасте 2+ и старше ерш переходит от питания бентосом к хищничеству. Поэтому смена пищевых объектов, возможно, отражается на его упитанности.
В ходе анализа была проведена проверка нормального характера распределения ершей Лахтинской губы Онежского озера для уровня значимости α = 0,05 (для 94 особей) с помощью W-теста по следующим признакам: возраст, масса, длина тела, длина тела без хвостового плавника и упитанность. Результаты проверок представлены в табл. 3. Заметим, что только для длины тела (D) извлеченное из теста значение p > 0,05 (p = 0,150215), следовательно, гипотеза о нормальном распределении длины тела ершей Лахтинской губы принимается с надежностью, равной 95 %. Для остальных признаков с 5 % вероятностью ошибки гипотеза о нормальном распределении отвергается, причем для упитан- ности значение p настолько мало, что оно обозначено как 0,00.
На диаграммах нормального распределения (также называемых Quantile-Quantile Plot) видно, как близко заданное распределение приближается к нормальному (рис. 3). На диаграммах каждое наблюдаемое значение (ось ординат) сравнивается со значением, ожидаемым при нормальном распределении (ось абцисс). Также оценка соответствия данного распределения нормальному распределению представлена на гистограммах (рис. 4), тонкой линией показана кривая нормального распределения.
Таблица 3
Тесты на нормальное распределение
Тест |
Статистика |
p (P-Value) |
Шапиро – Уилкс (W-тест) |
Для возраста 0,88831 1 1,59334E-9 Для массы тела 0,899556 2,66951E-8 Для длины тела (D) 0,969756 1 0,150215 Для длины тела без хвоста (d) 0,960445 0,0306247 Для упитанности 0,649041 0,00 |

Рис. 3. Диаграммы нормального распределения для возраста, массы (m), длины тела D и d, упитанности (u) ерша

Рис. 4. Гистограммы соответствия нормальному распределению для возраста, массы (m), длины тела D и d, упитанности (u) ерша
Таблица 4
Общая статистическая характеристика для массы, длины тела
и упитанности ерша (2+…5+) |
||||||
Статистики |
Масса (m) |
Длина тела (D) |
Упитанность (u) |
|||
самки |
самцы |
самки |
самцы |
самки |
самцы |
|
Объем |
57 |
19 |
57 |
19 |
57 |
19 |
Среднее |
31,614 |
27,0263 |
14,5404 |
14,1684 |
1,26601 |
1,28563 |
Медиана |
32,0 |
26,0 |
14,5 |
14,3 |
1,2444 |
1,2289 |
Дисперсия |
117,884 |
50,9576 |
4,16102 |
3,63228 |
0,117528 |
0,213995 |
Стандартное отклонение |
10,8574 |
7,13846 |
2,03986 |
1,90585 |
0,342824 |
0,462596 |
Коэффициент вариации |
34,3437 % |
26,413 % |
14,0289 % |
13,4514 % |
27,0791 % |
35,982 % |
Минимальный предел варьирования |
9,0 |
12,5 |
9,5 |
9,6 |
0,4395 |
0,3662 |
Максимальный предел варьирования |
68,0 |
45,0 |
19,3 |
17,3 |
2,1478 |
2,3324 |
Размах выборки |
59,0 |
32,5 |
9,8 |
7,7 |
1,7083 |
1,9662 |
Стандартизированная асимметрия |
1,44497 |
0,545528 |
-0,481319 |
-1,14782 |
1,29383 |
0,670024 |
Стандартизированный эксцесс |
2,1101 |
1,65011 |
0,895464 |
0,63802 |
0,967495 |
0,929577 |
С помощью статистического пакета StatGraphics был проведен анализ отличий между самцами и самками ерша Лахтинской губы Онежского озера по различным признакам – массе, длине тела, упитанности, а также по индексам органов. Так как самцы возраста старше 5+ в выборке не представлены, для исследования отличий между самцами и самками ерша по этим признакам рассмотрена группа особей возраста 2+…5+ (всего 76 особей). В табл. 4 приведена общая статистика обеих выборок для массы (m), длины тела (D) и упитанности (u) ерша.
Подробнее результаты табл. 4 представлены на рис. 5 с помощью диаграмм частот встречаемости разных значений двух выборок и «ящиков с усами». «Ящики» (окрашенные прямоугольники) демонстрируют доверительные интервалы для средних, «усы» (тонкие линии) - размах варьирования данных, вертикальные линии, делящие каждый ящик на две части, - медианы, а крестиками обозначены средние выборок.

Рис. 5. Диаграммы частот встречаемости разных значений двух выборок и «ящики с усами» (для m, D, u)
Наибольший интерес из статистических характеристик в рамках данного анализа представляют стандартизированные эксцесс и асимметрия, с помощью которых можно определить, относятся ли эти выборки к нормально распределенным. Так, для самцов выборочные значения стандартизированных асимметрии и эксцесса массы тела попадают в нормальный диапазон, для самок стандартизированная асимметрия попадает в этот диапазон, а стандартизированный эксцесс немного (на 0,1101) превысил правую границу нормального диапазона. При этом доверительные интервалы на уровне надежности 95 % для выборочных средних по массе для самок и самцов соответственно равны [28,7332; 34,4949] и [23,5857; 30,467].
Проверка с помощью t-теста подтвердила гипотезу о равенстве среднего значения массы тела самок и самцов на уровне значимости α = 0,05: для альтернативной гипотезы о неравенстве средних вычисленное значение вероятности p = 0,0899609 > 0,05, то есть альтернативная гипотеза отвергается на уровне значимости α = 0,05. Доверительные интервалы на уровне надежности 95 % для стандартных отклонений для самок и самцов соответственно равны [8,23998; 12,0144] и [5,39391; 10,5565]. Отметим, что при- менение критерия Стьюдента предполагает равенство выборочных дисперсий.
Проверка равенства дисперсий с помощью критерия Фишера показала, что гипотеза о равенстве дисперсий выборок по массе тела не отвергается на уровне значимости α = 0,05 (p = 0,0519313 > 0,05). Это подтверждается также тем, что 95 % доверительный интервал для отношения различий [0,992631; 4,6309] содержит 1. Заметим, что для данного признака среди самок исследуемой возрастной группы 2+…5+ одна особь возраста 3+ значительно отличается (68 г при длине тела 17,1 см) от особей этого же возраста, для которых максимальный предел варьирования равен 55. Вероятно, возникла ошибка в определении возраста (возможно, при внесении данных в среду Excel). Исключив из рассмотрения данную особь (выброс), получим следующие значения стандартизированных асимметрии и эксцесса по массе: -0,338909 и -0,120373 соответственно, находящиеся в нормальном диапазоне.
Аналогично были проанализированы отличия между самцами и самками ерша по длине тела. Поскольку вычисленные значения стандартизированных асимметрии и эксцесса для обеих выборок находятся в пределах диапазона от -2 до 2 (табл. 4), нет оснований считать, что рас- пределение выборок по длине тела существенно отличается от нормального. Доверительные интервалы на уровне надежности 95 % для выборочных средних длины тела для самок и самцов соответственно равны [13,9991; 15,0816] и [13,2498; 15,087]. Проверка с помощью t-теста подтвердила гипотезу о равенстве среднего значения по длине тела самок и самцов на уровне значимости α = 0,05: для альтернативной гипотезы о неравенстве средних вычисленное значение вероятности p = 0,486634 > 0,05, то есть альтернативная гипотеза отвергается на уровне значимости α = 0,05.
Поскольку предположение о нормальном распределении выборок (длина тела самок и самцов) не отвергается, для проверки равенства дисперсий можно воспользоваться F-тестом. Доверительные интервалы на уровне надежности 95 % для стандартных отклонений по длине тела для самок и самцов соответственно равны [1,72216; 2,50239] и [1,44009; 2,81842]. Поскольку вычисленная вероятность p = 0,778036 > 0,05, гипотеза о равенстве дисперсий (то есть гипотеза о том, что рассмотренные выборочные дисперсии оценивают одну и ту же генеральную дисперсию) не отвергается на уровне значимости α = 0,05. На это же указывает доверительный интервал для отношения различий, который простирается от 0,491544 до 2,29319. Так как этот доверительный интервал содержит значение 1, нет статистически существенного различия между стандартными отклонениями этих выборок на уровне значимости α = 0,05.
Для упитанности вычисленные значения стандартизированных асимметрии и эксцесса для обеих выборок также находятся в пределах нормального диапазона от -2 до 2 (табл. 4). На уровне значимости α = 0,05 t-тест подтвердил гипотезу о равенстве среднего значения упитанности самок и самцов (p = 0,844168 > 0,05), а F-тест - гипотезу о равенстве дисперсий обеих выборок (p = 0,0908117 > 0,05).
В качестве морфофизиологических индикаторов состояния водной среды и организмов используют показатели массы органов рыб. Для анализа был проведен расчет индексов органов (в ‰), которые позволяют устранить существующую прямую пропорциональную зависимость массы органа от массы тела (табл. 5). Индекс органа (i) определяется по следующей формуле:
i=mo , m где mo - масса органа (в мг), m - масса тела (в г).
Для ерша Лахтинской губы характерны высокие показатели относительной массы сердца (2,7–3,9 ‰). Например, индексы сердца ерша (2+…7+) из Кондопожской губы Онежского озера составляют 1,8–2,0 ‰, ерша района Куза-ранды (1+…5+) – 1,9–2,2 ‰ [9]. Высокая относительная масса сердца, вероятно, является следствием того, что в имеющихся условиях среды Лахтинской губы у ерша высокие затраты энергии на обеспечение жизнедеятельности. В ходе исследований была подтверждена закономерность уменьшения индекса сердца ерша с возрастом.
Для анализа были исследованы отличия между самцами и самками ерша возрастной группы (2+…5+) по индексам сердца. В табл. 6 представлена общая статистика обеих выборок для самок и самцов ерша.
Индексы органов (‰) ерша Лахтинской
Таблица 5 губы Онежского озера (июль, 2002 год)
Органы |
Половой состав |
Возраст |
||||||
2+ |
3+ |
4+ |
5+ |
6+ |
7+ |
8+ |
||
Сердце |
♀♂ |
3,9 ± 0,6 |
3,6 ± 0,2 |
3,0 ± 0,1 |
3,0 ± 0,2 |
2,9 ± 0,2 |
2,8 ± 0,3 |
2,7 ± 0,3 |
♀ |
3,8 ± 1,1 |
3,7 ± 0,2 |
3,0 ± 0,1 |
3,1 ± 0,2 |
2,9 ± 0,2 |
2,8 ± 0,3 |
2,7 ± 0,3 |
|
♂ |
4,0 ± 0,2 |
3,2 ± 0,4 |
3,2 ± 0,2 |
2,7 ± 0,3 |
– |
– |
– |
|
Жабры |
♀♂ |
30,8 ± 5,3 |
31,8 ± 2,4 |
29,9 ± 1,6 |
31,2 ± 2,6 |
32,2 ± 2,7 |
28,5 ± 2,1 |
21,6 ± 4,6 |
♀ |
78,3 ± 21,4 |
104,3 ± 12,6 |
109,5 ± 7,7 |
31,6 ± 3,1 |
32,2 ± 1,7 |
28,5 ± 2,1 |
21,6 ± 4,6 |
|
♂ |
26,4 ± 1,3 |
26,1 ± 8,2 |
26,2 ± 3,1 |
29,4 ± 3,5 |
– |
– |
– |
|
Печень |
♀♂ |
15,0 ± 3,3 |
17,0 ± 1,7 |
16,5 ± 1,4 |
19,4 ± 2,0 |
19,2 ± 2,0 |
23,0 ± 1,9 |
20,0 ± 4,3 |
♀ |
42,4 ± 13,3 |
51,6 ± 9,9 |
60,8 ± 6,1 |
20,0 ± 2,4 |
19,2 ± 2,0 |
23,0 ± 1,9 |
20,0 ± 4,3 |
|
♂ |
10,2 ± 1,1 |
24,4 ± 5,5 |
14,3 ± 3,2 |
16,9 ± 0,4 |
– |
– |
– |
|
Селезенка |
♀♂ |
2,4 ± 0,6 |
2,6 ± 0,2 |
2,0 ± 0,1 |
2,0 ± 0,3 |
2,3 ± 0,8 |
1,4 ± 0,4 |
1,2 ± 0,3 |
♀ |
5,0 ± 1,3 |
7,9 ± 1,0 |
7,2 ± 0,7 |
1,8 ± 0,2 |
2,3 ± 0,8 |
1,4 ± 0,4 |
1,2 ± 0,3 |
|
♂ |
2,3 ± 1,2 |
2,4 ± 0,2 |
1,9 ± 0,1 |
2,9 ± 1,4 |
– |
– |
– |
|
Желудок |
♀♂ |
11,2 ± 1,8 |
10,2 ± 0,5 |
9,2 ± 0,7 |
10,3 ± 0,7 |
11,3 ± 1,4 |
7,8 ± 0,7 |
8,9 ± 1,3 |
♀ |
26,6 ± 6,1 |
31,0 ± 3,4 |
32,9 ± 2,0 |
10,9 ± 0,7 |
11,3 ± 1,4 |
7,8 ± 0,7 |
8,9 ± 1,3 |
|
♂ |
10,1 ± 1,8 |
10,0 ± 1,1 |
7,4 ± 1,0 |
7,2 ± 1,5 |
– |
– |
– |
|
Кишечник |
♀♂ |
12,3 ± 4,5 |
11,2 ± 0,9 |
10,8 ± 1,1 |
12,1 ± 1,5 |
9,9 ± 0,8 |
10,7 ± 1,6 |
10,7 ± 1,9 |
♀ |
29,1 ± 6,0 |
35,5 ± 4,7 |
37,0 ± 3,6 |
12,8 ± 1,7 |
9,9 ± 0,8 |
10,7 ± 1,6 |
10,7 ± 1,9 |
|
♂ |
7,4 ± 1,3 |
9,0 ± 0,7 |
9,7 ± 1,5 |
8,1 ± 1,2 |
– |
– |
– |
Отметим, что вычисленные значения стандартизированных асимметрии и эксцесса для самок ершей превышают правую границу диапазона от -2 до 2, что указывает на существенные отклонения от нормального распределения. Среди самок исследуемой возрастной группы 2+^5+ у одной самки возраста 2+ индекс сердца намного больше (6,58 %о), чем у других самок изучаемой группы (2+...5+), для которых максимальный предел варьирования равен 5,19. Это указывает на то, что либо у этой особи была патология сердца, либо она отличалась высокими энергетическими затратами для обеспечения жизнедеятельно сти.
Исключив из рассмотрения данную особь, мы получили следующие значения стандартизированных асимметрии и эксцесса для самок: 0,762207 и 0,653316 соответственно, которые, как и значения данных статистик для самцов, принадлежат нормальному диапазону. При этом для уровня значимости а = 0,05 проверка с помощью t-теста подтвердила гипотезу о равенстве среднего значения по индексу сердца самок и самцов (так как р = 0,40498 > 0,05); критерий Фишера показал, что гипотеза о равенстве дисперсий выборок также подтверждается на данном уровне значимости (p = 0,387831 > 0,05).
Индекс жабр ерша был в пределах 21,631,8 %о, что близко к величине относительной массы жабр ерша Кондопожской губы (30,234,8 ^) и Кузаранды (28,1-33,1 ^) [9]. Индексы жабр имеют прямую связь с изменением уровня газообмена, а интенсивность газообмена и потребность в кислороде зависят от характера пищи и активности рыб. Не выявлена зависимость изменений этого показателя от возраста рыб, однако у ерша в возрасте 8+ был наименьший индекс жабр по сравнению с рыбами младших возрастов.
Были исследованы отличия между самцами и самками ерша возрастной группы 2+^5+ по индексам жабр. Для обеих выборок значения стандартизированных асимметрии и эксцесса попадают в нормальный диапазон (табл. 6). При этом доверительные интервалы на уровне надежности 95 % для выборочных средних для самок и самцов соответственно равны [29,4207; 34,8688] и [22,1326; 31,2674]. Проверка с помощью t-теста не подтвердила гипотезу о равенстве среднего значения по индексу жабр самок и самцов на уровне значимости а = 0,05 (р = 0,0450151 < 0,05). Доверительные интервалы на уровне надежности 95 % для стандартных отклонений для самок и самцов соответственно равны [8,66739; 12,5942] и [7,16031; 14,0136]. Проверка равенства дисперсий с помощью критерия Фишера показала, что гипотеза о равенстве дисперсий выборок не отвергается на уровне значимости а = 0,05 (р = 0,731504 > 0,05), что подтверждается также тем, что 95 % доверительный интервал для отношения различий [0,503623; 2,34954] содержит 1.
Изменения массы печени за счет накопления или расходования углеводов и отчасти белков и жиров позволяют судить о характере метаболизма рыб. Так, индекс печени ерша Лахтинской губы колеблется в широком диапазоне: 15,023,0 %о. Для сравнения: индекс печени ерша Кондопожской губы составляет 21,7-26,8 %о, района Кузаранды - 11,8-14,2 %о [9]. Наименьший индекс печени был выявлен у ерша в возрасте 2+. С возраста 3+ происходит увеличение этого показателя, что, вероятно, связано со сменой пищевых объектов (переход к хищничеству) и с изменением метаболизма.
Таблица 6
Общая статистическая характеристика для индексов органов ерша (2 + ...5 + )
Статистики |
Индекс сердца |
Индекс жабр |
Индекс печени |
Индекс селезенки |
||||
самки |
самцы |
самки |
самцы |
самки |
самцы |
самки |
самцы |
|
Объем |
57 |
19 |
57 |
19 |
57 |
19 |
57 |
19 |
Среднее |
3,24614 |
3,26263 |
32,1447 |
26,7 |
17,5647 |
16,1663 |
2,18175 |
2,22789 |
Медиана |
3,2 |
3,4 |
32,35 |
28,0 |
17,59 |
14,0 |
2,1 |
2,0 |
Дисперсия |
0,934878 |
0,469532 |
105,397 |
89,7978 |
55,7116 |
84,3439 |
0,882325 |
1,30587 |
Стандартное отклонение |
0,966891 |
0,685224 |
10,2663 |
9,47617 |
7,46402 |
9,18389 |
0,939322 |
1,14275 |
Коэффициент вариации |
29,7859 % |
21,0022 % |
31,9378 % |
35,4913 % |
42,4944 % |
56,8088 % |
43,0535 % |
51,2927 % |
Минимальный предел варьирования |
1,37 |
2,04 |
10,0 |
3,46 |
3,2 |
5,19 |
0,53 |
0,02 |
Максимальный предел варьирования |
6,58 |
4,25 |
61,54 |
41,94 |
46,4 |
36,67 |
4,0 |
5,56 |
Размах выборки |
5,21 |
2,21 |
51,54 |
38,48 |
43,2 |
31,48 |
3,47 |
5,54 |
Стандартизированная асимметрия |
2,21536 |
-0,53324 |
1,62126 |
-1,55477 |
3,09413 |
2,08921 |
0,705979 |
2,13091 |
Стандартизированный эксцесс |
2,61973 |
-1,0336 |
1,98678 |
0,741828 |
4,36626 |
0,705396 |
-0,948447 |
3,2522 |
Также были проведены исследования отличий между самцами и самками ерша возрастной группы (2+…5+) по индексам печени. Выборочные значения стандартизированной асимметрии для самок и самцов, а также выборочное значение стандартизированного эксцесса находятся за пределами нормального диапазона, что указывает на существенные отклонения исследуемых выборок от нормального распределения (табл. 6). Однако среди самок исследуемой возрастной группы 2+…5+ индекс печени у одной самки возраста 5+ (46,4 ‰) намного больше, чем у остальных самок этой выборки. Среди самцов исследуемой возрастной группы 2+…5+ у самца возраста 4+ и у самца возраста 3+ индекс печени (36,67 и 35,48 ‰ соответственно) намного больше, чем у других самцов изучаемой группы, для которых максимальный предел варьирования равен 30,56. Это указывает на то, что у данных особей могла быть патология печени .
Исключив из рассмотрения эти особи, мы получили следующие значения стандартизированных асимметрии и эксцесса для индекса печени: 0,455248 и -0,377295 для самок, 1,67233 и 1,62625 для самцов, которые принадлежат нормальному диапазону. При этом для уровня значимости α = 0,05 проверка с помощью t-теста подтвердила гипотезу о равенстве среднего значения индекса печени самок и самцов (так как p = 0,0728045 > 0,05), а критерий Фишера показал, что гипотеза о равенстве дисперсий выборок также подтверждается на данном уровне значимости (p = 0,968315 > 0,05).
Известно, что масса селезенки по сравнению с массой других органов варьирует в бóльших пределах. Это связано с ее разнообразной функциональной деятельностью: продуцирование форменных элементов крови (эритроциты и лейкоциты), депонирование крови, образование лимфоцитов. Индекс селезенки ерша Лахтин-ской губы был в пределах 1,2–2,6 ‰. Для ерша Кондопожской губы и Кузаранды этот показатель находится в диапазоне 1,0–1,3 ‰ и 1,4– 1,5 ‰ соответственно. Вероятно, в условиях Лахтинской губы интенсивность работы селезенки у ерша очень высока. Наибольшие показатели отмечены у рыб возраста 3+ (2,6 ‰), что может быть связано с перестройкой пищевой активности и метаболизма в этом возрасте.
Исследованы отличия между самцами и самками ерша возрастной группы 2+…5+ по индексам селезенки. Отметим, что выборочные значения стандартизированных асимметрии и эксцесса для самцов, в отличие от самок, находятся за пределами нормального диапазона, что указывает на существенные отклонения данной выборки от нормального распределения (табл. 6). Однако заметим, что среди самцов исследуемой возрастной группы 2+…5+ у одного самца возраста 5+ индекс селезенки (5,56 ‰) намного больше, чем у других самцов изучаемой группы (2+…5+), для которых максимальный предел варьирования равен 3,93. Это объясняется патологией или тем, что у этих особей интенсивность работы селезенки в силу каких-то причин была больше, чем у других.
Исключив из рассмотрения данную особь, мы получили следующие значения стандартизированных асимметрии и эксцесса для индекса селезенки самцов: -0,309174 и 1,69731, которые принадлежат нормальному диапазону. При этом для уровня значимости α = 0,05 проверка с помощью t-теста подтвердила гипотезу о равенстве среднего значения по индексу селезенки самок и самцов (так как p = 0,576234 > 0,05), а критерий Фишера показал, что гипотеза о равенстве дисперсий выборок также подтверждается на данном уровне значимости (p = 0,597917 > 0,05).
Весовые показатели пищеварительного тракта (индексы желудка и кишечника) могут быть использованы в качестве характеристики обмена веществ в организме рыб, так как пищеварительный тракт перерабатывает разные пищевые компоненты, которые могут иметь разную пищевую ценность и разную степень усвояемости. Поэтому нагрузка на него может существенно колебаться, что способствует изменению его весовых параметров [13], [14].
Исследованы отличия между самцами и самками ерша возрастной группы 2+…5+ по индексам желудка. Общая статистика обеих выборок для самок и самцов ерша представлена в табл. 7.
В отличие от самцов, выборочные значения стандартизированных асимметрии и эксцесса для индекса желудка самок находятся за пределами нормального диапазона, что указывает на существенные отклонения данной выборки от нормального распределения. Однако среди самок исследуемой возрастной группы у самки возраста 2+ и особей возраста 4+ индексы желудка (19,78 и 22,50 ‰ соответственно) гораздо больше, чем у других самок изучаемой группы, для которых максимальный предел варьирования равен 17,65.
Таблица 7
Общая статистическая характеристика для индексов желудка и кишечника ерша (2+…5+)
Статистики |
Индекс желудка |
Индекс кишечника |
||
самки |
самцы |
самки |
самцы |
|
Объем |
57 |
19 |
57 |
19 |
Среднее |
10,4218 |
8,33316 |
12,1174 |
8,94211 |
Медиана |
10,0 |
8,23 |
10,6 |
9,14 |
Дисперсия |
12,4935 |
8,12126 |
44,7042 |
10,9703 |
Стандартное отклонение |
3,53462 |
2,84978 |
6,68612 |
3,31214 |
Коэффициент вариации |
33,9158 % |
34,1981 % |
55,178 % |
37,0398 % |
Минимальный предел варьирования |
4,67 |
2,89 |
3,71 |
4,27 |
Максимальный предел варьирования |
22,5 |
13,53 |
38,89 |
18,0 |
Размах выборки |
17,83 |
10,64 |
35,18 |
13,73 |
Стандартизированная асимметрия |
3,31705 |
0,161026 |
5,99432 |
1,85111 |
Стандартизированный эксцесс |
2,83249 |
-0,229915 |
7,10964 |
1,51976 |
Исключив из рассмотрения данные особи, мы получили следующие значения стандартизированных асимметрии и эксцесса для индекса желудка самок: 1,23733 и -0,117339, которые принадлежат нормальному диапазону. При этом для уровня значимости α = 0,05 проверка с помощью t-теста гипотезы о равенстве среднего значения по индексу желудка самок и самцов также не подтвердилась (p = 0,0307113 < 0,05). Критерий Фишера показал, что гипотеза о равенстве дисперсий выборок подтверждается на данном уровне значимости (p = 0,963301 > 0,05).
Результаты исследования отличия между самцами и самками ерша возрастной группы 2+…5+ по индексам кишечника показали, что для самцов выборочные значения стандартизированных асимметрии и эксцесса для индекса кишечника принадлежат нормальному диапазону от -2 до 2 (табл. 7). Для самок значения данных статистик превышают правый предел нормального диапазона более чем в 2,98 и 3,55 раза соответственно, что указывает на существенные отклонения данной выборки от нормального распределения. Таким образом, по индексу кишечника самки и самцы имеют существенные различия, поскольку не являются однородными (заданы разными законами распределения).
Также в рамках данного исследования была решена задача оценки степени различия между самками (57 особей) и самцами (19 особей) половозрелых ершей Лахтинской губы, возраст которых варьируется от 2+ до 5+, по следующим 9 морфологическим признакам: массе тела (m), длине D и d, индексам сердца, печени, селезенки, желудка, кишечника, жабр. Даная задача решалась с помощью компонентного анализа в среде статистического пакета StatGraphics.
Отметим, что компонентный анализ (метод главных компонент) как один из многомерных методов является отражением двух основных тенденций развития современной биометрии. С одной стороны, это стремление к более полному (многоплановому, многомерному) изучению действительности, что требует количественной оценки большого числа свойств исследуемых объектов. С другой стороны, это формирование все более наглядного, интегрированного, обобщенного представления об огромных массивах информации [7]. Наблюдения за эколого-биологическими феноменами дают n объектов (в нашем случае особей ерша), охарактеризованных набором из s признаков (в нашем случае это масса тела, длина тела, длина тела без хвостового плавника, масса органов и др.).
Компонентный анализ стремится сконцентрировать информацию об отличиях объектов в гораздо меньшее число расчетных индексов (k < s), называемых главными компонентами. Новые показатели отображают в себе общие причины, в силу которых группы признаков изменяются параллельно, а группы объектов оказываются сходными. В качестве таких общих факторов, например, мо- гут выступать пол, возраст, состояние здоровья. Значение каждого признака xj (j = 1, $) особи распадается на несколько частей, каждая из которых связана с действием определенной причины (a, b, c ...): Xj = xj + xbj + xj + ... Тогда сумма вкладов, например фактора a , в значения всех s признаков вычисляется по формуле: ГКа = xa 1 + xa 2 + ••• + xa$, гДе xaj (j = 1, $) -вклад фактора a в значение варианты j-го признака. Выразим эффект действия фактора a на j-й признак как долю от общего значения варианты следующим образом: xaj = aj • xj, где aj - относительный вклад данного фактора в конечное значение варианты, xj - значение варианты j-го признака. Тогда уравнение первой главной компоненты имеет вид: FKa = a1 • x1 + a2 • x2 +. „ + a$ • x$. Аналогично определяются уравнения второй главной компоненты ГКb и всех последующих. Коэффициенты (aj и т. п.), стоящие перед значениями признаков Xj (j = 1, $) в уравнениях компонент, являются аналогами коэффициентов корреляции данных признаков с общим фактором и называются факторными нагрузками [8]. Низкое значение факторной нагрузки говорит о независимости признака от данного фактора. Компонентный анализ эффективно «концентрирует» информацию: дисперсии первых компонент выражают большую долю общей изменчивости. Так, первые 2-3 компоненты накапливают в себе основную содержательную информацию об отличиях объектов, заменяя собой s исходных признаков.
Применение общего фактора «пол особи» позволило дать оценку степени различия между самками и самцами половозрелых ершей Лах-тинской губы (2+…5+) по описанным выше 9 морфологическим признакам (то есть число признаков s = 9). Результаты расчетов представлены в табл. 8 (в расчетах используются нормированные значения исходных признаков).
Таблица 8
Результаты расчетов компонентного анализа для общего фактора «пол особи»
Номер компоненты |
Значение дисперсий главных компонент (Eigenvalue) |
Относительный вес главных компонент (Percent of variance) |
Суммативный вес главных компонент (Cumulative percentage) |
1 |
3,37857 |
37,540 |
37,540 |
2 |
2,09512 |
23,279 |
60,819 |
3 |
1,26672 |
14,075 |
74,893 |
4 |
0,719251 |
7,992 |
82,885 |
5 |
0,53556 |
5,951 |
88,836 |
6 |
0,454349 |
5,048 |
93,884 |
7 |
0,289387 |
3,215 |
97,100 |
8 |
0,237479 |
2,639 |
99,738 |
9 |
0,0235633 |
0,262 |
100,000 |
Заметим, что первая компонента учла в себе большую долю исходной изменчивости (ее дисперсия S ГК К 1 = 3,37857 из 9, относительный вес составляет 37,54 %). Информативность второй компоненты выражается дисперсией S ГКК 2 = 2,09512 и относительным весом, равным 23,279 %. Информативность третьей компоненты выражается дисперсией S ГК к 3 = 1,09512 и относительным весом, равным 14,075 %. Дисперсии остальных компонент меньше единицы, что указывает на их низкую информативность, поскольку по умолчанию минимальное значение дисперсии главной компоненты (Minimum Eigenvalue) принимается равным 1. График Scree Plot отражает изменение дисперсий компонент (сплошная линия) и минимальный уровень компонент (пунктир) (рис. 6).
Значения факторных нагрузок для первых трех главных компонент рассматриваемых морфологических признаков представлены в табл. 9.

Рис. 6. Изменение дисперсий компонент для общего фактора «пол особи»
Таблица 9
Таблица значений факторных нагрузок в главных компонентах
Признаки |
Компонента 1 |
Компонента 2 |
Компонента 3 |
Масса тела (m) (г) |
0,412956 |
0,223743 |
-0,0206958 |
Длина тела (D) (см) |
0,402107 |
0,370433 |
-0,233382 |
Длина тела без хвостового плавника (d) (см) |
0,415668 |
0,359758 |
-0,204312 |
Индекс сердца (m1) (мг) |
-0,394719 |
0,0904513 |
-0,460582 |
Индекс печени (m2) (мг) |
-0,01825 |
0,494377 |
0,270413 |
Индекс селезенки (m3) (мг) |
-0,31011 |
0,117718 |
-0,550079 |
Индекс желудка (m4) (мг) |
-0,315346 |
0,393559 |
0,279557 |
Индекс кишечника (m5) (мг) |
-0,30411 |
0,306402 |
0,419827 |
Индекс жабр (m7) (мг) |
-0,225314 |
0,409968 |
-0,2472 |
Уравнение расчета значений первой компоненты следующее:
Component 1 = 0,413 • m + 0,402 • D + 0,416 • d -- 0,395 • m 1 - 0,018 • m 2 + 0,31 • m 3 -- 0,315 • m 4 - 0,304 • m 5 - 0,225 • m 6.
Все факторные нагрузки первых трех признаков положительны и больше 0,3, остальные отрицательны, причем для индексов печени и жабр факторные нагрузки по модулю меньше 0,3, то есть их влияние наименее значимо. Первую компоненту можно назвать «Уменьшение индексов сердца, селезенки, желудка и кишечника при увеличении массы и размеров ерша». Для второй компоненты все факторные нагрузки положительны, причем для массы тела и индексов сердца и селезенки их величина меньше 0,3. Поэтому вторую компоненту можно назвать «Увеличение индексов печени, жабр, желудка и кишечника при увеличении размеров тела». Третью компоненту можно назвать «Уменьшение индекса селезенки и индекса сердца при увеличении индекса кишечника», поскольку она объясняется в основном за счет переменных «индекс селезенки», «индекс сердца» и «индекс кишечника», их факторные нагрузки равны 0,55, -0,461 и 0,42 соответственно.
Рассмотрим численные характеристики особей, выраженные в трех главных компонентах (табл. 10). Для всех главных компонент и самцов (f), и самок (с) эти значения принимают как положительные, так и отрицательные значения, то есть для этих компонент знак численной характеристики особи не зависит от пола.
ЗАКЛЮЧЕНИЕ
В условиях Лахтинской губы Онежского озера ерш имеет максимальные размеры, высокие приросты и упитанность, а также максимальные индексы сердца и пищеварительного тракта, широкий диапазон колебаний индекса селезенки и печени по сравнению с этими показателями у ерша в других районах Онежского озера (Кондопожская губа, Кузаранда). С большой вероятностью можно считать, что в 2002 году в Лахтин-ской губе сложились наиболее благоприятные пищевые условия для ерша, что, возможно, связано с работой форелевого хозяйства. По результатам отлова и опросным данным известно, что в период функционирования хозяйства в губе качественный и количественный состав рыб в уловах был высоким. Появился лещ, не встречавшийся в этом районе до организации хозяйства. В больших количествах стали вылавливать ряпушку, плотву и другие виды. Известно, что при переходе к хищничеству ерш часто охотится на рыбу, включая молодь того же вида. Он охотно поедает икру различных видов рыб (ряпушки, леща, плотвы). Сложившиеся в Лахтинской губе благоприятные кормовые условия для ерша, обусловленные работой хозяйства, отразились на его морфофизиологическом состоянии.
Для ершей возрастной группы 2+…8+ с 95 % достоверной вероятностью было выявлено, что длина его тела имеет нормальное распределение.
С 95 % достоверностью выявлено отличие самок и самцов одной возрастной группы 2+…5+ по весовым показателям пищеварительного тракта (индексам желудка и кишечника) и по индексу жабр. Для остальных морфофизиологических признаков (масса, длина тела, длина тела без хвостового плавника, упитанность, индексы сердца, печени, селезенки), за исключением особей с патоло- гией, существенных различий между самками и самцами на уровне значимости 5 % не обнаружено (при нормальном характере распределения выборок). Также для ершей возрастной группы 2+…5+ были замечены следующие закономерности: 1) уменьшение индексов сердца, селезенки, желудка и кишечника при увеличении массы и размеров ерша; 2) увеличение индексов печени, жабр, желудка и кишечника при увеличении размеров тела; 3) уменьшение индексов селезенки и сердца при увеличении индекса кишечника, которые не зависят от пола рыбы (ерша Лахтинской губы Онежского озера).
Таблица 10
Численные характеристики особей, выраженные в главных компонентах |
|||||||||
№ ряда |
Метка (c – самка, f – самец) |
Компонента 1 |
Компонента 2 |
Компонента 3 |
№ ряда |
Метка (c – самка, f – самец) |
Компонента 1 |
Компонента 2 |
Компонента 3 |
1 |
c |
-6,12308 |
1,41647 |
2,3042 |
39 |
c |
1,6869 |
-0,521527 |
0,801433 |
2 |
c |
-4,45018 |
1,58058 |
-2,11261 |
40 |
c |
1,65057 |
1,70079 |
2,3146 |
3 |
c |
2,28985 |
-2,10592 |
-0,454106 |
41 |
c |
2,07744 |
1,81038 |
-1,71901 |
4 |
c |
1,84703 |
-0,125954 |
1,44667 |
42 |
c |
0,589043 |
-1,37391 |
-0,317117 |
5 |
f |
-2,0432 |
-3,30558 |
1,5707 |
43 |
c |
1,34662 |
-1,86487 |
0,0435226 |
6 |
f |
-0,258481 |
0,599828 |
-1,92146 |
44 |
c |
1,45519 |
0,251025 |
0,0562624 |
7 |
f |
0,384944 |
-0,276337 |
-2,14482 |
45 |
c |
1,11879 |
-0,832046 |
-1,07978 |
8 |
c |
-4,33018 |
-1,00817 |
0,481406 |
46 |
c |
0,960517 |
0,0433192 |
-0,975247 |
9 |
c |
-2,27951 |
-3,06538 |
0,0450712 |
47 |
c |
2,00865 |
0,898473 |
-1,25265 |
10 |
c |
-2,76973 |
-1,88369 |
0,178356 |
48 |
c |
0,111796 |
0,696392 |
-0,491625 |
11 |
c |
-1,10599 |
-1,38561 |
0,579863 |
49 |
c |
3,09015 |
-0,753296 |
-0,0735608 |
12 |
c |
0,1922 |
0,130438 |
-1,96484 |
50 |
f |
-0,557063 |
-3,05861 |
-0,315866 |
13 |
c |
-0,876327 |
-1,71357 |
0,0678525 |
51 |
f |
-0,347011 |
-1,09079 |
-1,01857 |
14 |
c |
-0,433809 |
0,34987 |
-0,604301 |
52 |
f |
-1,26734 |
0,294449 |
0,582169 |
15 |
c |
-0,211629 |
0,831523 |
0,681087 |
53 |
f |
-0,489403 |
-0,864343 |
0,571373 |
16 |
c |
-1,97457 |
0,864856 |
-1,72525 |
54 |
f |
-0,625915 |
-0,796724 |
-0,453623 |
17 |
c |
-3,58448 |
1,11715 |
-1,20036 |
55 |
f |
0,399526 |
-0,376344 |
0,226167 |
18 |
c |
0,808437 |
1,00611 |
-0,949569 |
56 |
f |
0,0658157 |
-3,55223 |
0,366113 |
19 |
c |
-0,117693 |
-0,328563 |
-0,115229 |
57 |
f |
-0,157844 |
1,25817 |
0,675066 |
20 |
c |
-1,59095 |
1,14218 |
-0,0949976 |
58 |
f |
3,19381 |
-2,34388 |
0,0729681 |
21 |
c |
0,29709 |
-0,737892 |
-1,10293 |
59 |
c |
-0,196112 |
-0,504879 |
0,799256 |
22 |
c |
2,75963 |
2,47286 |
-1,58159 |
60 |
c |
0,0431036 |
-0,994454 |
-0,359415 |
23 |
c |
4,49462 |
0,704298 |
0,53886 |
61 |
c |
0,559504 |
0,361771 |
-0,124298 |
24 |
f |
1,04563 |
-0,828179 |
0,707224 |
62 |
c |
1,04332 |
1,36969 |
-1,14579 |
25 |
f |
-0,722265 |
-0,00961921 |
-0,736094 |
63 |
c |
0,315808 |
2,05829 |
1,38357 |
26 |
f |
-0,669869 |
1,26956 |
-0,273727 |
64 |
c |
-0,182978 |
0,571661 |
0,454876 |
27 |
f |
1,35892 |
0,560471 |
0,66861 |
65 |
c |
0,779205 |
0,444373 |
0,138653 |
28 |
c |
-1,99887 |
-0,469097 |
2,18703 |
66 |
c |
2,02667 |
-0,153904 |
1,40277 |
29 |
c |
-3,26901 |
2,36922 |
1,35051 |
67 |
c |
0,441921 |
1,6244 |
-0,497019 |
30 |
c |
-1,00794 |
0,33046 |
0,290611 |
68 |
c |
-2,84282 |
1,94485 |
-1,91719 |
31 |
c |
-1,23955 |
-1,18104 |
0,101608 |
69 |
c |
0,42474 |
1,25727 |
2,85178 |
32 |
c |
-1,79385 |
1,78507 |
-0,755595 |
70 |
c |
2,0433 |
0,064037 |
0,44329 |
33 |
c |
-0,0532005 |
-0,959125 |
0,337348 |
71 |
c |
0,913886 |
-0,685667 |
-0,645632 |
34 |
c |
-0,213046 |
-0,577952 |
0,222258 |
72 |
c |
2,41815 |
-0,0865143 |
0,359717 |
35 |
c |
0,266583 |
-0,841282 |
1,52479 |
73 |
c |
1,97936 |
4,88588 |
1,99735 |
36 |
c |
-0,761346 |
1,38452 |
1,22434 |
74 |
f |
-0,578833 |
-0,783703 |
-2,26556 |
37 |
c |
0,199767 |
0,890789 |
-0,611525 |
75 |
f |
0,606633 |
0,0867251 |
0,534439 |
38 |
c |
0,0990222 |
-0,559134 |
0,512173 |
76 |
f |
1,72993 |
-0,428431 |
-0,094992 |
Список литературы Математико-статистический анализ данных по результатам морфофизиологического исследования ерша Лахтинской губы Онежского озера
- Айвазян С. А., Енюков И. С., Мешалкин Л. Д. Основы моделирования и первичная обработка данных: Справочное издание. М.: Финансы и статистика, 1983. 471 с.
- Айвазян С. А., Енюков И. С., Мешалкин Л. Д. Прикладная статистика: Классификация и снижение размерности. М.: Финансы и статистика, 1989. 607 с.
- Боровиков В. П., Боровиков И. Д. Statistica. Статистический анализ и обработка данных в среде Windows. М.: Информационно-издательский дом «Филин», 1997. 608 с.
- Гублер Е. В., Генкин А. А. Применение непараметрических критериев статистики в медико-биологических исследованиях. Л.: Медицина, 1973. 142 с.
- Ивантер Д. Э., Рыжков Л. П. Рыбы. Петрозаводск: Изд-во ПетрГУ, 2004. 176 с.
- Ковалев П. М. Биологические особенности окуня и ерша озер Ильмень и Псковско-Чудского//Сб. науч. тр. ГосНИОРХ. Вып. № 236. Л., 1985. С. 117-127.
- Коросов А. В. Специальные методы биометрии: Учеб. пособие. Петрозаводск: Изд-во ПетрГУ, 2007. 364 с.
- Коросов А. В., Горбач В. В. Компьютерная обработка биологических данных: Метод. пособие. Петрозаводск: Изд-во ПетрГУ, 2010. 84 с.
- Крупень И. М. Морфофизиологические особенности ерша Gymnocephalus cernuus L. из различных условий обитания: Дис. … канд. биол. наук. Петрозаводск, 1999. 190 с.
- Малета Ю. С., Тарасов В. В. Непараметрические методы статистического анализа в биологии и медици-не. М.: Изд-во МГУ, 1982. 178 с.
- Мельянцев В. Г. Рыбы. Петрозаводск: Карелия, 1974. 120 с.
- Правдин И. Ф. Руководство по изучению рыб. М.: Пищ. пром-сть, 1966. 376 с.
- Смирнов В. С., Божко А. М., Рыжков Л. П., Добринская Л. А. Применение метода морфофизиологических индикаторов в экологии рыб//Тр. Т. 7. Петрозаводск, 1972. 190 с.
- Строганов Н. С. Экологическая физиология рыб. М.: Изд-во АН СССР, 1962. 444 с.