Indipendenza delle statistiche dalla distribuzione gamma


9

Sia un campione casuale dalla distribuzione gamma .G a m m a ( α , β )X1,...,XnGamma(α,β)

Sia X¯ e S2 la media e la varianza del campione.

Quindi prova o confuta che X¯ e S2/X¯2 sono indipendenti.


Il mio tentativo: poiché S2/X¯2=1n1i=1n(XiX¯1)2 , dobbiamo verificare l'indipendenza di X¯ e (XiX¯)i=1n , ma come dovrei stabilire l'indipendenza tra loro?


2
Si consideri il Laplace congiunta trasformata di somma U:=iXi e il vettore W di proporzioni Wi:=Xi/U . Questo è E{exp[tUzW]} ; si può dimostrare che questo è il prodotto di una funzione di t e una funzione di z .
Yves,

@Yves Potresti controllare la mia risposta postata di seguito?
Bellcircle,

Risposte:


4

C'è una dimostrazione carina, semplice, intuitivamente ovvia per integraleα. Si basa solo su proprietà ben note della distribuzione uniforme, distribuzione Gamma, processi di Poisson e variabili casuali e va così:

  1. Ogni è il tempo di attesa fino a quando si verificano i punti di un processo di Poisson. αXiα

  2. La somma è quindi il tempo di attesa fino a quando si verificano punti di quel processo. Chiamiamo questi punti n α Z 1 , Z 2 , , Z n α .Y=X1+X2++XnnαZ1,Z2,,Znα.

  3. In base a , i primi punti sono distribuiti in modo indipendente uniformemente tra en α - 1 0 Y .Ynα10Y.

  4. Pertanto i rapporti sono distribuiti in modo indipendente uniformemente tra e In particolare, le loro distribuzioni non dipendono da0 1. Y .Zi/Y, i=1,2,,nα101.Y.

  5. Di conseguenza, qualsiasi funzione (misurabile) di è indipendente daY .Zi/YY.

  6. Tra queste funzioni ci sono (dove le parentesi indicano le statistiche dell'ordine di ). []Zi

    X1/Y=Z[α]/YX2/Y=Z[2α]/YZ[α]/YXn1/Y=Z[(n1)α]/YZ[(n2)α]/YXn/Y=1Z[(n1)α]/Y
    []Zi

A questo punto, nota semplicemente che può essere scritto esplicitamente come una funzione (misurabile) di e quindi è indipendente daX i / Y ˉ X = Y / n .S2/X¯2Xi/YX¯=Y/n.


3

Volete dimostrare che la media e il rv.s sono indipendenti, o equivalentemente che la somma e i rapporti sono indipendente. Possiamo dimostrare un risultato leggermente più generale supponendo che abbia probabilmente diverse forme , ma la stessa scala che si può presumere essere . nXi/ ˉ X U:=XinWi:=XiX¯nXi/X¯U:=XinX i α i β > 0 β = 1Wi:=Xi/UXiαiβ>0β=1

Considera la trasformata congiunta di Laplace di e ie, Questo si esprime come un integrale dimensionale over dove la costante è relativa a . Se introduciamo nuove variabili sotto il segno integrale impostando U ψ (W=[Wi]i=1nn(0,)nCst

ψ(t,z):=E{exp[tUzW}=E{exp[tiXiiziXiU]}
n(0,)n
x y : = ( 1 + t )
Cstexp[-(1+t)(X1++Xn)-z1X1++znXnX1++Xn]X1α1-1...Xnαn-1dX
X t z U Wy: =(1+t)X , vediamo facilmente che l'integrale può essere scritto come prodotto di due funzioni, una seconda dall'altro a seconda del vettore . Ciò dimostra che e sono indipendenti.tzUW

Disclaimer . Questa domanda riguarda il teorema di Lukacs sull'indipendenza proporzionale , quindi l'articolo di Eugene Lukacs A Characterization of the Gamma Distribution . Ho appena estratto qui la parte rilevante di questo articolo (vale a dire p. 324), con alcune modifiche alle notazioni. Ho anche sostituito l'uso della funzione caratteristica con quella della trasformata di Laplace per evitare cambiamenti di variabili che coinvolgono numeri complessi.


1
(+1) Per l'articolo sulla caratterizzazione della distribuzione gamma.
Testardo:

1

Lascia . Si noti che è una statistica accessoria di , ovvero la sua distribuzione non dipende da . ( X i / U ) i β βU=ΣioXio(Xio/U)ioββ

Poiché è una statistica sufficientemente completa di , è indipendente da dal teorema di Basu, quindi la conclusione segue.β ( X i / U ) iUβ(Xio/U)io

Non sono sicuro della costruzione della statistica accessoria, poiché è indipendente solo da , non da .αβα


Buona. Il teorema può essere invocato con considerato come fisso, quindi considerando un modello statistico a un parametro. α
Yves,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.