Zufallszahlen manuell generieren

30

Wie kann ich aus einer gegebenen Verteilung manuell eine Zufallszahl generieren, beispielsweise 10 Realisierungen aus der Standardnormalverteilung?

Dan ish
quelle
10
Können Sie erklären, warum Sie dies tun möchten und welche anderen Einschränkungen Sie haben?
mdewey
2
Sie können eine Tabelle mit Zufallszahlen abrufen (Rand wird unter anderem zum Veröffentlichen verwendet).
Batman
4
@Batman: Ja, und das RAND-Buch mit 10⁶ Zufallszahlen erhielt 655 Kommentare bei Amazon. Alles natürlich vorhersehbar.
Xi'an
10
(Ich bin überrascht, dass niemand dies zuvor kommentiert hat.) Sofern es nicht unbedingt erforderlich ist, sollte man nicht versuchen, benutzerdefinierte Implementierungen zum Generieren von Zufallszahlen zu verwenden. Ja, es ist großartig zu wissen, wie man es macht, und es ist wahrscheinlich das erste, was Sie in einer Klasse mit " MC-Methoden " sehen werden (und das ist großartig!), Aber machen Sie es nicht in realen Projekten. Es gibt aus einem Grund spezialisierte Routinen zur Erzeugung von Zufallszahlen. Die Implementierung eines Kernalgorithmus von Grund auf ist ein Zeitverlust, eine Fehlerquelle und ein Mangel an Feldbewusstsein.
usεr11852 sagt Reinstate Monic
2
@ Xi'an: Ja, ich stimme zu, das ist eine sichere Vermutung. Wie bereits erwähnt, ist dies nur ein Kommentar, um die Leute davor zu warnen, ihr eigenes RNG zu verwenden, ohne zu bemerken, dass RNG-Design ein sehr ernstes Geschäft ist. Um zitiert von Neumann : „Jeden, der zur Herstellung von Zufallszahlen rechnerische Methoden betrachtet ist natürlich, in einem Zustand der Sünde.“
usεr11852 sagt Reinstate Monic

Antworten:

46

Wenn "manuell" "mechanisch" einschließt, stehen Ihnen viele Optionen zur Verfügung. Um eine Bernoulli-Variable mit der Hälfte der Wahrscheinlichkeit zu simulieren, können wir eine Münze werfen: für Schwänze, 1 für Köpfe. Um eine geometrische Verteilung zu simulieren, können wir zählen, wie viele Münzwürfe benötigt werden, bevor wir Köpfe erhalten. Um eine Binomialverteilung zu simulieren, können wir unsere Münze n- mal werfen (oder einfach n Münzen werfen ) und die Köpfe zählen. Die "quincunx" oder "bean machine" oder "Galton box" ist eine kinetischere Alternative - warum nicht eine in Aktion setzen und selbst sehen ? Es scheint, als gäbe es keine "gewichtete Münze"01nnWenn wir jedoch den Wahrscheinlichkeitsparameter unserer Bernoulli- oder Binomialvariablen auf andere Werte als ändern möchten, können wir dies mit der Nadel von Georges-Louis Leclerc, Comte de Buffon , tun. Um die diskrete Gleichverteilung auf { 1 , 2 , 3 , 4 , 5 , 6 } zu simulieren , werfen wir einen sechsseitigen Würfel. Fans von Rollenspielen werden auf exotischere Würfel gestoßen sein , zum Beispiel tetraedrische Würfel, die gleichmäßig aus { 1 , 2 , 3 , 4 } entnommen werden können.p=0.5{1,2,3,4,5,6}{1,2,3,4}Während mit einem Spinner oder Roulette-Rad kann man noch weiter gehen. ( Bildnachweis )

Vielzahl von Würfeln

Müssten wir heute verrückt sein, um auf diese Weise Zufallszahlen zu generieren, wenn es nur ein Befehl auf einer Computerkonsole ist - oder wenn wir eine geeignete Tabelle mit Zufallszahlen zur Verfügung haben, einen Streifzug zu den staubigeren Ecken des Bücherregals? Nun vielleicht, obwohl ein physikalisches Experiment etwas angenehm Tastbares hat. Für Menschen, die vor dem Computerzeitalter arbeiteten, war es jedoch praktischer, Zufallsvariablen manuell zu simulieren, und zwar vor weit verbreiteten Zufallszahltabellen (von denen später mehr zur Verfügung stehen). Als Buffon das St. Petersburg-Paradoxon untersuchte- das berühmte Münzwurfspiel, bei dem sich der Gewinn jedes Mal verdoppelt, wenn ein Kopf geworfen wird, der Spieler bei den ersten Schwänzen verliert und dessen erwartete Auszahlung kontraintuitiv unendlich ist - er musste die geometrische Verteilung mit simulieren = 0,5 . Zu diesem Zweck hat er anscheinend ein Kind angeheuert, das eine Münze wirft, um 2048 Spiele des St. Petersburg-Spiels zu simulieren. Diese simulierte geometrische Verteilung ist in Stigler (1991) wiedergegeben :p=0.5

Tosses Frequency
1      1061
2      494
3      232
4      137
5      56
6      29
7      25
8      8
9      6

In demselben Aufsatz, in dem er diese empirische Untersuchung des St. Petersburg-Paradoxons veröffentlichte, stellte Buffon auch die berühmte " Buffon's Needle " vor. Wenn eine Ebene in einem Abstand durch parallele Linien in Streifen unterteilt ist und eine Nadel der Länge l d darauf fällt, beträgt die Wahrscheinlichkeit, dass die Nadel eine der Linien kreuzt, 2 ldld .2lπd

Buffons Nadelexperiment

Mit Buffons Nadel kann daher eine Zufallsvariable X Bernoulli ( 2 l) simuliert werdenoderXBinomial(n,2lXBernoulli(2lπd), und wir können die Erfolgswahrscheinlichkeit anpassen, indem wir die Länge unserer Nadeln oder (vielleicht bequemer) den Abstand ändern, in dem wir die Linien anordnen. Eine alternative Verwendung von Buffon-Nadeln ist ein äußerst ineffizienter Weg, um eine probabilistische Approximation fürπ zu finden. Das Bild (Kredit) zeigt 17 Streichhölzer, von denen 11 eine Linie kreuzen. Wenn der Abstand zwischen den Linien gleich der Länge des Streichholzes eingestellt ist, wie hier, beträgt der erwartete Anteil an sich kreuzenden Streichhölzern2XBinomial(n,2lπd)π und somit können wir abschätzen , π als doppelt reziproken Wert der beobachteten Fraktion: hier erhält man π =2172ππ^. Verwendung2,5 cm Nadeln mit Linien 3 cm voneinander entfernt 1901 beanspruchte Mario Lazzarinium das Experiment durchgeführt wird, und nach 3408 Würfen erhielten π =355π^=217113.1 . Dies ist ein bekanntes Rational fürπ, das auf sechs Dezimalstellen genau ist. Badger (1994) liefert überzeugende Beweise dafür, dass dies betrügerisch war, nicht zuletzt, dass, um mit Lazzarinis Apparat eine Genauigkeit von sechs Dezimalstellen zu 95% zu erreichen, ein Geduldsspiel von 134 Billionen Nadeln geworfen werden muss! Mit Sicherheit ist Buffons Nadel als Zufallszahlengenerator nützlicher als als Methode zur Schätzung vonπ.π^=355113ππ


Unsere Generatoren waren bisher enttäuschend diskret. Was ist, wenn wir eine Normalverteilung simulieren wollen? Eine Möglichkeit besteht darin, zufällige Ziffern zu erhalten und daraus gute diskrete Näherungen für eine gleichmäßige Verteilung auf bilden. Anschließend werden einige Berechnungen durchgeführt, um diese in zufällige normale Abweichungen umzuwandeln. Ein Dreh- oder Rouletterad kann Dezimalstellen von null bis neun angeben. ein Würfel kann binäre Ziffern erzeugen; Wenn unsere arithmetischen Fähigkeiten mit einer funkigeren Basis fertig werden, würde sogar ein Standardsatz Würfel ausreichen. Andere Antworten haben diese Art von transformationsbasiertem Ansatz ausführlicher behandelt. Ich verschiebe jede weitere Diskussion bis zum Ende.[0,1]

Im späten neunzehnten Jahrhundert war die Nützlichkeit der Normalverteilung bekannt, und so gab es Statistiker, die zufällige normale Abweichungen simulieren wollten. Es erübrigt sich zu erwähnen, dass langwierige Handberechnungen nicht geeignet gewesen wären, außer den Simulationsprozess an erster Stelle einzurichten. Sobald dies festgestellt wurde, musste die Erzeugung der Zufallszahlen relativ schnell und einfach sein. Stigler (1991) listet die Methoden auf, die drei Statistiker dieser Ära anwenden. Alle untersuchten Glättungstechniken: Zufällige normale Abweichungen waren von offensichtlichem Interesse, z. B. um Messfehler zu simulieren, die geglättet werden mussten.

Der bemerkenswerte amerikanische Statistiker Erastus Lyman De Forest war an der Glättung von Lebenstabellen interessiert und stieß auf ein Problem, das die Simulation der absoluten Werte normaler Abweichungen erforderte. In einem Thema, das sich als laufend erweisen wird, hat De Forest tatsächlich eine Stichprobe aus einer Halbnormalverteilung gezogen . Außerdem, anstatt eine Standardabweichung von einem mit (dem wir Aufruf „Standard“ verwendet), wollte De Forest einen „wahrscheinlichen Fehler“ (mittlere Abweichung) von einem. Dies war die Form, die in der Tabelle "Wahrscheinlichkeit von Fehlern" in den Anhängen von "Ein Handbuch der sphärischen und praktischen Astronomie, Band II" von angegeben istZN(0,12)William Chauvenet . Aus dieser Tabelle interpolierte De Forest die Quantile einer Halbnormalverteilung von bis p = 0,995 , die er als "Fehler gleicher Häufigkeit" ansah .p=0,005p=0,995

De Forest-Fehlertabelle gleicher Häufigkeit

Wenn Sie die Normalverteilung nach De Forest simulieren möchten, können Sie diese Tabelle ausdrucken und ausschneiden. De Forest (1876) schrieb, dass die Fehler "auf 100 gleichgroße Pappbits geschrieben wurden, die in einer Schachtel geschüttelt und alle einzeln herausgezogen wurden".

Der Astronom und Meteorologe Sir George Howard Darwin (Sohn des Naturforschers Charles) veränderte die Dinge, indem er ein sogenanntes "Roulette" zur Erzeugung zufälliger normaler Abweichungen entwickelte. Darwin (1877) beschreibt, wie:

x720π0xe-x2dx+-+-

"Index" sollte hier als "Zeiger" oder "Indikator" gelesen werden (vgl. "Zeigefinger"). Stigler weist darauf hin, dass Darwin wie De Forest eine halbnormale kumulative Verteilung um die Festplatte verwendete. Die anschließende Verwendung einer Münze zum Anbringen eines zufälligen Zeichens führt zu einer vollständigen Normalverteilung. Stigler merkt an, dass es unklar ist, wie fein die Skala abgestuft wurde, geht aber davon aus, dass die Anweisung zum manuellen Anhalten der Scheibe während der Drehung "die potentielle Vorspannung in Richtung eines Abschnitts der Scheibe zu verringern und den Vorgang zu beschleunigen" lautete.

Sir Francis Galton , übrigens ein Halbcousin von Charles Darwin, wurde bereits im Zusammenhang mit seinem Quincunx erwähnt. Während dies mechanisch eine Binomialverteilung simuliert, die nach dem De Moivre-Laplace-Theorem eine auffallende Ähnlichkeit mit der Normalverteilung aufweist (und gelegentlich als Lehrmittel für dieses Thema verwendet wird), hat Galton tatsächlich ein weitaus aufwändigeres Schema erstellt, als er es wünschte Probe aus einer Normalverteilung. Noch außergewöhnlicher als die unkonventionellen Beispiele oben in dieser Antwort, entwickelte Galton normalverteilte Würfel- oder genauer gesagt, ein Würfelsatz, der eine ausgezeichnete diskrete Annäherung an eine Normalverteilung mit einer mittleren Abweichung von eins ergibt. Diese Würfel aus dem Jahr 1890 werden in der Galton Collection am University College London aufbewahrt.

Galton normale Würfel

In einem Artikel von 1890 in Nature schrieb Galton:

Als Instrument zur zufälligen Auswahl habe ich nichts Besseres als Würfel gefunden. Es ist am mühsamsten, die Karten zwischen den aufeinanderfolgenden Zügen gründlich zu mischen, und die Methode, markierte Bälle in einem Beutel zu mischen und aufzurühren, ist noch mühsamer. Ein Teetotum oder irgendeine Form von Roulette ist diesen vorzuziehen, aber Würfel sind besser als alle. Wenn sie geschüttelt und in einen Korb geworfen werden, rasen sie so unterschiedlich gegeneinander und gegen die Rippen der Korbarbeit, dass sie wild umherstolpern, und ihre Positionen am Anfang geben keinen erkennbaren Hinweis darauf, was sie nach einer Weile sein werden einmal gut schütteln und werfen. Die Chancen, die ein Würfel bietet, sind vielfältiger als allgemein angenommen. Es gibt 24 gleiche Möglichkeiten und nicht nur 6, da jede Fläche vier Kanten hat, die verwendet werden können, wie ich zeigen werde.

Für Galton war es wichtig, in der Lage zu sein, schnell eine Folge normaler Abweichungen zu erzeugen . Nach jedem Wurf richtete Galton die Würfel nur durch Berühren aus und notierte die Punkte an den Vorderkanten. Er würfelte anfangs mehrere Würfel vom Typ I, deren Kanten halbnormale Abweichungen aufwiesen, ähnlich wie bei De Forest, aber mit 24 nicht 100 Quantilen. Für die größten Abweichungen (die auf den Würfeln vom Typ I als Leerzeichen markiert waren) würfelte er so viele der empfindlicheren Würfel vom Typ II (die nur bei einer feineren Abstufung große Abweichungen aufwiesen), wie er brauchte, um die Felder in seiner Sequenz auszufüllen . Um von halbnormalen zu normalen Abweichungen zu konvertieren, würfelt er mit Würfel III, der zuweist.+-114

Galton normales Würfeldesign

Das Labor für mathematisch-statistische Experimente von Raazesh Sainudiin umfasst ein Studentenprojekt der Universität von Canterbury, Neuseeland, das Galtons Würfel reproduziert . Das Projekt umfasst empirische Untersuchungen, bei denen die Würfel viele Male gewürfelt wurden (einschließlich einer empirischen CDF, die beruhigend "normal" aussieht), und eine Anpassung der Würfelbewertungen, damit sie der Standardnormalverteilung folgen. Unter Verwendung von Galtons Originalwerten gibt es auch ein Diagramm der diskretisierten Normalverteilung, der die Würfelwerte tatsächlich folgen.

Galton würfelt diskrete Verteilung


Wenn Sie im großen Stil bereit sind, das "Mechanische" zum Elektrischen zu dehnen, beachten Sie, dass RANDs Epos " A Million Random Digits" mit 100.000 normalen Abweichungen auf einer Art elektronischer Simulation eines Rouletterads basierte. Aus dem technischen Bericht (von George W. Brown, ursprünglich Juni 1949) entnehmen wir:

So motiviert entwickelten die RAND-Leute mit Unterstützung des Ingenieurs der Douglas Aircraft Company ein Elektro-Roulette-Rad, das auf einer Variation eines Vorschlags von Cecil Hastings basierte. Für die Zwecke dieses Vortrags ist eine kurze Beschreibung ausreichend. Eine Zufallsfrequenz-Impulsquelle wurde etwa einmal pro Sekunde durch einen Impuls mit konstanter Frequenz gesteuert, was durchschnittlich etwa 100.000 Impulse in einer Sekunde ergab. Impulsstandardisierungsschaltungen leiteten die Impulse an einen fünfstelligen Binärzähler weiter, so dass die Maschine im Prinzip wie ein Rouletterad mit 32 Positionen aussieht und durchschnittlich etwa 3000 Umdrehungen pro Umdrehung ausführt. Eine Binär-zu-Dezimal-Umwandlung wurde verwendet, wobei 12 der 32 Positionen weggeworfen wurden, und die resultierende Zufallsziffer wurde in einen IBM-Durchschlag eingespeist, wodurch Lochkartentabellen mit Zufallsziffern erhalten wurden.

χ2Tests der Häufigkeiten von geraden und ungeraden Ziffern ergaben, dass einige Chargen eine leichte Unwucht aufwiesen. Dies war in einigen Chargen schlimmer als in anderen, was darauf hindeutet, dass "die Maschine in dem Monat seit ihrer Einstellung ausgefallen ist ... Die Angaben zu dieser Maschine erfordern eine übermäßige Wartung, um sie in einwandfreiem Zustand zu halten". Es wurde jedoch ein statistischer Weg zur Lösung dieser Probleme gefunden:

Zu diesem Zeitpunkt hatten wir unsere ursprünglichen Millionen Ziffern, 20.000 IBM-Karten mit 50 Ziffern auf einer Karte, wobei die statistische Analyse eine geringe, aber wahrnehmbare ungerade Tendenz aufzeigte. Es wurde nun beschlossen, den Tisch neu zu ordnen oder zumindest zu verändern, indem mit ihm ein kleines Roulette gespielt wurde, um die ungerade-gerade-Tendenz zu beseitigen. Wir haben (Mod 10) die Ziffern jeder Karte ziffernweise zu den entsprechenden Ziffern der vorherigen Karte hinzugefügt. Die abgeleitete Tabelle mit einer Million Stellen wurde dann den verschiedenen Standardtests, Frequenztests, Serientests, Pokertests usw. unterzogen. Diese Millionen Stellen weisen eine saubere Gesundheitsbewertung auf und wurden als moderne Tabelle mit zufälligen Stellen von RAND übernommen.

Es gab natürlich guten Grund zu der Annahme, dass das Hinzufügen etwas Gutes bewirken würde. Im Allgemeinen ist der zugrunde liegende Mechanismus der begrenzende Ansatz von Summen von Zufallsvariablen, die das Einheitsintervall in der Rechteckverteilung auf die gleiche Weise modulieren, wie sich uneingeschränkte Summen von Zufallsvariablen der Normalität nähern. Diese Methode wurde von Horton und Smith von der Interstate Commerce Commission verwendet, um einige gute Chargen von scheinbar zufälligen Zahlen aus größeren Chargen von schlecht nicht zufälligen Zahlen zu erhalten.

[0,1]u[0,1]FF-1(u)

Ziggurat für halbnormal

Verweise


()"Wer arithmetische Methoden zur Erzeugung von Zufallszahlen in Betracht zieht, befindet sich natürlich in einem Zustand der Sünde. Denn es gibt, wie bereits mehrfach erwähnt, keine Zufallszahl - es gibt nur Methoden zur Erzeugung von Zufallszahlen und ein strenges Rechenverfahren ist natürlich keine solche Methode. "

Silberfisch
quelle
5
Sie müssen ein Scherz sein ... (aber, ja, das ist wirklich manuell +1)
Nalzok
8
Wie viel Arbeit steckt dahinter!
Richard Hardy
3
Beachten Sie, dass Würfel, die nicht dem Standard entsprechen, unfair sein können. Daher sollten Sie sie zuerst testen, z. B. youtube.com/watch?v=VI3N4Qg-JZM
Tim
2
@RichardHardy Paradoxerweise ist es für mich tatsächlich einfacher, all diese Dinge aufzuschreiben, während sie vor mir liegen und zumindest vorübergehend in meinem Gedächtnis bleiben, als zu versuchen, sie in meinem Kopf herumzutragen!
Silverfish
2
Auf jeden Fall finde ich es beeindruckend!
Richard Hardy
44

X=-2LogU1cos(2πU2)
Y.=-2LogU1Sünde(2πU2)X

Zum Beispiel kann ich auf meinem Linux-Betriebssystem überprüfen

$ date +%s.%N
1479733744.077762986
$ date +%s.%N
1479733980.615056616

U1=.077762986, U2=.615056616
X
> sqrt(-2*log(.077762986))*cos(2*pi*.615056616)
[1] -1.694815

Nachtrag: Da die Berechnung von Logarithmen und Cosinus als nicht manuell genug angesehen werden kann, gibt es für Box-Müller eine Variante, die die Verwendung dieser transzendentalen Funktionen vermeidet (siehe Aufgabe 2.9 in unserem Buch Monte Carlo Statistical Methods ):

eine Variante von Box-Muller

Nun kann man wegen der Exponentialvariablen gegen diese Version argumentieren. Es gibt jedoch auch eine sehr clevere Möglichkeit, diese Variablen ohne Aufruf transzendenter Funktionen zu simulieren , wie in diesem Algorithmus aus Luc Devroyes uneinheitlichen Zufallsvarianten zusammengefasst :

Bildbeschreibung hier eingeben

Zugegeben, es ist die Berechnung von 1 / e erforderlich, jedoch nur einmal.

(0,1)2(0,1)(0,2π)

Die Verwendung der CLT zur Annäherung an die Normalität ist sicherlich keine Methode, die ich jemals empfehlen würde, da (1) Sie immer noch andere Variablen benötigen, um den Durchschnitt zu ermitteln. Daher können Sie auch Uniformen im Box-Müller-Algorithmus verwenden, und (2) die Genauigkeit wächst erheblich langsam mit der Anzahl der Simulationen. Vor allem, wenn Sie eine diskrete Zufallsvariable wie das Ergebnis eines Würfels verwenden, auch bei mehr als sechs Gesichtern . Zitat von Thomas et al. (2007), eine Umfrage zum Für und Wider von Gaußschen Zufallsgeneratoren:

Der zentrale Grenzwertsatz ist natürlich ein Beispiel für eine "Näherungsmethode" - selbst wenn eine perfekte Arithmetik verwendet wird, ist die Ausgabe für endliches K nicht Gaußsch.

Hier ist ein kurzer Versuch, um das Problem zu veranschaulichen: Ich habe den 100-fachen Durchschnitt von 30 Würfelergebnissen generiert:

dies=apply(matrix(sample(1:6,30*100,rep=TRUE),ncol=30),1,mean)

dann normalisierte man diese Mittelwerte zu einem Mittelwert von Null - Varianz, die man variiert

stdies=(dies-3.5)/sqrt(35/12/30)

und betrachtete den normalen Sitz (oder das Fehlen davon) dieser Probe:

qqnorm (stdies, col = "gold2", pch = 19); abline (a = 0, b = 1, col = "steelblue", lwd = 2, lty = 2)

dies76/30122/30Dich

U=ich=1kDich-16ich
k=15
dies=matrix(apply(matrix(sample(0:5,15*200,rep=TRUE),nrow=15)/6^(1:15),2,sum),ncol=2) 
norma=sqrt(-2*log(dies[,1]))*c(cos(2*pi*dies[,2]),sin(2*pi*dies[,2]))

Die Passform ist so gut, wie es für eine normale Stichprobe der Größe 200 zu erwarten ist (zeichnen Sie eine weitere Stichprobe für eine echte normale Stichprobe auf norma=rnorm(100)):

Bildbeschreibung hier eingeben

wie ein Kolmogorov-Smirnov-Test weiter zeigt:

> ks.test(norma,pnorm)

        One-sample Kolmogorov-Smirnov test

data:  norma
D = 0.06439, p-value = 0.3783
alternative hypothesis: two-sided
Xi'an
quelle
3

Dies ist nicht genau zufällig, aber es sollte nahe genug sein, wie Sie scheinen ein grobes Experiment zu wollen.

Verwenden Sie Ihr Telefon, um einen Chronometer einzurichten. Nach gut 10 Sekunden hören Sie auf (je länger Sie warten, desto eher nähern Sie sich einem wirklich "zufälligen" Ergebnis, aber 10 Sekunden sind in Ordnung). Nehmen Sie die letzten Ziffern (z. B. 10,67 Sekunden ergeben 67). Wenden Sie die Perzentiltabelle für die Normalverteilung an. In diesem Beispiel müssen Sie nur nach 0,67 suchen, und Sie finden die Nummer. In diesem Fall beträgt Ihr Wert ungefähr 0,45. Dies ist nicht ganz genau, aber es gibt Ihnen eine solide Schätzung.

Wenn Sie unter 50 rollen, tun Sie einfach 100- [Ihr Ergebnis] und verwenden Sie die Tabelle. Ihr Ergebnis ist aufgrund der Symmetrie von N (0,1) dasselbe mit einem Minuszeichen.

Thalantas
quelle
3

n+1-1nnnN(0,1)


Warum? Lassen

Xk: ={+1 ob kMünzwurf ist Köpfe-1 ob k-te Münzwurf ist Schwänze

P(Xk=±1)=12

E(Xk)=0Var(Xk)=1

Y.: =X1+X2++Xn

E(Y.)=0Var(Y.)=n

Normalisieren,

Z: =Y.n

Wir erhalten eine Zufallsvariable mit Einheitsvarianz

E(Z)=0Var(Z)=1
Rodrigo de Azevedo
quelle
1
Wie bereits in meiner Antwort erwähnt, befürchte ich, dass die CLT die Zufälligkeit der Münzwürfe nicht effizient nutzt. Sie sollten besser ausgenutzt werden, wenn sich die Binärziffern eines pseudozufälligen U (0,1) mit einfacher oder doppelter Genauigkeit ändern.
Xi'an
@ Xi'an Ich habe Ihre Antwort gelesen, bevor ich meine poste, und Ihr Einwand gegen das CLT schien auf einer langsamen Konvergenz zu beruhen. Da dies ein Gedankenexperiment ist, kostet es nichts, 10 Milliarden Mal eine Münze zu werfen. Und es ist wirklich ein manuelles Verfahren, das keine Computer, keine Berechnung von Logarithmen, Quadratwurzeln oder Cosinus erfordert. Sicher, man könnte einen Rechenschieber verwenden, aber das könnte zu weit gehen.
Rodrigo de Azevedo
1
:}}: 10 Milliarden Münzwürfe klingen für mich nicht sehr manuell ...!
Xi'an
Manuell = von Hand. Das manuelle Berechnen von Logarithmen und Cosinus benötigt ebenfalls Zeit.
Rodrigo de Azevedo
0

Es ist erwähnenswert, dass Sie, sobald Sie eine Uniform (0,1) generieren können, eine beliebige Zufallsvariable generieren können, für die die inverse cdf berechnet werden kann, indem Sie einfach die einheitliche Zufallsvariable in die inverse CDF einfügen.

Wie kann man also eine Uniform (0,1) manuell berechnen? Nun, wie von @Silverfish erwähnt, gibt es eine Vielzahl von Würfeln, die von traditionellen RPG-Spielern verwendet werden. Eines davon ist ein zehnseitiger Würfel. Unter der Annahme, dass dies ein fairer Würfel ist, können wir jetzt eine diskrete Uniform erzeugen (0, 9).

0,01,0,02,...,0,99,1.000,001,0,002,...,0,999,1.000

Wir können also einer durchgehenden Uniform (0,1) sehr nahe kommen, indem wir sie mit einer fein gitterten diskreten Uniformverteilung mit einigen 10-seitigen Würfeln approximieren. Dieses kann dann in eine inverse CDF eingesteckt werden, um die interessierende Zufallsvariable zu erzeugen.

Cliff AB
quelle