Лекция 5
Експоненциално семейство

Това семейство разпределения заема изключително положение в математическата статистика по две причини. Първо, в него влизат почти всички разпределения с аналитична форма на запис и, второ, както ще видим в тази лекция , то притежава всички необходими теоретични качества.

Към експоненциалното семейство спадат гаусовото, бета, гама, стюдент, фишер, поасоново, биномно и редица други разпределения.

Определение 1 Ще наричаме експоненциално семейство E семейството от разпределения с плътности от вида

f(x,q) = h(x) exp((a(q),U(x)) +V(q))
(5.1)
където q О Q М Rk, a(q), U(x) О Rk

Преобразованието a:Rk® Rk, определено в някаква област Q, обикновено се предполага взаимно еднозначно на Q. Ние ще предположим, че функциите a0(q) = 1,a1(q),a2(q),...,ak(q) са линейно независими.

Функцията на правдоподобие на експоненциалното разпределение ще има вида

L(x,q) = n
Х
i = 1 
f(xi,q) = ( n
Х
i = 1 
h(xi) )exp((a(q), n
е
i = 1 
U(xi)) + nV(q)).
(5.2)

5.1  Минималност на S

Лесно се вижда, че разпределението породено от n наблюдения в има проста форма, която може да се изрази по следния начин. Да означим с S = еni = 1 U(xi). Да разгледаме функцията
r(x,q) = fq(x)
fq0(x)
= exp( (a(q)- a(q0),S)+n(V(qj) - V(q0))).

Значи можем да разгледаме m = fq0(x). От теорема 4.2 следва, че породената от нея s-алгебра s(h(x,q) е минимална достатъчна.

Теорема 1 Статистиката S е минимална достатъчна.

Доказателство:  От линейната независмост на 1,a1(q),a2(q),...,ak(q)

следва линейна независимост на a1(q)-a1(q0),a2(q)-a2(q0),...,ak(q)-ak(q0) Значи ще се намерят в Rk точки q1,q2,...,qk, такива че матрицата A = (ai,j = ai(qj)-ai(q0) е неизродена. Следователно системата:


е
i 
si ai,j = lnr(x,qj)-n(V(qj) - V(q0))
е разрешима еднозначно относно S. Това означава, че s(S) М s(h(x,q)) и значи е минимална достатъчна. Q.E.D.

5.2  Пълнота на S

Теорема 2 Ако образът на преобразованието a(Q) съдържа непразен паралелопипед, то статистиката S е пълна.

Доказателство: (Теорема 5.2) В нашия случай теоремата на Нейман - Фишер за факторизацията на правдоподобието по лебеговата мярка води до представянето:

g(s,q) = exp(a(q),s)+nV(q))
h(x) = n
Х
i = 1 
h(xi)

Лема 1 Нека G1 и G2 са две s-крайни мярки в Rk такива, че интегралите съществуват и са равни:

у
х
e(a,x) d G1 = у
х
e(a,x) d G2
за "a от някой непразен куб в Rk. Тогава G1 = G2.

Доказателство: Нека първо кубът има вида |aj| < a. Тъй като функциите са аналитични по а следва, че за всяко -Ґ < b < Ґ ще бъде изпълнено:

у
х
ei(b,x) d G1 = у
х
ei(b,x) d G2
От теоремата за взаимно еднозначното съответствие на х.ф. и ф.р. в Rk следва твърдението на лемата. Ако кубът има друг център |aj-cj| < a, то следва да разгледаме мерките G* = ei(c,x)G . Q.E.D.

Да се върнем към доказателството на теорема 5.2. Можем да разгледаме в Rk не зависещата от q мярка:

n(B) = у
х


S(-1(B) 
h(x) dl.
Трябва да покажем, че от
у
х
f(x) d Gq = 0     "q О Q
следва f(x) за всяка измерима функция f. Първо f = f+ - f-. Значи

у
х
f+(x) d Gq = у
х
f-(x) d Gq
у
х
g(s,q)f+(x) dn = у
х
g(s,q)f-(x) dn
у
х
g(s,q) dn+ = у
х
g(s,q)(x) dn-
у
х
exp(a(q),s)dn+ = у
х
exp(a(q),s)dn-
Сега твърдението следва от лемата.Q.E.D.

Пример 1 Гама разпределение

Сега имаме два параметъра: q = {a,l}.

f(x,q) = la
G(a)
xa-1 e-lx,     x > 0.
(5.3)

Получаваме: h(x) = x-1, S(x) = {U1(x) = lnx, U2(x) = x}, V(a,l) = [(la)/( G(a))], a1(q) = l, a2(q) = - a.


Начало на лекцията | Съдържание | Индекс


File translated from TEX by TTH, version 2.10.
On 4 Jun 1999, 15:57.