Теория эксперимента. Учебное пособие
.pdf
В таблице 4.4. обозначено:
i=1,k |
– число уровней фактора ХА (число столбцов); |
|
J=1,m |
– |
число уровней фактора ХВ (число строк); |
l=1,n |
– |
число наблюдений ячейке (число параллельных опытов); |
y111, y112,…, y11n,…, уkmn – наблюдавшиеся значения выходного параметра;
yij – среднее значение в ячейке;
|
|
|
1 |
n |
|
|
y |
ij |
|
yijl ; |
(4.42) |
||
|
||||||
|
|
|
n l 1 |
|
||
y j |
– среднее значение по уровням фактора ХВ (средние по |
строкам);
|
|
|
1 |
k |
|
|
|
|
y |
j |
|
|
y |
ij*; |
(4.43) |
||
k |
||||||||
|
|
|
i 1 |
|
||||
yi |
– среднее значение по уровням фактора ХА (средние по |
столбцам);
|
|
|
1 |
m |
|
|
|
|
y |
i* |
|
|
y |
ij*; |
(4.44) |
||
|
||||||||
|
|
|
m j 1 |
|
||||
y – общее среднее;
|
|
1 |
k |
m |
|
|
|
|
|
|
|
||||||
y |
|
|
y |
ij* . |
(4.45) |
|||
|
||||||||
|
|
km i 1 |
j 1 |
|
||||
|
|
|
|
|
|
|
181 |
|
|
|
|
yijl di j |
di i |
ijl |
, |
(4.46) |
|
где |
|
– общая средняя (математическое ожидание для среднего |
||||||
|
|
|
по всей таблице 4.4); |
|
|
|
||
|
di |
– |
эффект |
фактора |
ХА |
на i-ом уровне (отклонение |
||
|
|
|
математического ожидания выходного параметра при i-ом |
|||||
|
|
|
уровне фактора ХА от общего математического ожидания); |
|||||
|
j |
– |
эффект |
фактора |
ХВ |
на j-ом уровне (отклонение |
||
|
|
|
математического ожидания выходного параметра при j-ом |
|||||
|
|
|
уровне |
фактора |
ХB |
от |
общего математического |
|
|
|
|
ожидания); |
|
|
|
|
|
|
di i |
– |
эффект взаимодействия факторов ХА и ХВ; |
|
||||
|
ijl |
– случайный остаток или вариация результатов |
внутри |
|||||
|
|
|
отдельной ячейки (ошибка воспроизводимости). |
|
||||
Эффект взаимодействия факторов di i |
представляет собой отклонение |
|||||||
среднего по наблюдениям в ij-ой серии от суммы первых трех членов модели
(4.46). Если предположить, что между факторами нет взаимодействия, то можно построить линейную модель вида:
yij di j ijl . (4.47)
Эта модель, как правило, применяется лишь в случае отсутствия параллельных опытов, например, как показано в таблице 4.5. (при анализе трех
182
и более факторов отдельные эффекты взаимодействия удается оценить и без параллельных наблюдений).
Таблица 4.5 – Исходные данные для двухфакторного ДА с одним наблюдением в ячейке
|
|
|
|
|
|
|
|
|
|
XA |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
ХВ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y j |
|||||||||
|
|
xA1 |
|
xA2 |
… |
|
|
xAi |
… |
|
|
xAK |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
xB1 |
|
y11 |
|
y21 |
… |
|
|
yi1 |
… |
|
|
yK1 |
|
|
|
|
|
1 |
|||||||
|
|
y |
||||||||||||||||||||||||
|
xB2 |
|
y12 |
|
y22 |
… |
|
|
y |
… |
|
|
yK2 |
|
|
|
|
|
|
|
||||||
|
|
y |
|
2 |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i2 |
|
|
|
|
|
|
|
|
|
|
|
||
|
… |
|
… |
|
… |
… |
|
|
… |
… |
|
|
… |
|
|
… |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xBj |
|
y1j |
|
y2 j |
… |
|
|
yij |
… |
|
|
yKj |
|
|
|
|
j |
||||||||
|
|
|
|
y |
||||||||||||||||||||||
|
… |
|
… |
|
… |
… |
|
|
… |
… |
|
|
… |
|
|
… |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xBm |
y1m |
|
y2m |
… |
|
|
yim |
… |
|
|
|
|
|
|
|
|
|
m |
||||||||
y |
|
|
||||||||||||||||||||||||
|
Km |
|
|
y |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
… |
|
|
|
|
… |
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
y |
|
|
y |
|
|
|
y |
|
|
y |
|
|
|
|
|
|
||||||
|
i |
|
1 |
|
2 |
|
|
i |
|
K |
|
|
|
|
y |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В таблице 4.5 обозначено: |
||||
i=1,k |
– число уровней фактора ХА (число столбцов); |
|||
J=1,m |
– число уровней фактора ХВ (число строк); |
|||
|
y11 ,…, уkmn |
– |
наблюдавшиеся значения выходного параметра; |
|
|
|
j |
– |
среднее значение по уровням фактора ХB (средние по |
|
y |
|||
столбцам);
|
|
|
1 |
k |
|
|
y |
*j |
|
yij ; |
(4.48) |
||
k |
||||||
|
|
|
i 1 |
|
yi |
– среднее значение по уровням фактора ХА (средние по |
столбцам);
183
|
|
1 |
m |
|
|
y |
i * |
yij ; |
(4.49) |
||
|
|||||
|
|
m j 1 |
|
||
y |
– общее |
среднее;
|
|
1 |
k |
m |
|
|
|
|
|
||||
y |
|
yij . |
(4.50) |
|||
|
||||||
|
|
km i 1 |
j 1 |
|
||
Врассматриваемом случае (модель (4.47)) оценки общей дисперсии, как
иранее, можно получить из основного тождества ДА. Однако следует подчеркнуть, что в двухфакторном ДА, в отличие от однофакторного, общая
сумма квадратов отклонений наблюдений yij от общего среднего y (числитель
(4.9)) раскладывается согласно (4.47) уже не на две, а на три части: часть,
обусловленную влиянием фактора ХА, часть, обусловленную влиянием фактора
ХВ, и часть, обусловленную влиянием всех неучтенных факторов.
|
|
|
k |
|
m |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
SSобщ |
|
|
yji |
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
i 1 |
|
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
k |
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
yij |
|
y |
i* |
y |
j* |
y |
|
y |
i* |
y |
|
y |
*j |
|
y |
|
|
|
||||||||||||||||||
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.51) |
m |
|
|
|
|
|
|
|
2 |
|
k |
|
|
2 |
|
k |
m |
|
|
|
|
|
|
|
2 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
k |
y |
*j |
y |
|
m y |
i* |
y |
|
|
yij |
|
y |
i* |
y |
j* |
y |
|
|
|||||||||||||||||||
j 1 |
|
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
|
|
||||||
SSXB SSXA S .
184
Слагаемое представляет собой сумму квадратов разностей между
средними по строкам и общим средним и характеризует изменение выходного параметра по фактору ХB:
m |
|
|
|
|
2 |
|
|
|
|
|
|||
SSXB k |
y |
*j |
|
y |
. |
(4.52) |
j 1 |
|
|
|
|
|
|
Слагаемое SSX A представляет собой сумму квадратов разностей между средними по строкам и общим средним и характеризует изменение выходного
параметра по фактору ХA:
k |
|
|
|
|
2 |
|
|
|
|
|
|||
SSXA m |
y |
i* |
|
y |
. |
(4.53) |
i 1
Слагаемое SSЕ называется остаточной суммой квадратов и характеризует влияние неучтенных моделью (4.47) факторов:
|
k |
m |
|
|
|
2 |
|
||||
|
|
|
|
||||||||
S |
|
|
yij |
|
y |
i* |
y |
j* |
y |
. |
(4.54) |
|
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
Зная суммы квадратов отклонений, определим оценки соответствующих дисперсий:
S2 |
общ |
SSобщ |
|
|
SSобщ |
; |
(4.55) |
|||
mk 1 |
|
|
||||||||
|
|
|
N 1 |
|
||||||
|
S2XB |
|
SSXB |
; |
(4.56) |
|||||
|
|
|||||||||
|
|
|
|
m 1 |
|
|||||
|
S2XA |
|
SSX A |
; |
(4.57) |
|||||
|
|
|||||||||
|
|
|
|
K 1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
185 |
S |
2 |
|
|
SS |
|
|
|
|
|
|
. |
(4.58) |
|
|
m 1 K 1 |
|||||
Нулевые гипотезы о незначимости факторов ХА и ХВ проверяют по критерию Фишера. Если дисперсионное отношение:
Fрасч |
S2XB |
Fтабл ( , f1, f2 ), |
(4.59) |
|
S2 |
|
|||
где
f1 m 1;
f2 m 1 k 1,
принимается гипотеза H0 : j 0.
Если дисперсионное отношение:
S2
Fрасч X2B Fтабл ( , f1, f2), (4.60)
S
нулевая гипотеза отвергается, и влияние фактора ХВ считается значимым.
Аналогично, если дисперсионное отношение
|
S2 |
|
|
|
Fрасч |
XA |
Fтабл ( , f1, f2 ), |
(4.61) |
|
|
||||
|
S2 |
|
|
|
где
186
f1 k 1;
f2 k 1 m 1,
принимается гипотеза H0 :di 0. При несправедливости неравенства:
Fрасч |
|
S2XA |
Fтабл ( , f1, f2 ), |
(4.62) |
|
S |
2 |
||||
|
|
|
|
|
|
влияние фактора ХА считается значимым. При проверке нулевых гипотез,
также принимается односторонний F-критерий (приложение Б).
При проведении ДА в условиях линейной модели (4.47) удобно использовать следующий алгоритм расчета. Находят:
- итоги по столбцам:
m
XAi yij ; j 1
- итоги по строкам:
k
XBj yij ;
i 1
-сумму квадратов всех наблюдений:
k m
Q1 yij2 ;
i 1 j 1
- сумму квадратов итогов по столбцам, деленную на число наблюдений в столбце:
187
|
1 |
k |
|
Q2 |
XA2i ; |
||
|
|||
|
m i 1 |
||
- сумму квадратов итогов по строкам, деленную на число наблюдений в строке:
|
|
1 |
m |
|
Q3 |
|
XB2j ; |
||
k |
||||
|
|
j 1 |
- квадрат общего итога, деленный на число всех наблюдений
(корректирующий член):
1 |
k |
|
2 |
1 |
m |
2 |
|
||
|
|
|
|
|
|
|
|
|
|
Q4 m k |
XAi |
|
m k |
XBj |
; |
||||
|
|
i 1 |
|
|
|
|
j 1 |
|
|
-сумму квадратов для столбца:
SSXA Q2 Q4 ;
-сумму квадратов для строки:
SSXB Q3 Q4 ;
-общую сумму квадратов:
SSобщ. Q1 Q4;
-остаточную сумму квадратов:
188
SS SSобщ SSXB SSXA
Q1 Q4 Q3 Q4 Q2 Q4
Q1 Q2 Q3 Q4;
-оценку дисперсии S2XA :
S2XA |
|
SSX A |
; |
|
k 1 |
||||
|
|
|
- оценку дисперсии S2XB :
S2XB |
|
SSXB |
|
; |
|
m 1 |
|||||
|
|
|
|||
- оценку дисперсии S2 :
S |
2 |
|
SS |
|
|
|
. |
||
|
m 1 k 1 |
|||
Результаты расчета удобно представлять в виде таблицы 4.6.
Таким образом, установив при помощи ДА значимость влияния факторов, выясняют затем, какие именно средние значения Y различны.
Как было отмечено ранее, линейная модель вида (4.47) справедлива лишь при отсутствии взаимодействия факторов. В противном случае этому взаимодействию как самостоятельному фактору присуща своя дисперсия
X2 A XB .
Вприведенном выше алгоритме при наличии взаимодействия между
факторами ХА и ХВ дисперсия X2 A XB , как составная часть, входит в дисперсию
189
ошибки 2 . Выделить и оценить эту дисперсию можно только при наличии параллельных наблюдений (таблица 4.4). При указанном расположении наблюдений их рассеяние в каждой ячейке относительно среднего той же ячейки обусловлено действием только случайных причин с дисперсией X2 A XB .
Рассеяние же самих средних в ячейках по всем возможным сочетаниям уровней факторов ХА и ХВ вокруг общего среднего, помимо фактора случайности, вызывается влиянием фактора взаимодействия ХА ХВ с
дисперсией X2 A XB . Кроме этих факторов на рассеяние средних по строкам
оказывает влияние только один фактор ХВ с дисперсией X2B , а на рассеяние
средних по столбцам – только один фактор ХА с дисперсией X2 A , так как все уровни другого фактора в каждом из этих случаев осреднены.
В этом случае для получения модели вида (4.46) общую сумму квадратов отклонений, наблюдений от общего среднего SSобщ. необходимо разложить уже на четыре компонента, отвечающие перечисленным выше факторам.
Таблица 4.6 – Двухфакторный ДА с одним наблюдением в ячейке
Компонен |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Математи |
|
Число |
|
|
|
|
|
|
Средний квадрат |
ческое |
||||||||||
ты |
степеней |
|
Сумма квадратов |
|
|
(оценка |
|
|
|
ожидание |
||||||||
дисперси |
|
|
|
|
|
|
||||||||||||
и |
свободы |
|
|
|
|
|
|
дисперсии) |
среднего |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
квадрата |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ХА |
k-1 |
|
|
SSXA |
Q2 Q4 |
|
S2XA |
|
SSX |
A |
|
|
2 |
2 |
||||
|
|
|
|
|
|
|
m XA |
|
||||||||||
|
|
|
k 1 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ХВ |
m-1 |
|
|
SSXB |
Q3 Q4 |
|
S2XB |
|
SSX |
B |
|
|
2 |
2 |
||||
|
|
|
|
|
|
|
k XB |
|
||||||||||
|
|
|
m 1 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Остаточн |
|
SS |
|
Q Q Q Q |
|
S2 |
|
|
|
|
SS |
|
|
|||||
(k-1)(m-1) |
|
4 |
m 1 k 1 |
|
2 |
|||||||||||||
ая |
|
1 |
2 |
3 |
|
|
|
|||||||||||
(ошибки) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Общая |
km-1 |
|
SSобщ. Q1 |
Q4; |
|
S |
2 |
|
|
SSобщ |
– |
|||||||
(полная) |
|
|
y |
|
|
|
|
|
|
|||||||||
|
mk 1 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
190
