- •4 Обработка данных методами регрессионного анализа
- •4.1 Теоретическое введение
- •4.1.1 Оценка коэффициентов регрессии
- •4.1.2 Построение доверительных интервалов для коэффициентов регрессии
- •4.1.3 Проверка гипотезы об адекватности регрессионной модели
- •4.2 Содержание типового расчета
- •4.3 Пример выполнения типового расчета
- •4.3.1 Первичная обработка результатов экспериментов
- •4.3.3 Расчет линейной и квадратичной регрессионных моделей
- •4.3.2 Подготовка данных для расчета моделей регрессии. Построение ортогональных многочленов
- •4.3.4 Графический анализ результатов расчета
- •4.3.5 Проверка адекватности регрессионных моделей и принятие решения о выборе модели регрессии
- •4.3.6 Построение доверительных интервалов
- •4.3.7 Выводы по результатам типового расчета
- •Литература
- •5 Обработка данных методами линейного корреляционного анализа
- •5.1 Теоретическое введение
- •5.1.1 Двумерный случайный вектор. Линейная корреляция
- •5.2 Содержание типового расчета
- •5.3 Порядок выполнения типового расчета. Примеры
- •5.4 Оформление отчета
- •Литература
4.1.2 Построение доверительных интервалов для коэффициентов регрессии
Оценки параметров регрессии , определяемые формулами (4.5) и (4.11), являются точечными оценками истинных значений параметров Bj . Если результаты экспериментов независимы и подчиняются нормальному закону распределения с дисперсией σ2, то доверительный интервал с доверительной вероятностью P = 1 – α для каждого параметра Bjможно определить неравенством
| Bj – | < εj или Bj = ± εj , |
(4.15) |
где
εj = t1-α
/2(k) SB ;
SB = S / ||Tj||;
||Tj|| = |
(4.16) |
здесь число ||Tj|| называется нормой базисной функции Tj ; S 2 – несмещенная оценка дисперсии σ2 с числом степеней свободы k, t1-α /2(k) – квантиль распределения Стьюдента; Wi = 1 в случае равноточных измерений.
4.1.3 Проверка гипотезы об адекватности регрессионной модели
Регрессионная модель называется адекватной, если предсказанные по ней значения переменной Y согласуются с результатами эксперимента. Если модель адекватна, то отклонения результатов эксперимента от полученной функции регрессии ΔYi = Yi – (xi) являются реализациями случайных ошибок эксперимента Zi, которые, в силу предположений (4.1), (4.9) должны быть независимыми нормально распределенными случайными величинами с нулевыми средними и дисперсиями, равными σ2 / Wi. Проверка выполнения этих предположений осуществляются статистическими методами и лежит в основе оценки адекватности модели регрессии. Для проверки адекватности регрессионной модели вычисляют остаточную дисперсию (так называемую дисперсию адекватности) по формуле
|
(4.17) |
где ΔYi – отклонения средних Yi от проверяемой модели регрессии; kад – число степеней свободы дисперсии адекватности; n – число точек, в которых проводился эксперимент; m – число оцениваемых параметров Bj в проверяемой модели. Если истинная функция регрессии имеет тот же вид, что и рассматриваемая модель (например, так же, как и модель, представляет собой квадратичную функцию), то дисперсия адекватности служит несмещенной оценкой истинной дисперсии эксперимента и ее можно сравнивать с другими подобными оценками. В частности, может быть проведена независимая серия измерений для получения оценки дисперсии эксперимента S2экс. Если в каждой точке проводилось несколько измерений, то за независимую оценку можно взять сводную оценку дисперсии S2св :
|
(4.18) |
с
числом степеней свободы
.
Здесь L число
серий измерений, kj = nj -
1, ..., nj -
число измерений в каждой
серии.
Следовательно, S2экс оценивает
дисперсию эксперимента Dэкс,
S2ад характеризует
степень отклонения экспериментальных
точек от регрессионной модели, т.е.
оценивает некую дисперсию адекватности Dад.
Проверка адекватности модели заключается
в проверке гипотезы H0:
Dад = Dэкс при
альтернативной гипотезе H1:
Dад > Dэкс (если
модель неадекватна, отклонения
экспериментальных точек от модели будут
больше погрешностей эксперимента).
Таким образом, задача сводится к проверке
гипотезы о равенстве дисперсий, которая
решается с помощью критерия Фишера.
Вычисляем отношение
F = S2ад / S2экс . |
(4.19) |
Если при заданном уровне значимости α отношение F окажется меньше квантили F1-α(k1, k2), где k1 = kад, k2 = kэкс, то рассматриваемая модель не противоречит результатам эксперимента и принимается; в противоположном случае модель отвергается с уровнем значимости α, как противоречащая результатам эксперимента. В построенной регрессионной модели (4.4) некоторые коэффициенты могут быть незначимы, т.е. может выполняться гипотеза H0 : Bj = 0. Для проверки этой гипотезы можно найти доверительный интервал для коэффициента Bj с уровнем значимости α. Если этот интервал «накрывает» значение Bj = 0, гипотеза H0 принимается и коэффициент Bj признается незначимым, в противоположном случае коэффициент Bj значим.

( j = 1, 2, 3, ...),
,