Zusammenfassung meiner Antwort. Ich mag die Markov-Kettenmodellierung, aber sie vermisst den "zeitlichen" Aspekt. Wenn Sie sich hingegen auf den zeitlichen Aspekt konzentrieren (z. B. durchschnittliche Zeit beiDies ist eine Zwischenstufe aus dem Fall, in dem Sie nur die Übergangswahrscheinlichkeit schätzen und aus dem Fall, in dem Sie nur die in einem bestimmten Zustand verbrachte Zeit messen. Ich hoffe das hilft.−1
(VDi)i≥1(Si)i≥1
Yt=Y+t−Y−t
Y+t=∑i=0∞1VDi≤t,Si=1 and Y−t=∑i=0∞1VDi≤t,Si=−1
ϵ
λϵt=limdt→01dtP(Yϵt+dt−Yϵt=1|Ft)
ϵ−+FtFt=σ(Y+t,Y−t,VD1,…,VDY+t+Y−t,S1,…,SY+t+Y−t)
Aber in Anlehnung an Ihre Frage nehmen Sie implizit an, dass
Dies bedeutet, dass für eine deterministische Folge existiert so dass .
P(Yϵt+dt−Yϵt=1|Ft)=P(Yϵt+dt−Yϵt=1|Yt)
ϵ=+,−(μϵi)i∈Zλϵt=μϵYt
Innerhalb dieses Formalismus kann Ihre Frage wie folgt umformuliert werden: "Es ist wahrscheinlich, dass " (oder zumindest ist der Unterschied größer als a gegebene Schwelle).μ+−1−μ+0>0
Unter dieser Annahme ist es einfach zu zeigen, dass ein [homogener Markov-Prozess] [3] auf wobei der Generator durch gegeben istYtZQ
∀i,j∈ZQi,i+1=μ+iQi,i−1=μ−iQii=1−(μ+i+μ−i)Qij=0 if |i−j|>1
Beantwortung der Frage (durch Vorschlagen eines Maximum-Likelihood-Schätzwerts für das statistische Problem) Nach
dieser Neuformulierung wird das Problem gelöst, indem geschätzt und ein Test auf dessen Werten erstellt wird. Lassen Sie uns den Index ohne Verlust der Allgemeinheit reparieren und vergessen . Die Schätzung von (und ) kann vor der Beobachtung von erfolgen(μ+i)iμ+μ−
(T1,η1),…,(Tp,ηp) wobei die Länge der der im Zustand verbrachten Perioden sind (dh aufeinanderfolgende Zeiten mit ) und ist wenn die Frage positiv bewertet wurde, wenn sie negativ bewertet wurde, und wenn sie der letzte Beobachtungszustand war.TjjthpiYt=iηj+1−10
Wenn Sie den Fall mit dem letzten Beobachtungsstand vergessen, stammen die erwähnten Paare aus einer Verteilung, die von und abhängt : Sie wird verteilt als (wobei Exp eine Zufallsvariable aus einer Exponentialverteilung ist und + oder -1 ist, je nachdem, wer das Maximum realisiert). Dann können Sie das folgende einfache Lemma verwenden (der Beweis ist einfach):μ+iμ−i(min(Exp(μ+i),Exp(μ−i)),η)η
Lemma Wenn und dann ist und . X+⇝Exp(μ+)X−⇝Exp(μ−)T=min(X+,X−)⇝Exp(μ++μ−)P(X+1<X−)=μ+μ++μ−
Dies impliziert, dass die Dichte von gegeben ist durch:
wobei für die Dichtefunktion einer exponentiellen Zufallsvariablen ist mit Parameter . Aus diesem Ausdruck lässt sich leicht der Maximum-Likelihood-Schätzer von und :f(t,ϵ)(T,η)
f(t,ϵ)=gμ++μ−(1(ϵ=+1)∗μ++1(ϵ=−1)∗μ−μ++μ−)
gaa>0aμ+μ−
(μ^+,μ^−)=argminln(μ−+μ+)((μ−+μ+)∑i=1pTi+p)−p−ln(μ−)−p+ln(μ+)
wobeiund.
p + = | i : δ i = + 1 |p−=|i:δi=−1|p+=|i:δi=+1|
Kommentare für fortgeschrittenere Ansätze
Wenn Sie Fälle berücksichtigen möchten, bei denen der letzte beobachtete Zustand ist (sicherlich klüger, weil es bei oft Ihre letzte Punktzahl ist ...), müssen Sie die Neuzuordnung ein wenig ändern. Die entsprechende Zensur ist relativ klassisch ...- 1i−1
Mögliche andere Ansätze können die Möglichkeit von
- Eine Intensität haben, die mit der Zeit abnimmt
- Eine Intensität, die mit der Zeit abnimmt, die seit der letzten Abstimmung vergangen ist (ich bevorzuge diese. In diesem Fall gibt es klassische Methoden, um zu modellieren, wie die Dichte abnimmt ...
- Sie können annehmen, dass eine glatte Funktion von iμ+ii
- .... Sie können andere Ideen vorschlagen!