Электронная библиотека диссертаций и авторефератов России
dslib.net
Библиотека диссертаций
Навигация
Каталог диссертаций России
Англоязычные диссертации
Диссертации бесплатно
Предстоящие защиты
Рецензии на автореферат
Отчисления авторам
Мой кабинет
Заказы: забрать, оплатить
Мой личный счет
Мой профиль
Мой авторский профиль
Подписки на рассылки



расширенный поиск

Непараметрическое оценивание плотности и функции регрессии для слабо зависимых случайных полей Миллионщиков Николай Владимирович

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

Диссертация - 480 руб., доставка 10 минут, круглосуточно, без выходных и праздников

Автореферат - бесплатно, доставка 10 минут, круглосуточно, без выходных и праздников

Миллионщиков Николай Владимирович. Непараметрическое оценивание плотности и функции регрессии для слабо зависимых случайных полей : диссертация ... кандидата физико-математических наук : 01.01.05.- Москва, 2005.- 82 с.: ил. РГБ ОД, 61 06-1/512

Содержание к диссертации

Введение

1 Асимптотическая нормальность ядерных оценок плотности 19

1.1 Применение техники секционирования к слабо зависимым случайным полям 19

1.2 Скорость нормального приближения ядерных оценок плотности 22

2 Асимптотическая нормальность ядерных оценок функции регрессии 34

2.1 Применение метода Стейна к слабо зависимым случайным полям 34

2.2 Скорость нормального приближения ядерных оценок функции регрессии 39

3 Состоятельность ядерных оценок плотности 54

3.1 Экспоненциальные неравенства 54

3.2 Неравенства для моментов сумм зависимых случайных величин 60

3.3 Сходимость почти наверное ядерных оценок плотности 64

Список литературы 76

Введение к работе

Задача оценивания плотности распределения случайных величин является традиционной задачей статистического анализа данных, представляющей не только теоретический интерес, но и допускающей приложения на практике. Напомним, что имеются два основных типа статистических оценок: параметрические и непараметрические. Непараметрическая процедура обычно определяется как процедура, которая справедлива независимо от типа распределения, которому принадлежит выборка. Привлекательность непараметрических стохастических моделей обусловлена весьма общими условиями, в терминах которых они описываются, например, независимость и одинаковая распределенность изучаемых случайных величин. Именно поэтому основное внимание в диссертации уделяется непараметрическому подходу к оцениванию плотностей. При этом акцент делается на получении результатов, справедливых для зависимых наблюдений. Разумеется, параметрические модели также имеют свои достоинства (часто в большей простоте и законченности статистических выводов). Поэтому целесообразным является взаимодополняющее развитие как параметрических, так и непараметрических методов математической статистики. Имеется множество работ, посвященных изучению свойств иепараметрических оценок плотностей и функции регрессии (см., например, [2], [14], [23], [40], [75]-[72]).

Значительное внимание специалистов в области математической статистики уделяется введенным Розеблаттом ([70]) ядерным оценкам плотности, основные свойства которых были изучены Парзеном ([68]).

Далее все используемые случайные величины и векторы считаются определенными на некотором вероятностном пространстве (П,^7, Р). Как обычно, Е - усреднение по мере Р, var(Y) - дисперсия случайной величины Y. В качестве расстояния между точками і = (ц,... ,) и j = (ji,..., jd) из Zd всегда будет выступать ||г — j\\ = maxi<g\iq — jq\.

Пусть Xi,...,Xn - одинаково распределенные случайные векторы со значениями в Rs, имеющие общую плотность распределения /. Ядерная оценка, или оценка Парзена-Розенблатта, плотности / задается формулой

/»(*);=;^Е

К (^г^) , х Є Ж3, п Є N, (0.0.1)

где К - плотность некоторого вероятностного распределения (ядро),

{hn}n>i - последовательность положительных чисел, стремящаяся к нулю при п —> со.

Различные примеры ядер можно найти в работах [48], [68]. Выбор функции К в формуле (0.0.1) играет большое значение. На практике этот выбор может осуществляться по-разному, исходя из конкретных условий. При работе с ядерными оценками используются разнообразные численные методы, например, быстрое преобразование Фурье (см. [78]).

Существует много подходов к выбору последовательности {hn}n>1 в (0.0.1) (см., например, [2], [14], [53], [77]).

Большое количество работ, посвящено исследованию близости /„ к / в различных метриках, например, метриках пространств Lp, С. Изучению свойств непараметрических, в том числе ядерных, оценок плотности как элементов пространства L\ посвящена книга [14]. Результаты о равномерной в W почти наверное сходимости /п к / были получены, например, в статьях [49], [50].

Ядерные оценки плотности применяются и для процессов с непрерывным временем. Уделяется внимание построению доверительных интервалов и асимптотике дисперсий соответствующих оценок (см., например, [52]).

Наряду с ядерными используются и другие непараметрические оценки плотности, такие как проекционные, гистограммные (см. [14], [23], [26], [81]). Применяются также различные модификации ядерных оценок: адаптивные ядерные оценки, рекуррентные ядерные оценки, преобразованные ядерные оценки и другие (см., например, [14] ). Идеи Парзена и Розенблатта нашли свое применение и при построении семипараметрических оценок, которые в определенных условиях обладают одновременно преимуществами непараметрической и параметрической процедур (см. [32], [58]).

Оценки Парзена-Розенблатта являются классическим примером непараметрических статистических оценок, благодаря их появлению в непараметрической статистике стали активно использоваться многие идеи гармонического анализа и теории аппроксимации функций. При работе с ними широко применяются результаты и методы теории суммирования случайных величин. Их удобство и универсальность привели к появлению и последующему детальному изучению ядерных оценок функций распределения, квантилей (см. [39]), моды (см. [63], [67]), функции регрессии (см. [23])'.

Ядерные оценки функции регрессии были введены в работах Надарая и

Ватсона (см. [22], [82]) и с тех пор использовались и изучались в различных моделях и ситуациях (см., например, [72], [59], [20]).

Имеются две основные постановки задачи оценки функции регрессии. В первой из них рассматривается выборка (Xi,Y\),... ,(Xn,Yn) из наблюдений (s + т)-мерного случайного вектора (X, Y), где X и Y принимают значения в Rs и Rm соответственно, s,m Є N; требуется оценить функцию регрессии У по X, то есть

(п(х),... ,rm{x)) := E(Y\X = х), хЄ W.

Во второй постановке задачи условное распределение У при условии X = х не известно, но точки Xi,... п, в которых производятся наблюдения Y[,...,Yn, задает статистик. Ядерные оценки функции регрессии были предложены Надарая ([22]) в первой постановке задачи, и определяются следующим образом

(f„i(a0,...,rnm(*)):=2=l ; г-» *GRs,nN, (0.0.2)

где, как и в случае ядерных оценок плотности, К - некоторое ядро, {hn}n>1 - последовательность положительных чисел, стремящаяся к нулю при п —» со. Если знаменатель в (0.0.2) равен нулю, то считаем саму оценку также равной нулю.

Оценки Надарая-Ватсона широко используются во многих исследованиях. В некоторых случаях для учета специфики конкретных моделей их нужным образом модифицируют (см. [83]).

Ряд работ посвящен вопросам оптимальной скорости сходимости оценок функции регрессии в различных метриках, например, метриках пространств Lpu С (см. [80], [19], [18]).

Для выборок, состоящих из независимых случайных величин, основные статистические свойства, например, состоятельность и асимптотическая нормальность ядерных оценок как плотности так и функции регрессии хорошо известны (см. [23], [14], [50]).

Имеется большое количество важных стохастических моделей, в которых нужно учитывать зависимость между случайными элементами. Например, модели, в которых рассматриваются марковские процессы и поля (см. [15]), мартингалы (см. [30]), процессы и поля с премешиванием (см. [17], [44]), гиббсовские поля (см. [25]).

Ряд работ посвящен оцениванию переходных вероятностей для цепей Маркова (см., например, [12], [71]). Во многих статьях изучаются свойства ядерных оценок плотности (см. [40]) и функции регрессии (см. [72], [59]) для процессов и полей, удовлетворяющих тем или иным условиям перемешивания.

Одним из наиболее распространенных способов описания стохастической зависимости является задание ограничений на ковариации некоторых функций от конечных наборов случайных величин. В частности, если рассмотреть покоординатно неубывающие функции и потребовать, чтобы соответствующие ковариации были неотрицательны, то мы придем к определению ассоциированных случайных величин, введенному в работе [46].

Определение 0.1([46]). Семейство X = {Xt,t Є Т} случайных величин, где Т-некоторое параметрическое множество, называется ассоциированным, если для любых конечных множеств индексов I,JcT и ограниченных покоординатно неубывающих функций / : R'7' —» R и д : RlJl —> R справедливо неравенство

cov(f(Xs,s Є /),g(Xt,t Є J)) > 0.

Здесь и далее |/| обозначает число элементов конечного множества /.

В [46] показано, что произвольное семейство взаимно независимых случайных величин ассоциировано. Питт доказал ([69]), что гауссовское случайное поле ассоциировано тогда и только тогда, когда его ковариационная функция всюду неотрицательна. Примеры ассоциированных случайных полей возникают в статистической физике (см. [51], [25]).

В статье [47] Фортуин, Кастелейн и Жинибр использовали условия зависимости случайных полей, задаваемые в виде неравенств для определенных мер на решетках (FKG-неравенства). Эти неравенства влекут ассоциированность и применяются в статистической физике, квантовой механике, теории гиббсовских полей, теории перколяции (см. [25], [47], [65]).

Имеется следующее обобщение определения 0.1 на многомерный случай.

Определение 0.2([38], [46]). Семейство X = {Xt,t Є Т} случайных векторов со значениями в R5, называется слабо ассоциированным (или положительно ассоциированным), если для любых конечных непересекающихся множеств индексов I,J С Ти ограниченных покоординатно неубывающих функций / : R5'7' —» R и д : RS'J' — R

справедливо неравенство

cov{f(Xs,s Є I),g{Yut Є J)) > 0.

Йаг-Дев и Прошан ввели понятие отрицательно ассоциированных случайных векторов.

Определение 0.3((57]). Семейство X = {Xt,t Є Т} случайных векторов со значениями в Rs, называется отрицательно ассоциированным, если для любых конечных непересекающихся множеств индексов I,JgT и ограниченных покоординатно неубывающих функций / : RS'7I —» R и д : RSIJI —» R справедливо неравенство

cov(f(Xs,s Є I),g{Yt,t Є J)) < 0.

Со свойствами и примерами отрицательно ассоциированных случайных величин можно ознакомится в работе [57]. Приведем некоторые из них.

Независимые случайные величины удовлетворяют определению 0.3. Отрицательно коррелированные гауссовские случайные величины являются отрицательно ассоциированными. Примеры' отрицательно ассоциированных случайных величин возникают в математической статистике. Многие известные распределения обладают свойством отрицательной ассоциированности, например, мультиномиальное, многомерное гипергеометрическое и другие.

Случайные процессы и поля, удовлетворяющие условиям зависимости 0.1-0.3, являются предметом активного изучения многих исследователей. В течение нескольких последних десятилетий были доказаны различные варианты центральной предельной теоремы, принципы инвариантности, законы повторного логарифма, законы больших чисел, экспоненциальные неравенства и другие результаты для сумм положительно и отрицательно ассоциированных случайных величин. При этом предполагалось, что зависимость между случайными величинами, выражаемая в терминах ковариаций, убывает с определенной скоростью, когда расстояния между индексами этих величин растут (см., например, [64], [33], [34], [5], [66], [28]).

Много работ посвящено изучению свойств непараметрических оценок плотности, функции регрессии и их производных, построенных по выборкам, образующим положительно или отрицательно ассоциированные процессы и поля. Для этих оценок доказана асимптотическая нормальность слабая и сильная состоятельность (см. [74], [61], [60], [54]), получены

оценки скорости нормального приближения (см. [10], [37], [20]). При доказательстве этих результатов, активно используются методы теории суммирования случайных величин.

Для формулировки некоторых определений и результатов, которые нам потребуются в дальнейшем, введем ряд понятий и обозначений.

Рассмотрим метрические пространства S и М соответственно с метриками бит. Функция / : S —> М называется липшицевой, если

ЫрШ= ІМШШ<00.

x,yeS,x^y 0\Х,У) t

Класс ограниченных липшицевых функций / : Sm —> R, где m Є N, обозначим BL(Sm); в пространстве Sm используется метрика 5m(x,y) = ХГ=15(хк, Ук), Для х = (яь..., хт) Є Sm и у = (г/ь..., ут) Є 5m, а в R берется евклидово расстояние.

Для положительно и отрицательно ассоциированных случайных величин справедливы некоторые неравенства, которые связывают ковариацию функций от наборов случайных величин с суммой ковариаций самих величин (см. [65], [33], [35]). Выполнения этих неравенств часто бывает достаточно для доказательства предельных теорем, неравенств и других результатов. В связи с этим появились новые определения, обобщающие свойства положительной и отрицательной ассоциированности.

Определение 0.4([35]). Семейство X — {Xut Є Т} случайных величин, имеющих конечный второй момент, называется квазиассоциир-ованным, если для любых непересекающихся конечных множеств I, J С Т и любых / Є L(Rl7l),д Є BL{m)J\)

\cov{f{Xuie I),g(Xj,j є J))\ < Lip(f)Lip(g) \cov(XuXj)\.

iEljZJ

(0.0.3) В [35] использовались липшицевы функции, для которых учитывалась возможность разной скорости роста по различным аргументам, что несущественно видоизменяло (0.0.3). В статье [9] показано, что соотношение (0.0.3) выполнено и для положительно и для отрицательно ассоциированных случайных величин с конечными вторыми моментами. Любое семейство гауссовских случайных величин удовлетворяет неравенству (0.0.3) (см. [43]). Определение 0.4 несложно обобщить для случайных векторов со значениями в Ш8.

Определение 0.5([6]). Семейство X = {Xt, t ЄТ} случайных векторов со значениями в Rs, имеющих конечный второй момент, называется

квазиассоциированным, если для любых непересекающихся конечных множеств I, J С Т и любых / Є BL(RSW), д Є BL(RSW)

\cov{f{Xi,ieI),g{Xhj Є J))| < Lip(f)Lip(g) ± \cov(Xik,Xjv)\.

il,j&J k,v=l

(0.0.4)

А.П.Шашкин ([27]) показал, что любая гауссовская система случайных векторов, принимающих значения в Rs, обладает свойством квазиассоциированности.

Аналог неравенства (0.0.3) для гладких функций fug был впервые доказан в работе [33], примеры использования ковариационных неравенств можно найти в статье [43].

В данной диссертации наряду с понятием квазиассоциированности используются более общие подходы к описанию стохастической зависимости, которые применялись в работах П.Дукана, С.Луиши, А.В.Булинского, Ш.Сюке, Г.Ланга, А.П.Шашкина, Н.В.Миллионщикова (см. [43], [45], [35], [И], [37], [7], [20], [21]).

Пусть Т - метрическое пространство с метрикой р, и пусть в = {0Г}Г>1 - такая числовая последовательность, что вг —» 0 при г —> со.

Определение 0.6([35]). Случайное поле X — {Xut Є Т}, принимающее при каждом t значения в метрическом пространстве S, называется (BL,9) - зависимым, если для любого s Є Т и произвольного непустого конечного множества J С Т такого, что і ф J,

|«w(/(X,),0«,J Є J))\ < erLip(f)Lip(g) (0.0.5)

при всех / Є BL{S),g Є BL{S^), где г = inf{p(z, j) : j Є J}.

Для описания зависимости поля X широко используется коэффициент Кокса-Гримметта ([41])

0r* = sup Yl \cov(Xj,Xq)\, reN. (0.0.6)

jeZ" qZZ*:\\q-j\\>r

(\\x\\ = maxi\хі\ для x = {х\,...,хд) Є Rd). Легко видеть, что квазиассоциированное поле X = {Xj,j Є Zd} является (ВЬ,в) -зависимым, если в{ < со и в* —> 0 при г —> со, причем можно положить вг = в*.

В диссертации [29] построены новые примеры случайных полей скользящего среднего, гауссовских и марковских, являющихся (ВЬ,в)-зависимыми.

В течение последних десяти лет для описания зависимости между элементами случайных полей активно используется следующая схема, которая охватывает определение 0.6. Пусть L - некоторый класс действительнозначных функций конечного числа аргументов, а неотрицательный функционал и : С х С х 2Т х 2Т —> Е. Условие зависимости поля {Xt, t Є Т} задается в виде неравенства

\cov{f(Xi,i є I),g(XjJ Є J))\ < u{f,gJ,J), (0.0.7)

справедливого для всех f,g Є и множеств 7, J С Т определенного вида (например, состоящих из одного элемента). Требуется подобрать подходящим образом класс С, вид множеств I, J и функцию и, чтобы неравенство (0.0.7) обеспечивало решение поставленной задачи, например, получения оценки скорости сходимости в центральной предельной теореме или законе больших чисел. В разных задачах и моделях этот выбор может осуществляться по-разному. Если в качестве С взять класс BL(S) и потребовать выполнения (0.0.7) с u(f,g,I,J) = 9rLip(f)Lip(g) для всех непустых 7, J, таких что |/| = 1, | J\ < со и I П J = 0, то получим определение (BL, ^-зависимости. Естественно считать, что функция и учитывает расстояние между множествами I и J, а также мощности этих множеств. "Увеличение расстояния" приводит к уменьшению зависимости, а "рост" хотя бы одного из множеств, вообще говоря, может увеличить зависимость между f(Xi,i Є I) и g(Xj,j Є J).

В работах [45], [37] использовались неравенства (0.0.7), в которых u(f,g,I,J) = ф(1,д, \I\,\J\Wr, где ф : С2 х N2 -> R, а г = dist(I,J) = min{||i-j|| :ielj Є J}.

Чаще всего применяются два типа функционалов ф:

Фі(ІдАіШ\)=тт(\ІІ\А)ЬІР(/)Ьір(д), (0.0.8)

ф2(/,д, |/|, |J|) = ЩЬірі/Ш^ + \J\Lip(g)\\f\L- (0.0.9)

В качестве класса С "пробных функций" кроме ограниченных липшицевых функций берутся и другие классы функций. Например, в работе [1] для получения оценок абсолютных моментов частичных сумм мультииндексированных зависимых случайных величин использовались функции степенного типа. Рассматриваются также ковариационные неравенства для линейных функций, функций с ограниченной вариацией; в контексте коэффициентов перемешивания возникают индикаторные функции (см. [44]).

Заметим, что в работах [37], [45] для описания удаленности множеств I и J использовалось расстояние dist(I, J). Это не давало простой возможности учитывать разную скорость убывания зависимости поля {Xt,t Є 7Ld,} по различным направлениям. В данной работе будет использоваться модификация определения 0.6.

Пусть заданы некоторый класс С функций конечного числа действительных аргументов со значениями в R, неотрицательный функционал ф :XxxNxN—> R, и массив в = {0r}reNd, состоящий из неотрицательных чисел, причем 9Г —> 0, когда все компоненты гг-, г = 1,...,(1, мультииндекса г стремятся к бесконечности. Определим для I С Zd, г ENd множество Н(1,г) С Rd следующим образом:

H{I,r) := Un^ - Гк>ік + г*)-

Определение 0.7 Случайное поле X = {Xj,j Є U}, U С Zd, d Є N со значениями в Rs, s Є N, называется (, ф, 0)-зависимым, если для любых непустых непересекающихся конечных множеств I,J a U таких, что Н(1,г) П J = 0, и всех функций / : Rsl;l -* R и д : RSIJI -* R, принадлежащих классу С, справедливо неравенство

\cov(f(Xi,ie I),g(Xj,j Є J))\ < ф(/,д, \I\,\J\)9r.

В данной работе, как правило, С = BL, а ф = ф\, то есть рассматриваются (5Ь,'0ь^)-зависимые случайные поля, и под "слабой зависимостью" будет пониматься (ВЬ,фі, ^-зависимость.

В отличие от семейств независимых случайных величин, для зависимых случайных полей большое значение имеет расположение индексов элементов поля в пространстве, так как оно существенным образом учитывается большинством условий зависимости. В многомерном случае (т.е. при изучении случайных полей, заданных на подмножествах Ъл, где d > 1) уже не существует понятий "прошлое" и "будущее", традиционно используемых в анализе случайных процессов. Кроме того, суммирование, взятие максимума и другие действия с элементами поля могут производиться по множествам более сложной структуры, чем для случайных процессов. В результате методы, применяемые для случайных процессов, удовлетворяющих тем же условиям зависимости, и тем более независимых, требуют изменений. Эти изменения часто приводят к появлению различных ограничений на множества индексов случайных

элементов. Например рассматривают только кубы или параллелепипеды, используют последовательности конечных множеств, стремящихся к бесконечности в смысле Ван Хова (см. [8], [36]).

Для выборок, образованных случайными полями {Х{,і Є Un}, {{Xi,Yj),i Є Un} аналогами (0.0.1) и (0.0.2) служат формулы

іє[/„

(rnl{x)... ,rnm(x)) := — /' ч , xer.neN. (0.0.11)

Как уже было отмечено, на асимптотические свойства оценок (0.0.10) и (0.0.11) существенным образом влияют характер зависимости изучаемого семейства величин Xj, поведение ядра К, скорость стремления к нулю последовательности {hn}n>l, вид множеств Un. В случае /п, являющихся "регулярно растущими целочисленными прямоугольниками", состоятельность и скорость нормального приближения оценок (0.0.10), (0.0.11) получены в работах [10], [40], [37], [59]. В этих работах использовались условия сильного перемешивания и (BL, #)-зависимости. Отметим, что в статье [20] при доказательстве асимптотической нормальности оценок (0.0.11) рассматривались множества Un произвольной конфигурации.

Асимптотическая нормальность оценок (0.0.1) для положительно и отрицательно ассоциированных выборок была доказана в работе [74]. А именно, для положительно или отрицательно ассоциированного строго стационарного случайного процесса X = {Хі,і Є N} при iER таких, что f(x) > 0 и / непрерывна в х, установлено соотношение

Vnh~n{?n{x) - f(x))/cr(x) N(0,1) при п -> со,

где cr2(x) = f(x)fRK2(x)dx, " —> " обозначает сходимость по распределению, iV(0,1) - стандартный нормальный закон. При этом предполагается выполнение следующих условий.

Пусть EXf < со и J2Zi \cov(^i^i+i)\ < со. Производная /' существует и ограничена на R. Для любых і > 0 существует совместная плотность /і,і величин ХъХ{+1 и |/i,i(w,u) — f(u)f(v)\ < С при всех u,v Є Ми некоторой константе С > 0. Ядро К является ограниченной плотностью

вероятностного распределения, имеет ограниченную производную и fR\u\K(u)du < со. Кроме того, \imn^oo nh^ = 0, lim^oon/in = со. Существуют последовательности натуральных чисел 0?п}пЄм> {Яп}пеи такие, что Ишп^ооРпкп/п = 1, lim,,.^ pnhn = 0, Ііт^^ р2п/nhn - О, limrj^oo/i-^SgJc^^b^+i)! = 0, где кп = [п/(рп + qn)] (здесь [] обозначает целую часть числа).

В первой главе данной диссертационной работы получен аналог этого результата для слабо зависимых случайных полей, при этом множества Un, фигурирующие в формуле (0.0.10), представляют собой "регулярно растущие целочисленные прямоугольники". Доказательство этого результата основано на применении следующей теоремы, установленной в работе [35].

Теорема 0.1. Пусть X = {Xt,t Є Zd} есть (ВL,9)-зависимое случайное поле такое, что EXj = 0 и EXj < со для каждого j Є Zd. Тогда для любого конечного U С Zd, для которого В2 := YLjeu ^-Щ > ^> всех х Є Ш и произвольных положительных є, 7 верна оценка

в~1 Xj > х ) - ?{Z > х)

(0.0.12)

+ (2 + f) (f + (4 + є + 2e~l)Ct + 3|І7|Б-2^) ,

здесь Сє = В~2 X^jef/ ^Xjl{\Xj\>eB} * функция Линдеберга, Z ~ JV(0,1).

Заметим, что для независимых слагаемых 9\ = 0, и тем самым оценка (0.0.12) упрощается. В частности, получается справедливость классической теоремы Линдеберга. Для зависимых слагаемых оценка скорости нормального приближения, как правило, хуже, чем для сумм независимых случайных величин. Достаточно напомнить эталонные оценки Берри-Эссеена (см., например, [24] стр. 154), а также неклассические варианты этих оценок (см., например, [16], [31], [34]).

Теорема 0.1 получена в [35] с помощью метода Стейна. Напомним, что идея метода Стейна (см. [79]) состоит в следующем. Рассматривается дифференциальное уравнение

f(x)-xf(x)=g(x)-Eg(Z), xeR, (0.0.13)

где g - заданная ограниченная непрерывная функция, a Z - стандартная нормальная величина. Данное уравнение имеет единственное ограниченное

решение, которое определяется формулой

f(x) = е*2/2 J e-t2'2{g{t) - Eg{Z))dt, xeR.

Заменив х в (0.0.13) на случайную величину W и взяв математическое ожидание от каждой части этого уравнения, получим

Ef(W) - EWf(W) = Eg(W) - Eg(Z). (0.0.14)

Таким образом, чтобы оценить близость распределения W к нормальному, выражаемому величиной Eg(W)—Eg(Z), достаточно оценить левую часть уравнения (0.0.14). Заметим, что такая оценка, вообще говоря, не является тривиальной.

Для выборок, состоящих из независимых случайных векторов (Xi,Yi), г = 1,...,п, хорошо известен следующий результат (см. [23] стр. 124). Пусть в (0.0.11) 5 = 1, К - ограниченная симметричная функция такая, что lim^i^oo \х\К(х) = 0, функция х2К(х) принадлежит пространству Li(R). Функция V(х) = fRy2f(x,y)dy непрерывна. Кроме того, плотность / распределения величины X и функция ср(х) = fRyf(x,y)dy имеют ограниченные производные первого порядка, |У| < С п.н., lim^oo nhn = со, lim^oo nhn — 0- Тогда

\Znhn(rn(x) — r(x))/a(x) > N(Q, 1) при n —> oo,

o2{x) := var(Y\X = x)/f(x) f K2{u)du.

Возникает вопрос, при каких условиях ядерные оценки функции

регрессии обладают свойством асимптотической нормальности (точнее

говоря, величина л/пкпп(х) г(х))/сг(х) распределена асимптотически

нормально), если поле {(Xi,Yi),i Є Un} является (BL,ipi,^-зависимым,

a Un С Ъл. Можно ли в качестве множеств Un брать не только

"целочисленные прямоугольники", а множества более общей структуры?

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

убывания величин вг, достаточные для получения оценки скорости

нормального приближения совместного распределения ядерных оценок

функции регрессии (0.0.11), вычисленных в конечном числе точек X Є Rs.

При этом никаких ограничений на форму множеств Un не налагается. Для

доказательства этого результата используется предложенный в работе [37]

модифицированный метод Стейна, позволяющий оценивать приближение к нормальному закону случайных векторов.

Сильная состоятельность ядерных оценок плотности для положительно и отрицательно ассоциированных выборок изучалась в недавних работах Масри и Оливейра, см. [61], [54].

Пусть {-^і}ієК - строго стационарная последовательность положительно (отрицательно) ассоциированных случайных величин. В статье [61] изучаются ядерные оценки плотности распределения случайных векторов Xj = (Xj+i,... ,Xj+m), т Є N, а также оценки ее производных. Предполагается, что плотность / случайных величин X,- ограничена. Кроме того, для всех j > 0 существует совместная плотность fj векторов Хо и Xj (для j < т плотность распределения вектора (Xq,. .. ,Xj)), причем fj ограничены равномерно по j. Также предполагается, что ядро К имеет ограниченные частные производные; В этой работе, в частности, доказано, что если в* = J2\j\>rcov(Xo>Xj) < Саг для некоторых С > О, а Є (0,1), то

sup \fn(x) - Efn(x)\ = 0{(nhy/2) п.н. при п -> со

для любого компакта Q С Rm и всех 7 > 0.

При доказательстве этого результата использовались моментные неравенства для сумм положительно и отрицательно ассоциированных случайных величин из статьи [76].

В [54] для получения оценки скорости сходимости почти наверное к нулю величин \fn{x) — Е/П(х)| использовались экспоненциальные неравенства для сумм ассоциированных случайных величин. При этом требовалась экспоненциальная скорость убывания коэффициента Кокса-Гримметта в* = Х)ь|>гсог,(^о)-^і)) когда г, стремится к бесконечности. Полученная

скорость сходимости имеет вид 0((lnn/n(1-7)/in) ) для некоторого 7 ^ (0,1).

Состоятельность ядерных оценок плотности изучалась и для случайных полей. Так в работе [40] рассматривались оценки Парзена-Розенблатта плотности, построенные по наблюдениям, образующим векторное случайное поле, удовлетворяющее условию сильного перемешивания. При этом в качестве множеств Un для построения оценок (0.0.10) брались "регулярно растущие целочисленные прямоугольники".

В третьей главе диссертации с помощью новых экспоненциальных и моментных неравенств для сумм элементов (BL,ipi, ^-зависимых

случайных полей доказывается сильная состоятельность ядерных оценок плотности. При этом никаких ограничений на форму множеств Un не налагается. Используемое здесь определение 0.7, позволяет учитывать различный характер зависимости поля по разным направлениям. Это важно для некоторых моделей, применяемых, например в геостатистике, где для описания зависимости гауссовского случайного поля {Z,-, і Є Zd} используют вариограмму (см., например, [42])

7(r) = -var(Z0 - Zr), г Є Zd.

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

Из определения 0.4 видно, что, как и в случае гауссовского поля, зависимость стационарного квазиассоциированного случайного поля определяется его ковариационной функцией. Поэтому естественно возникает возможность учета анизотропии в близком к квазиассоциированности понятии (,^, ^-зависимости.

Сделаем несколько замечаний. Далее буква С будет обозначать, вообще говоря, различные для разных формул положительные множители, не зависящие от п, а в некоторых случаях и от других параметров, что будет оговариваться отдельно. Поскольку все оценки в главах 1 и 2 доказываются с точностью до множителей, которые могут зависеть от размерности рассматриваемых конечномерных пространств, то выбор нормы в рассматриваемом пространстве не играет роли.

Структура работы

Работа, объемом 82 страницы, состоит из введения, трех глав и списка литературы, насчитывающего 83 наименования.

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

Вторая глава также состоит из двух параграфов и посвящена изучению асимптотики совместных распределений ядерных оценок функции регрессии для слабо зависимых случайных полей. В 2.1 излагается метод

Стейна и доказываются результаты, которые играют ключевую роль в 2.2 при получении оценки скорости нормального приближения ядерных оценок функции регрессии (теоремы 2.2.2 и 2.2.3).

Третья глава содержит три параграфа. В первом параграфе устанавливается новое экспоненциальное неравенство (лемма 3.1.1) для сумм элементов (BL, -0і,0)-зависимого случайного поля, во втором доказано неравенство для моментов сумм элементов (ВЬ+фі, #)-зависимого случайного поля (лемма 3.2.1), а в 3.3.3 на основе результатов 3.3.1 и 3.3.2 получены оценки скорости сходимости почти наверное ядерных оценок плотности (см. теоремы 3.3.2 и 3.3.3).

Основные результаты диссертации опубликованы в работах автора [10], [20], [62], [21]. В работе [10] А.В.Булинскому принадлежит постановка задачи и разработка подхода, основанного на использовании ковариационных неравенств, Н.В.Миллионщикову принадлежит асимптотический анализ вспомогательных функций и оптимизация используемых параметров. Все остальные результаты получены Н.В.Миллионщиковым самостоятельно. Результаты диссертации докладывались на международной конференции "Колмогоров и современная математика" (Москва, 2003), в 2005 году на Большом семинаре кафедры теории вероятностей механико-математического факультета МГУ (руков. чл.-корр. РАН, проф. А.Н.Ширяев), а также в 2004 и 2005 годах на семинаре "Асимптотический анализ случайных процессов и полей" (мех-мат МГУ, руков. проф. А.В.Булинский и проф. В.И.Питербарг). Автор благодарен своему научному руководителю профессору А.В.Булинскому за постановку задач, постоянное внимание к работе и ценные советы.

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

Различные примеры ядер можно найти в работах [48], [68]. Выбор функции К в формуле (0.0.1) играет большое значение. На практике этот выбор может осуществляться по-разному, исходя из конкретных условий. При работе с ядерными оценками используются разнообразные численные методы, например, быстрое преобразование Фурье (см. [78]). Существует много подходов к выбору последовательности {hn}n 1 в (0.0.1) (см., например, [2], [14], [53], [77]). Большое количество работ, посвящено исследованию близости /„ к / в различных метриках, например, метриках пространств Lp, С. Изучению свойств непараметрических, в том числе ядерных, оценок плотности как элементов пространства L\ посвящена книга [14]. Результаты о равномерной в W почти наверное сходимости /п к / были получены, например, в статьях [49], [50]. Ядерные оценки плотности применяются и для процессов с непрерывным временем. Уделяется внимание построению доверительных интервалов и асимптотике дисперсий соответствующих оценок (см., например, [52]). Наряду с ядерными используются и другие непараметрические оценки плотности, такие как проекционные, гистограммные (см. [14], [23], [26], [81]). Применяются также различные модификации ядерных оценок: адаптивные ядерные оценки, рекуррентные ядерные оценки, преобразованные ядерные оценки и другие (см., например, [14] ). Идеи Парзена и Розенблатта нашли свое применение и при построении семипараметрических оценок, которые в определенных условиях обладают одновременно преимуществами непараметрической и параметрической процедур (см. [32], [58]). Оценки Парзена-Розенблатта являются классическим примером непараметрических статистических оценок, благодаря их появлению в непараметрической статистике стали активно использоваться многие идеи гармонического анализа и теории аппроксимации функций. При работе с ними широко применяются результаты и методы теории суммирования случайных величин. Их удобство и универсальность привели к появлению и последующему детальному изучению ядерных оценок функций распределения, квантилей (см. [39]), моды (см. [63], [67]), функции регрессии (см. [23]) . Ядерные оценки функции регрессии были введены в работах Надарая и Ватсона (см. [22], [82]) и с тех пор использовались и изучались в различных моделях и ситуациях (см., например, [72], [59], [20]). Имеются две основные постановки задачи оценки функции регрессии. В первой из них рассматривается выборка (Xi,Y\),... ,(Xn,Yn) из наблюдений (s + т)-мерного случайного вектора (X, Y), где X и Y принимают значения в Rs и Rm соответственно, s,m Є N; требуется оценить функцию регрессии У по X, то есть Во второй постановке задачи условное распределение У при условии X = х не известно, но точки Xi,... ,Хп, в которых производятся наблюдения Y[,...,Yn, задает статистик. Ядерные оценки функции регрессии были предложены Надарая ([22]) в первой постановке задачи, и определяются следующим образом где, как и в случае ядерных оценок плотности, К - некоторое ядро, {hn}n 1 - последовательность положительных чисел, стремящаяся к нулю при п —» со. Если знаменатель в (0.0.2) равен нулю, то считаем саму оценку также равной нулю. Оценки Надарая-Ватсона широко используются во многих исследованиях. В некоторых случаях для учета специфики конкретных моделей их нужным образом модифицируют (см. [83]). Ряд работ посвящен вопросам оптимальной скорости сходимости оценок функции регрессии в различных метриках, например, метриках пространств Lpu С (см. [80], [19], [18]). Для выборок, состоящих из независимых случайных величин, основные статистические свойства, например, состоятельность и асимптотическая нормальность ядерных оценок как плотности так и функции регрессии хорошо известны (см. [23], [14], [50]). Имеется большое количество важных стохастических моделей, в которых нужно учитывать зависимость между случайными элементами. Например, модели, в которых рассматриваются марковские процессы и поля (см. [15]), мартингалы (см. [30]), процессы и поля с премешиванием (см. [17], [44]), гиббсовские поля (см. [25]). Ряд работ посвящен оцениванию переходных вероятностей для цепей Маркова (см., например, [12], [71]). Во многих статьях изучаются свойства ядерных оценок плотности (см. [40]) и функции регрессии (см. [72], [59]) для процессов и полей, удовлетворяющих тем или иным условиям перемешивания. Одним из наиболее распространенных способов описания стохастической зависимости является задание ограничений на ковариации некоторых функций от конечных наборов случайных величин. В частности, если рассмотреть покоординатно неубывающие функции и потребовать, чтобы соответствующие ковариации были неотрицательны, то мы придем к определению ассоциированных случайных величин, введенному в работе [46]. Определение 0.1([46]). Семейство X = {Xt,t Є Т} случайных величин, где Т-некоторое параметрическое множество, называется ассоциированным, если для любых конечных множеств индексов I,JcT и ограниченных покоординатно неубывающих функций / : R 7 —» R и д : RlJl — R справедливо неравенство Здесь и далее / обозначает число элементов конечного множества /. В [46] показано, что произвольное семейство взаимно независимых случайных величин ассоциировано. Питт доказал ([69]), что гауссовское случайное поле ассоциировано тогда и только тогда, когда его ковариационная функция всюду неотрицательна. Примеры ассоциированных случайных полей возникают в статистической физике (см. [51], [25]). В статье [47] Фортуин, Кастелейн и Жинибр использовали условия зависимости случайных полей, задаваемые в виде неравенств для определенных мер на решетках (FKG-неравенства). Эти неравенства влекут ассоциированность и применяются в статистической физике, квантовой механике, теории гиббсовских полей, теории перколяции (см. [25], [47], [65]). Имеется следующее обобщение определения 0.1 на многомерный случай. Определение 0.2([38], [46]). Семейство X = {Xt,t Є Т} случайных векторов со значениями в R5, называется слабо ассоциированным (или положительно ассоциированным), если для любых конечных непересекающихся множеств индексов I,J С Ти ограниченных покоординатно неубывающих функций / : R5 7 —» R и д : RS J —справедливо неравенство

Применение метода Стейна к слабо зависимым случайным полям

Отметим, что в следствии 3.2.1 учитывается анизотропное изменение зависимости. Если этого не делать, то согласно результатам работы [13], вместо условия (3.2.27) надо было бы использовать более сильное предположение вг Сгпг , где mr = minj=iv..idri. Лемма 3.2.1 позволяет также рассматривать поля, характер зависимости которых неодинаков в различных частях индексного множества. В данном параграфе рассматриваются серии случайных полей Х = {Х ,і Є Un} Un С Zd, n Є N, со значениями в Rs, удовлетворяющих условию (BL, i, -зависимости, которое обобщает понятия положительной (отрицательной) ассоциированности. Все случайные векторы Xj,j Є Un,n Є N, имеют общую плотность распределения /. Для каждого п Є N ядерные оценки /„(ж) плотности / строятся по элементам полей в соответствии с формулой (0.0.10). При определенных условиях, налагаемых на распределения полей Х и ядро К, удается получить более точную оценку скорости сходимости fn(x) к /(ж), чем в работах [54] и [61], в которых выборки состоят из элементов ассоциированных последовательностей. Для формулировки основных результатов нам потребуется ряд условий, которые удобно разделить на группы. В первой из них представлены свойства распределений элементов полей (CI) {Un}n&N - последовательность конечных подмножеств Zd. При всех п Є N случайные поля Х = {Xj,j Є Un} принимают значения в Rs,s Є N и являются (ВЬ,-фі,в )-зависииьіми, где 9 = {#r )veNd- Для всех п Є N и любых j Є Un случайные величины Xj имеют общую плотность распределения /, являющуюся ограниченной липшицевой функцией, причем suPa.gR» f(x) В\. Для любых п Є N, hj Є Un, і j, случайные векторы pQ , Xj ) имеют плотность распределения Д- и существует В2 0 такое, что sup g » Д- (х) Bi для всех п Є N и любых несовпадающих i,j ЄІ/п. Во второй группе приводятся ограничения на ядро К. (С2) Ядро К(х),х Є Rs, является липшицевой функцией и существуют положительные константы Яі,Яг такие, что supieR, К(х) #i и / \\x\\K(x)dx = Я2. В третьей группе собраны предположения о поведении множеств Un и последовательности { n}n6N а также скорости убывания зависимости между элементами полей Х п\ выражаемой величинами R(Un,pn,a) (см. (3.1.12)), когда расстояние между их индексами растет. Последовательность {Pn}nN» которой здесь идет речь, задает размеры блоков Ak(U,pn, а), участвующих в используемом методе секционирования (см. формулы (3.1.9) и (3.1.10)). А последовательность { 7п}пЄн возникает при получении асимптотики дисперсий случайных величин YljeUn К{(х — Xnj)/hn). Формулировка условий ((73) представляется довольно сложной. Упрощенные аналоги этих условий будут приведены ниже. Как обычно, (а) - произведение Yii=iа» компонент вектора а Є Rd. (СЗ) Существуют последовательности {pn}neN и { n}neN элементов Nd такие, что при разбиении (3.1.10) множества Un на блоки Ak(Un, рп) размера рп, для некоторого є 0, любых к Є D(Un,pn) (см. (3.1.11)) и всех достаточно больших п Є N и, кроме того, для некоторого С 1 и всех а Є {0,1} справедливо неравенство Теорема 3.3.1. Пусть выполняются условия (С1),(С2) и (СЗ). Тогда при всехп Є N, для которых выполняются соотношения (3.3.29), (3.3.31), и всех х Є Ж8 таких, что f(x) О, (3.3.32) где a"j, Рп є, и С фигурируют в условии (СЗ). Кроме того, Нам потребуются некоторые дополнительные обозначения. Определим случайные величины Knj := К((х - Xf)/K) - ЕЩх - X}n))/hn),x Є Rs,j Є Un,n Є N. Для каждого V CUn обозначим jev Докажем несколько вспомогательных утверждений. Лемма 3.3.1. Пусть выполнены условия (С1),(С2) uVnC Un, qn Є Nd, пвК Тогда где Pn задается формулой (3.3.30). Доказательство. Введем Q{i,q) := Vn\H({i}tq), где H({i},q) фигурирует в определении (BL, ір, 0)-зависимости, і Є Zd, g Є Nd. Представим var(Kn(Vn)) в виде.

Скорость нормального приближения ядерных оценок функции регрессии

Ряд работ посвящен оцениванию переходных вероятностей для цепей Маркова (см., например, [12], [71]). Во многих статьях изучаются свойства ядерных оценок плотности (см. [40]) и функции регрессии (см. [72], [59]) для процессов и полей, удовлетворяющих тем или иным условиям перемешивания.

Одним из наиболее распространенных способов описания стохастической зависимости является задание ограничений на ковариации некоторых функций от конечных наборов случайных величин. В частности, если рассмотреть покоординатно неубывающие функции и потребовать, чтобы соответствующие ковариации были неотрицательны, то мы придем к определению ассоциированных случайных величин, введенному в работе [46].

Определение 0.1([46]). Семейство X = {Xt,t Є Т} случайных величин, где Т-некоторое параметрическое множество, называется ассоциированным, если для любых конечных множеств индексов I,JcT и ограниченных покоординатно неубывающих функций / : R 7 —» R и д : RlJl — R справедливо неравенство

Здесь и далее / обозначает число элементов конечного множества /. В [46] показано, что произвольное семейство взаимно независимых случайных величин ассоциировано. Питт доказал ([69]), что гауссовское случайное поле ассоциировано тогда и только тогда, когда его ковариационная функция всюду неотрицательна. Примеры ассоциированных случайных полей возникают в статистической физике (см. [51], [25]). В статье [47] Фортуин, Кастелейн и Жинибр использовали условия зависимости случайных полей, задаваемые в виде неравенств для определенных мер на решетках (FKG-неравенства). Эти неравенства влекут ассоциированность и применяются в статистической физике, квантовой механике, теории гиббсовских полей, теории перколяции (см. [25], [47], [65]). Имеется следующее обобщение определения 0.1 на многомерный случай. Определение 0.2([38], [46]). Семейство X = {Xt,t Є Т} случайных векторов со значениями в R5, называется слабо ассоциированным (или положительно ассоциированным), если для любых конечных непересекающихся множеств индексов I,J С Ти ограниченных покоординатно неубывающих функций / : R5 7 —» R и д : RS J —справедливо неравенство. Йаг-Дев и Прошан ввели понятие отрицательно ассоциированных случайных векторов. Определение 0.3((57]). Семейство X = {Xt,t Є Т} случайных векторов со значениями в Rs, называется отрицательно ассоциированным, если для любых конечных непересекающихся множеств индексов I,JGT и ограниченных покоординатно неубывающих функций / : RS 7I —» R и д : RSIJI —» R справедливо неравенство cov(f(Xs,s Є I),g{Yt,t Є J)) 0. Со свойствами и примерами отрицательно ассоциированных случайных величин можно ознакомится в работе [57]. Приведем некоторые из них. Независимые случайные величины удовлетворяют определению 0.3. Отрицательно коррелированные гауссовские случайные величины являются отрицательно ассоциированными. Примеры отрицательно ассоциированных случайных величин возникают в математической статистике. Многие известные распределения обладают свойством отрицательной ассоциированности, например, мультиномиальное, многомерное гипергеометрическое и другие. Случайные процессы и поля, удовлетворяющие условиям зависимости 0.1-0.3, являются предметом активного изучения многих исследователей. В течение нескольких последних десятилетий были доказаны различные варианты центральной предельной теоремы, принципы инвариантности, законы повторного логарифма, законы больших чисел, экспоненциальные неравенства и другие результаты для сумм положительно и отрицательно ассоциированных случайных величин. При этом предполагалось, что зависимость между случайными величинами, выражаемая в терминах ковариаций, убывает с определенной скоростью, когда расстояния между индексами этих величин растут (см., например, [64], [33], [34], [5], [66], [28]). Много работ посвящено изучению свойств непараметрических оценок плотности, функции регрессии и их производных, построенных по выборкам, образующим положительно или отрицательно ассоциированные процессы и поля. Для этих оценок доказана асимптотическая нормальность слабая и сильная состоятельность (см. [74], [61], [60], [54]), получены оценки скорости нормального приближения (см. [10], [37], [20]). При доказательстве этих результатов, активно используются методы теории суммирования случайных величин. Для формулировки некоторых определений и результатов, которые нам потребуются в дальнейшем, введем ряд понятий и обозначений. Рассмотрим метрические пространства S и М соответственно с метриками бит. Функция / : S — М называется липшицевой, если Класс ограниченных липшицевых функций / : Sm — R, где m Є N, обозначим BL(Sm); в пространстве Sm используется метрика 5m(x,y) = ХГ=15(хк, Ук), Для х = (яь..., хт) Є Sm и у = (г/ь..., ут) Є 5m, а в R берется евклидово расстояние.

Для положительно и отрицательно ассоциированных случайных величин справедливы некоторые неравенства, которые связывают ковариацию функций от наборов случайных величин с суммой ковариаций самих величин (см. [65], [33], [35]). Выполнения этих неравенств часто бывает достаточно для доказательства предельных теорем, неравенств и других результатов. В связи с этим появились новые определения, обобщающие свойства положительной и отрицательной ассоциированности.

Неравенства для моментов сумм зависимых случайных величин

Возникает вопрос, при каких условиях ядерные оценки функции регрессии обладают свойством асимптотической нормальности (точнее говоря, величина л/пкп(гп(х) — г(х))/сг(х) распределена асимптотически нормально), если поле {(Xi,Yi),i Є Un} является (BL,ipi, -зависимым, a Un С Ъл. Можно ли в качестве множеств Un брать не только "целочисленные прямоугольники", а множества более общей структуры? Во второй главе данной работы даны условия, в частности, на скорость убывания величин вг, достаточные для получения оценки скорости нормального приближения совместного распределения ядерных оценок функции регрессии (0.0.11), вычисленных в конечном числе точек X Є Rs. При этом никаких ограничений на форму множеств Un не налагается. Для доказательства этого результата используется предложенный в работе [37] модифицированный метод Стейна, позволяющий оценивать приближение к нормальному закону случайных векторов.

Сильная состоятельность ядерных оценок плотности для положительно и отрицательно ассоциированных выборок изучалась в недавних работах Масри и Оливейра, см. [61], [54].

Пусть {- І}ІЄК - строго стационарная последовательность положительно (отрицательно) ассоциированных случайных величин. В статье [61] изучаются ядерные оценки плотности распределения случайных векторов Xj = (Xj+i,... ,Xj+m), т Є N, а также оценки ее производных. Предполагается, что плотность / случайных величин X,- ограничена. Кроме того, для всех j 0 существует совместная плотность fj векторов Хо и Xj (для j т плотность распределения вектора (XQ,. .. ,Xj)), причем fj ограничены равномерно по j. Также предполагается, что ядро К имеет ограниченные частные производные; В этой работе, в частности, доказано, что если в для любого компакта Q С Rm и всех 7 0.

При доказательстве этого результата использовались моментные неравенства для сумм положительно и отрицательно ассоциированных случайных величин из статьи [76]. В [54] для получения оценки скорости сходимости почти наверное к нулю величин \fn{x) — Е/П(х) использовались экспоненциальные неравенства для сумм ассоциированных случайных величин. При этом требовалась экспоненциальная скорость убывания коэффициента Кокса-Гримметта в = Х)ь гсог,( о)- і)) когда г, стремится к бесконечности. Полученная скорость сходимости имеет вид 0((lnn/n(1-7)/in) ) для некоторого 7 (0,1).

Состоятельность ядерных оценок плотности изучалась и для случайных полей. Так в работе [40] рассматривались оценки Парзена-Розенблатта плотности, построенные по наблюдениям, образующим векторное случайное поле, удовлетворяющее условию сильного перемешивания. При этом в качестве множеств Un для построения оценок (0.0.10) брались "регулярно растущие целочисленные прямоугольники".

В третьей главе диссертации с помощью новых экспоненциальных и моментных неравенств для сумм элементов (BL,ipi, -зависимых случайных полей доказывается сильная состоятельность ядерных оценок плотности. При этом никаких ограничений на форму множеств Un не налагается. Используемое здесь определение 0.7, позволяет учитывать различный характер зависимости поля по разным направлениям. Это важно для некоторых моделей, применяемых, например в геостатистике, где для описания зависимости гауссовского случайного поля {Z,-, і Є Zd} используют вариограмму (см., например, [42])

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

Из определения 0.4 видно, что, как и в случае гауссовского поля, зависимость стационарного квазиассоциированного случайного поля определяется его ковариационной функцией. Поэтому естественно возникает возможность учета анизотропии в близком к квазиассоциированности понятии (, , -зависимости.

Сделаем несколько замечаний. Далее буква С будет обозначать, вообще говоря, различные для разных формул положительные множители, не зависящие от п, а в некоторых случаях и от других параметров, что будет оговариваться отдельно. Поскольку все оценки в главах 1 и 2 доказываются с точностью до множителей, которые могут зависеть от размерности рассматриваемых конечномерных пространств, то выбор нормы в рассматриваемом пространстве не играет роли. Работа, объемом 82 страницы, состоит из введения, трех глав и списка литературы, насчитывающего 83 наименования.

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

Вторая глава также состоит из двух параграфов и посвящена изучению асимптотики совместных распределений ядерных оценок функции регрессии для слабо зависимых случайных полей. В 2.1 излагается метод Стейна и доказываются результаты, которые играют ключевую роль в 2.2 при получении оценки скорости нормального приближения ядерных оценок функции регрессии (теоремы 2.2.2 и 2.2.3).

Третья глава содержит три параграфа. В первом параграфе устанавливается новое экспоненциальное неравенство (лемма 3.1.1) для сумм элементов (BL, -0і,0)-зависимого случайного поля, во втором доказано неравенство для моментов сумм элементов (ВЬ+фі, #)-зависимого случайного поля (лемма 3.2.1), а в 3.3.3 на основе результатов 3.3.1 и 3.3.2 получены оценки скорости сходимости почти наверное ядерных оценок плотности (см. теоремы 3.3.2 и 3.3.3).

Основные результаты диссертации опубликованы в работах автора [10], [20], [62], [21]. В работе [10] А.В.Булинскому принадлежит постановка задачи и разработка подхода, основанного на использовании ковариационных неравенств, Н.В.Миллионщикову принадлежит асимптотический анализ вспомогательных функций и оптимизация используемых параметров. Все остальные результаты получены Н.В.Миллионщиковым самостоятельно. Результаты диссертации докладывались на международной конференции "Колмогоров и современная математика" (Москва, 2003), в 2005 году на Большом семинаре кафедры теории вероятностей механико-математического факультета МГУ (руков. чл.-корр. РАН, проф. А.Н.Ширяев), а также в 2004 и 2005 годах на семинаре "Асимптотический анализ случайных процессов и полей" (мех-мат МГУ, руков. проф. А.В.Булинский и проф. В.И.Питербарг). Автор благодарен своему научному руководителю профессору А.В.Булинскому за постановку задач, постоянное внимание к работе и ценные советы.

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