ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 20.12.2024
Просмотров: 423
Скачиваний: 1
4.1. Определение коэффициентов регрессии При определении коэффициентов уравнения регрессии будем использовать следующие функции для транспонирования, перемножения и обращения матриц:
ТРАНСП(матрица) – возвращает транспонированный массив;
МОБР(матрица) – возвращает обратную матрицу для матрицы, хранящейся в массиве;
МУМНОЖ(матрица1; матрица2) - возвращает произведение матриц (матрицы хранятся в массивах). Результатом является матрица с таким же числом строк, как матрица1 и с таким же числом столбцов, как матрица2.
При вычислении коэффициентов регрессии (п.I табл.2) необходимо провести вычисление промежуточной матрицы (XTX)-1, которая будет использована для вычисления матрицы ковариаций (п.IV табл.2). Формула для вычисления такой матрицы будет выглядеть:
“=МОБР(МУМНОЖ(ТРАНСП(X_1);X_1))” , где X_1 - ссылка на поименованный массив для линейного уравнения регрессии.
Для вычисления коэффициентов регрессии используем формулу (33), для этого вычислим матрицу коэффициентов регрессии, предварительно выделив соответствующий диапазон ячеек, и в строке формул введем:
“=МУМНОЖ(МУМНОЖ(<диапазон>;ТРАНСП(X_1));Y)” ** , где <диапазон> – диапазон ячеек вычисленной матрицы (XTX)-1.
При вычислении коэффициентов регрессии для нелинейных видов уравнений в качестве матрицы X будем использовать ссылки на поименованные области X_2 и X_3.
4.2. Определение коэффициента детерминации
При определении коэффициента детерминации R2 воспользуемся формулой “=СРЗНАЧ(Y)^2” для вычисления Yср2 . Данное значение определим в ячейке F14.
Для упрощения вычислений предлагается отдельно вычислить знаменатель формулы (35) и занести результат в ячейку F15 :
“=МУМНОЖ(ТРАНСП(Y);Y)-n*F14”, где n – ссылка на количество наблюдений выборки и F14 - ссылка на вычисленное значение Yср2.
Далее по формуле (35) вычислим коэффициент детерминации R2 :
=(МУМНОЖ(МУМНОЖ(ТРАНСП(F9:F10);ТРАНСП(X_1));Y) -n*F14)/F15, где F15 - вычисленный знаменатель формулы (30).
|
|
|
|
|
|
Линейная регрессия |
Нелинейная регрессия |
Нелинейная регрессия |
||||||||
|
|
|
n= |
156 |
|
y = b0+b1*x |
y = b0+b1*x+b2*x2 |
y = b0+b1*x+b2*x2+b3*x3 |
||||||||
|
|
|
k= |
|
|
1 |
|
|
2 |
|
|
|
3 |
|
|
|
|
I. |
Определение коэффициентов регрессии |
|
|
|
|
|
|
|
|
||||||
|
|
Матрица (XTX)-1 |
0.03543 |
-0.00675 |
|
0.17023 |
-0.08257 |
0.00847 |
|
0.82886 |
-0.67643 |
0.15823 |
-0.01112 |
|||
|
|
|
|
|
|
-0.00675 |
0.00157 |
|
-0.08257 |
0.04421 |
-0.00476 |
|
-0.67643 |
0.57967 |
-0.13979 |
0.01002 |
|
|
|
|
|
|
|
|
|
0.00847 |
-0.00476 |
0.00053 |
|
0.15823 |
-0.13979 |
0.03458 |
-0.00253 |
|
|
|
|
|
|
|
|
|
|
|
|
|
-0.01112 |
0.01002 |
-0.00253 |
0.00019 |
|
|
Кэффициенты регрессии b |
1.37 |
|
|
2.82650 |
|
|
|
2.72380 |
|
|
|
|||
|
|
|
|
|
|
4.69429 |
|
|
3.87314 |
|
|
|
3.96574 |
|
|
|
|
|
|
|
|
|
|
|
|
0.09174 |
|
|
|
0.06839 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.00173 |
|
|
|
|
II. |
Определение коэффициента детерминации |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
Yср2 |
464.27 |
|
|
464.27 |
|
|
|
464.27 |
|
|
|
|
|
yTy=YTY-nYср2 |
14617 |
|
|
14617 |
|
|
|
14617 |
|
|
|
|||
|
|
|
|
|
R2 |
0.960177 |
|
|
0.961258 |
|
|
|
0.961260 |
|
|
|
|
III. |
Оценка значений дисперсии |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
3.780 |
|
|
3.701 |
|
|
|
3.726 |
|
|
|
|
IV. |
Определение дисперсии оценок i |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
S2(i) |
0.133909 |
|
|
0.630084 |
|
|
|
3.08796 |
|
|
|
|
|
|
|
|
|
0.005935 |
|
|
0.163651 |
|
|
|
2.15958 |
|
|
|
|
|
|
|
|
|
|
|
|
0.00197 |
|
|
|
0.12884 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6.99E-04 |
|
|
|
|
|
Среднеквадратические отклонения i |
0.36594 |
|
|
0.79378 |
|
|
|
1.75726 |
|
|
|
|||
|
|
|
|
|
|
0.07704 |
|
|
0.40454 |
|
|
|
1.46955 |
|
|
|
|
|
|
|
|
|
|
|
|
0.04439 |
|
|
|
0.35894 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.02644 |
|
|
|
|
V. |
Проверка значимости коэффициентов i |
|
|
|
|
|
|
|
|
||||||
|
|
Доверительная вероятность |
0.05 |
|
|
0.05 |
|
|
|
0.05 |
|
|
|
|||
|
|
Табличное значение t-статистики |
1.97549 |
|
|
1.97559 |
|
|
|
1.97570 |
|
|
|
|||
|
|
Оценочные значения для i |
3.73 |
|
|
3.56 |
|
|
|
1.55 |
|
|
|
|||
|
|
|
|
|
|
60.94 |
|
|
9.57 |
|
|
|
2.70 |
|
|
|
|
|
|
|
|
|
|
|
|
2.07 |
|
|
|
0.19 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.07 |
|
|
|
|
VI. |
Построение доверительных интервалов для коэффициентов регрессии |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
0 |
1 |
|
0 |
1 |
2 |
|
0 |
1 |
2 |
3 |
|
|
Верхняя граница |
2.089421 |
4.84647 |
|
4.394595 |
4.6723 |
0.17943 |
|
6.195243 |
6.86882 |
0.77748 |
0.05397 |
|||
|
|
Нижняя граница |
0.643621 |
4.5421 |
|
1.258399 |
3.07398 |
0.00406 |
|
-0.74764 |
1.06266 |
-0.6407 |
-0.0505 |
|||
|
VI. |
Проверка значимости уравнения регрессии |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
Fф |
3713.091 |
|
|
1898.126 |
|
|
|
1257.184 |
|
|
|
|
|
|
|
|
F |
3.902557 |
|
|
3.055163 |
|
|
|
2.664109 |
|
|
|
4.3. Оценка значений дисперсии
Для вычисления значений дисперсии используем формулу (34). Числитель приведенной формулы определим через вычисленные значения:
.
Определяем дисперсию:
=F15*(1-R2)/(n-k-1), где R2- ссылка на вычисленное значение коэффициента детерминации, F15 - вычисленный знаменатель формулы (30).
Результаты вычислений значений дисперсии для каждого вида уравнения необходимо занести в соответствующие ячейки F18, I18, M18.
4.4. Определение дисперсии оценок i
Для определения дисперсии оценок i по формуле (37) используем результаты вычисления предыдущего расчета и диагональные элементы матрицы (XTX)-1 вычисленные в п.4.1.
Для вычисления
среднеквадратических отклонений
используйте функциюКОРЕНЬ(число).
4.5. Проверка значимости коэффициентовi
Для проверки значимости коэффициентов регрессии используем t-статистику. Для проверки неравенства (41) в ячейках F30, I30, M30 (для каждого уравнения регрессии) вычисляем критическую значение распределения Стьюдента по уровню значимости a = 0,05 и n-(k+1) степеням свободы с помощью встроенной функции:
СТЬЮДРАСПОБР(0.05; <степеньсвободы>).
Вычисляем соотношение (41), сравниваем с критическим значением tкр и делаем вывод о значимости коэффициентов регрессии для каждого уравнения.
4.6. Построение доверительных интервалов для коэффициентов регрессии
Для выбранного уравнения регрессии, у которого в предыдущем пункте определили значимые коэффициенты регрессии, используя формулу (39), определяем доверительные интервалы для коэффициентов регрессии.
Например для линейного уравнения в ячейки F38,G38,F39,G39 внесены формулы:
=F$9+$F$30*F24;
=F$10+$F$30*F25;
=F$9-$F$30*F24;
=F$10-$F$30*F25, где $F$30 - ссылка на критическое значение распределения Стьюдента tкр.
4.7. Проверка значимости уравнения регрессии
Для оценки значимости уравнения регрессии используется критерий Фишера. Наблюдаемое значение Fнабл в формуле (43) вычислим по следующему соотношению:
Критическое значение F-распределения Фишера Fкр с (k;n-k-1) степенями свободы определяем с помощью встроенной функции:
=FРАСПОБР(0.05;k;n-k-1) - где n - число наблюдений, k - число независимых параметров уравнения.
Проводим сравнение Fнабл и Fкр.
4.8. Построение диаграммы
Для графической интерпретации наблюдений построим диаграмму.
5. ЗАДАНИЕ К ЛАБОРАТОРНОЙ РАБОТЕ
Для каждого варианта лабораторной работы необходимо выполнить следующие расчеты по полученной с помощью программы выборке:
1) Сгенерировать выборку наблюдений.
2) Сформировать матрицы X и Y для проведения регрессионного анализа.
3) Подготовить таблицу для анализа трех видов уравнения регрессии.
y = b0+b1x
y = b0+b1x+b2x2
y = b0+b1x+b2x2+b3x3
4) Используя функции действия с матрицами вычислить коэффициенты регрессии.
5) Вычислить коэффициенты детерминации.
6) Определить значения дисперсии.
7) Определить дисперсии оценок i.
8) Провести проверку значимости коэффициентов i при уровне значимости g = 0,95 и выбрать одно из трех рассчитанных уравнений регрессии, у которого значимы все коэффициенты и наибольший коэффициент детерминации.
9) Рассчитать доверительные интервалы для коэффициентов выбранного уравнения регрессии при доверительной вероятности g = 0,95.
10) Провести оценку значимости уравнения регрессии.
11) Построить диаграмму наблюдений выборки.
Отчет по лабораторной работе должен содержать результаты всех промежуточных расчетов и сделанные на их основе выводы. Дополнительно к отчету должны быть представлены контрольные результаты, представленные в рабочем листе «Контрольный_результат».