Sind unverzerrte effiziente Schätzer stochastisch dominant gegenüber anderen (mittleren) unverzerrten Schätzern?

Allgemeine Beschreibung

Maximiert ein effizienter Schätzer (dessen Stichprobenvarianz gleich der Cramér-Rao-Grenze ist) die Wahrscheinlichkeit, nahe am wahren Parameter ?$\theta$

Angenommen, wir vergleichen die Differenz oder die absolute Differenz zwischen der Schätzung und dem wahren Parameter

Ist die Verteilung von für einen effizienten Schätzer stochastisch dominant gegenüber der Verteilung von für einen anderen unverzerrten Schätzer?$\hat\Delta$$\tilde\Delta$

Motivation

Ich denke darüber nach wegen der Frage Schätzer, der unter allen sinnvollen Verlustfunktionen (Bewertungsfunktionen) optimal ist, wobei wir sagen können, dass der unverzerrte beste Schätzer in Bezug auf eine konvexe Verlustfunktion auch ein unvoreingenommener bester Schätzer in Bezug auf eine andere Verlustfunktion ist (From Iosif Pinelis, 2015, Eine Charakterisierung der besten unverzerrten Schätzer. ArXiv-Vorabdruck arXiv: 1508.07636 ). Die stochastische Dominanz für die Nähe zum wahren Parameter scheint mir ähnlich zu sein (dies ist eine ausreichende Bedingung und eine stärkere Aussage).

Genauere Ausdrücke

Die obige Fragestellung ist weit gefasst, z. B. welche Art von Unparteilichkeit wird berücksichtigt, und haben wir dieselbe Abstandsmetrik für negative und positive Unterschiede?

Betrachten wir die folgenden zwei Fälle: , um die Frage weniger umfassend zu gestalten:$^\dagger$

Vermutung 1: Wenn ein effizienter Mittelwert und ein Median-unverzerrter Schätzer ist. Dann für jeden Mittelwert und Median-unverzerrten Schätzer wobei und$\hat \theta$$\tilde \theta$

Vermutung 2: Wenn ein effizienter Mittelwert-unverzerrter Schätzer ist. Dann gilt für jeden mittleren unverzerrten Schätzer und$\hat \theta$$\tilde \theta$$x>0$

• Sind die obigen Vermutungen wahr?
• Wenn die Sätze zu stark sind, können wir sie anpassen, damit sie funktionieren?

^{$\dagger$ Der zweite bezieht sich auf den ersten, hebt jedoch die Einschränkung für die Median-Unvoreingenommenheit auf (und dann müssen wir beide Seiten zusammen nehmen, sonst wäre der Satz für jeden Schätzer falsch, der einen anderen Median als der effiziente Schätzer hat).}

Beispiel, Abbildung:

Betrachten sie die Schätzung des mittleren der Verteilung der Bevölkerung (die angenommen wird , normale verteilt) durch (1) die Probe Median und (2) die Probe Mittelwert.$\mu$

Im Fall einer Stichprobe der Größe 5 und wenn die wahre Verteilung der Population beträgt, sieht dies so aus$N(0,1)$

Im Bild sehen wir, dass der gefaltete CDF des Stichprobenmittelwerts (der ein effizienter Schätzer für ) unter dem gefalteten CDF des Stichprobenmedians liegt. Die Frage ist, ob der gefaltete CDF des Stichprobenmittelwerts auch unter dem gefalteten CDF eines anderen unverzerrten Schätzers liegt.$\mu$

Alternativ können wir bei Verwendung der CDF anstelle von gefalteten CDFs die Frage stellen, ob die CDF eines Mittelwerts den Abstand von 0,5 an jedem Punkt maximiert. Wir wissen, dass

Haben wir dies auch, wenn wir für die Verteilung eines anderen Mittelwerts und eines Median-unverzerrten Schätzers ersetzen ?$F_{median}(\hat \theta)$

$(X_1,\ldots,X_N)$$\mathcal{C}(\mu,1)$$\mu$

1. $\hat{\mu}_1= \text{median}(X_1,\ldots,X_N)=X_{(N/2)}$
2. $\hat{\mu}_2= \text{mean}(X_{(N/4)},\ldots,X_{(3N/4)})=\frac{2}{N}(X_{(N/4)}+\ldots+X_{(3N/4)})$
3. $\hat{\mu}_3=\mu^\text{MLE}$
4. $\hat{\mu}_4=\hat{\mu}_1+\frac{2}{N}\frac{\partial \ell}{\partial \mu}(\hat{\mu}_1)$

$\hat{\mu}_3$$\hat{\mu}_4$$\hat{\mu}_1$$\hat{\mu}_2$

Eine Darstellung der Unterschiede zum empirischen cdf der MLE macht es klarer:

Hier ist der entsprechende R-Code:

T=1e4
N=11
mlechy=function(x){
  return(optimize(function(theta) -sum(dcauchy(x, 
    location=theta, log=TRUE)),  c(-100,100))$minimum)
}
est=matrix(0,T,4)
for (t in 1:T){
cauc=sort(rcauchy(N))
est[t,1]=median(cauc)
est[t,2]=mean(cauc[4:8])
est[t,3]=mlechy(cauc)
est[t,4]=est[t,1]+(4/N)*sum((cauc-est[t,1])/(1+(cauc-est[t,1])^2))
}

plot(ecdf(est[,1]),col="steelblue",cex=.4,xlim=c(-1,1),main="",ylab="F(x)")
plot(ecdf(est[,2]),add=TRUE,col="sienna",cex=.4)
plot(ecdf(est[,3]),add=TRUE,col="gold",cex=.4)
plot(ecdf(est[,4]),add=TRUE,col="tomato",cex=.4)