Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Статистический анализ временных рядов

..pdf
Скачиваний:
4
Добавлен:
19.11.2023
Размер:
47.02 Mб
Скачать

6.3. РАВНОМЕРНО НАИБОЛЕЕ МОЩНЫЕ КРИТЕРИИ 293

скую область

f9)

h

MQtlQo. . . . .

Qi- й vl2>)

 

 

 

 

u

^

MQ/ I Qo - . . . ,

Q/_ 1; y<.‘>)

 

 

 

 

 

 

= exp

 

,<2)

mj(Q0............ Q t_ й

yP )

 

 

 

m i (Qo-

• • •

> Qi—1> Y/2*)

 

 

 

 

 

Константа

k, зависящая

от

значений

Q0........ Qi-\,

определяется

так,

чтобы условная вероятность события (8 ) при

у, <= уР

равня­

лась

желаемому уровню

значимости. Если у*1*>

y f \ то

крити­

ческую область (9) можно записать в виде

 

 

 

(Ю)

 

0,>с,(<г0.......

Qt-i:

У('\У12\

 

 

 

В частности, такой вид будет иметь критерий для проверки гипоте­

зы у, = 0 против альтернативы у, = у‘2) <

0 .

 

 

Т е о р е м а 6 . 3 . 3 . Наилучший подобный

критерий

для

проверки

нулевой гипотезы у(- = у!1’ против альтернативы

уг =

у'-2) < у*1*

с уровнем значимости е; имеет критическую область (10), где кон­ станта ct (Q0, ..., Qt-i', y i\ y f’) определяется так, чтобы при у, = у*1' вероятность события (10), соответствующая плотности (8 ), была равна

Поскольку сс (Q0, ..., Qi_i; у'0, у‘2))

определяется из условной

плотности (8 ) при у( = y V \ то она не

зависит от значения у*2) и

поэтому ее можно записать в виде с( (Q0, ...,

Qi~\\ у;1}). Таким об­

разом, критическая область (10) не зависит

от у|2). Поэтому она

определяет критерий, являющийся наилучшим против любой аль­

тернативы

у!2»<

уГ.

 

 

Т е о р е м а

6.3.4.

Равномерно наиболее мощный подобный критерий

для проверки нулевой гипотезы у(- =

у,-1* против альтернатив уf <;

< у*1*с уровнем значимости 8,. имеет критическую область

(11)

 

Q,>c,(Qo.

•••.

Qi-б У^)-

Константа ct (Q0......... Qi-i', у!1*)

определяется так, чтобы при

у, = у*-1* вероятность события

(11), вычисленная в соответствии

с плотностью (8 ), равнялась е(-.

 

 

Нас будет особенно интересовать случай yll) = 0. В этом случае условная плотность Q, при нулевой гипотезе равна

< 1 2 >

......................%му

294

 

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Гл. 6.

Следствие 6.3.1.

Равномерно наиболее мощный подобный крите­

рий

для проверки

нулевой гипотезы yL= 0 против альтернатив

Т/ <

0 с уровнем значимости ei имеет критическую область

(13)

 

Q,>c,(Qo, . . . , Qi-i).

 

Константа ct (Q0,

..., Q,_i) определяется так,

чтобы интеграл

от плотности (12) помножеству (13) был равен е,-.

 

Когда мы перейдем к рассмотрению конкретных моделей, то уви­ дим, что обычно критическая область (13) имеет вид

(14)

 

г(- > С,(гх,

.... r<_i),

 

где /*! == Qi/Q0,

rt

= Q//Q0 — сериальные

коэффициенты кор­

реляции. В частности,

при i = 1

критерий для

проверки гипотезы

о независимости против альтернатив о положительной зависимости имеет критическую область гг > Сг.

6.3.2. Равномерно наиболее мощные несмещенные критерии

При проверке гипотезы о том, что порядок зависимости меньше i (т. е. y-L = 0 ), нас может интересовать альтернатива, состоящая в том, что порядок зависимости равен /(т. е. у-ьФ 0 ). Эта альтер­ натива двусторонняя, а для двусторонних альтернатив не существу­ ет равномерно наиболее мощного критерия. Например, наилучший

подобный критерий для проверки гипотезы у,- = у^ против альтер­

нативы у( = у|-2) < ур с уровнем значимости е, имеет критическую область (11). В то же время наилучший подобный критерий для про­

верки той же гипотезы, но при альтернативе у£. = у*2) > у/1* будет иметь критическую область

(15)

 

Q i < C i ( Q 0.......... Q i - 1; у;1’),

где с] (Q0, ....

Qf_ i; yl1’) определяется так, чтобы вероятность собы-

тия

(15), подсчитываемая согласно

плотности (8 ), была равна г1

при

уг = у (Д

При рассмотрении

двусторонних альтернатив по­

требуем, чтобы критерий был несмещенным, т. е. чтобы его функция

мощности достигала минимума при yt = у (/ ). В классе несмещен­ ных критериев с данным уровнем значимости найдем равномерно наиболее мощный критерий. Такой критерий действительно будет существовать, поскольку семейство плотностей экспоненциально, и область принятия гипотезы будет интервалом. Опишем вывод рав­ номерно наиболее мощного несмещенного критерия только в общих чертах. Полное и строгое доказательство приведено у Лемана (1959, гл. 4, разд. 4).

6.3.

РАВНОМЕРНО НАИБОЛЕЕ МОЩНЫЕ КРИТЕРИИ

295

Для заданных Q0, Q/_i пусть R будет критической областью с уровнем значимости et, т. е.

(16)

^ h{(Qi |Qo, • • •»Q*— у,- dQi = Е/.

Пусть соответствующей функцией мощности будет р (у,). При этом

р (у< ‘>) = е , и

<17>

Для того чтобы критерий был локально несмещенным, надо прирав­ нять нулю первую производную функции мощности в точке yt =

= у!-1*. Дифференцируя под знаком интеграла, получаем

(18) О - ^ Р Ы !

 

 

 

 

kj (Q0>

. .

• , Qi)

 

exp [—

) dQi +

 

 

 

 

mi (Q0, . . .

,

Y,(I))

 

 

dyi

[ mi (Q0, . ..,

 

Qe-—!;

yt)

 

Yt*=V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X ^

 

(Qo»

• * •

f

Qi)

f

 

 

 

 

 

 

R

 

 

 

 

'

 

 

— —

2~

QA (Q/1 Qo»

• • • » Q*«-i»Y* *) dQi +

 

+

(Q0,

... ,

Qi—i;

у!1*)

 

 

mi(Q0,

, Q(_,; yi)

v ^ 1»-

Дифференцируя

подобным же образом

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

(19)

 

 

1 = 1

У 0’

Q<)

 

e~y‘Q‘/2dQ,

 

v

 

 

 

J

(Qo..........Qi_i; уi)

 

 

 

 

 

 

 

 

— oo

 

 

 

 

 

 

 

 

в точке yi = y il)

и

сравнивая

результат

с (18),

получаем

соотно­

шение

 

 

 

 

 

 

 

 

 

 

 

 

(20) JQA(Q,|Q0, .... Q,_i;

y!1,)dQ/ =

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

oo

= et j QA (Qi I Qo. — Qi—it Y*1’) dQi-

296

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Гл. 6.

Таким образом, (2 0 ) является условием локальной несмещенности критерия.

Мы хотим теперь найти такую критическую область R, чтобы она удовлетворяла соотношениям (16) и (2 0 ) и, кроме того, макси­ мизировала функцию р (у;) для некоторого частного значения у19

например для yt =yi-2) > у\1К В силу обобщенной фундаменталь­ ной леммы Неймана — Пирсона (упр. 12) такая область определя­ ется неравенством

(2 1 ) MQilQo, . . . .

Qi-n v'2)) >

> cthi (Qt | Q0,

. . . , Qi—j', yl -f- bQfti (Qi I Qo> • • • > Qi—ii Y<' *)>

если существуют функции а и b от Q0, .... Q,-i, такие, что указан­

ная область удовлетворяет (16) и (2 0 ). Неравенство (2 1 ) равносильно

неравенству g (Qf) >

0 , где

 

 

 

 

 

 

(22)

 

 

 

 

 

+

BQi),

 

 

и

 

 

 

 

 

 

 

 

 

(23) А «

 

. Qi-i> v‘2))

в

mi (Q0.

• • • .

Q i-v

Vi2>)

mi(Q0, .

. Qi-ь

у!11) ’

= b

 

 

 

 

 

 

mi (Q0,

. . . .

Q t_ t;

yj'b

Поэтому

 

 

 

 

 

 

 

 

(24)

 

g(Qi) = —

k~(Y*2>

»t*h

 

,<2 > .

-Y(/>

 

в.

 

dQi

у*1*) exP

- T W

 

 

 

 

 

 

 

) Q ] ~

 

Если В > 0 (и Y;2) >

уГ*), т 0

последняя производная отрицательна

и неравенство (2 1 ) определяет некоторый полуинтервал (— оо, с). Однако такая критическая область не удовлетворяет соотношениям (16) и (20). (См. упр. 13.) Поэтому величину В следует брать отри­ цательной. Указанная производная будет непрерывно возрастающей функцией, принимающей нулевое значение только в одной точке. Эта точка будет точкой минимума функции g (Qt). Если этот мини­ мум отрицателен (что имеет место, когда А достаточно велико), то (21) определяет два полуинтервала (— сю, с) и (d, сю). Область принятия в таком случае равна (с, d). Отметим, что значения c u d

определяются соотношениями (16) и (2 0 ), которые не зависят от у*21.

Подобным же образом при у}2> < у*'* можно доказать, что констан­ та В должна быть положительной и что область принятия будет ин­

тервалом, удовлетворяющим (16) и (2 0 ) независимо от значения у*2). Таким образом, мы получили равномерно наиболее мощный несме­ щенный критерий.

Теорема 6.3.5. Равномерно наиболее мощный критерий для про­ верки нулевой гипотезы у, = у!-1* против альтернативы уг =£ у!1’

6.3. РАВНОМЕРНО НАИБОЛЕЕ МОЩНЫЕ КРИТЕРИИ 297

с уровнем значимости et имеет критическую область

(25)

 

Qi ^

(Q0> • • *> Qi—li

У( ^)

 

 

Qi ^

СLi (Q0» • • • » Qi—lf

Тi ^)*

 

 

 

 

где Су. (Qo,

•••,

Qi-i;

Т;0) и cLi (Q0, ....

Qi_i; •у|1))

определяются

так, чтобы

 

 

 

 

 

.... Qi—V Г]1*)

 

 

 

(26)

J

A, (Qt | Qo. • • •, Qi—и Y.!1))dQt =

l —e,

CLSQO.... Qt—i; vj1))

 

 

 

и

 

 

 

 

 

cu^Qo,...»Qt—i; vp)

 

 

 

(27)

j

 

QA (Q, I Qo. . . . . Qi—i; yi{))dQt =

е/.г(в».... v*l))

oo

=( i — 8,) J QA (Q/1 Qo. . . . . Qi—i: TiVQi-

Следствие 6.3.2. Равномерно наиболее мощный критерий для про­

верки гипотезы у, = 0 против альтернатив yt Ф 0 с уровнем зна­ чимости г, имеет критическую область

 

 

Qi

сиi (Qo. . . . .

Qi—О»

(28)

 

Q, <

CLt (Q0, . . . .

Qi—I),

гЗе су. (Qo,

..., Qi_i) и cLi (Q0,

....

Q.-i)

определяются соотноше­

ниями

cu.iQo

Qi—1)

 

 

 

 

 

 

 

 

 

 

 

^OQ\

f

 

^ ^Q> • • •

Q^

dQi = 1 — 8;

 

J

^-i (Qo»

...»

Qt-1)

a/Qo....Qf-i)

и

cu^Qo Qi—l)

<3o> J

oi^Qo Qt—l )

Если условная плотность распределения величины Qt при за­

данных Qq, ..., Qi_ 1 для Y/ = 0 симметрична, то (30) автоматически

298

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Гл. 6.

будет выполняться для

 

(31)

Ci. (Q0, ..., Qi—i) = Cut (Qo,

• **»Qi—1)>

поскольку тогда обе части (30) будут равны нулю. При этом (29) можно записать в виде

(32)

(Qpi

Qf)

—i (Q0.

i)

s

q/((Qo....Qt—i)

 

 

Равномерно наиболее мощный односторонний подобный критерий можно называть также равномерно наиболее мощным односторон­ ним несмещенным критерием, поскольку из несмещенности крите­ рия вытекает подобие критической области.

Следует отметить, что условное распределение величины Q, при заданных Q0, ..., Q,_i не зависит от значений параметров у0, ...

..., Yi_i. Ввиду этого последним можно приписать любые значения, которые окажутся более удобными. Например, можно взять у0 =

= 1, Vi = -

= Т«—1 = 0 . т

Обычно Q0 будет равно 2 У/- Отношение

 

<=I

(33)

г , - %

можно рассматривать как сериальный коэффициент корреляции.

Если Yi = ... =Y(

= 0 , то совместное распределение гь ..., г,

не зависит от Y0 (>

0) и Q0. (См. теорему 6.7.2.) Оптимальные кри­

терии (13) и (28) могут быть определены с использованием условного распределения г, = QJQa при заданных гъ ..., г,_\.

Частную сериальную корреляцию порядка i при неизменных сериальных корреляциях порядков 1, ..., t — 1 можно определить обычным способом. Мы уже говорили о том, что этот частный се­ риальный коэффициент корреляции может быть использован для проверки гипотезы у* = 0. Поскольку он является линейной функ­

цией аргумента г, (с коэффициентами, зависящими от г1( ...,

r,_i),

то использование его при заданных значениях гъ

..., rt_i

экви­

валентно использованию самого гь при тех же гъ

..., п_ь

(См.

разд. 6.9.2.)

 

 

6.3.3. Случаи, когда равномерно наиболее мощного критерия не существует

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

(34)

= y P = Q

6 .3 . РАВНОМЕРНО НАИБОЛЕЕ МОЩНЫЕ КРИТЕРИИ 299

против альтернативы уг = уи ..., ур = у°Р> где уи • ••, уР задан­ ные числа, то равномерно наиболее мощный подобный критерий мог бы иметь критическую область

(3 5 )

2 Y / Q / < C ( Q O),

 

/=1

где константа с (Q0) подбирается так, чтобы условная вероятность события (35) при нулевой гипотезе была равна е. Если альтернатив­

ные значения параметров Yi,

ур имеют вид

(3 6 )

y, = ky),

/ = 1, . . . . р,

где у] фиксированы, то критическая область (35), с у0/ замененными

на у/, будет наилучшей для всех альтернатив k > 0 . Однако, если альтернативные значения не лежат на прямой (36) в пространстве параметров, область (35) не будет одной и той же для всех допусти­

мых альтернатив (у?, ..., у£). Но это и означает, что в данном слу­ чае не существует равномерно наиболее мощного критерия.

Рассмотрим стационарный случайный процесс, определяемый соотношением

(37)

yt + P*/*--i =

t — ...

f — 1, 0, 1, »* *,

в котором | Р | < 1, а каждое ut распределено нормально, с нулевым средним и дисперсией а2, и не зависит от y t-1, yt- 2, ... и щ -\у щ^2, ... . При этом маргинальное распределение выборки уъ ...

..., ут задается выражением (1), в котором

(38)

А0 I,

0 1 0 . . . 0 0 -

10 1 ... о о

!0 1 о ... о о

Ах F

 

:

• :

 

: •

 

 

0

 

0 0

 

. 0 1

 

 

0

 

0 0

... 1 0

1

о

0

0

.,..

о

1 0

 

0

1

0 . ..

0

0

 

о о

1 ..

о

0

 

. •

 

 

 

......

 

..

 

о

оо ..

i

0

1

о

0

0 . ..

0

I 0

300

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Гл. 6.

Множество альтернатив для ух = 2 р/о2 и у2 = р22

не лежит на

прямой. Поэтому здесь не существует равномерно наиболее мощно­ го критерия даже для односторонних альтернатив.

Если | Р | мало (т. е. близко к нулю), этот случай подобен моде­ ли разд. 6.5.4. Если р близко к —1, такой случай подобен мо­ дели разд. 6.5.3. Следовательно, можно предположить, что для проверки гипотезы р = 0 в этой модели (когда среднее значение известно) полезной статистикой должны оказаться одна или две

соответствующие корреляции. Положительная

зависимость соответ­

ствует р < 0 .

t = 1, ..., Т при ус­

Пусть (37) определяет распределение для

ловии у 0 == 0. Тогда А„ и Ах будут такими же, как и в предыдущем

случае, а

 

 

 

 

_

“ 1

0

0

. , . .

0

I о

 

0

1

0

. , . .

0

0

О

О

1 . . . .

0

0

.

 

 

 

. .

 

 

 

 

 

0

0

0 . . .

1

_ 0

0

0 . . .

0

о о

1

Множество альтернатив для ух = 2р/а2 и у2 = Р22 не лежит на прямой, так что и здесь для проверки гипотезы Р = 0 не суще­ ствует равномерно наиболее мощного критерия.

Модель, определяемая разностным уравнением второго порядка

(40)

ijt + рху*—

+ Рг№—2 =

Щ, t — 1, •••> Т,

с циклическим определением у0 =

ут, у~\ =

у т - 1, приводит к плот-

ности (1) с

А0 =

I,

 

 

 

 

 

 

 

 

 

 

 

 

(

( .

 

 

 

 

 

“ 0

1

0

0

г

 

 

 

 

1

0

1

• • .

0

0

 

 

А, =

1

0

1

0

 

.

0

0

 

 

 

 

 

 

 

;

 

 

 

 

0

0

0

• 1. 0

1

 

 

 

_

1

0

0

. .

1

0 _

(41)

 

0

0

1

0 .. . .

1

0

 

 

 

0

0

0

1

. . . 0

 

1

 

 

 

1

0

0

0

.. . . 0

0

 

а , =

- г

0

1

0

0

.. . . 0

0

 

 

 

 

 

 

 

 

 

 

 

 

1

0

0

0

. . . 0

0

 

 

 

0

1

0

0

.. . . 0

0

6 .4 .

ВЫБОР ПОРЯДКА ЗАВИСИМОСТИ

301

Если альтернативные гипотезы не связывают специальным обра­ зом Рх и р2, то равномерно наиболее мощного критерия для проверки гипотезы рх = р2 = 0 не существует.

Теория, изложенная в этом параграфе, была развита Т. Андер­ соном (1948), (1963).

6.4.ВЫБОР ПОРЯДКА ЗАВИСИМОСТИ КАК ЗАДАЧА СО МНОГИМИ РЕШЕНИЯМИ

Займемся теперь вопросом об определении порядка зависимости на основании наблюдений уъ ..., ут. Будем предполагать, что ис­ следователь может определить некоторый наименьший пг и неко­ торый наибольший р — возможные порядки зависимости. Во мно­ гих случаях m оказывается равным нулю (указывая тем самым на независимость). В других случаях исследователь может в силу ка­ ких-нибудь причин ожидать, что порядок зависимости не меньше некоторого положительного числа т. Ввиду того что имеющиеся наблюдения удалены друг от друга не более чем на Т — 1 единиц времени, ясно, что максимальный порядок р не может быть больше чем Т — 1. Однако ясно также, что для того, чтобы процедура об­ ладала хорошими свойствами в отношении мощности, максималь­ ный порядок р должен быть выбран по возможности много мень­ шим Т — 1.

Исследователь имеет, таким образом, задачу со многими реше­ ниями. Именно, ему необходимо выбрать порядок m, т + 1, ...

..., р — 1 или р. Эта задача подобна задаче выбора степени поли­ нома, представляющего тренд, изучавшейся в разд. 3.2.2. Хотя распределения здесь уже другие, общая структура задачи и ее ре­ шение будут теми же, что и в разд. 3.2.2. Мы формализуем задачу, стоящую перед исследователем, считая, что ему необходимо решить, какому из перечисленных ниже взаимно не пересекающихся мно­ жеств принадлежит параметрическая точка (ym+ i, • • • , ур):

 

Нр'- УРФ0.

Ур—1Ф0,

 

 

Hp—i Ур —0,

 

(0

:

1• Ур— • ' '

— Ym+2 — 0,

Тт+1

 

 

 

Нт '■ УР = ' ' '

= Y»«-H = 0-

 

Принадлежность

параметрической

точки

множеству #,■ означа­

ет, что

порядок

зависимости

равен

i. Альтернативная формули­

ровка состоит в том, что исследователю необходимо решить, вер­ на ли хотя бы одна и (если да) какая именно из следующих нулевых

302

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Г л . 6.

гипотез:

Н'р’-Ъ =

Нр.р—1 • Yр — Ур^-1 — 0 *

(2)

;

Нр,р—1 т+ 1 '• Ур — • • • = Y^+i = 0*

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

(3)

Нр,р—1, .... m+l CZ

• • • CZ Нр.

 

Эти два семейства множеств связаны соотношениями

 

я ; =

Яр-1

и

•••

U я т ,

(4)

Нр , р 1 =

Я р —2

U

• • •

U Я т ,

Яр.р—1...... т+1 — Я т .

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

рр= Рг (отвергнуть Н*р| Н*р),

рр+ рр- 1 = Рг (отвергнуть Яр,р_!

(5)

:

 

 

Рр +

• • • + рт+1 =

Рг (отвергнуть

Я р ,р _ ,......т+1 | Яр,р_ 1 ........т+1 } =

 

=

Рг (отвергнуть

Нт \Нт),

причем Pi > 0 и Рр + ... + рт+1 < 1. Поскольку здесь каждая нулевая гипотеза включает в себя последующую (т. е. каждая по­ следующая гипотеза является более сильной), вероятности откло­ нения берутся монотонно неубывающими (т. е. вероятность откло­ нения более сильной нулевой гипотезы, когда она верна, не меньше аналогичной вероятности для менее сильной гипотезы). С помощью взаимно непересекающихся множеств, указанное разнесение по уровням значимости записывается в виде