Erwarteter Wert der Statistik der Mindestbestellmenge aus einer normalen Stichprobe

9

UPDATE 25. Januar 2014: Der Fehler wurde behoben. Bitte ignorieren Sie die berechneten Werte des erwarteten Werts im hochgeladenen Bild - sie sind falsch - ich lösche das Bild nicht, weil es eine Antwort auf diese Frage generiert hat.

UPDATE 10. Januar 2014: Der Fehler wurde gefunden - ein Tippfehler in einer der verwendeten Quellen. Korrektur vorbereiten ...

Die Dichte der Statistik minimaler Ordnung aus einer Sammlung von iid kontinuierlichen Zufallsvariablen mit cdf und pdf beträgt F X ( x ) f X ( x ) f X ( 1 ) ( x ( 1 ) ) = n f X ( x ( 1 ) ) [ 1 - F X ( x ( 1 ) ) ] n - 1nFX(x)fX(x)

fX(1)(x(1))=nfX(x(1))[1FX(x(1))]n1[1]

Wenn diese Zufallsvariablen Standardnormal sind, dann

fX(1)(x(1))=nϕ(x(1))[1Φ(x(1))]n1=nϕ(x(1))[Φ(x(1))]n1[2]
und daher ist sein erwarteter Wert
E(X(1))=nx(1)ϕ(x(1))[Φ(x(1))]n1dx(1)[3]

wo wir die symmetrischen Eigenschaften der Standardnormalen verwendet haben. In Owen 1980 , S.402, Gleichung [ n, 011 ] finden wir, dass

zϕ(z)[Φ(az)]mdz=am(a2+1)(2π)ϕ(z)[Φ(aza2+1)]m1dz[4]

Übereinstimmende Parameter zwischen den Gleichungen und ( , ) erhalten wir[3][4]a=1m=n1

E(X(1))=n(n1)2πϕ(x(1))[Φ(x(1)2)]n2dx(1)[5]

Wieder in Owen 1980, p. 409, Gleichung [ n0,010.2 ] finden wir das

[i=1mΦ(hidiz1di2)]ϕ(z)dz=Zm(h1,...,hm;{ρij})[6]

wobei die multivariate Standardnormalen ist, sind die paarweisen Korrelationskoeffizienten und .Zm()ρij=didj,ij1di1

Wenn wir und übereinstimmen, haben wir , und [5][6]m=n2hi=0,i

di1di2=12di=±13iρij=ρ=1/3

Mit Hilfe dieser Ergebnisse, eq wird[5]

E(X(1))=n(n1)2πZn2(0,...,0;ρ=1/3)[7]

Dieses multivaririate Standard-Normalwahrscheinlichkeitsintegral von äquikorrelierten Variablen, die alle bei Null bewertet wurden , wurde ausreichend untersucht, und es wurden verschiedene Methoden zur Approximation und Berechnung abgeleitet. Eine ausführliche Übersicht (im Zusammenhang mit der Berechnung multivariater Normalwahrscheinlichkeitsintegrale im Allgemeinen) ist Gupta (1963) . Gupta bietet explizite Werte für verschiedene Korrelationskoeffizienten und für bis zu 12 Variablen (es deckt also eine Sammlung von 14 Variablen ab). Die Ergebnisse sind (DIE LETZTE SPALTE IST FALSCH) :

Geben Sie hier die Bildbeschreibung ein

Wenn wir nun grafisch darstellen, wie sich der Wert von mit ändert , erhalten wirZn2(0,...,0;ρ=1/3)n

Geben Sie hier die Bildbeschreibung ein

Ich komme also zu meinen drei Fragen / Wünschen:
1) Könnte jemand analytisch prüfen und / oder durch Simulation überprüfen, ob die Ergebnisse für den erwarteten Wert korrekt sind (dh die Gültigkeit von Gleichung prüfen )?[7]

2) Könnte jemand unter der Annahme, dass der Ansatz korrekt ist, die Lösung für Normalen mit einem Mittelwert ungleich Null und einer nicht einheitlichen Varianz geben? Bei all den Transformationen ist mir wirklich schwindelig.

3) Der Wert des Wahrscheinlichkeitsintegrals scheint sich reibungslos zu entwickeln. Wie wäre es mit einer Funktion von approximieren ?n

Alecos Papadopoulos
quelle

Antworten:

6

Ihre Ergebnisse erscheinen nicht korrekt. Dies ist ohne Berechnung leicht zu erkennen, da in Ihrer Tabelle Ihr mit der Stichprobengröße zunimmt . Es ist klar, dass der erwartete Wert des Stichprobenminimums kleiner werden muss (dh negativer werden muss), wenn die Stichprobengröße größer wird.E[X(1)] nn

Das Problem ist konzeptionell recht einfach.

Kurz gesagt: wenn ~ mit pdf :XN(0,1)f(x)

Geben Sie hier die Bildbeschreibung ein

... dann lautet das PDF der Statistik 1. Ordnung (in einer Stichprobe der Größe ):n

Geben Sie hier die Bildbeschreibung ein

... hier mit der OrderStatFunktion in erhalten mathStatica, mit Bereich der Unterstützung:

Geben Sie hier die Bildbeschreibung ein

Dann kann für leicht genau erhalten werden als:E[X(1)]n=1,2,3

Geben Sie hier die Bildbeschreibung ein

Der genaue Fall ist ungefähr , was sich offensichtlich von Ihrer Arbeitsweise von -1,06 (Zeile 1 Ihrer Tabelle) unterscheidet. Es scheint also klar zu sein, dass etwas mit Ihrer Arbeitsweise nicht stimmt (oder vielleicht mein Verständnis dessen, was Sie suchen). .n=30.846284

Für ist es schwieriger, Lösungen in geschlossener Form zu erhalten, aber selbst wenn sich die symbolische Integration als schwierig erweist, können wir die numerische Integration immer verwenden (auf Wunsch mit beliebiger Genauigkeit). Das ist wirklich sehr einfach ... hier ist zum Beispiel für die Stichprobengröße bis 14 unter Verwendung von Mathematica :n4E[X(1)]n=1

 sol = Table[NIntegrate[x g, {x, -Infinity, Infinity}], {n, 1, 14}]

{0, -0,56419, -0,846284, -1,02938, -1,16296, -1,26721, -1,35218, -1,4236, -1,48501, -1,53875, -1,58644, -1,62923, -1,66799, -1,70338}

Alles erledigt. Diese Werte unterscheiden sich offensichtlich stark von denen in Ihrer Tabelle (rechte Spalte).

Um den allgemeineren Fall eines Elternteils , gehen Sie genau wie oben vor, beginnend mit dem allgemeinen Normal-PDF.N(μ,σ2)

Wölfe
quelle
Danke für die Antwort. In der Tat habe ich auch festgestellt, dass etwas mit den numerischen Ergebnissen nicht stimmt - schließlich sollte der erwartete Wert in absoluter Größe zunehmen und nicht abnehmen, wenn zunimmt. Ich habe die Antwort unverändert gelassen, um zu sehen, ob ich aus einer Antwort einen Einblick gewinnen kann. Ich suche immer noch auf der theoretischen Ebene, wo genau der Fehler liegt. Der Verdächtige ist die erste Gleichung, die ich von Owen verwende (weil die zweite von anderen Quellen verifiziert wurde). Könnten Sie übrigens überprüfen, ob diese Gleichung in Mein Beitrag (als eigenständige Transformation) ist korrekt? Ich wäre dankbar. 4n4
Alecos Papadopoulos