Главная страница

Элементарная биометрия. Книга служит элементарным пособием для практического применения вариационной статистики в биологических исследованиях


Скачать 3.04 Mb.
НазваниеКнига служит элементарным пособием для практического применения вариационной статистики в биологических исследованиях
АнкорЭлементарная биометрия.doc
Дата04.01.2018
Размер3.04 Mb.
Формат файлаdoc
Имя файлаЭлементарная биометрия.doc
ТипКнига
#13667
страница20 из 32
1   ...   16   17   18   19   20   21   22   23   ...   32

Регрессионный анализ


Коэффициент корреляции указывает лишь на степень (тес­ноту) связи в изменчивости двух переменных величин, но не позволяет судить о том, как меняется одна величина по мере изменения другой. Ответ на этот вопрос дает вычисление коэф­фициента регрессии, показывающего, на какую величину в среднем изменяется один признак при изменении другого на единицу измерения. Регрессионный анализ, в отличие от корреляционного, изучает эффект влияния одного признака на другой, зависимость признака от фактора, характер влияния фактора на признак. Его основные результаты таковы:

1. Таблица дисперсионного анализа, в которой показана сила и достоверность влияния на признак изучаемого фактора или другого признака.

2. Уравнение регрессии, выражающее пропорциональность сопряженного изменения признаков, тенденции их взаимосвязанной изменчивости или динамики.

3. Оценки значимости коэффициентов уравнения регрессии.

Регрессионный анализ методически ориентирован односторонне – на изучение зависимости одного признака от другого (зависимость y от x или, напротив, зависимость xотy), хотя может применяться к случаям, когда фактически имеется взаимозависимость двух переменных.

Основную тенденцию взаимосвязанного изменения двух признаков можно отобразить с помощью простого графического приема. Разобьем ось x на несколько интервалов. Найдем для каждого из них частные средние значения признака y (My). Теперь проведем через эти средние точки ломаную линию. Это будет линия регрессии Y по x. Регрессия изменение среднего уровня одного признака при изменении другого (рис. 12).

Линейная регрессия


К сожалению, ход ломаной линии нельзя передать простым уравнением, к тому же на нем сказывается способ интервального разбиения оси абсцисс, а также уровень репрезентативности в разных областях распределения. В этом смысле предпочтительнее единственная прямая линия регрессии, подчеркивающая основные тенденции зависимости признаков, которая может быть выражена простым уравнением линии: y = ax + b.

Судить о том, как меняется одна величина по мере изменения другой, позволяет коэффициент регрессии (a), показывающий, на какую величину в среднем изменяется один признак (y) при изменении другого (x) на единицу измерения (точнее, на какую величину один признак отклоняется от своей средней при некотором отклонении другого признака от своей средней):

 My = a ∙ (Mx).

Простые преобразования:

y = a ∙ x + My − a ∙ Mx, b = My − a ∙ Mx

и приводят к уравнению линии: y = ax + b.



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

Вычисления коэффициентов линейной регрессии y = ax + b ведутся по следующему алгоритму. Сначала найдем вспомогательные величины:

Cx = Σx² − (Σx)² / n,

Cy = Σy² − (Σy)² / n,

Cxy = Σ(x ∙ y− (Σx) ∙ (Σy) / n,

M= Σn, Mx = Σn.

Затем рассчитаем коэффициенты: a = CxyCx, b= M− a ∙ Mx.

Оценить значимость коэффициента регрессии позволяет критерий tСтьюдента, проверяющий нулевую гипотезу Но: а = 0, коэффициент регрессии значимо от нуля не отличается. С этой целью рассчитывается ошибка коэффициента регрессии ma:

, где mr – ошибка коэффициента корреляции (см. с. 62),

и вычисляется значение критерия:

t = (a − 0) / ma = a / mat(0.05, n− 2).
Смысл этого критерия состоит в следующем. Коэффициент регрессии a характеризует сопряженность пропорционального изменения двух признаков, т. е. отвечает за то, что линия регрессии имеет некоторый угол относительно оси абсцисс. Значение a = 0 означает, что линия регрессии идет параллельно оси ОХ, что при изменении признака xпризнак yне меняется, т. е. что yне зависит от x. Значения коэффициента, отличные от нуля, говорят о том, что взаимосвязь признаков имеет место, при a > 0 зависимость положительная, при a < 0 – отрицательная.

Вернемся к примеру с описанием зависимости меж­ду живым весом коров и их приплода (стр. 61). Расчеты для построения уравнения регрессии показаны в таблице 16. Сначала вычисляются квадраты вариант и их произведения, а также суммы вариант, квадратов и произведений. Вычисления ведут­ся по точным рабочим формулам. Проще всего это делать в среде Excel, с помощью команды Сервис \ Анализ данных \ Регрессия.

Таблица 16

i

у

х

у²

х²

ху

Y

(yYi

tmY

minY

maxY

1

25

352

625

123904

8800

25.6

0.31

2.0

23.6

27.5

2

26

376

676

141376

9776

27.1

1.29

1.7

25.5

28.8

3

31

402

961

161604

12462

28.8

4.65

1.4

27.4

30.2

4

32

453

1024

205208

14496

32.2

0.04

1.2

31.0

33.4

5

34

484

1156

234256

16456

34.2

0.06

1.3

32.9

35.5

6

38

528

1444

278784

20064

37.1

0.76

1.7

35.4

38.9

7

38

555

1444

308025

21090

38.9

0.81

2.1

36.8

41.0

Σ

224

3150

7330

1453158

103144




7.92











Проведем последовательные расчеты вручную. Сначала определим вспомогательные величины:

n = 7,

Cxy = Σ(xy)−(Σx)∙(Σy)/n = 1031443150224 / 7 = 2344,

Cy = Σy² − (Σy)² / n = 7330 − 224² / 7 = 162,

Cx = Σx² − (Σx)² / n = 1453158 − 3150² / 7 = 35658,

затем – параметры:

My= Σn = 224 / 7 = 32,

Mx = Σn = 3150 / 7 = 450,

= 5.2,

= 77.1,

= 0.975,

,

b = M− a∙Mx = 32 − 0.0657∙450 = 2.419.

Получено уравнение линейной регрессии Y= 0.0657x+2.419, которое позволяет рассчитать теоретические значения Y (табл. 16, графа 7).

Далее найдем ошибку коэффициента регрессии:

,

,

и, наконец, критерий t Стьюдента для проверки значимости коэффициента регрессии: ta= a / ma = 0.0657 / 0.00667 = 9.84.

Для уровня значимости α = 0.05 и числа степеней свободы df=n−2=5 находим табличное значение критерия Стьюдента t(0.05,5) = 2.57. Полученная величина (9.84) превышает табличную (2.57), что говорит о статистической значимости коэффициента регрессии (a), о достоверности его отличия от нуля. Масса тела теленка действительно возрастает вслед за ростом массы тела коровы.

Рассчитаем доверительную зону (интервал), в которой с той или иной вероятностью заключены теоретические средние значения веса новорожденных. Критерий Стьюдента (нормированное отклонение) для уровня значимости α = 0.05, и числа степеней свободы df = п − 1 = 6 составит 2.45. Далее находим границы. Так, для значения x= 352 кг прогноз по уравнению регрессии равен Y= 25.56, а возможное отклонение средней составит:

tmY =  =

= 2.45∙0.81 = 1.98.

Отсюда находим границу доверительного интервала (табл. 16):

верхнюю: maxY = YitmY = 25.56 + 1.98 = 27.54

и нижнюю: minY = Yi− tmY = 25.56 − 1.98 = 23.58.

Средняя масса новорожденного теленка для коров весом 352 кг с вероятностью P = 0.95 должна находиться в диапазоне от 23.6 до 27.5 кг (рис. 14).

Регрессионный анализ позволяет проверить значимость и второго коэффициента уравнения регрессии, свободного члена b. Математический смысл свободного члена уравнения линии состоит в том, что этому значению равна функция (y) при условии, что аргумент равен нулю (x = 0):

y = ax + b= a0+ b = b.

В рамках регрессионного анализа рассматривается именно эта гипотеза Но: b = 0, т. е. что линия регрессии проходит через начало осей координат, точку пересечения осей координат, через нуль. Если гипотеза опровергается, значит, линия регрессии не пересекает ось ординат. Если гипотеза не опровергается, мы можем считать, что между признаками существует простая пропорция (Y = ax) и расчет коэффициента регрессии a упрощается: a= Σ(xy) / Σx². Нулевая гипотеза Но: b = 0 проверяется по критерию Стьюдента: t= (b − 0/mb = b /mbt(0.05, n−2), где mb – ошибка коэффициента b.



Рис. 14. Линия регрессии Y = 0.0657∙x+ 2.1347 и ее доверительный интервал
Ошибка второго коэффициента регрессии рассчитывается в два этапа. Сначала находим общую ошибку регрессионной средней (или остаточное стандартное отклонение), которая может вычисляться по-разному.

Точная формула для небольших выборок дает величину:

.

Общая точная формула показывает практически такой же результат:

 = 1.2582

(величина Cостат.– это сумма квадратов разности между расчетными и реальными значениями признака, она найдена в табл. 16, внизу 7 графы, Cостат.= 7.92).

Теперь вычисляем ошибку коэффициента b:

3.0359

и критерий tСтьюдента: tb = b / mb = 2.419 / 3.0359 = 0.797.

Для уровня значимости α = 0.05 и числа степеней свободы dfn− 2 = 5 табличное значение составляет t(0.05, 5) = 2.57. Анализ показал, что критерий Стьюдента для свободного члена уравнения (0.797) оказался ниже табличного значения (2.57), т. е. коэффициент b значимо от нуля не отличается (при данном объеме собранных материалов). Это позволяет пересчитать коэффициент регрессии: a= Σ(xy) / Σx² = 0.071. Теперь можно пользоваться уравнением регрессии вида: Y= 0.071∙x.

Оценить достоверности взаимодействия признаков можно и с помощью дисперсионного анализа (табл. 17). В этом случае общая дисперсия зависимого признака y (Cобщ.) разлагается на две составляющие – регрессионную дисперсию (изменчивость признака y, связанная с влиянием признака x, Срегр.) и случайную, или остаточную, дисперсию (изменчивость признака y, связанная с влиянием неучтенных случайных факторов, Состат. (рис. 14, табл. 17, 18).

Общую сумму квадратов (Собщ. = Cy= Σ(y− My)2 = Σyi− (Σyi)n) находят непосредственно как сумму квадратов отличий между значением yi для каждой варианты и общей средней признака y. Остаточную сумму квадратов (Состат. = Σ(y− Yi)2) находят также непосредственно как сумму квадратов отличий между значением yi для каждой варианты и значением, предварительно рассчитанным по уравнению регрессии Yi= axi + b(для соответствующих значений xi). Модельную сумму квадратов (Смод. = Σ(Y− My)2) рассчитывают как разность между общей и остаточной (Смод. = Cобщ− Cостат.).

Рис. 15. Модель варианты в регрессионном анализе
Таблица 17

Составляющие дисперсии

Суммы квадратов, С

Формулы расчета сумм квадратов


df


S²


F

Регрессия

Срегр.

Σ(YiMy)2

Cобщ.  Cостат.

1

S2 регр.

 = 




Отклонения вариант от линии

регрессии

Состат.

 = Σ(yi− Yi)2




n− 2

S2остат.

 = 

F(0.05, 1, n−2)

Общая

(всего)

Собщ.

 = Σ(yiMy)2

yi2 − Σyi)2 n= 

 = Cy










Таблица 18

Составляющие дисперсии

С




df

S²

F

Регрессия

Срегр.

 = Σ (YiY)2

154.08

1

S2 регр.

= 154.08

F

 =  = 

 = 97.3

Отклонения

вариант от линии регрессии

Состат.

 = Σ (yiYxi)2

7.92

5

S2остат. =

= 1.58

F(0.05, 1, 5) = 6.6

Общая

(всего)

Собщ.

 = Σ (yi − Y)2

162











Показателем «силы влияния признака на признак» служит коэффициент детерминации, отношение регрессионной суммы квадратов к общей сумме квадратов (принимает значения от 0 до 1): 0.95. Между коэффициентом детерминации и коэффициентом корреляции существует простое соответствие: r = = 0.975.

Построив таблицу дисперсионного анализа с помощью критерия Фишера можно проверить нулевую гипотезу Но: предсказания регрессионной модели в целом неадекватно описывают исходные данные, зависимости между признаками нет. Конструкция критерия исследует вопрос, превышает ли варьирование, учтенное моделью, случайное (остаточное) варьирование? Критерий Фишера вычисляется как отношение модельной и остаточной дисперсии:

F = S2модS2остат. = 154.08 / 1.58 = 97.3.

Табличное значение F(0.05, 1, 5) = 6.6. Поскольку полученное значение критерия оказалось выше табличного, дисперсия реального признака y приближается по величине к дисперсии расчетных значений признака Y, т. е. существенно превышает (случайные) отличия между ними. Регрессионная модель в целом адекватно описывает исходные данные.

1   ...   16   17   18   19   20   21   22   23   ...   32


написать администратору сайта