, Simulation über Prognosezeitraum

18

Ich habe Zeitreihendaten und ich habe ein als Modell verwendet, um die Daten . Das ist eine Indikator-Zufallsvariable, die entweder 0 (wenn ich kein seltenes Ereignis sehe) oder 1 (wenn ich das seltene Ereignis sehe) ist. Basierend auf früheren Beobachtungen, die ich für , kann ich ein Modell für Verwendung der Markov-Ketten-Methode mit variabler Länge entwickeln. Dadurch kann ich das über den Prognosezeitraum simulieren und eine Folge von Nullen und Einsen angeben. Da dies ein seltenes Ereignis ist, sehe ich nicht oft . Ich kann die Vorhersageintervalle basierend auf den simulierten Werten für vorhersagen und erhalten . ARIMA(p,d,q)+XtXtXtXtXtXt=1Xt

Frage:

Wie kann ich ein effizientes Simulationsverfahren entwickeln, um das Auftreten von im simulierten während des Prognosezeitraums zu berücksichtigen ? Ich muss den Mittelwert und die Prognoseintervalle ermitteln. Xt

Die Wahrscheinlichkeit, 1 zu beobachten, ist zu gering, als dass ich annehmen könnte, dass die reguläre Monte-Carlo-Simulation in diesem Fall gut funktioniert. Vielleicht kann ich "Wichtigkeitsproben" verwenden, aber ich bin nicht sicher, wie genau.

Vielen Dank.

Stat
quelle
6
Leute, bitte ändert den Titel und den Text meiner Frage nicht zu sehr! "Mischen" und "Markov-Kette variabler Länge" ist nicht meine Frage. Die Frage dreht sich um Vorhersage und Simulation. Bitte lassen Sie mich entscheiden, wie ich die Frage stellen soll ...
26.08.12
Welche Bedeutung hat die Arima-Komponente in Ihrer Frage? Es scheint, dass es überhaupt nicht mit der Frage zusammenhängt?
mpiktas
Ein anderer Gedanke ist, dass, wenn die Wahrscheinlichkeit von sehr niedrig ist, im Vergleich zu das Vorhersageintervall von die Abdeckungswahrscheinlichkeit . Vielleicht sind Vorhersageintervalle in Ihrem Fall nicht so nützlich? Wenn für Ihr Modell ist, dominiert der das . P(Xt=1)=pXt=0[0,0]1pd>0ARIMA(p,d,q)ARIMA(p,d,q)Xt
mpiktas
@mpiktas: danke für die kommentare. Arima ist in der Tat wichtig für meine Frage, da dies das Hauptmodell ist, zu dem ich gepasst habe. Was meinst du mit "Vorhersageintervall von [0,0]"? Ich denke, die Vorhersageintervalle sind auch in diesem Fall nützlich. Ich habe , aber die Auswirkung von auf die angepassten Werte ist auffällig. Auch im prognostizierten Zeitraum hat einen eigenen Effekt. X t A R I M A ( p , d , q ) X td>0XtARIMA(p,d,q)Xt
Stat

Antworten:

0

Zunächst betrachten wir einen allgemeineren Fall. Sei , wobei A f A ( ) und X f X ( ) . Dann, unter der Annahme , die Unterstützung von g x ( ) dominiert die eine von f X ( ) und alle die Integrale unter exist, haben wir: P ( Y y ) = E f A , fY=Y(A,X)AfA()XfX()gx()fX()

P(Yy)=EfA,fX[I(Yy)]=EfX[EfA[I(Yy)X]]=supp(fX)EfA[I(Yy)X=x]fX(x)dx=supp(fX)EfA[I(Yy)X=x]fX(x)gX(x)gX(x)dx=supp(gX)EfA[I(Yy)fX(X)gX(X)X=x]gX(x)dx=EgX[EfA[I(Yy)fX(X)gX(X)X]]=EfA,gX[I(Yy)fX(X)gX(X)]

In Ihrem Fall ist und g X ( ) kann wie folgt definiert werden: g X ( x ) = { 0,5 x = 1 0,5 x = 0 Deshalb sind Sie kann X über die Verteilung g X ( ) simulieren , aber alle Beobachtungen mit X = 1

fX(x)={px=11px=0
gX()
gX(x)={0.5x=10.5x=0
XgX()X=1wird das Gewicht und alle Beobachtungen mitX=0haben das Gewicht1-pp0.5=2pX=0. Die Simulation des ARIMA-Prozesses wird nicht beeinflusst.1p0.5=2(1p)
LeonM
quelle