2.6.        Регрессионный анализ результатов однофакторного эксперимента

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

1) подбор вида уравнения регрессии. Осуществляется с помощью, например, функциональных шкал;

2) определение постоянных коэффициентов в уравнении регрессии. Осуществляется с помощью МНК;

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

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

1) По результатам повторных изменений в каждом опыте вычисляются дисперсия функции отклика (σу2), среднеквадратичная погрешность (), доверительный интервал 9), математическое ожидание ().

2) Осуществляется проверка однородности дисперсии , где j = 1,2,…, N – число опытов. Дисперсии  называются однородными, если при неограниченном увеличении числа повторных измерений п в каждом из опытов они стремятся к общему пределу, т.е.

.

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

Отношение / сравнивают с табличным значением критерия Фишера (F):

· если /< F – дисперсии  и  однородны, следовательно, однородны и все дисперсии, лежащие в промежутке <<;

· если /> F – дисперсии неоднородны.

Значения критерия Фишера берутся из таблицы 2.3, составленной при доверительной вероятности Р = 0,95, для числа степеней свободы числителя (дисперсии ) и числа степеней свободы знаменателя (дисперсии ).

Таблица 2.3

fзнам

fчисл

1

2

3

4

5

6

12

24

1

2

3

4

5

6

7

8

9

10

12

24

164,4

18,5

10,1

7,7

6,6

6,0

5,6

5,3

5,1

5,0

4,8

4,3

3,8

199,5

19,2

9,6

6,9

5,8

5,1

4,7

4,5

4,3

4,7

3,9

3,4

3,0

215,7

19,2

9,3

6,6

5,4

4,8

4,4

4,1

3,9

3,7

3,5

3,0

2,6

224,6

19,3

9,1

6,4

5,2

4,5

4,1

3,8

3,6

3,5

3,3

2,8

2,4

230,2

19,3

9,0

6,3

5,1

4,4

4,0

3,7

3,5

3,3

3,1

2,6

2,2

234,9

19,3

8,9

6,2

5,0

4,3

3,9

3,6

3,4

3,2

3,0

2,5

2,1

244,9

19,4

8,7

5,9

4,7

4,0

3,6

3,3

3,1

2,9

2,7

2,2

1,8

249,0

19,5

8,6

5,8

4,5

3,8

3,4

3,1

2,9

2,7

2,5

2,0

1,5

254,3

19,5

8,5

5,6

4,4

3,7

3,2

2,9

2,7

2,5

2,3

1,7

1,0

Если в каждом из опытов число повторных измерений одинаково и равно п, то fчисл = fзнам = п – 1. Если дисперсии  и неоднородны, то это означает, что число повторных измерений недостаточно или что среди результатов повторных измерений опыта с дисперсией находится промах.

Следует очистить результаты повторных измерений от промахов, увеличить число повторных измерений, снова определить  и проверить их однородность. Если все дисперсии  и однородны, можно переходить к следующему действию.

3) Вычисляется дисперсия воспроизводимости (S2воспр). Дисперсия воспроизводимости – это среднее из дисперсий () всех опытов. Дисперсия воспроизводимости характеризует средний разброс результатов повторных измерений во всех опытах относительно своих математических ожиданий.

Если в каждом опыте число повторных измерений одинаково и равно п, то S2воспр определяется по формуле:

S2воспр =

где N – число опытов,  – математическое ожидание в j-м опыте. Если в опытах число повторных измерений различно, то S2воспр определяется как средневзвешенная величина

S2воспр =

где - число повторных измерений в j-м опыте.

Число степеней свободы дисперсии воспроизводимости равно сумме чисел степеней свободы дисперсий опытов:

.

4) Вычисляется дисперсия адекватности (S2ад). Дисперсия адекватности – это сумма квадратов отклонений расчётных и экспериментальных значений функции отклика в каждом опыте, отнесённая к числу степеней свободы. Она характеризует разброс экспериментальных результатов относительно расчётных и определяется по формуле:

S2ад =

где т – число постоянных коэффициентов в уравнении регрессии;  — математическое ожидание в j-м опыте; N – m = fад – число степеней свободы дисперсии адекватности;

3) Проверяется однородность дисперсий адекватности и воспроизводимости. Дисперсии адекватности и воспроизводимости будут однородными, если их отношение будет меньше табличного значения критерия Фишера:

S2ад / S2воспр < F,

где значения F берутся (см. табл. 2.3)для числа степеней свободы числителя (дисперсии адекватности) и числа степеней свободы знаменателя (дисперсия воспроизводимости).

Если дисперсии S2ад и S2воспр однородны , то с доверительной вероятностью Р = 0,95 можно утверждать, что составленное уравнение регрессии адекватно. Это следует из того, что при и  разброс экспериментальных значений функции отклика относительно её расчётных значений равен среднему разбросу результатов повторных измерений в каждом опыте относительно своих математических ожиданий.

Если S2ад / S2воспр > F, то выбранное уравнение регрессии неадекватно. Следует перейти к уравнению регрессии более высокого порядка или выбрать уравнение регрессии другого вида, определить значения коэффициента и снова проверить адекватность.

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