Genaue Probenahme aus ungeeigneten Gemischen

10

Angenommen, ich möchte aus einer kontinuierlichen Verteilung . Wenn ich einen Ausdruck von in der Form habepp(x)p

p(x)=i=1aifi(x)

wobei und f_i Verteilungen sind, aus denen leicht abgetastet werden kann, dann kann ich leicht Abtastwerte aus p erzeugen durch:ai0,iai=1fip

  1. Abtasten eines Labels i mit der Wahrscheinlichkeit ai
  2. Abtastung Xfi

Ist es möglich, diese Prozedur zu verallgemeinern, wenn die ai gelegentlich negativ sind? Ich vermute, ich habe dies irgendwo gesehen - möglicherweise in einem Buch, möglicherweise für die Kolmogorov-Distribution -, daher würde ich gerne eine Referenz als Antwort akzeptieren.

Wenn ein konkretes Spielzeug Beispiel nützlich ist, sagen wir , ich auf Probe möchten

p(x,y)exp(xyαxy)x,y>0
dann werde ich Nehmen Sie α(0,2) aus technischen Gründen, die im großen Schema der Dinge nicht allzu wichtig sein sollten.

Im Prinzip könnte ich dies dann als folgende Summe erweitern:

p(x,y)n=0(1)nαn(n2)!(n2)!n!(xn/2ex(n2)!)(yn/2ey(n2)!).

Die (x,y) -Terms innerhalb der Summe können dann unabhängig voneinander als Gamma-Zufallsvariablen abgetastet werden. Mein Problem ist offensichtlich, dass die Koeffizienten "gelegentlich" negativ sind.

Bearbeiten 1 : Ich stelle klar, dass ich aus p exakte Stichproben generieren möchte , anstatt die Erwartungen unter p zu berechnen . Für Interessenten wird in den Kommentaren auf einige Verfahren hingewiesen.pp

Edit 2 : Ich habe die Referenz, die einen bestimmten Ansatz für dieses Problem enthält, in Devroyes 'Non-Uniform Random Variate Generation' gefunden . Der Algorithmus stammt aus 'A Note on Sampling from Combinations of Distributions' von Bignami und de Matteis . Das Verfahren besteht effektiv darin, die Dichte von oben durch die positiven Terme der Summe zu begrenzen und dann eine Ablehnungsabtastung basierend auf dieser Hüllkurve zu verwenden. Dies entspricht der in der Antwort von @ Xi'an beschriebenen Methode.

πr8
quelle
1
Warum können Sie nicht einfach den absoluten Wert von und dann Ihre Stichprobe negieren ? Mit anderen Worten definieren Sie(vorausgesetzt , es ist endlich) und dann renormieren Ihre Summe von . aiXfiZ:=i=1|ai|Z
Alex R.
2
@AlexR. Wenn ich Sie verstehe, wäre eine Version davon praktisch, um die Erwartungen unter zu berechnen , aber immer noch nicht, um genaue Stichproben aus . Dies ist sicherlich eine Antwort auf ein relevantes Problem, obwohl es nicht ganz das ist, wonach ich suche. pp
πr8
4
Dies hängt davon ab, was Sie mit diesem Beispiel tun möchten. Zum Beispiel zum Berechnen von Momenten erscheint es einfach, die Abtastung aus Dichtemischungen zu verallgemeinern, indem zusätzlich jeder Punkt, der aus einer Komponente mit negativem Koeffizienten ausgewählt wurde, als "negativer" Punkt markiert und sein Beitrag in der Momentschätzung negativ gewichtet wird. Ebenso könnten Sie eine KDE mit solchen negativen Gewichten konstruieren, vorausgesetzt, Sie können die Möglichkeit akzeptieren, dass einige ihrer Werte negativ sind! (cc @ Xi'an)
whuber
1
Was wäre eine "genaue" Stichprobe einer Verteilung? Ob und wie Sie eine Mischung mit negativen Gewichten nutzen können, hängt wiederum davon ab, wie Sie die Probe verwenden möchten.
whuber
1
Dies beantwortet Ihre Frage nicht, aber Sie könnten daran interessiert sein, Informationen
Tim

Antworten:

5

Ich habe über diese Frage gerätselt, aber nie eine zufriedenstellende Lösung gefunden.

Eine mögliche Eigenschaft ist, dass wenn eine Dichte schreibt, wobei ist Dichte, so dass , Simulation von und Zurückweisen dieser Simulationen mit der Wahrscheinlichkeit Simulationen von liefert . Im aktuellen Fall ist die normalisierte Version der positiven Gewichtskomponenten und ist der Rest

f(x)=g(x)ωh(x)1ωω>0
gg(x)ωh(x)gωh(x)/g(x)fg
g(x)=αi>0αifi(x)/αi>0αi
ωh
h(x)=αi<0αifi(x)/αi<0αi
Dies findet sich zwar in der Simulationsbibel von Devroye, Ungleichmäßige Erzeugung zufälliger Variablen , Abschnitt II.7.4, folgt jedoch aus einer einfachen Annahme-Ablehnungs-Argumentation.

Ein erster rechnerischer Nachteil dieses Ansatzes besteht darin, dass trotz der ersten Simulation aus einer ausgewählten Komponente die Summen in und für den Zurückweisungsschritt berechnet werden müssen. Wenn die Summen ohne geschlossene Form unendlich sind, kann die Methode "Akzeptieren / Ablehnen" nicht implementiert werden .figh

Eine zweite Schwierigkeit besteht darin, dass die Ablehnungsrate ist, da beide Gewichtssummen in derselben Reihenfolge sindhat keine Obergrenze. Wenn die mit verknüpfte Reihe nicht absolut konvergiert, ist die Akzeptanzwahrscheinlichkeit tatsächlich Null! Und die Methode kann in dieser Situation nicht implementiert werden.

αi>0αi=1αi<0αi
1ϱaccept=αi<0|αi|/i|αi|
αi

Wenn im Fall einer Mischungsdarstellung geschrieben werden kann als Die Komponente kann zuerst ausgewählt und dann die auf die Komponente angewendete Methode. Dies kann jedoch schwierig zu implementieren sein, da die Identifizierung von Paaren , die zu passen, aus der möglicherweise unendlichen Summe nicht unbedingt machbar ist.f

f(x)=i=1αigi(x)ωih(xi)1ωiωi>0
(gi,hi)gi(x)ωih(xi)>0

Ich denke, eine effizientere Auflösung könnte von der Seriendarstellung selbst kommen. Devroye, Uneinheitliche Erzeugung zufälliger Variablen , Abschnitt IV.5, enthält eine Vielzahl von Serienmethoden. Wie zum Beispiel der folgende Algorithmus für eine alternative Seriendarstellung des Ziels wenn ' s konvergieren mit gegen Null und ist eine Dichte:

f(x)=κh(x){1a1(x)+a2(x)}
ai(x)nhDevroyes alternative Serienmethode

Das Problem wurde kürzlich im Zusammenhang mit dem Debiasing von voreingenommenen Schätzern für MCMC betrachtet, wie beispielsweise beim Glynn-Rhee-Ansatz . Und der russische Roulette- Schätzer (mit einem Zusammenhang mit dem Bernoulli-Fabrikproblem). Und die unvoreingenommene MCMC-Methodik . Aber es gibt kein Entkommen vor dem Vorzeichenproblem ... Dies macht seine Verwendung bei der Schätzung von Dichten wie bei pseudo-marginalen Methoden schwierig.

Bei weiterem Denken, ist meine Schlussfolgerung , dass es keine generische Methode ist eine tatsächliche Simulation von dieser zu erzeugen Serie [anstatten Mischung , die eine falsche Bezeichnung erweist], ohne die Einführung weitere> Struktur mit den Elementen der Serie, wie die in der obige Algorithmus aus Devroyes Bibel . Da die meisten (?) Dichten eine Serienerweiterung der oben genannten Art ermöglichen, würde dies ansonsten die Existenz einer Art universeller Simulationsmaschine implizieren ...

Xi'an
quelle
Vielen Dank! Ich schätze auch die zusätzlichen Referenzen.
πr8
1
Zusätzlicher Dank für die sehr gründliche Antwort und Referenzen. Ich bin froh, diese Antwort zu akzeptieren, da es gelingt, exakte Stichproben aus in endlicher Zeit zu generieren . Ich werde wahrscheinlich bis zu einem gewissen Grad weiter über das Problem nachdenken. Die einzige zusätzliche Idee, die ich hatte und die vielversprechend erscheint, besteht darin, die Abtastung von als Abtastung , abhängig von , und dass es möglicherweise geometrische gibt Einsicht, die für diese Charakterisierung nützlich ist (ich denke wie ein Slice-Sampler auf ). Prost! pp=λgμhXgλgμh{(x,y):μh(x)<y<λg(x)}
8.
1
Ich habe den bedingten Sampler ziemlich schlecht erklärt; Die satzbasierte Charakterisierung ist (meiner Meinung nach) etwas klarer. Mein entscheidender Punkt ist, dass, wenn Sie gleichmäßig aus der zweidimensionalen Menge in der letzten Zeile abtasten können , die Koordinate die richtige Verteilung hat. Ob diese Charakterisierung für längere summenbasierte unsachgemäße Gemische nützlich sein kann, bleibt abzuwarten. (x,y)x
πr8
1
Ich habe auch an einen Slice-Sampler gedacht, aber dies ist im Sinne einer Simulation nicht "genau".
Xi'an
1

Ich habe den Entwurf einer Idee, die funktionieren könnte. Es ist nicht genau , aber hoffentlich asymptotisch genau. Um daraus eine wirklich strenge Methode zu machen, bei der die Annäherung kontrolliert wird oder etwas daran bewiesen werden kann, ist wahrscheinlich viel Arbeit erforderlich.

Erstens können Sie, wie von Xi'an erwähnt, die positiven Gewichte einerseits und die negativen Gewichte andererseits gruppieren, so dass das Problem schließlich nur zwei Verteilungen und :gh

p=λgμh

mit . Beachten Sie, dass Sie .λμ=1λ1

Meine Idee ist die folgende. Sie möchten Beispiel Beobachtungen von . Machen:Np

  • Probieren Sie Werte aus und speichern Sie sie in einer ListeλNg
  • Entfernen Sie für jeden von abgetasteten -Wert den nächsten (verbleibenden) Nachbarn aus der Liste.μNh

Am Ende erhalten Sie Punkte. Es muss nicht genau der nächste Nachbar sein, sondern nur ein Punkt, der "nah genug" ist. Der erste Schritt ist wie das Erzeugen von Materie. Der zweite Schritt ist wie die Erzeugung von Antimaterie und deren Kollision und Aufhebung mit Materie. Diese Methode ist nicht genau, aber ich glaube, unter bestimmten Bedingungen ist sie für großes asymptotisch genau (um sie für kleines fast genau zu machen, müssen Sie zuerst ein großes und dann einen kleinen zufälligen Teil der endgültigen Liste nehmen). . Ich gebe ein sehr informelles Argument, das eher eine Erklärung als ein Beweis ist.(λμ)N=NNnN

Betrachten Sie im Beobachtungsraum und ein kleines Volumen um mit Lebesgue-Volumen . Nach dem Abtasten von beträgt die Anzahl der Elemente in der Liste, die sich ebenfalls in befinden, ungefähr . Nach dem zweiten Schritt wird ungefähr daraus entfernt, und Sie haben ungefähr die gewünschte Anzahl . Dazu müssen Sie davon ausgehen, dass die Anzahl der Punkte im Volumen ausreichend groß ist.xvxϵgvλNg(x)ϵμNh(x)ϵNp(x)ϵ

Es ist sehr unwahrscheinlich, dass dieses Verfahren einer großen Dimension oder einigen Pathologien von und widersteht, es kann jedoch in einer kleinen Dimension und ausreichend glatten, "ausreichend gleichmäßigen" Verteilungen arbeiten.gh

Hinweis zu einer genauen Methode:

Ich habe dies zuerst für diskrete Verteilungen gedacht, und in diesem Fall ist die Methode eindeutig nicht genau, da sie Stichproben mit der Wahrscheinlichkeit 0 erzeugen kann. Ich habe die starke Intuition, dass eine genaue Methode mit endlicher Verarbeitungszeit nicht möglich ist und dass dies der Fall ist Unmöglichkeit konnte zumindest für diskrete Verteilungen nachgewiesen werden. Die Regel des Spiels ist, dass Sie nur exakte "Orakel" -Sampler für und aber und als Funktionen von . Der Einfachheit halber beschränken Sie sich auf Bernoulli-Verteilungen. Die Nichtexistenz einer exakten Methode hängt mit der Bernoulli-Factory- Theorie zusammen: Wenn Sie aus einem eine -Münze erstellen könntenghghx(λpμq)p-coin und eine -coin, dann könnten Sie eine -coin aus einer -coin erstellen, von der bekannt ist, dass sie für unmöglich ist .qλppλ>1

Benoit Sanchez
quelle
1
Ich dachte darüber nach, lehnte es jedoch ab, weil meine anfänglichen Bemühungen, zu demonstrieren, dass es funktionieren könnte, zu der Erkenntnis führten, dass es bestenfalls eine Annäherung und möglicherweise eine schlechte sein wird. Ja, asymptotisch könnte es funktionieren, aber es wird die Anforderung des OP nach "exakten" Stichproben aus der Verteilung nicht erfüllen.
whuber
Die Effizienz dieser Methode liegt genau in der gleichen Größenordnung wie die exakte Annahme-Zurückweisungs-Methode.
Xi'an
1
Einverstanden. Sie sind jedoch ganz anders. Die Accept-Reject-Methode muss und als Funktionen von berechnen . Ich konzentrierte mich darauf, nur Sampling von und als "Orakel" -Sampler zu verwenden, wie in einer echten Mischung. Je mehr ich darüber nachdenke, desto mehr bin ich davon überzeugt, dass eine exakte Methode, die auf der Probenahme von Orakeln basiert, nicht existieren kann. ghxgh
Benoit Sanchez
1
Ich denke , das ist grundsätzlich richtig, aber es können nützliche Klassen von Sonderfällen, in denen eine solche genaue Methode tut exist. Dies liegt daran, dass (1) in einigen Fällen die Berechnung von einfach ist und (2) Sie nicht sowohl als auch berechnen müssen - Sie müssen nur dieses Verhältnis berechnen. g/(g+h)gh
whuber
@ BenoitSanchez Vielen Dank für Ihre ausführliche Antwort; Ich schätze besonders die Kommentare am Ende über die (mögliche) Unmöglichkeit der Genauigkeit. Ich bin in der Vergangenheit auf Bernoulli-Fabriken gestoßen und fand sie ziemlich herausfordernd. Ich werde versuchen, das Thema erneut zu betrachten und zu prüfen, ob es Erkenntnisse liefert.
8.