Sull'esistenza di UMVUE e la scelta dello stimatore della popolazione


10

Let essere un campione casuale tratto da della popolazione in cui .(X1,X2,,Xn)N(θ,θ2)θR

Sto cercando l'UMVUE di .θ

La densità congiunta di è(X1,X2,,Xn)

fθ(x1,x2,,xn)=i=1n1θ2πexp[12θ2(xiθ)2]=1(θ2π)nexp[12θ2i=1n(xiθ)2]=1(θ2π)nexp[1θi=1nxi12θ2i=1nxi2n2]=g(θ,T(x))h(x)(x1,,xn)Rn,θR

, dove e .h(x)=1g(θ,T(x))=1(θ2π)nexp[1θi=1nxi12θ2i=1nxi2n2]h(x)=1

Qui, dipende da \ theta e da x_1, \ cdots, da x_n a T (\ mathbf x) = \ left (\ sum_ {i = 1} ^ nx_i, \ sum_ {i = 1} ^ nx_i ^ 2 \ right) e h è indipendente da \ theta . Quindi, secondo il teorema di fattorizzazione di Fisher-Neyman, la statistica bidimensionale T (\ mathbf X) = \ left (\ sum_ {i = 1} ^ nX_i, \ sum_ {i = 1} ^ nX_i ^ 2 \ right) è sufficiente per \ theta .θ x 1 , , x n T ( x ) = ( n i = 1 x i , n i = 1 x 2 i ) h θ T ( X ) = ( n i = 1 X i , n i = 1 X 2 i ) θgθx1,,xnT(x)=(i=1nxi,i=1nxi2)hθT(X)=(i=1nXi,i=1nXi2)θ

Tuttavia, T non è una statistica completa. Questo perché

Eθ[2(i=1nXi)2(n+1)i=1nXi2]=2n(1+n)θ2(n+1)2nθ2=0θ

e la funzione non è identicamente zero.g(T(X))=2(i=1nXi)2(n+1)i=1nXi2

Ma so che è una statistica minima sufficiente.T

Non ne sono certo, ma penso che potrebbe non esistere una statistica completa per questa famiglia esponenziale curva. Quindi, come dovrei ottenere l'UMVUE? Se non esiste una statistica completa, può essere uno UMVUE uno stimatore imparziale (come in questo caso) che è una funzione di statistica minima sufficiente? ( Discussione correlata: qual è la condizione necessaria affinché uno stimatore imparziale sia UMVUE? )X¯

Cosa succede se considero il miglior stimatore lineare imparziale (BLU) di ? Il BLU può essere l'UMVUE?θ

Supponiamo che io consideri lo stimatore lineare imparziale di dove e . Poiché sappiamo che . La mia idea è di minimizzare modo che sia il BLU di . Sarebbe essere poi l'UMVUE di ?θ c ( n ) = T(X)=aX¯+(1a)cSθS2=1c(n)=n12Γ(n12)Γ(n2)Eθ(cS)=θVar(T)TθTθS2=1n1i=1n(XiX¯)2Eθ(cS)=θVar(T)TθTθ

Ho preso uno stimatore imparziale lineare basato su e poiché è anche sufficiente per . S( ˉ X ,S2)θX¯S(X¯,S2)θ

Modificare:

Molto lavoro è stato fatto nella stima di nella più generale famiglia dove è noto . Di seguito sono riportati alcuni dei riferimenti più rilevanti:θN(θ,aθ2)a>0

Ho trovato il primo di questi riferimenti in questo esercizio di Statistical Inference di Casella / Berger:

inserisci qui la descrizione dell'immagine

inserisci qui la descrizione dell'immagine

La mia domanda non riguarda questo esercizio però.

La nota finale (l'estratto del capitolo) afferma che l' UMVUE di non esiste inθ quanto la statistica minima sufficiente non è completa. Vorrei sapere cosa ci consente di concludere che un UMVUE non esiste semplicemente perché non è possibile trovare una statistica completa completa? C'è qualche risultato correlato in merito? Vedo l'esistenza di UMVUE anche quando non esistono statistiche complete complete nel thread collegato.

Ora supponendo che non esista uno stimatore imparziale di varianza uniformemente minima, quali dovrebbero essere i nostri prossimi criteri per scegliere lo stimatore "migliore"? Cerchiamo l'MSE minimo, la varianza minima o l'MLE? O la scelta dei criteri dipenderà dal nostro scopo di stima?

Ad esempio, supponiamo di avere uno stimatore imparziale e un altro stimatore distorto di . Supponiamo che l'MSE di (che è la sua varianza) sia maggiore di quello di . Poiché la minimizzazione di MSE significa minimizzare la distorsione e la varianza contemporaneamente, penso che dovrebbe essere una scelta 'migliore' di stimatore rispetto a sebbene il primo sia distorto.T1T2θT1T2T2T1

Le possibili scelte degli stimatori di sono elencate dalla pagina 4 dell'ultima nota.θ

Il seguente estratto è tratto da Theory of Point Estimation di Lehmann / Casella (seconda edizione, pagina 87-88):

inserisci qui la descrizione dell'immagine

inserisci qui la descrizione dell'immagine

È molto probabile che io abbia frainteso tutto, ma l'ultima frase dice che in determinate condizioni, l'esistenza di una statistica completa è necessaria per l'esistenza di UMVUE? Se è così, è questo il risultato che dovrei guardare?

L'ultimo risultato dovuto a RR Bahadur, menzionato alla fine, fa riferimento a questa nota.

Dopo ulteriori ricerche, ho trovato un risultato che afferma che se la statistica minima sufficiente non è completa, allora non esiste una statistica completa. Quindi almeno sono quasi convinto che qui non esista una statistica completa.

Un altro risultato che ho dimenticato di considerare è quello che dice approssimativamente che una condizione necessaria e sufficiente affinché uno stimatore imparziale sia l'UMVUE è che deve essere non correlato con ogni stimatore imparziale pari a zero. Ho provato a usare questo teorema per dimostrare che qui non esiste un UMVUE, e anche il fatto che uno stimatore imparziale come non è l'UMVUE. Ma questo non funziona come semplice come fatto, ad esempio qui , nell'illustrazione finale.X¯

Risposte:


3

Aggiornare:

Si consideri lo stimatore 0 = ˉ X - c S dove c è dato nel tuo post. Questo è uno stimatore imparziale di 0 e sarà chiaramente correlato con lo stimatore indicato di seguito (per qualsiasi valore di a ).

0^=X¯cS
c0a

Il teorema 6.2.25 di C&B mostra come trovare statistiche complete complete per la famiglia esponenziale fintanto che contiene un set aperto in R k . Sfortunatamente questa distribuzione produce w 1 ( θ ) = θ - 2 e w 2 ( θ ) = θ - 1 che NON forma un insieme aperto in R 2 (poiché

{(w1(θ),wk(θ)}
Rkw1(θ)=θ2w2(θ)=θ1R2 ). È per questo che la statistica ( ˉ X , S 2 ) non è completa per θ , ed è per lo stesso motivo che possiamo costruire uno stimatore imparziale di 0 che sarà correlato con qualsiasi stimatore imparziale di θ che si basa sulle statistiche sufficienti.w1(θ)=w2(θ)2(X¯,S2)θ0θ

Un altro aggiornamento:

θ~Var(θ~)<Var(θ^)θΘ

E(θ^)=θE(0^)=0Cov(θ^,0^)<0θ

θ~=θ^+b0^
Var(θ~)=Var(θ^)+b2Var(0^)+2bCov(θ^,0^)
M(θ)=2Cov(θ^,0^)Var(0^)

θ0M(θ0)>0B(0,M(θ0))Vun'r(θ~)<Vun'r(θ^) θ0θ^

θ^0^aθ^ θ0θ^


Diamo un'occhiata più da vicino alla tua idea di combinazioni lineari.

θ^=aX¯+(1a)cS

θ^

MSE(θ^)=a2Var(X¯)+(1a)2c2Var(S)=a2θ2n+(1a)2c2[E(S2)E(S)2]=a2θ2n+(1a)2c2[θ2θ2/c2]=θ2[a2n+(1a)2(c21)]

an

aopt(n)=c211/n+c21
c2=n12(Γ((n1)/2)Γ(n/2))2

ainserisci qui la descrizione dell'immagine

naopt13

Sebbene non vi sia alcuna garanzia che si tratti dell'UMVUE, questo stimatore è lo stimatore di varianza minima di tutte le combinazioni lineari imparziali delle statistiche sufficienti.


Grazie per l'aggiornamento. Non ho seguito C&B come un libro di testo, ho solo guardato gli esercizi.
Testardo

1
θ^
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.