Ich versuche, das Innenleben des Hamiltonian Monte Carlo (HMC) zu verstehen, kann aber den Teil nicht vollständig verstehen, wenn wir die deterministische Zeitintegration durch einen Vorschlag von Metropolis-Hasting ersetzen. Ich lese das großartige Einführungspapier A Conceptual Introduction to Hamiltonian Monte Carlo von Michael Betancourt, daher werde ich der darin verwendeten …
Angenommen, wir möchten einige Erwartungen berechnen: EYEX|Y[f(X,Y)]EYEX|Y[f(X,Y)]E_YE_{X|Y}[f(X,Y)] Angenommen, wir möchten dies mithilfe der Monte-Carlo-Simulation approximieren. EYEX|Y[f(X,Y)]≈1RS∑r=1R∑s=1Sf(xr,s,yr)EYEX|Y[f(X,Y)]≈1RS∑r=1R∑s=1Sf(xr,s,yr)E_YE_{X|Y}[f(X,Y)] \approx \frac1{RS}\sum_{r=1}^R\sum_{s=1}^Sf(x^{r,s},y^r) ABER nehmen wir an, es ist teuer, Proben aus beiden Verteilungen zu ziehen, so dass wir es uns nur leisten können, eine feste Zahl . KKK Wie sollen wir zuordnen ? Beispiele …
Ich hatte nie einen formalen Statistikkurs, aber aufgrund meiner Forschungsrichtung stoße ich ständig auf Artikel, die verschiedene statistische Konzepte anwenden. Oft sehe ich eine Beschreibung eines Monte-Carlo- Prozesses, der auf eine bestimmte Situation angewendet wird, und für das, was ich 9 von 10 Mal erfassen kann, kommt es auf eine …
Markov-Kette Monte Carlo ist eine auf Markov-Ketten basierende Methode, mit der wir Proben (in einer Monte-Carlo-Einstellung) aus nicht standardmäßigen Verteilungen erhalten können, aus denen wir keine Proben direkt ziehen können. Meine Frage ist, warum die Markov-Kette für die Monte-Carlo-Probenahme "auf dem neuesten Stand" ist. Eine alternative Frage könnte sein, ob …
Ich bin mir ziemlich sicher, dass ich die Funktionsweise der Monte-Carlo-Integration verstehe, aber ich verstehe nicht die Formulierung, wie sie zur Schätzung von Pi verwendet wird. Ich gehe nach dem in der 5. Folie dieser Präsentation beschriebenen Verfahren vor: http://homepages.inf.ed.ac.uk/imurray2/teaching/09mlss/slides.pdf Ich verstehe die vorbereitenden Schritte. Pi entspricht der vierfachen Fläche …
Ich möchte die Verwendung der Monte-Carlo-Simulation in der chisq.test()Funktion in R verstehen . Ich habe eine qualitative Variable mit 128 Stufen / Klassen. Meine Stichprobengröße beträgt 26 (ich konnte nicht mehr "Einzelpersonen" befragen). Also werde ich natürlich einige Level mit 0 "Individuen" haben. Tatsache ist jedoch, dass ich nur eine …
Ich hoffe, dies ist der richtige Ort, um zu fragen, ob Sie es nicht in ein geeigneteres Forum verschieben können. Ich habe mich schon eine ganze Weile gefragt, wie man nicht quadratisch integrierbare Funktionen mit Monte Carlo Integration behandelt. Ich weiß, dass MC immer noch eine korrekte Schätzung liefert, aber …
Ich versuche, die Grenzwahrscheinlichkeit für ein statistisches Modell mit Monte-Carlo-Methoden zu berechnen: f( x ) = ∫f( x ∣ θ ) π( θ )dθf(x)=∫f(x∣θ)π(θ)dθf(x) = \int f(x\mid\theta) \pi(\theta)\, d\theta Die Wahrscheinlichkeit ist gut verhalten - glatt, logarithmisch konkav - aber hochdimensional. Ich habe versucht, wichtige Stichproben zu erstellen, aber die …
Ich habe versucht, mit Metropolis-Algorithmen in R aus einer bivariaten Dichte zu simulieren und hatte kein Glück. Die Dichte kann ausgedrückt werden als , wobei die Singh-Maddala-Verteilung istp ( y | x ) p ( x ) p ( x )p ( x , y)p(x,y)p(x,y)p ( y| x)p(x)p(y|x)p(x)p(y|x)p(x)p ( x …
Ich schätze derzeit ein stochastisches Volatilitätsmodell mit Markov-Ketten-Monte-Carlo-Methoden. Dabei implementiere ich Gibbs- und Metropolis-Stichprobenverfahren. Angenommen, ich nehme eher den Mittelwert der posterioren Verteilung als eine Zufallsstichprobe daraus. Wird dies allgemein als Rao-Blackwellization bezeichnet ? Insgesamt würde dies dazu führen, dass der Mittelwert über die Mittelwerte der posterioren Verteilungen als Parameterschätzung …
Ich nehme an einem Kurs über Monte-Carlo-Methoden teil und wir haben in der letzten Vorlesung die Rejection Sampling-Methode (oder Accept-Reject Sampling-Methode) gelernt. Es gibt viele Ressourcen im Web, die den Beweis dieser Methode zeigen, aber irgendwie bin ich nicht davon überzeugt. In der Ablehnungsabtastung haben wir also eine Verteilung der …
Nach dem, was ich gelesen habe, ist Hamiltonian Monte Carlo die "goto" MCMC-Methode, wenn Ihr Problem hochdimensional ist. Wie viele Dimensionen 10, 100, 1000, 10.000, 100.000, ... sind praktisch zu viele? Die Rechenkosten werden zweifellos zu einem Problem, und ich nehme an, dass das verwendete Modell wichtig ist, aber abgesehen …
Ich sehe mir an, wie sich der erwartete minimale euklidische Abstand zwischen zufällig einheitlichen Punkten und dem Ursprung ändert, wenn wir die Dichte zufälliger Punkte ( Punkte pro Quadrateinheit ) um den Ursprung erhöhen . Ich habe es geschafft, eine Beziehung zwischen den beiden als solche zu finden: Expected Min …
Ich lese gerade ein großartiges HMC-Einführungspapier von Prof. Michael Betancourt, aber ich verstehe nicht, wie wir die Verteilung des Impulses wählen sollen. Zusammenfassung Die Grundidee von HMC besteht darin, eine Impulsvariable in Verbindung mit der Zielvariablen einzuführen . Sie bilden gemeinsam einen Phasenraum .pppqqq Die Gesamtenergie eines konservativen Systems ist …
Kann einer unserer Monte-Carlo-Experten die "unerwartete" Erwartung am Ende dieser Antwort erklären ? Ex-post- Zusammenfassung der anderen Frage / Antwort: Wenn IID-Zufallsvariablen sind und die Erwartungen existieren, zeigt ein einfaches Symmetrieargument , dass , aber ein Monte-Carlo-Experiment mit scheint diesem Satz zu widersprechen.X1,…,XnX1,…,XnX_1,\dots,X_nE[Xi/X¯]E[Xi/X¯]\mathrm{E}[X_i/\bar{X}]E[Xi/X¯]=1E[Xi/X¯]=1\mathrm{E}[X_i/\bar{X}]=1Xi∼N(0,1)Xi∼N(0,1)X_i\sim\mathrm{N}(0,1) x <- matrix(rnorm(10^6), nrow = 10^5) mean(x[,2]/rowMeans(x)) …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.