Файл: Химмельблау Д. Анализ процессов статистическими методами.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). Отметим, что главные оси эллипса параллельны координатным осям, а оценки


21

20

 

 

19

 

 

18

 

 

17

-

 

16

 

 

15

 

 

14

 

 

13

 

'////

12 ~

11 "

/

/

10 /

 

,

 

 

 

/ /

 

2

/

V

 

/

//

У

 

/ У

 

/

7

V

A

 

/

У

/

 

/

/•Значения f£

/г/

20

1

!

1

 

30

АО

50

60

Ф и г . П И . З . І а . Модель,

оценка

линии регрессии

и

значения

\ ~т»™1£^™™

М , 0 Д е л ь 1 1 =

+ ß.oc;

2 - оценка линии

регрессии

У = Ъ

3 - геометрическое

место точек

верхних

доверительных

пределов

для

тг 4 —

рическое место точек нижних

доверительных

пределов для т).

8,0

8,5

9,0

9,5

 

 

Щ5

 

110

w\\ '

' ' ч • ' І ІТ'

' ' ч ' ' ' ' ч

 

i <

T '

»

 

• • м . .

030

 

 

 

 

 

 

 

 

020

 

 

 

 

 

 

 

 

О,10

 

Модель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

? " Л ' М а

 

 

 

 

 

 

О

J

I L

15,0

-I

1

L

 

 

 

 

 

 

 

 

 

 

 

16,0

 

 

 

 

 

A

 

 

 

 

 

Ф л г. П . 4 . 3 . 1 6 . ГПЯНИТГЫ пиа

оно/


Линейные

модели с одной

переменной

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 можно

проверить,

составляя

отношение дис-

персий

 

 

 

 

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