Дисперсионный анализ

Содержание.                                                   ст   

Введение…………………………………………………………………………....3

1.Задача дисперсионного анализа……………………………………………...4

2. Однофакторный дисперсионный анализ ………………….………………..6

3.Двухфакторный  дисперсионный анализ…………………………………….11

4. Планирование  эксперимента при дисперсионном  анализе. 

Латинские и гипер-греко-латинские квадраты.………………..…………......22

Заключение...............................................................................................................30

Литература................................................................................................................32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введение.

Результаты  эксперимента в химии и химической промышленности были и остаются главным  критерием при решении практических задач и при проверке теоретических  гипотез. Изучение сложных технологических  процессов сопряжено с трудоемким и длительным экспериментом. Для  увеличения эффективности научных исследований, сокращения сроков разработки новых технологических процессов необходима оптимизация экспериментальных исследований на всех стадиях разработки, исследования, внедрения и эксплуатации химико-технологических процессов.

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

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

Математической  моделью служит функция отклика, связывающая параметр оптимизации, характеризующий результаты эксперимента, с переменными параметрами, которыми варьируют при проведении опытов.

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

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

1. Задача дисперсионного анализа.

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

Дисперсионный анализ использует свойство аддитивности дисперсии изучаемой случайной величины, обусловленной действием независимых факторов. Р.А. Фишер в 1938 г. впервые определил дисперсионный анализ как «отделение дисперсии, приписываемой одной группе причин от дисперсии, приписываемой другими группами». В зависимости от числа источников дисперсии различают однофакторный  и многофакторный дисперсионный анализ.

Дисперсионный анализ особенно эффективен при изучении нескольких факторов. При классическом методе исследования варьируют только один фактор, а остальные оставляют постоянными. при этом для каждого фактора проводится серия наблюдений, не используемая при изучении других факторов. Кроме того, при таком методе исследования не удается определить взаимодействие факторов при одновременном их изменении. При дисперсионном анализе каждое наблюдение служит для одновременной оценки всех факторов и их взаимодействий.

Дисперсионный анализ состоит в выделении и  оценке отдельных факторов, вызывающих изменчивость изучаемой случайной  величины. Для этого производится разложение суммарной выборочной дисперсии на составляющие, обусловленные независимыми факторами. Каждая из этих составляющих представляет собой оценку дисперсии генеральной совокупности. Чтобы решить, значимо ли влияние данного фактора, необходимо оценить значимость соответствующей выборочной дисперсии в сравнении с дисперсией воспроизводимости, обусловленной случайными факторами. Проверка значимости оценок дисперсий проводится по критерию Фишера. Если рассчитанное значение критерия Фишера окажется меньше табличного, то влияние рассматриваемого фактора нет оснований считать значимым. Если же рассчитанное значение критерия Фишера окажется больше табличного, то рассматриваемый фактор влияет на изменчивость средних. В дальнейшем будем полагать, что выполняются следующие допущения:

- случайные  ошибки наблюдений имеют нормальное  распределение;

- факторы  влияют только на изменение  средних значений, а дисперсия  наблюдений остается постоянной;

- эксперименты  равноточны.

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

Факторы рассматриваемые в дисперсионном анализе, бывают двух родов:

- со случайными  уровнями;

- с фиксированными.

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

Дисперсионный анализ может применяться в различных  формах в зависимости от структуры исследуемого процесса; выбор соответствующей формы является обычно одной из главных трудностей в практическом применении анализа.

2. Однофакторный дисперсионный  анализ.

Рассмотрим  действие единичного фактора А (количественного или качественного), который принимает k различных значений (уровней фактора). На i-м уровне производится n_i наблюдений, результаты которых можно записать следующим образом:

 

Будем предполагать, что результат любого наблюдения можно представить в виде модели:

                                              (1.1)

где,

μ – суммарный эффект во всех опытах;

- эффект фактора на i-м уровне (i=1,2…,k);

ε_ij – ошибка измерения на i-м уровне. Предположим также, что наблюдения на фиксированном уровне фактора нормально распределены относительно среднего значения μ+ с общей дисперсией . Общее число опытов равно N:

                                       (1.2)

проверяется нулевая гипотеза равенства средних  значений на различных уровнях А:

 

Наиболее  простые расчеты получаются при  равном числе опытов на каждом уровне фактора А: n₁=n₂=…=n_k=n.

При этом общее число наблюдений N равно kn. Обозначим через , среднее значение наблюдений на i-м уровне.

                                                  (1.3) 

а общее среднее  значение для всей выборки из N наблюдений:

                                  (1.4)

Для проведения дисперсионного анализа необходимо общую выборочную дисперсию s².

                  (1.5)

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

 

i=1,2,…,k.                                          (1.6)

Если  нет уверенности в равноточности экспериментов, однородность дисперсий s₁², s₂²,…, s_k², можно проверить по критерию Кохнера.

Если  между дисперсиями нет значимых различий, для оценки генеральной  дисперсии σ², характеризующнй фактор случайности. используют выборочную дисперсию s²_общ:

      (1.7)

Число степеней свободы дисперсии s²_общ равно k(n-1)=N-k. Приближенную оценку для дисперсии фактора А можно получить следующим образом:

                                                (1.8)

Более точную оценку для можно получить, рассматривая отклонение средних на отдельных уровнях от общего среднего всей выборки . Действительно,

                    (1.9)

Отсюда

                                       (1.10)

Дисперсия фактора А для модели с фиксированными уровнями не связана ни с какой случайной, это условное название для математического ожидания среднего квадрата отклонений, обусловленного влиянием фактора А. Такое обозначение удобно, так как определяет рассеяние, вызванное влиянием фактора А аналогично показателю влияния случайного фактора. что позволяет непосредственно сравнивать фактор А с эффектом случайности. Введем также следующее обозначение:

                         (1.11)

Эта дисперсия  имеет k-1степеней свободы. Если дисперсия значимо отличается от , нулевая гипотеза m₁=m₂=…=m_k=m отвергается, влияние фактора А считается существенным. Проверяется нулевая гипотеза по критерию Фишера. Так как альтернативой является равенство , для проверки гипотезы применяется односторонний критерий Фишера. Влияние фактора А считается значимым:

                                                 (1.12)

 

 

Дисперсионный анализ можно провести по следующему алгоритму: 1. итоги по столбцам

                                             (1.13)

2. сумму  квадратов всех наблюдений

                                         (1.14)

3. сумму квадратов по столбцам, деленную на число наблюдений в столбце,

                                               (1.15)

4. квадрат  общего итога, деленный на число  всех наблюдений (корректирующей член),

                                        (1.16)

5. сумму  квадратов для столбца

                                        (1.17)

6. - общую сумму квадратов, равную разнице между суммой квадратов всех наблюдений и корректирующим членом,

                                     (1.18)

7. - остаточную сумму квадратов для оценки ошибки эксперимента

                                     (1.19)

8. дисперсию

                                                    (1.20)

9. дисперсию 

                                             (1.21)

Результаты  расчета представляются в виде таблицы  дисперсионного анализа:

Таблицы №1.

Однофакторный дисперсионный анализ (с равным числом повторений опытов)

 

Источник дисперсии	Число степеней свободы	Сумма квадратов	Средний квадрат	Математическое ожидание среднего квадрата
А	k-1
Остаток	k(n-1)
Общая сумма	kn-1			-

 

Если  отношение , то влияние фактора А следует считать незначимым. При этом общая дисперсия s² связана только с фактором случайности и может служить оценкой для дисперсии воспроизводимости. Такая оценка лучше, чем, так как имеет большее число степеней свободы, равное kn-1. при интерпретации результатов дисперсионного анализа необходимо иметь в виду, что очень низкое значение дисперсионного отношения может быть связано с тем, что влияние какого-то важного неконтролируемого в ходе эксперимента фактора не было рандомизированно. Это может увеличить дисперсию внутри уровней, а дисперсию между уровнями оставить неизменной, что уменьшает дисперсионное отношение.