Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 09.04.2024
Просмотров: 688
Скачиваний: 2
Линейные модели с |
одной |
переменной |
261 |
и, следовательно, |
|
|
|
x = с + |
dY*, |
|
(4.3.34) |
где x служит оценкой величины х, соответствующей заданному измерению величины У. Крючкоф считает, что метод обращения Эйзенхарта для оценивания величины х по данному значению Y*
. является более удовлетворительным.
Пример 4.3.1 (продолжение)
Доверительный интервал для ß^ вычисляется с помощью выра жения (4.3.18), и соотношение для Ь'0, аналогичное неравенству (4.3.20) (t = 2,306 для объединенной дисперсии с 3 + 5 = 8 степенями свободы), имеет вид
10,26 - 2,306-(0,825)1/2<ß;< ю,26 + 2,306-(0,82s)1 /2 ,
8 , 1 6 < ß ; < 12,36.
Подобным образом из неравенства (4.3.20) получаем доверитель ный интервал для ß 0
14,99 < ß 0 < 16,77.
Из неравенства (4.3.22) следует доверительный интервал для ßj: 0,1874 - 2,306 (7,5 - Ю - 4 ) 1 / 2 < ß t < 0,1874 + 2,306 (7,5-Ю"4 )*/2 , 0,124 < ß f < 0,251.
Заметим, что истинные значения ß^ и ß t (которые в данном случае известны) попадают внутрь доверительного интервала. Наконец, доверительный интервал для ті получается из неравенства (4.3.24):
Y - 2,306 sY < Tj < Y + 2,306 sf.
Дл я построения графика можно вычислить доверительные пределы
ввыбранных точках:
X = |
10 |
12,13 |
- |
1,55 |
< |
т| < |
12,13 |
+ |
1,55, |
|
X = |
20 |
14,01 |
- |
1,09 |
< |
т] < |
14,01 + |
1,09, |
||
X = |
30 |
15,88 |
- |
0,89 |
< |
т | < |
15,88 |
+ |
0,89, |
|
X |
= |
40 |
17,76 |
- |
1,09 |
< |
т] < |
17,76 |
+ |
1,09, |
X |
= |
50 |
19,63 |
- |
1,55 |
< |
ц < |
19,63 |
+ |
1,55. |
На фиг. П.4.3.1а нанесены значения У, и сравниваются между собой предполагаемая модель и линия регрессии, которая эту модель оценивает.
На фиг. П.4.3.16 изображена совместная доверительная область для ß 0 и ßi, определяемая уравнением (4.3.28). Отметим, что главные оси эллипса параллельны координатным осям, а оценки
Линейные |
модели с одной |
переменной |
263 |
Ь0 и &! попадают внутрь эллипса. Напротив, эллипс, ограничиваю щий совместную доверительную область для ßo и ßt , повернут отно сительно осей координат; кроме того, очевидно различие между совместной доверительной областью и отдельными доверитель ными интервалами, вычисленными для ß0 , ßt и ß0 .
Дисперсионный анализ по образцу табл. 4.3.3 приведен в табл. П.4.3.1в.
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица |
ПА.З.Ів |
||
|
|
|
|
|
|
|
|
Сумма |
Число |
|
|
|
|
|
||
|
Источник рассеяния |
|
|
степеней |
|
Средний квадрат |
||||||||||
|
|
|
квадратов |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
свободы V |
|
|
|
|
||
Обусловленный |
регрессией |
|
|
|
|
|
|
|
|
|
|
|||||
V |
Y^Pi |
|
|
|
|
|
252,81 |
|
1 |
|
252,81 |
|
|
|
||
ЬІ |
с поправкой на |
Ь 0 : |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
70,28 |
|
1 |
|
70,28 |
|
|
|
|
Отклонения относительно ли |
|
|
|
|
|
|
|
|
|
|||||||
нии регрессии: |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
5,02 |
|
3 |
|
1,67ч |
|
|
|
|
Отклонения |
внутри2 |
серий |
|
|
|
|
|
1 |
объединен- |
|||||||
2 |
>г(Уг-Уг) |
|
|
|
|
|
|
|
|
|
Г |
ный |
1,50 |
|||
(ошибка |
эксперимента): |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
6,95 |
|
5 |
|
1,39.1 |
|
|
|
|
22 |
Pti-Yt)* |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Общий |
|
|
|
|
|
|
335,06 |
10 |
|
|
|
|
|
|||
|
Гипотезу |
ßi |
0 можно |
проверить, |
составляя |
отношение дис- |
||||||||||
персий |
|
|
|
|
s§ |
70,28 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
= 46,8, |
|
|
|
|
|
|
||||
|
|
|
|
|
|
"Y, |
1,50 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
которое |
значительно |
|
превышает |
значение |
|
^о.эв (1 >5) = |
6,61 |
|||||||||
из табл. В 4; следовательно, гипотеза |
ß t = |
0 отвергается. |
Ясно, |
|||||||||||||
что величина |
|
отношения |
s^/syi =252,81/1,50 |
указывает |
на то, что |
|||||||||||
ß o |
является |
|
значимым |
параметром |
модели. |
Двусторонний |
t- |
|||||||||
критерий для нулевой гипотезы Н0 |
о том, что ß t |
= |
0 с 8 степенями |
|||||||||||||
свободы |
(из табл. В . З при а |
= 0,05, t = 2,306); |
|
|
|
|
|
|||||||||
|
|
|
|
|
, = |
і ^ 0 _ _ У 8 7 4 _ |
= |
|
|
|
|
|
|
|||
|
|
|
|
|
|
«Ы |
"1/7,5-10-* |
|
|
|
|
|
|
|||
также показывает, |
что гипотеза ß i = 0 |
должна |
быть отвергнута, |
|||||||||||||
так |
как |
6,85 > 2,306. |
|
|
|
|
|
|
|
|
|
|
||||
|
Интересно, что /^-критерий для гипотезы ßi = |
0 точно совпадает |
||||||||||||||
с ^-критерием |
для ß 4 |
= |
0, если s2 |
= |
ибо |
|
|
|
|
|
264 |
Глава |
4 |
а из оценки |
уравнения регрессии |
Yt = b0-j- bt (xt — Х), где b0 = Y, |
следует |
|
|
|
Y.PiÔ'i-Yf^b^iXi-xf, |
|
так что |
|
42 (ж*-*)2]1/2 |
|
|
|
|
|
•У, |
Вводя оценку дисперсии Ь4 , основанную на выражении (4.3.11), получим
4 V «ы '
Так как переменная F [1, (и — 2)] точно равна квадрату t [п — 2],
использование |
или f-критерия означает фактически одну и ту |
же проверку. |
|
В качестве иллюстрации можно также проверить гипотезу о том, что ßj = 0,150, применяя ^-критерий:
0,1874 — 0,150 |
_ 0,0374 |
„ 7 |
1/7,5-10-4 |
- 2 , 7 4 - 1 0 - 2 |
• |
Эту гипотезу следует принять, так как 1,37 < 2,306. Можно было бы выдвинуть и проверить и другие гипотезы как для ßg и ßi по отдельности, так и для обоих параметров вместе.
Пример 4.3.2. Простая линейная регрессия
В этом примере проводится анализ реальных эксперименталь ных данных (приведенных в табл. П.4.3.2а), которые собирались для обнаружения связи между показаниями х поплавкового рас ходомера и скоростью потока Y. Показания расходомера можно было отмечать с большой точностью и стабильно поддерживать на заданном уровне, так что х вполне можно считать детерминиро ванной переменной.
|
|
|
|
|
Таблица |
П.4.3.2а |
Набор |
У, |
смз/мин |
X, см |
Набор |
У, смЗ/мин |
X, см |
данных |
данных |
|||||
1 |
|
112 |
1,14 |
9 |
233 |
3,19 |
2 |
|
115 |
1,37 |
10 |
259 |
3,09 |
3 |
• |
152 |
1,89 |
11 |
287 |
3,05 |
4 |
|
199 |
2,09 |
12 |
240 |
3,10 |
5 |
|
161 |
2,45 |
13 |
281 |
3,34 |
6 |
|
209 |
2,04 |
14 |
311 |
3,75 |
7 |
|
237 |
2,73 |
15 |
392 |
4,19 |
8 |
|
231 |
3,04 |
16 |
357 |
4,59 |
|
|
Линейные |
модели |
с |
одной |
переменной |
265 |
|
|
Требуется: 1) оценить |
значения ß^ и ßi в предполагаемой моде |
||||||
ли У] = |
ßo + ßi^; 2) |
найти доверительные интервалы (для довери |
||||||
тельной |
вероятности |
1 — а = |
0,95) для |
ß^, ßt и и; 3) изобразить |
||||
на |
чертеже оценку линии регрессии, доверительные пределы для |
|||||||
и |
и экспериментальные |
данные; |
4) проверить гипотезу, |
отличен |
||||
ли от нуля угловой коэффициент |
линии |
регрессии. |
|
Решение
Этот пример отличается от предыдущего тем, что основан на подлинных экспериментальных данных, и тем, что никаких повторных измерений не производилось. По этой причине нельзя
450, |
. |
SO I |
1 |
1 |
1 |
1 |
I |
I |
0 |
1 |
2 |
3 |
4 |
5 |
|
|
|
Показания |
х,см |
|
|
|
Ф и г . П . 4 . 3 . 2 . Оценка линии регрессии, доверительные пределы и экспе риментальные данные.
" • экспериментальные данные", оценка линии регрессии;
— — — геометрическое место доверительных пределов для г\ (а = 0,05).
использовать F-критерий, чтобы определить, является ли прямая линия подходящей моделью"для подгонки экспериментальных дан ных. Однако эти данные можно нанести на график и визуально оце нить пригодность линейной модели, как показано на фиг. П.4.3.2