51443

Выбор модели методом шаговой регрессии

Лабораторная работа

Математика и математический анализ

Составить множество регрессоров, включив в него факторы, квадраты факторов и их взаимные произведения Проверить условие Fk1k2 Fpk1k2 выполнение которого свидетельствует о целесообразности произведенного усложнения модели что обусловило существенное увеличение точности аппроксимации моделью исходных данных.

Русский

2014-02-11

33 KB

6 чел.

Лабораторная работа

Выбор модели методом шаговой регрессии

1. Составить множество регрессоров, включив в него факторы, квадраты факторов и их взаимные произведения:

X1, X2, X3, X4, X12, X22, X32, X42, X1X2, X1X3, X1X4, X2X3, X2X4, X3X4.

2. С помощью статистической функции КОРРЕЛ (массив1, массив2) найти коэффициенты корреляции Y со всеми регрессорами (строка 1, табл.1)

Табл.1

X1

X2

X3

X4

X12

...

X2X4

X3X4

MAX

Y

max

Y-Y1

max

Y-Y2

...

3.Выбрать регрессор, имеющий максимальное по модулю значение коэффициента парной корреляции с зависимой переменной Y. Включить этот регрессор в модель

Y1 = a0 + a1*X2Х4. (1)

4.Методом наименьших квадратов (МНК) с помощью статистической функции ЛИНЕЙН найти коэффициенты  аi   и рассчитать модельные значения Y1 по формуле (1).

5.Найти невязку (Y-Y1) и определить коэффициенты корреляции невязки со всеми регрессорами. Занести данные в табл.1.

6. Опять выбрать регрессор, имеющий максимальное значение коэффициента парной корреляции с переменной (Y-Y1). Включить этот регрессор в усложненную модель

Y2 = a0 + a1*X2Х4 + а2*Х3. (2)

7. Для модели (2) с помощью МНК рассчитать оценки а0, а1, а2.

8.Рассчитать для двух моделей (1) и (2)   F-статистику

F12=k*(S - S')/S', (3)

где k  - степень свободы для сложной модели,  S' - остаточная сумма квадратов для сложной модели,  S  - остаточная сумма квадратов для простой модели.

  1.  Проверить условие  Fk1,k2 > F(p,k1,k2), выполнение которого свидетельствует о целесообразности произведенного усложнения модели, что обусловило существенное увеличение точности аппроксимации моделью исходных данных.

Как только процесс пошагового усложнения модели окажется неэффективным, необходимо выполнить по t-критерию Стьюдента

(4)

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

11. П.п.3-5 повторяются, пока пошаговое усложнение модели будет эффективным.

12. При появлении сомнений в необходимости последующего усложнения модели производится комплексный анализ качества модели по F-, t-статистикам, по показателям коэффициента детерминации R2, по значениям стандартной ошибки оценивания Y.

13. Построить график Yp и Y с помощью Мастера диаграмм.

Табличные значения статистик F(p,k1,k2),    t(p,k):

F(0,95;1;28)=4,20

t(0,95;28)=2,05

F(0,95;1;27)=4,21

t(0,95;27)=2,05

F(0,95;1;26)=4,22

t(0,95;26)=2,06

F(0,95;1;25)=4,24

t(0,95;25)=2,06

F(0,95;1;24)=4,26

t(0,95;24)=2,06


 

А также другие работы, которые могут Вас заинтересовать

14620. АВТОМАТИЧЕСКОЕ ВКЛЮЧЕНИЕ РЕЗЕРВНОГО ИСТОЧНИКА ПИТАНИЯ 107 KB
  Отчет по лабораторному практикуму Автоматическое включение резервного источника питания Цель работы Работа предназначена для практического ознакомления с принципом выполнения схемы автоматического включения резервного источника питания АВР. Исход
14621. Усталостные испытания 2.57 MB
  Усталостные испытания Методические указания к лабораторным практическим работам и КНИРС по специальным дисциплинам для студентов всех металловедческих и материаловедческих специальностей Усталостные испытания: Методические указания к лабор
14622. Решение обратной задачи кинематики трехзвенного манипулятора 96 KB
  Лабораторная работа №5: Вариант 1 Решение обратной задачи кинематики трехзвенного манипулятора. Цель работы: изучение алгоритмов решения обратной задачи кинематики. Решение ПЗП для трехзвенного манипулятора с вращательными парами: Дано: ...
14623. Решение обратной задачи кинематики двухзвенного манипулятора 176.5 KB
  Лабораторная работа №4: Вариант 2 Решение обратной задачи кинематики двухзвенного манипулятора. Цель работы: изучение алгоритмов решения обратной задачи кинематики Решение ПЗП для двухзвенного манипулятора с вращательными парами: Дано: Получим р...
14624. Решение прямой задачи кинематики манипулятора 294 KB
  Лабораторная работа №3: Вариант 1 Решение прямой задачи кинематики манипулятора. Цель работы: решение прямой задачи о положении манипуляционной системы на ЭВМ на основе формализованного описания кинематических цепей Геометрические характеристики звеньев: ...
14625. Изучение метода преобразования систем координат промышленных роботов и кодирования кинематических цепей «иркутским методом» на примере робота МП-9С (Ритм -01-02) 196.5 KB
  Лабораторная работа №2: Вариант 1 Изучение метода преобразования систем координат промышленных роботов и кодирования кинематических цепей иркутским методом на примере робота МП9С Ритм 0102 . ЦЕЛЬ РАБОТЫ: – выбор абсолютной и связанных систем координат ПР; –
14626. Изучение методики разработки программ в системе MATLAB при изучении кинематического управления роботами 156.5 KB
  Лабораторная работа №1: Вариант 1 Изучение методики разработки программ в системе MATLAB при изучении кинематического управления роботами. Изучение методики разработки программ в системе MATLAB при изучении кинематического управления роботами Цель работы: Изу...
14627. ИЗУЧЕНИЕ ДЕЙСТВИЯ УЛЬТРАЗВУКОВЫХ КОЛЕБАНИЙ НА ВЕЩЕСТВО И ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ УЛЬТРАЗВУКА 54.5 KB
  2 ЛАБОРАТОРНАЯ РАБОТА ИЗУЧЕНИЕ ДЕЙСТВИЯ УЛЬТРАЗВУКОВЫХ КОЛЕБАНИЙ НА ВЕЩЕСТВО И ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ УЛЬТРАЗВУКА ЦЕЛЬ РАБОТЫ: изучить свойства ультразвука его взаимодействие с веществом; ознакомиться с устройством и работой ультразву...
14628. ДИНАМИЧЕСКИЕ ИСПЫТАНИЯ НА УДАРНЫЙ ИЗГИБ 452 KB
  ДИНАМИЧЕСКИЕ ИСПЫТАНИЯ НА УДАРНЫЙ ИЗГИБ Методические указания к лабораторным и практическим работам по специальным дисциплинам для студентов металловедческих специальностей Данные методические указания включают в себя понятие ударной вязкости и методы ее опред