- •Самарский государственный университет.
- •Глава 5. Анализ временных рядов. 29
- •Глава 6. Дисперсионный анализ. 34
- •Глава 7. Метод главных компонент 47
- •Глава 8. Факторный анализ 53
- •Глава 1. Математические модели.
- •1.1. Классификация моделей.
- •1.2. Этапы построения математических моделей
- •2.1.Коэффициент парной корреляции
- •Глава 2. Корреляционный анализ: вид связи и проверка гипотез.
- •2.2. Коэффициент детерминации
- •2.3. Статистическая оценка значимости коэффициентов корреляции.
- •2.4. Интервальные оценки для значимых парных коэффициентов корреляции.
- •2.5. Коэффициент корреляции рангов.
- •Решение
- •Глава 3. Регрессионный анализ: парная линейная регрессия.
- •3.1. Общие сведения. Парная линейная регрессия.
- •3.1.1. Функция регрессии
- •3.1.2. Последовательность проведения регрессионного анализа
- •3.1.3. Предпосылки к проведению регрессионного анализа
- •3.2. Метод наименьших квадратов
- •3.2.1. Проверка адекватности модели
- •3.2.2. Проверка значимости параметров модели
- •3.3. Описание типового примера. Алгоритм построения регрессионной модели.
- •3.4. Использование оцененной модели для прогнозирования
- •Глава 4. Множественная линейная регрессия
- •4.1. Корреляционно - регрессионный анализ.
- •4.2. Некоторые нелинейные модели, сводящиеся к линейным
- •4.3. Проверка предпосылок регрессионного анализа
- •1. Проверка нормальности закона распределения ошибок
- •3. Проверка на автокорреляцию случайных ошибок
- •Замечание. Если в модели много незначимых коэффициентов - ее надо упростить, т.Е. Уменьшить число искомых величин b. Если модель неадекватна, ее надо усложнить, т.Е. Добавить новые члены.
- •4.4. Пример решения задачи моделирования с использованием метода наименьших квадратов.
- •4.5. Метод наименьших квадратов (мнк) – общий случай
- •Глава 5. Анализ временных рядов.
- •5.1. Проверка законов распределения.
- •1. Проверка подчинения эмпирического распределения нормальному закону распределения
- •2. Проверка подчинения эмпирического распределения закону распределения Пуассона.
- •5.2. Изучение динамики
- •5.2.1. Средние показатели тенденции динамики
- •5.2.2. Выявление типа тенденции динамики тренда
- •5.2.3. Определение оптимального значения тренда
- •Глава 6. Дисперсионный анализ.
- •Однофакторный дисперсионный анализ.
- •6.2. Применение однофакторного анализа.
- •6.3. Проверка гипотезы об отсутствии влияния факторов на параметр
- •Можно показать [5], что, как и в случае однофакторного анализа, разность
- •Трехфакторный дисперсионный анализ.
- •6.5. Отсеивающие эксперименты и дисперсионный анализ
- •6.6. Рандомизация эксперимента
- •6.7. Определение аналитического вида зависимости параметра от факторов.
- •6.7.1. Матрица планирования
- •Глава 7. Метод главных компонент
- •7.1. Математическая модель метода главных компонент.
- •7.2. Алгоритм метода главных компонент
- •7.3. Нахождение коэффициентов характеристического уравнения
- •Глава 8. Факторный анализ
- •8.1. Модель факторного анализа.
- •8.2. Вращение факторов.
- •8.3. Применение факторного анализа.
- •8.4. Некоторые результаты факторного анализа
- •Приложение 1. Статистические таблицы
- •При числе степеней свободы
- •Критические точки распределения критерия g
- •Значения верхнего предела n2 в зависимости от уровня значимости
Можно показать [5], что, как и в случае однофакторного анализа, разность
Q3= Q – Q1- Q2(6.3.9.)
характеризует отклонение от среднего не из-за влияния факторов, а из-за случайных погрешностей эксперимента (к ним следует отнести и влияние факторов, не контролируемых в данном эксперименте). Вычисляем Q3:
Q3= 42.02 - 20.37 - 19.16 = 2.49.
Если гипотеза об отсутствии влияния факторов верна, то отношение всех Qк общей дисперсии распределено по закону2, причем число степеней свободы для Q/2 равно (rv-1), для Q1/2– (r-1), для Q2/2- (v-1), тогда на «долю» Q3/2 остается (rv-1)-(r-1)-(v-1) = (r-1)(v-1) степеней свободы. По этим значениям находим соответствующие оценки дисперсий:
![]()
![]()
![]()
Если факторы не влияют на результат, то все это - оценки одной и той же дисперсии. Если влияют, то в S12отразилось влияние первого фактора, в S22- второго. Проверка влияния факторов ведется по критерию Фишера.
,m=1, 2.
Если Fmэ>Fq, тоm-й фактор влияет на результат. Выбрав= 0.05, имеем:
для первого фактора при k1=6,k2=12Fq=2.996,F1э=3.395/0.207=16.38.
Гипотеза о том, что первый фактор (конструкция) не влияет на светоотдачу, отвергается.Отсюда следует важный практический вывод: все горелки соответствовали действующим техническим условиям на их изготовление, разброс в конструктивных параметрах не превышал допусков на соответствующие элементы конструкций. Следовательно, действующие технические условия нуждались в корректировке с целью закрепления конструктивных параметров, обеспечивающих максимальную светоотдачу. Для второго фактора приk1=2,k2=12 получим:
Fq=3.88,F2э=9.58/0.207=46.2. ,
тогда гипотеза о том, что второй фактор (тепловой режим) не влияет на светоотдачу, такжеотвергается. Следовательно, тепловой режим заметно влияет на светоотдачу.
Как видно, результат статистической обработки даже небольшой серии измерений приводит к важным практическим результатам, показывает необходимость дальнейшей работы по оптимизации разрабатываемого изделия и позволяет наметить пути достижения стабильных параметров.
Трехфакторный дисперсионный анализ.
Рассмотрим пример, иллюстрирующий трехфакторный дисперсионный анализ.
Предположим, мы хотим выяснить влияние различных компонент наполнения на выход излучения металлогалоидных ламп. Изготовим 16 образцов горелок, в которых варьируются:
дозировка NaJ (первый фактор);
давление Ar (второй фактор);
дозировка TlJ (третий фактор).
Выполним эксперимент по плану, представленному в таблице 3.8. Условные обозначения уровней факторов здесь примут конкретные значения. Для первого фактора (NaJ):
а - соответствует 5 мг;
b - соответствует 20 мг;
с - соответствует 30 мг;
d - соответствует 40 мг.
Для второго фактора (давление Ar):
А - 10 мм рт. ст.;
В - 15 мм рт. ст.;
С - 20 мм рт. ст.;
D – 25 мм рт. ст.
Для третьего фактора (TlJ):
- 2 мг;
- 4 мг;
- 6 мг;
- 8 мг.
Результат измерения светоотдачи (в относительных единицах) приведен в той же таблице.
Светоотдача ламп в условных единицах : Таблица 6.4.1
-
Фактор 1
j
Фактор 2 i n-уровней
n
j
A
B
C
D

у
р
a
14
17
22
39
23
о
в
b
45
11
15
18
22.25
н
е
c
19
41
12
14
21.5
й
d
13
23
36
10
20.5

22.7
23
21.2
20.25
21.8
Результаты
надо понимать, например, так: горелка с
5 мг NaJ, 5 мм рт.ст. Ar и 8 мг TlJ дала 45 ед.
светоотдачи. Как и раньше, вычислим
средние по строкам и по столбцам и
занесем эти данные в таблицу. Кроме
того, вычислим средние значения по
уровням третьего фактора
(p=,,,).
![]()
![]()
![]()
![]()
Найдем
также генеральное среднее всех значений
в таблице
и поместим его в правом нижнем углу
таблицы.
Оценка влияния факторов опять выполняется с помощью сравнения соответствующих оценок дисперсий, которые вычисляются из соответствующих сумм квадратов:
![]()
![]()
,
(6.4.1)
характеризующих влияние первого (Q1), второго (Q2)и третьего (Q3) факторов соответственно и «остаточной суммы квадратов»
Q4= Q – Q1- Q2- Q3, (6.4.2)
где Q - сумма квадратов отклонений всех значений в таблице от генерального среднего.
Число степеней свободы для каждой оценки, характеризующей влияние факторов, равно n-1, дляQоно, очевидно, равноn2-1, поэтому для «остаточной» оценки остаетсяn2-1-3(n-1) степеней свободы.
Проверка гипотезы об отсутствии влияния фактора с номером m(m=1,2,3) проводится по критерию Фишера, экспериментальное значение которого определяется отношением (6.19)
,
где
(m=1,2,3),
(6.4.3)
Критическое значение Fqдля всех факторов при выбранном уровне значимости одинаково и берется при числе степеней свободыk1=n-1,k2=n2-1-3(n-1).
В нашем примере:
Q1= 4[(23-21.8)2+ (22.25-21.8)2+(21.5-21.8)2+(20.5-21.3)2] = 13,
Q2= 4[(22.7-21.8)2+ (23-21.8)2+(21.2-21.8)2+ (20.25-21.8)2] = 17,
Q3 = 4[(11.7-21.8)2+ (14.7-21.8)2+ (20.5-21.8)2+(40.25-21.8)2] = 1968,
Q = (14-21.8)2+ (17-21.8)2+(22-21.8)2+(39-21.8)2+(45-21.8)2+ (11-21.8)2+(15-21.8)2+
(18-21.8)2 +(19-21.8)2 +(41-21.8)2 +(12-21.8)2 +(14-21.8)2+(13-21.8)2+ (23-21.8)2+ (36-21.8)2+ (10-21.8)2= 2047
Q4= 2047 - 13 - 17 - 1968 = 49.
k1=3,k2=6, соответственно:
S12=Q1/3=4.33, S22=Q2/3=5.66, S32=Q3/3=656, S42=Q4/6=8.16.
F1э=0.53, F2э=0.69, F3э=80.4 ; Fq =4.77 (при k1=3, k2=6, =0.05)
Делаем вывод, что факторы первый и второй не влияют на светоотдачу, т. е. дозировки NaJ и аргона при вариации их в изученных пределах не сказываются на светоотдаче (это, разумеется, не исключает возможности заметить влияние этих факторов, варьируя их в более широких пределах).
В то же время дозировка TlJ при изменении ее в диапазоне от 2 до 8 мг заметно влияет на светоотдачу.
Следовательно, в дальнейших экспериментах интервал варьирования количества NaJ и буферного газа следует расширить, если необходимо найти, когда они станут влиять на параметр, а характер влияния TlJ на светоотдачу требуется проследить более детально.
