Mittelwert der Bootstrap-Stichprobe vs. Statistik der Stichprobe

18

Angenommen, ich habe ein Beispiel und das Bootstrap-Beispiel aus diesem Beispiel für ein stastitisches (z. B. den Mittelwert). Wie wir alle wissen, schätzt dieses Bootstrap-Beispiel die Stichprobenverteilung des Schätzers der Statistik.χ

Ist der Mittelwert dieser Bootstrap-Stichprobe eine bessere Schätzung der Bevölkerungsstatistik als die Statistik der ursprünglichen Stichprobe ? Unter welchen Bedingungen wäre das der Fall?

Amelio Vazquez-Reina
quelle
2
Der Mittelwert des Bootstrap-Beispiels ist der Mittelwert des Beispiels. In diesem Fall benötigen Sie kein Bootstrap-Beispiel.
Xi'an,
1
Danke @ Xi'an Ich bin nicht sicher, ob ich folge. Der Mittelwert des Bootstrap-Beispiels kann sich numerisch vom Mittelwert des Beispiels unterscheiden. Versuchen Sie zu sagen, dass die beiden theoretisch immer noch gleichwertig sind? Können Sie an beiden Enden bestätigen?
Amelio Vazquez-Reina
2
Lassen Sie uns unsere Terminologie klarstellen: "Bootstrap-Beispiel" könnte sich entweder auf ein bestimmtes Beispiel mit Ersetzung aus den Daten beziehen oder auf eine (multivariate) Zufallsvariable, von der ein solches Beispiel als eine Realisierung angesehen würde. Sie haben Recht, dass der Mittelwert einer Realisierung vom Mittelwert der Daten abweichen kann, aber @ Xi'an liefert die relevantere Beobachtung, dass der Mittelwert der Zufallsvariablen (der per Definition die Bootstrap-Schätzung des Populationsmittelwerts ist ) übereinstimmen muss mit dem Mittelwert der Daten.
Whuber
1
Dann ist Ihre Frage fast identisch mit stats.stackexchange.com/questions/126633/… ; Der einzige Unterschied besteht darin, dass sich die Bootstrap-Beispielrealisierungen überschneiden können, die in der Antwort angegebene Analyse jedoch leicht auf die Bootstrap-Situation übertragen werden kann, mit demselben Ergebnis.
Whuber
1
Ich sehe die Verbindung @whuber, obwohl man in Bootstrap "Subsets with Replacement" hat und die Realisierungen sich überschneiden können, wie Sie sagten. Ich würde mir vorstellen, dass die Verteilung (z. B. Pseudozufälligkeit), die zum Abrufen der Neuabtastungen im Bootstrap verwendet wird, auch die Abweichung der Schätzung vom Bootstrap-Beispiel beeinflussen kann. Vielleicht ist die Antwort, dass der Unterschied in allen praktischen Angelegenheiten vernachlässigbar ist. Darauf kommt es an: Bedingungen, Feinheiten und der Unterschied in der Praxis.
Amelio Vazquez-Reina

Antworten:

19

Lassen Sie uns verallgemeinern, um uns auf den Kern der Sache zu konzentrieren. Ich werde die kleinsten Details darlegen, um keine Zweifel zu hinterlassen. Die Analyse erfordert nur Folgendes:

  1. Das arithmetische Mittel einer Menge von Zahlen ist definiert alsz1,,zm

    1m(z1++zm).
  2. Erwartung ist ein linearer Operator. Das heißt, wenn Zufallsvariablen und α i Zahlen sind, dann ist die Erwartung einer linearen Kombination die lineare Kombination der Erwartungen.Zi,i=1,,mαi

    E(α1Z1++αmZm)=α1E(Z1)++αmE(Zm).

Sei eine Stichprobe ( B 1 , , B k ), die aus einem Datensatz x = ( x 1 , , x n ) erhalten wird, indem k Elemente einheitlich aus x mit Ersetzung genommen werden. Lassen m ( B ) das arithmetische Mittel sein B . Dies ist eine Zufallsvariable. DannB(B1,,Bk)x=(x1,,xn)kxm(B)B

E(m(B))=E(1k(B1++Bk))=1k(E(B1)++E(Bk))

folgt aus der Linearität der Erwartung. Da die Elemente von alle auf die gleiche Weise erhalten werden, haben sie alle die gleiche Erwartung, b sagen:Bb

E(B1)==E(Bk)=b.

Dies vereinfacht das Vorhergehende

E(m(B))=1k(b+b++b)=1k(kb)=b.

Per Definition ist die Erwartung die wahrscheinlichkeitsgewichtete Summe von Werten. Da angenommen wird, dass jeder Wert von die gleiche Chance hat, dass 1 / n ausgewählt wird,X1/n

E(m(B))=b=E(B1)=1nx1++1nxn=1n(x1++xn)=x¯,

das arithmetische Mittel der Daten.

Zur Beantwortung der Frage, ob man das , um das Populationsmittel zu schätzen, ist das Bootstrap-Mittel (was der Fall ist, k = n ) ebenfalls gleich ˉ x und ist daher identisch mit einem Schätzer des Populationsmittels.x¯k=nx¯


Für Statistiken, die keine linearen Funktionen der Daten sind, gilt nicht unbedingt dasselbe Ergebnis. Es wäre jedoch falsch, einfach das Bootstrap-Mittel für den Statistikwert in den Daten zu ersetzen: So funktioniert Bootstrapping nicht. Stattdessen erhalten wir durch Vergleichen des Bootstrap-Mittels mit der Datenstatistik Informationen über die Abweichung der Statistik. Hiermit kann die ursprüngliche Statistik angepasst werden, um die Abweichung zu beseitigen. Somit wird die vorspannungskorrigierte Schätzung eine algebraische Kombination der ursprünglichen Statistik und des Bootstrap-Mittels. Weitere Informationen finden Sie unter "BCa" (Bias-korrigierter und beschleunigter Bootstrap) und "ABC". Wikipedia bietet einige Referenzen.

whuber
quelle
Sie meinen, dass die Erwartung des Bootstrap-Mittelwerts gleich dem Datenmittelwert ist, nein? Der Bootstrap-Mittelwert selbst wird nicht durch das (Original-) Datenmuster bestimmt.
Capybaralet
@ user2429920 Der Bootstrap-Mittelwert ist eine Statistik, die von der Stichprobe bestimmt wird. In diesem Sinne ist es identisch mit dem Stichprobenmittelwert. Ihre Erwartung wird im Sinne der Stichprobenverteilung genommen. Ich vermute, Sie verwenden "Erwartung" in einem anderen Sinne als den Vorgang der Berechnung des Bootstrap-Mittels durch wiederholtes Unterabtasten mit Ersetzen.
whuber
1
Ich denke, der letzte Absatz ist die eigentliche Antwort auf diese Frage, da er allgemein gehalten ist und sich nicht nur auf die Durchschnittsstatistik konzentriert. Ich hatte den gleichen Zweifel wie das OP und war mir der Existenz von BCa nicht bewusst. Obwohl die Demonstration in dieser Antwort mir nicht viel half (ich verwende den Mittelwert nicht als meine Statistik), war der letzte Absatz sehr klar über den Kern der Sache. Ich glaube, dass Xi'ans Antwort auch den Fall anspricht, in dem die Durchschnittsstatistik verwendet wird, also dasselbe Problem. Vielen Dank!
Gabriel
1
@ Gabriel gute Punkte. Ich habe den Datensatz überprüft: Vor der Bearbeitung wurde bei dieser Frage ursprünglich nur der Mittelwert abgefragt . Aus diesem Grund scheinen sich die Antworten auf diese Statistik zu konzentrieren.
whuber
9

F^n(x)=1ni=1nIXixXiiidF(x),
EF^n[X]=1ni=1nXi=X¯n
EF^n[X]X¯n
Xi'an
quelle
2
+1 Dies ist die Antwort, die ich ursprünglich schreiben wollte, befürchtete jedoch, dass sie für einige Leser zu undurchsichtig sein könnte. Ich bin trotzdem froh, dass es so elegant präsentiert wird. Ich bin mir jedoch nicht sicher, was Sie in Ihrem letzten Satz meinen, in dem Sie die "Erwartung" der simulierten Näherung anscheinend vom "Grenzwert" unterscheiden: Da die Erwartung konstant ist (sie variiert nicht mit der Simulationsgröße) ), es gibt wirklich keine Grenzen.
whuber
@whuber: Danke für den Kommentar und sorry, dass du meine knappe Antwort genau zur gleichen Zeit wie deine geschrieben hast! Ihre Erklärungen sind für Anfänger in Bootstrap sicherlich besser lesbar. Ich habe den letzten Satz korrigiert, dessen einschränkender Teil das Gesetz der großen Zahlen ist.
Xi'an,
3
Ihre Verwendung von "mean" in diesem letzten Satz ist ziemlich mehrdeutig! Ich habe es anhand Ihres LLN-Hinweises herausgefunden. Für jede endliche Simulation der Bootstrap-Verteilung erzeugt jede Stichprobe in der Simulation ihren eigenen Mittelwert (es gibt eine Bedeutung von "Mittelwert"). Der Durchschnitt aller dieser Proben in einer bestimmten Simulation einer Simulation erzeugt Mittelwert (es gibt eine andere Bedeutung). Die Simulation Mittelwert konvergiert auf eine Konstante als Simulationsgröße wächst groß, was die Bootstrap ist Mittelwert (eine dritte Bedeutung), und dies entspricht die Probe Mittelwert (die vierte Bedeutung). (Und das schätzt die Bevölkerung bedeutet paradoxerweise gerade fünfte Bedeutung!)
whuber