Finite-State-Maschine mit Gamma-verteilten Wartezeiten

7

Ich habe eine Zustandsmaschine mit positiven und negativen Eingängen. Die Zeit zwischen positiven Eingaben folgt einer Gammaverteilung (X+Γ(k+,θ+)) und die Zeit zwischen negativen Eingaben folgt einer anderen Gammaverteilung (XΓ(k,θ)). Daher die EmpfangswahrscheinlichkeitenK positive und negative Eingaben über einen bestimmten Zeitraum sind für alle genau bekannt K. Die Zustandsmaschine ist unten dargestellt:

Zustandsmaschine mit vier Zuständen.

Die blauen Kästchen repräsentieren erreichbare Zustände in der Maschine, und die durchgezogenen und gestrichelten Linien repräsentieren einen positiven bzw. negativen Eingang. Wenn sich die Maschine beispielsweise im Zustand 3 befindet und ein positiver Eingang eintrifft, liefert die Maschine einen positiven Ausgang und wird auf den Zustand 2 zurückgesetzt. Wenn die Maschine dann einen negativen Eingang empfängt, wechselt sie in den Zustand 1, ohne einen Ausgang zu erzeugen.

Ist es möglich, den PMF für positive Ausgänge zu finden? Das heißt, wie hoch ist die Wahrscheinlichkeit zu bekommenK positive Ausgänge über das gleiche Zeitintervall für alle K?

Wälder
quelle
4
Das Problem bei der Präsentation einer Spielzeugversion Ihres eigentlichen Problems besteht darin, dass Sie es möglicherweise nicht so abstrahieren, wie es für Sie am hilfreichsten ist. (Nicht-Statistiker sind oft überrascht, welche Details statistisch wichtig sind.) Beschreiben Sie stattdessen das eigentliche Problem .
Kodiologe
Vielen Dank. Ich weiß nicht viel über diese Domain, aber es sieht so aus, als könnten Sie das System wahrscheinlich als Markov-Kette darstellen.
Kodiologe
Ich habe mir Markov-Ketten mit kontinuierlicher Zeit angesehen, und es sieht nicht so aus, als ob meine Zustandsmaschine die Markov-Eigenschaft erfüllt, da die in jedem Zustand verbrachte Zeit keiner Exponentialverteilung folgt. Bin ich tot im Wasser oder gibt es noch etwas, das ich versuchen könnte?
Mwoods
1
(+1) Sind alle Wartezeiten unabhängig voneinander, dh bilden die positiven und negativen Eingänge zwei unabhängige Erneuerungsprozesse? (Dies ist der Fall, wenn zur Zeit ein negativer Eingang auftritttIst die Zeit bis zum nächsten negativen Eingang unabhängig von allen vorherigen Vorkommen positiver und negativer Eingänge? und der Staat?) Dies ist höchstwahrscheinlich beabsichtigt, da keine Abhängigkeit angegeben ist (muss aber angenommen werden, da die Beschreibung dies nicht garantiert)
Juho Kokkala
@JuhoKokkala, die Antwort ist ja, sie sind unabhängig voneinander. Wenn ein negativer Eingang zum Zeitpunkt t eintrifft, hängt die Zeit bis zum nächsten negativen Eingang nicht von vorherigen Eingängen oder dem aktuellen Zustand ab.
Mwoods

Antworten:

2

Beachten Sie, dass dies KEIN Versuch ist, das Problem vollständig zu beantworten, sondern zu zeigen, wie das Fehlen der Markov-Eigenschaft für einen Sonderfall behoben werden kann, der möglicherweise nicht zutrifft - einer, der viel zu lang ist, um Kommentare abzugeben.

Wie Sie festgestellt haben, handelt es sich leider nicht um einen Markov-Prozess, sondern um einen Semi-Markov-Prozess. Wenn Sie a) eine ganze Zahl habenk+ und kund b) bereit sind, Ihren Zustandsraum zu erweitern, können Sie dies in einen Markov-Prozess umwandeln, indem Sie die Gamma-Verteilungen zu Erlang-Verteilungen machen und die Erlang-Variablen eine Summe von iid-Exponentialvariablen mit demselben Skalierungsparameter wie der ursprüngliche Erlang sind variieren.

Wir können den Zustandsraum um zwei neue Variablen erweitern, den "+ Zustand" und den "- Zustand", die aufzeichnen, "wie weit" wir bei der Erzeugung der nächsten positiven oder negativen Ankunft sind. Nehmen Sie zur Konkretheit ank+=5;; Die nächste positive Ankunft erfolgt, wenn die fünfte von fünf aufeinanderfolgenden exponentiellen Ankünften aufgetreten ist. Der "+ -Zustand" zeichnet also auf, wie viele positive Ankünfte seit der letzten positiven Eingabe aufgetreten sind. Die Reihenfolge der "+ state" -Werte ist{0,1,2,3,4,0,1,...};; Zustand0 kann nur zu wechseln 0 oder 1, Zustand 4 kann nur zu wechseln 4 oder 0, und so weiter.

Ihr Zustandsraum wird [BoxID,+,] Aufzeichnen, in welcher Box sich der Prozess befindet, wie viele positive Ankünfte im Modul aufgetreten sind k+und wie viele negative Ankünfte sind im Modul aufgetreten k.

Wir haben jetzt zwei zufällige Variablen - die Zeit bis zum nächsten Übergang "+ Zustand" und die Zeit bis zum nächsten Übergang "Zustand" - die beide exponentiell verteilt sind. Da das Minimum von zwei unabhängigen Exponentialvariablen selbst Exponential ist, ist die Zeit bis zum nächsten Übergang (eines beliebigen Typs) Exponential mit einer Rate, die der Summe der beiden Komponentenraten entspricht (θ++θ oder 1/θ++1/θabhängig davon, wie Ihre Gamma-Verteilungen parametrisiert sind). Die Wahrscheinlichkeit, dass der nächste Übergang ein "+ Zustand" -Übergang ist, ist gerechtθ+/(θ++θ), oder 1/θ+/(1/θ++1/θ), wiederum abhängig davon, wie Ihre Gamma-Verteilungen parametrisiert sind. Da die Zeit bis zum nächsten Übergang jetzt exponentiell verteilt ist, verfügen Sie über eine CTMC (Continuous Time Markov Chain), die auf standardmäßige Weise analysiert werden kann.

Für ein konkretes Beispiel wird angenommen, dass positive Ankünfte mit einer Rate von auftreten 0.5/Zeiteinheit und negative Ankünfte treten mit der Rate von auf 0.25/Zeiteinheit. Die Zeit bis zum nächsten Übergang ist exponentiell mit einer Rate von0.75/Zeiteinheit und die Wahrscheinlichkeit, dass der Übergang durch eine positive Ankunft ausgelöst wird, ist 0.5/(0.5+0.25)=2/3.

Jetzt haben Sie einen stark vergrößerten Zustandsraum, wobei jedes Feld in Ihrem ursprünglichen Diagramm vorhanden ist k+k gibt das Innere an, aber zumindest haben Sie die Markov-Eigenschaft und können die stationären Wahrscheinlichkeiten finden, in Box 3 zu sein, und die Zustände mit "+ state" = k+1Dies ist einer der Zustände, aus denen Sie einen Übergang erleben können, der zu einer positiven Ausgabe führt. Wenn Sie diese stationären Wahrscheinlichkeiten mit der Übergangsmatrix und der mittleren Zeit zwischen den Übergängen kombinieren, erhalten Sie die langfristige Durchschnittsrate für eine positive Ausgabe. Sie können die gewünschte Wahrscheinlichkeitsverteilung auch anhand der stationären Wahrscheinlichkeiten, der Übergangsmatrix und der Tatsache berechnen, dass die Zeit zwischen Übergängen eine Exponentialverteilung mit bekannter Rate aufweist.

Jbowman
quelle
Genau das brauchte ich. Ich habe sehr großk+ und k(in der Größenordnung von zehntausend), so dass die Übergangsmatrix ziemlich groß sein wird. Die Rettung ist, dass die Matrix auch sehr spärlich ist. Vielen Dank!
Mwoods
Nun, bei so großen Formparametern ist die Unterscheidung zwischen ganzzahligen und nicht ganzzahligen Werten umstritten! Ja, Sparsamkeit wird in der Tat Ihre rettende Gnade sein ...
Jbowman