Научная статья на тему 'Методы решения задачи непараметрической проверки гипотез Лемана при испытаниях параллельных систем'

Методы решения задачи непараметрической проверки гипотез Лемана при испытаниях параллельных систем Текст научной статьи по специальности «Математика»

CC BY
71
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕПАРАМЕТРИЧЕСКАЯ СТАТИСТИКА / ГИПОТЕЗА ЛЕМАНА / КРИТЕРИЙ КОЛМОГОРОВА-СМИРНОВА / ОЦЕНКА КАПЛАНА-МЕЙЕРА / NONPARAMETRIC STATISTICS / LEHMANN'S HYPOTHESIS / KOLMOGOROV-SMIRNOV TYPE CRITERION / KAPLAN MEIER ESTIMATE

Аннотация научной статьи по математике, автор научной работы — Тимонин Владимир Иванович, Тянникова Нина Дмитриевна

Рассмотрена задача проверки степенной гипотезы Лемана для двух цензурированных выборок. Для проверки степенной гипотезы разработан критерий типа Колмогорова Смирнова, основанный на сравнении оценок типа Каплана Мейера функций распределения по каждой цензурированной выборке. На основе модели случайного блуждания частицы по целочисленной решетке описан метод вычисления точных распределений статистики. Вычислены значения вероятностей для широкого набора возможных объемов выборок. Доказана сходимость распределения данной статистики к стандартному распределению Колмогорова Смирнова при условии справедливости проверяемой гипотезы. Методами статистического моделирования исследованы свойства оценки степенного параметра, получаемой минимизацией статистики.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Тимонин Владимир Иванович, Тянникова Нина Дмитриевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Methods for solving the problem of non-parametric testing of Lehmann's hypotheses when testing parallel systems

The article considers the problem of testing the Lehmann power hypothesis for two censored samples. The Kolmogorov Smirnov type criterion based on a comparison of the Kaplan Meier type estimates of the distribution functions for each censored sample is developed to test the power hypothesis. A method for calculating the exact statistics distributions is described on the basis of the model of particle random walk over an integer lattice. The probability values are calculated for a wide range of possible sample sizes. The convergence of this statistics distribution to the standard Kolmogorov Smirnov distribution is proved provided that the hypothesis being tested is valid. The properties of a power parameter estimate obtained by minimizing statistics are investigated by statistical modeling methods.

Текст научной работы на тему «Методы решения задачи непараметрической проверки гипотез Лемана при испытаниях параллельных систем»

УДК 519.248

DOI: 10.18698/2309-3684-2018-0-000

Методы решения задачи непараметрической проверки гипотез Лемана при испытаниях параллельных систем

© В.И. Тимонин, Н.Д. Тянникова МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

Рассмотрена задача проверки степенной гипотезы Лемана для двух цензурированных выборок. Для проверки степенной гипотезы разработан критерий типа Колмогорова - Смирнова, основанный на сравнении оценок типа Каплана -Мейера функций распределения по каждой цензурированной выборке. На основе модели случайного блуждания частицы по целочисленной решетке описан метод вычисления точных распределений статистики. Вычислены значения вероятностей для широкого набора возможных объемов выборок. Доказана сходимость распределения данной статистики к стандартному распределению Колмогорова - Смирнова при условии справедливости проверяемой гипотезы. Методами статистического моделирования исследованы свойства оценки степенного параметра, получаемой минимизацией статистики.

Ключевые слова: непараметрическая статистика, гипотеза Лемана, критерий Колмогорова - Смирнова, оценка Каплана - Мейера

Введение При испытаниях технических систем часто возникает задача проверки степенной гипотезы

Д0:ЖО = №(0)*, к> 1, (1)

где ^(0, ^2(0 — функции распределения для первой и второй независимых выборок соответственно [1-3]. Однако не всегда при испытаниях наблюдаются все отказы элементов, составляющих сложные технические системы. В настоящей работе решена задача проверки степенной гипотезы (1) в случае, когда данные являются цензуриро-ванными [4, 5].

Постановка задачи На испытания ставятся 1\ систем, состоящих из Щ параллельно соединенных элементов, которые работают в режиме 81, и /?2 систем, состоящих из /И2 параллельно соединенных элементов, работающих в режиме 82. Испытания проходят таким образом, что наблюдается только отказ последнего входящего в систему элемента, а наработки до отказа остальных элементов неизвестны. В результате имеются две цензурированные справа выборки ©1 =(е1,...,е?), ©2 =(02,...,022), где е{ = шах{&£2,...,Хщ},

/= 1,11, и 02 = шах< X/, Х2,.., £¿,1, -/= 1,л2, — максимумы наработок

до отказа элементов систем, работающих в режимах в! и 82 соответственно.

Рассмотрим задачу проверки степенной зависимости двух функций распределения наработок до отказа элементов на основе данных

об отказах систем, составленных из этих элементов.

Пусть ^(¿) — функция распределения наработок до отказа элементов в режиме 81з а ^2(0 — в режиме 82. Тогда проверяемая гипотеза имеет вид (1).

В настоящей работе предлагается критерий типа Колмогорова -Смирнова, основанный на сравнении оценок типа Каплана - Мейера

[6] /^(0, ^э2(0 функций распределения ^2(0 [7] по цензури-

рованным выборкам ©1 и ©2:

0,

пгЩ(

п

1=1

1,

1 -

1

<0(0 = о;

1 ^(О^у-1^=1, 2,

^(0=ПР

где — количество отказов систем выборки ©J к моменту вре-

мени £ Построение и свойства таких оценок по аналогии с классическими оценками Каплана — Мейера функций надежности были рассмотрены в работе [7].

При допущении, что ©1, ©2 — независимые полные выборки из отказов систем, их функции распределения можно оценить эмпири-

с! (/) б/ (/)

ческими функциями распределения г(1) = ——, Р2(() = ——. Для проверки гипотезы (1) предлагается статистика вида

Т=

_тах_

Ш2

Щ

к2 +Цг2)

КС И Ш

рл^

А г

где р = 1; л2 П2

п2

, ч ¿/ль 2

-т 1

+ к1

(2)

-л =

рп^к2 +п2' рп^к2 +п2

"2 к

щ

При этом для случая — -1 < 0 и к2уР } + ) = 0 по-

А^Р +к2{р2)кт2 Лк'1

ложим-----= 0.

щ

Статистика (2) является обобщением предложенной в работе [7] статистики типа Колмогорова - Смирнова, которая применяется для проверки гипотезы однородности, что эквивалентно случаю к = 1. Отметим также, что при щ = Щ = 1, к = 1 вид данной статистики

совпадает с видом классической статистики Смирнова.

Асимптотическое распределение Т. Покажем, что вид статистики Г связан с тем, что ее асимптотическое распределение может быть приближено классическим распределением Колмогорова -Смирнова [8].

Без ограничения общности будем считать, что 1\(() = .^(О = = 0<¿<1. Рассмотрим процесс Л^) = . В рабо-

те [7] для асимптотической ковариации процесса была доказа-

на следующая теорема.

Теорема 1. Пусть 0 < А < 5< ^ < 1, тогда асимптотическая кова-

риация процесса Х1^)

1 -Г1

К^) = ЕХ^) Х^) -ЕХ^) ЩЩ & -ф-

Рассмотрим вспомогательный процесс Х2 (¿) = -

0 <*<1.

Вывод асимптотической ковариации процесса Х2 полностью аналогичен выводу асимптотической ковариации Л^ (■£, процесса -^1(0 [7], поэтому приведем без доказательства вид асимптотической ковариации этого процесса:

К2(^) = ЕХ2(З) Х2(1) -ЕХ2(З) ЕХ2(().

Теорема 2. Пусть 0 < А < 5< ^< 1, тогда асимптотическая кова-риация процесса А^(^)

А2 (л, О = ЕХ2^) Х2(Ъ - ¿ЭД ЕХ2 (?)

^/к^/к

П2

1 -Г2

¡к

П.

Для вывода асимптотической ковариации процесса -^2(0 = = ^ ^Ь2(?)) ,0 <?<1, воспользуемся теоремой, доказанной в работе [9].

Теорема Зь Пусть Д,(лт) = Гп

л/л.

^ Л

1—

Л(л-) = А? * л,

лге Д0,1], Д0,1] — пространство функций без разрывов второго рода на [0,1]. Тогда последовательность Лл(л) сходится в метрике Скорохода к Ь^х) равномерно на ограниченных множествах из Д0,1].

Поскольку асимптотически ^2 (0 = -^(-^2(0), нетрудно получить вид функции асимптотической ковариации процесса Л^ (?) :

к2 (5Д) = ¿^ М Л2 (0 - (О-

¡¡2

^(1 -Г2

Рассмотрим случайный процесс 2п(р) = ^-(-^О) 0 <А<?< 1, который определяет статистику (2). Очевидно, что

ад ((^(о -*)-(( -'Т)=*1(о -ТР^о, 0 <А<?< 1.

Теорема 4. Пусть 0 < А < 5< ? < 1, тогда асимптотическая кова-риация процесса

а

П2

= £5

1 ^ р^(1 -Г* ^

щ2Л + /Ир^^

= 5

1 1

Л-1

+

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

р^(1 -Л

®2

ш2м

-1

Введем (?) =

ЩЩ

-2д(?). При стандартных ограничени-

ях 0 <А<5<?< 1 процесс !^(?) сходится к гауссовскому процессу с нулевым математическим ожиданием и ковариацией

^ -1 + щ -1

ф-

щ

щ

щ

= 5

+к1

щ

Г

Для перехода к асимптотическому распределению статистики (2) предварительно докажем две леммы.

Лемма 1. Существует монотонно убывающее преобразование ¿=у(т): [0,1] ^ [1,0], при котором процесс Ш(т) = К(у(т))х

щ. -1

х--является броуновским мостом.

Доказательство. Рассмотрим преобразование

щ

-1щ

Щ

х(1)= ^ +к1 : (0,1]^(1,0].

щ

щ

к21к +гкх

При необходимости доопределим функцию т(^) в точке 0: т(0) = 1. Нетрудно показать, что т'(^)< 0, т(0) = 1, т(1) = 0. Тогда для т(^) существует обратное преобразование £ = у(т).

Введем в рассмотрение процесс

щ -1

У(т) к

^(т)=У(у(т))

Щ

щ

= ^МФМт)).

При некотором А1: 0 <н<г<А1 < 1, у(и) = у(к) = £ имеем

Е[ Ш(т) ] = 0;

Щ -1 Щ -1 Щ Щ

5а М кЛк +к -гк

-5-2 1

Щ Щ Щ 1

-щ , - / -1

к^к + к1к^к +к1 Г

■ = и(1 - V). <

щ щ щ

5* к^к +к1 -¡к

ЛЬ ЛЬ

—2-Щ —2-щ

к2Бк +к1 к^к +к1

Рассмотрим выборки ^б^..., б^ и ^б2,..., б2^). Если считать их независимыми полными выборками из совокупностей с функциями распределения !?(() = ?Ш1, ^2(?) = ^ к соответственно, то справедливо следующее утверждение.

Лемма 2. Для эмпирических функций распределения =

= 4(0 р2ц) = ^

А П2

Г

7 sup V (

(

(РщУщ ^ 0 ] = 1; ^ sup (Р2(1))к/т2 -/

V '

^ 0 1 = 1.

п2 )

Доказательство. Утверждение леммы очевидно в силу теоремы Гливенко и равномерной непрерывности на [0, 1] степенных функций с положительным показателем степени. <

Асимптотическое распределение (2) приближенно может быть получено следующим образом. Для функции 1 (у(т)) = 1 (?), определенной в лемме 1, множители ? должны быть заменены на сходящиеся к ним множители Р1({)1/т1 или Р2^)кт2. Тогда приближенно можно считать, что распределение статистики (2) совпадает с распределением максимума модуля броуновского моста, которое является классическим распределением Колмогорова - Смирнова [8].

Точные распределения статистики Т. Для вычисления точных распределений статистики (2) составим объединенную выборку из

отказов систем в обоих режимах: 0= |б1,..., б!1,б2,..., б2^|. Расположим элементы выборки ф в порядке убывания и составим из них ряд Г = {У1 >У2 >Уз ^ .. ^Ул^}.

Введем вектор 2 = ^^2,..., £Л[+/)2), состоящий из 1\ единиц и щ нулей, где

1, если у/ = бъ i = 1, /1, ,

Щ=\ 1 J — 1= 1, ...,А

[0, если у1 = q2,1 = 1,12,

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

Предварительно запишем статистику (2) в виде Т = max /(f), где

t

/(f) определяется видом статистики.

Пусть (i,j),0< i< ij; 0<j<i2, — целочисленная решетка на плоскости (рис. 1). Частица на первом шаге выходит из точки (0, 0) и на (1\ + л^-м шаге заканчивает блуждание в точке (т^,/^), совершая ]\ скачков «вправо» и 12 скачков «вверх». Траектории ш частицы будут находиться во взаимно однозначном соответствии с векторами Z. Равенство Zj = 1, 1= 1, + в векторе Z соответствует

скачку «вправо» на 7-м шаге, а Zy = 0, /= 1,...+12), — скачку «вверх». При прохождении траектории блуждающей точки через точку (i,J) функция /(f) принимает значение Т^ :

К

(

щ-1 А у

М/Л!

(

+ к2

*h - J

_1

.\frm2 Л к

У

У

1J ^ к2рлг2 +Ш2

К

г „ л1/л* (

+ к2

щ-1

v А у

Ъ-J

щ.

лк/щ ^ук-101

-Л,

V п2 У

+¿1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где Лу =

П

-%=1

1-

щ(п1 - S1 + ^

ль

Л ( i i

п

V %=1

1-

1

(^2 + ^

Шг

"2\"2 т -,jj ( .Л1^ ( л^щ

При этом для случая —— -1 < 0 и ^

hZ1

V А У

+ к

Щ -J

=0

К

Г« ,Л1/й1 (

V "1 У

+ к2

П2 -J

V л2 У

V п2 У щ 1

сохраняется условие

К

Г „ Л1/л1 (

+ к2

v Ч у

V п2 у

Щ

лк/щ ^Т""1

= 0.

2

1

2

У "2

О щ г

Рис 1. Случайное блуждание частицы по целочисленной решетке

Отметим, что значение Ту одинаково для любой траектории, проходящей через точку (/,/).

Теорема 5. Вероятности векторов Z = ^гъ г^, • ••, ^ +а^ не зависят от распределения отказов и определяются следующим выражением:

К ) ы (4-Ц-^щк^щ-и-^щ

где = (Р0 = 0) — количество единиц в векторе Ъ до 1-го места ы

1

включительно; и1 = (и0 = 0) — количество нулей в векторе Ъ

1=1

до 1-го места включительно.

Доказательство. Утверждение теоремы следует из вида условных вероятностей перехода «вверх» или «вправо» в схеме случайного блуждания частицы [10, 11]. <

Теорема 6. Вероятность Р(Т <Н) равна величине р^ ^ (А), которую можно получить повторным применением соотношения

Мь)=

1= 0;У = 0;

щк^ц -1+1)

щк^ -1+1) +щ(п2-]) С щк^ц -/+1)

п,-

^ЛЬНАЯ 1<1<П1^=0;

щ^ц -/+1) +щ (п2 -у)

+

щк(щ-/)+тв,(л.

) +

Х9(Ь),1<1<П1;1^<щ.

(3)

Здесь гу(И) =

1, Ту<И, 0, Ту>И.

Доказательство. Вероятность каждой траектории ш можно записать в следующем виде:

Щ+П2

Р(2) = Р(ш)= П

ы

((д -У-^щ^Цъ -и-^щ) ^-^Щк+ ("2 -и1-1)щ>

1-г}Л

Ц+П2

П мш.

1=1

Пусть Шу — множество «частичных» траекторий, начинающихся в точке (0, 0) и оканчивающихся в точке (/, у) (соответствующие векторы 7. имеют / единиц и ] нулей на 7-м месте, где 1=1+ /).

Обозначим Ру(ш) = ^ Х5 (ш).

5=1

Вероятность любой траектории, совершающей скачок ( 1-1, у) ^ /,у), что соответствует г; = 1, имеет множитель ^(Ш:

„ щк(ц -1+ ^

-77-Г^-7-л.

4 7 щк(п1 -1+1) +1Щ{П2-у)

Если же происходит скачок (/.у-1) 1, у), что соответствует z, = 0, то

1](&) =

Щ(П2 + ^

щк(п1 -1)+Ш2(п2 -у' + ^

1

Пусть %у = ^ Ру- Тогда соотношение (3) следует из того, что

ту

в (7, у) за один скачок можно попасть только из точки (1-1, у) или из

точки (/.у-1). Множители обеспечивают обращение в нуль

вероятностей траекторий, для которых Ту>Ь. <

Выбор вида статистики (2) обусловлен тем, что ее предельное распределение при 1\ щ не зависит от количества щ, щ

элементов систем и является классическим распределением Колмогорова - Смирнова.

В таблице приведены вероятности точного распределения статистики Т, рассчитанные для квантилей /? = 1,22, И = 1,36 и А = 1,63, которые являются соответственно квантилями уровней 0,8981, 0,9505 и 0,9901 асимптотического распределения Колмогорова - Смирнова.

Точные вероятности Р(Т < й) для различных квантилей h в случае равных объемов выборок при ш\ = ш% = 2

Л1 -П2 А = 1,22 А = 1,36 А = 1,63

к = 1,5 к = 3 к = 1,5 к = 3 к = 1,5 к= 3

10 0,9040 0,7991 0,9492 0,9070 0,9852 0,9650

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

50 0,9138 0,8716 0,9563 0,9231 0,9916 0,9811

100 0,9108 0,8916 0,9572 0,9437 0,9913 0,9864

300 0,9060 0,9014 0,9551 0,9518 0,9911 0,9901

500 0,9046 0,9025 0,9542 0,9530 0,9909 0,9906

700 0,9041 0,9028 0,9536 0,9530 0,9908 0,9906

900 0,9033 0,9024 0,9531 0,9529 0,9908 0,9906

1100 0,9029 0,9021 0,9530 0,9528 0,9907 0,9906

1300 0,9023 0,9023 0,9528 0,9526 0,9907 0,9906

1500 0,9020 0,9020 0,9527 0,9525 0,9906 0,9906

ж 0,8981 0,8981 0,9505 0,9505 0,9901 0,9901

Как показывает сравнение точных и асимптотических вероятностей распределения статистики (2), разница между ними существенна, поэтому на практике при небольших объемах выборок следует использовать точные вероятности.

Оценка степенного параметра Полученный критерий типа Колмогорова - Смирнова позволяет проверять гипотезы о значениях степенного параметра. Однако существует также задача оценки данного параметра. Обычно в параметрическом случае при известном виде функций распределения ^(О оценка осуществляется ме-

тодом максимизации функции частного правдоподобия [12, 13]. Для

непараметрического случая в работе предлагается оценка, получаемая путем минимизации статистики Колмогорова - Смирнова:

k = arg min Т(к).

Исследование точности оценки степенного параметра проводилось методом статистического моделирования по следующему алгоритму.

1. Моделируем Щ1\ одинаково распределенных случайных величин (x1, ..., XL^ ) с функцией распределения F0(t). В качестве F0(t)

рассмотрим следующие функции распределения: экспоненциальную (с параметром ß = 0,001) и Вейбулла (с параметрами ß = 0,001,

р = 1,5). Величины (Х1,..., Хлг^ ) в дальнейшем будем называть

наработками.

2. Наработки (Х1, ..., Х^ ) случайным образом разбиваем на i\ групп по щ элементов в каждой. Элементы i-й группы обозначим (x1,y,..., Х^' = 1,-Д. Определим максимумы в каждой группе:

01 = max(XV,..,

3. Аналогичным образом моделируем ЩЩ одинаково распреде-

2 2

ленных случайных величин (Х1 , ..., Хщщ, ) с функцией распределения

FQ(tj1k, где к — заранее заданное значение степенного параметра.

22

Наработки (Х1 , ..., Хл^ ) случайным образом разбиваем на Щ групп по Щ элементов в каждой. Элементы i-й группы обозначим . ., Х2щ ), = 1,-°2. Определим максимумы в каждой группе:

02 = max(Х2^,..., Х22').

4. Для некоторого изменяемого значения к, 1 < к < К, по двум полученным выборкам ©1 = ..., 0^ j и ©2 = |02,..., 02^) вычислим

значение статистики Колмогорова - Смирнова Т(к).

5. Методом перебора к определим оценку к, минимизирующую значение статистики Т(к):

к = arg min Т(к).

Для определения статистических свойств оценки к пункты 1-5 повторялись 500 раз. По полученным значениям оценок была построена гистограмма этой выборки, а также вычислены эмпирические

средние Ек и дисперсии 52.

На рис. 2 изображены гистограммы полученных оценок к при щ = 2, щ = 3, д =12 = 100 для экспоненциального распределения с параметром Ь = 0,001 и для распределения Вейбулла с параметрами Р = 0,001, р = 1,5.

ЛГГ N

120

100

50

1,0

1,5

2,0 2,5 а

3,0

100 80 60 40 20 0

1,0 1,5

2,0 2,5 б

3,0

Рис. 2. Гистограммы оценок к при к = 2 (ТУ— число наблюдений):

а — экспоненциальное распределение; б — распределение Вейбулла

По результатам моделирования для экспоненциального распределения математическое ожидание полученной оценки Ек = 2,035 и среднеквадратическое отклонение а = 0,34, а для распределения

Вейбулла Ек = 2,055, а = 0,35, что свидетельствует о возможности применения на практике данных оценок.

Заключение. В работе предложен метод проверки гипотезы о степенной зависимости функций распределения для двух цензури-рованных выборок. Для проверки этой гипотезы разработан критерий типа Колмогорова - Смирнова, основанный на сравнении оценок типа Каплана - Мейера функций распределения по цензурированным данным. Показана сходимость точных распределений статистики к ее асимптотическим распределениям. Методом Монте-Карло проведено моделирование оценки степенного параметра, получаемой минимизацией статистики.

ЛИТЕРАТУРА

[1] Гнеденко Б.В., Беляев Ю.К., Соловьев А.Д. Математические методы в теории надежности. Основные характеристики надежности и их статистический анализ. Москва, Либроком, 2013, 584 с.

[2] Гнеденко Б.В. Вопросы математической теории надежности. Москва, Радио и связь, 1983, 376 с.

[3] Balakrishnan N., Tripathi R.C., Kannan N., Ng Н.К.Т. Some nonparametric precedence-type tests based on progressively censored samples and evaluation of power. J. Stat. Plan. Infer., 2010, vol. 140, no.10, pp. 559-573.

[4] Bagdonavicius V., Kruopis J. Nikulin M.S. Nonparametric tests for censored data. London, Wiley-ISTE, 2013, 233 p.

[5] Balakrishnan N., Cramer E. The Art of Progressive Censoring. Applications to Reliability and Quality. New York, Springer, 2014, 645 p.

[6] Kaplan E.L., Meier P. Nonparametric estimation from incomplete observations. J. of the American Statistical Association, 1958, vol. 53, no. 282, pp. 457-481.

[7] Тимонин В.И., Тянникова Н.Д. Проверка однородности распределений наработок до отказа элементов в двух различных режимах по результатам испытаний резервированных систем, составленных из этих элементов. Физические основы приборостроения [в печати].

[8] Большев Л.Н., Смирнов Н.В. Таблицы математической статистики. Москва, Наука, 1983, 416 с.

[9] Тимонин В.И. О предельном распределении статистики одного непараметрического критерия. Теория вероятностей и ее применение, 1987, т. 32, № 4, с. 790-792.

[10] Тянникова Н.Д. Разработка непараметрических методов анализа цензу-рированных данных при оценке надежности сложных технических систем в различных режимах испытаний. Автореф. дис. ... канд. физ.-мат. наук, Москва, 2014, 16 с.

[11] Тимонин В.И., Тянникова Н.Д. Сравнение прогрессивно цензурированных выборок — численные методы табулирования распределений статистик однородности и исследование оценки параметров связи их распределений методом Монте-Карло. Математическое моделирование и численные методы, 2015, № 3, с. 89-100.

[12] Hajek J., Sidak Z., Sen P.K. Theory of rank tests. London, Academic Press, 1999, 438 p.

[13] Садыхов Г.С., Крапоткин В.Г., Казакова О.И. Расчет и оценка показателей ресурса изделий с использованием модели аддитивного накопления повреждений. Математическое моделирование и численные методы, 2014, № 1, с. 82-98.

Статья поступила в редакцию 26.11.2018

Ссылку на эту статью просим оформлять следующим образом:

Тимонин В.И., Тянникова Н.Д. Методы решения задачи непараметрической проверки гипотез Лемана при испытаниях параллельных систем. Математическое моделирование и численные методы, 2018, № 2, с.

Тянникова Нина Дмитриевна — канд. физ.-мат. наук, доцент кафедры «Высшая математика» МГТУ им. Н.Э. Баумана. Область научных интересов: теория надежности, математическая статистика. e-mail: [email protected]

Тимонин Владимир Иванович — д-р физ.-мат. наук, профессор кафедры «Высшая математика» МГТУ им. Н.Э. Баумана. Автор более 40 работ в теории надежности и математической статистике. e-mail: [email protected]

Methods for solving the problem of non-parametric testing of Lehmann's hypotheses when testing parallel systems

© V.I. Timonin, N.D. Tyannikova Bauman Moscow State Technical University, Moscow, 105005, Russia

The article considers the problem of testing the Lehmann power hypothesis for two censored samples. The Kolmogorov — Smirnov type criterion based on a comparison of the Kaplan — Meier type estimates of the distribution functions for each censored sample is developed to test the power hypothesis. A method for calculating the exact statistics distributions is described on the basis of the model ofparticle random walk over an integer lattice. The probability values are calculated for a wide range of possible sample sizes. The convergence of this statistics distribution to the standard Kolmogorov — Smirnov distribution is proved provided that the hypothesis being tested is valid. The properties of a power parameter estimate obtained by minimizing statistics are investigated by statistical modeling methods.

Keywords: non-parametric statistics, Lehmann's hypothesis, Kolmogorov — Smirnov type criterion, Kaplan — Meier estimate

REFERENCES

[1] Gnedenko B.V., Belyaev Y.K., Soloviev A.D. Matematicheskie metody v teorii nadezhnosti. Osnovnye kharakteristiki nadezhnosti i ikh statisticheskiy analiz [Mathematical methods in reliability theory. The main characteristics of reliability and statistical analysis]. Moscow, Librokom Publ., 2013, 584 p.

[2] Gnedenko B.V. Voprosy matematicheskoy teorii nadezhnosti [Problems of mathematical theory of reliability]. Moscow, Radio i svyaz Publ., 1983, 376 p.

[3] Balakrishnan N., Tripathi R.C., Kannan N. Journal of Statistical Planning and Inference, 2010, no. 140, pp. 559-573.

[4] Bagdanovichus V., Kruopis J. Nikulin M.S. Nonparametric tests for censored data. London, ISTE Ltd Publ., 2011, 233 p.

[5] Balakrishnan N., Cramer E. The Art of Progressive Censoring. Applications to Reliability and Quality. New York, Springer Publ., 2014, 645 p.

[6] Kaplan E.L., Meier P. JAM Statistics Association, 1958, no. 53, pp. 57-481.

[7] Timonin V.I., Tyannikova N.D. Fizicheskie osnovypriborostroeniya - Physical Bases of Instrumentation, [in print].

[8] Bolshev L.N., Smirnov N.V. Tablitsy matematicheskoy statistiki [Mathematical Statistics Tables]. Moscow, Nauka Publ., 1983, 416 p.

[9] Timonin V.I. Teoriya veroyatnostey i ее primenenie - Theory of Probability and its Applications, 1987, vol. 32, no. 4, pp. 790-792.

[10] Tyannikova N.D. Razrabotka neparametricheskih metodov analiza tsenzuriro-vannykh dannykh pri otsenke nadezhnosti slozhnykh tekhnicheskikh sistem v razlichnykh rezhimakh ispytaniy. Avtoreferat diss. cand. fiz.-mat. nauk [Development of nonparametric methods for analyzing censored data in assessing the reliability of complex technical systems in various test modes. Cand. phys. and math. sc. diss. Abstract]. Moscow, 2014, 16 p.

[11] Timonin V.I., Tyannikova N.D. Matematicheskoe modelirovanie i chislennye menody - Mathematical Modeling and Computational Methods, 2015, no. 3 (7), pp. 89-100.

[12] Hajek J., Sidak Z. Theory of rank tests. London, Academic Press Publ., 2004, 438 p.

E.H. TUMOHUH, TnmuKoea

[13] Sadykhov G.S, Krapotkin V.G., Kazakov O.I. Matematicheskoe modelirovanie i chlslennye menody - Mathematical Modeling and Computational Methods,

2014, no. 1, pp. 82-98.

'tyaiuukova N.D., Cand. Sc. (Phys. & Math.) Assoc. Professor, Department of Higher Mathematics, Bauman Moscow State Technical University. Scientific interests: reliability theory, mathematical statistics. e-mail: [email protected]

Timonin V. I., Dr. Sc. (Phys. & Math.), Professor, Department of Higher Mathematics, Bauman Moscow State Technical University, author of over 40 research publications in the field of reliability theory and mathematical statistics. e-mail: [email protected]

i Надоели баннеры? Вы всегда можете отключить рекламу.