Идентификация параметров модели кинетики сложной химической реакции, страница 3

                                                       Sост ,                                              (12)

Таким образом, опыт, в котором  |ri|  > 2 Sост  является выбросом.

Проверка значимости  влияния факторов на отклик.

При  постановке  экспериментального  исследования,  как  правило, заранее неизвестна степень влияния  отдельных  факторов  на  отклик. Может возникнуть ситуация, когда результаты эксперимента не  зависят или слабо зависят от некоторого фактора, который, тем не менее, включен в уравнение регрессии. Для проверки значимости влияния факторов проверяется   гипотеза   о   равенстве   нулю математического ожидания коэффициента bi, H0: M[bi]=0.Проверка H0 выполняется по критерию Стьюдента

                                          t i =| bi | /  .                                                        (13)

По критерию Стьюдента для всех коэффициентов регрессии с помощью функции распределения можно рассчитать уровень значимости, назовём его plevel.  Если для коэффициента при каком-то факторе значение plevelпревысит 0.05 (уровень значимости при вероятности 0.95), то гипотеза H0 принимается и  предполагается, что  этот фактор не оказывает влияние на отклик, т.е. коэффициент незначим. Для проверки этого предположения следует исключить фактор из уравнения регрессии, пересчитать коэффициенты и проанализировать, как повлияло исключение фактора на точность регрессии. Если точность увеличилась ( уменьшилась) или не изменилась, фактор действительно следует исключить из регрессии. Если точность уменьшилась ( увеличилась), фактор исключать не следует.

Проверка адекватности регрессии.

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

·  В эксперименте выполнялись параллельные опыты. В этом случае   дисперсия воспроизводимости S2воспр характеризует усредненную погрешность измерения отклика и, если отклики измерены с  большой  погрешностью, то большая величина     вряд ли связана с плохим  выбором  уравнения  регрессии.   Таким  образом,  для   определения  соответствия   (адекватности) регрессии нужно сравнить остаточную дисперсию  с величиной дисперсии воспроизводимости. Если различие дисперсий незначимо, модель следует признать адекватной. Проверка значимости различия производится на основе статистической гипотезы о равенстве соответствующих генеральных дисперсий H0: .  Для сравнения двух дисперсий используют критерий Фишера 

                                  F = S2ОСТ/ S2воспр ,                                                               (14)

где в числителе и знаменателе стоят выборочные оценки генеральных дисперсий.

Проверка выполняется по стандартной схеме:

    - задается уровень значимости p=0.05;

    - по числу степеней свободы остаточной дисперсии f1=m-n-1  и числу степе         ней свободы дисперсии воспроизводимости  f2=m(L-1) (см. Примечание ниже)для заданного p по таблице находят квантиль распределения Фишера F1-p ;

    - если значение  F  в (14) меньше или равно F1-pто гипотезу H0 следует принять с доверительной вероятностью 1 - p (95%);

    - если F > F1-p , гипотезу следует отклонить и признать модель неадекватной.

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