, จำลองสถานการณ์ในช่วงการพยากรณ์


18

ฉันมีข้อมูลอนุกรมเวลาและฉันใช้เป็นโมเดลเพื่อให้พอดีกับข้อมูล เอ็กซ์ทีเป็นตัวบ่งชี้ตัวแปรสุ่มที่เป็นทั้ง 0 (เมื่อฉันไม่เห็นเหตุการณ์ที่ยาก) หรือ 1 (เมื่อฉันเห็นเหตุการณ์ที่หายาก) จากการสังเกตก่อนหน้านี้ที่ฉันมีสำหรับX tฉันสามารถพัฒนาแบบจำลองสำหรับX tโดยใช้วิธีการแบบ Variable Length Markov Chain สิ่งนี้ทำให้ฉันสามารถจำลองX tตลอดช่วงเวลาการพยากรณ์และให้ลำดับของศูนย์และอัน เนื่องจากนี่เป็นเหตุการณ์ที่หายากฉันจะไม่เห็นARIMA(p,d,q)+XtXtXtXtXtXt=1 บ่อยครั้ง ฉันสามารถคาดการณ์และได้รับการคาดการณ์ช่วงเวลาที่อยู่บนพื้นฐานของค่าจำลองสำหรับที Xt

คำถาม:

ฉันจะพัฒนาขั้นตอนการจำลองที่มีประสิทธิภาพที่จะคำนึงถึงการเกิดขึ้นของ 1 ในจำลองในช่วงคาดการณ์หรือไม่ ฉันต้องได้รับค่าเฉลี่ยและช่วงการพยากรณ์ Xt

ความน่าจะเป็นของการสังเกต 1 นั้นน้อยเกินไปสำหรับฉันที่จะคิดว่าการจำลองแบบมอนติคาร์โลปกติจะทำงานได้ดีในกรณีนี้ บางทีฉันสามารถใช้“ การสุ่มตัวอย่างที่สำคัญ” แต่ฉันไม่แน่ใจอย่างแน่นอน

ขอขอบคุณ.


6
พวกโปรดอย่าเปลี่ยนชื่อและเนื้อหาของคำถามของฉันมากเกินไป! "การผสม" และ "ลูกโซ่มาร์คอฟที่มีความยาวผันแปร" ไม่ใช่คำถามของฉัน คำถามเกี่ยวกับการพยากรณ์และการจำลอง กรุณาแจ้งให้เราตัดสินใจว่าจะถามคำถาม ...
สถิติ

อะไรคือความสำคัญขององค์ประกอบ Arima ในคำถามของคุณ ดูเหมือนว่ามันจะไม่เกี่ยวข้องกับคำถามเลย?
mpiktas

คิดอีกถ้าน่าจะเป็นของอยู่ในระดับต่ำมากเมื่อเทียบกับX T = 0ช่วงคาดการณ์ของ[ 0 , 0 ]จะมีความคุ้มครองความน่าจะเป็น1 -พี ดังนั้นช่วงการทำนายอาจจะไม่เป็นประโยชน์ในกรณีของคุณ? นอกจากนี้หากd > 0สำหรับโมเดลA R I M A ( p , d , q ) ของคุณดังนั้นA RP(Xt=1)=pXt=0[0,0]1pd>0ARIMA(p,d,q)ARIMA(p,d,q)ส่วนจะครองที Xt
mpiktas

@mpiktas: ขอบคุณสำหรับความคิดเห็น Arima มีความสำคัญในคำถามของฉันเนื่องจากนี่เป็นรุ่นหลักที่ฉันใช้ให้พอดี คุณหมายถึงอะไรโดย“ ช่วงเวลาการคาดคะเนของ [0,0]”? ฉันคิดว่าช่วงเวลาการพยากรณ์มีประโยชน์แม้ในกรณีนี้ ฉันมีแต่ผลของX t ที่มีต่อค่าที่ติดตั้งA R I M A ( p , d , q )นั้นโดดเด่น แม้ในช่วงเวลาที่คาดการณ์ไว้X tมีผลของมันเอง d>0XtARIMA(p,d,q)Xt
สถิติ

คำตอบ:


0

ประการแรกเราพิจารณากรณีทั่วไปมากขึ้น Let ที่~ ( )และX ~ X ( ) จากนั้นสมมติว่าการสนับสนุนของg x ( )ครอบงำหนึ่งในf X ( )และอินทิกรัลด้านล่างทั้งหมดที่มีอยู่เรามี: 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 your case

fX(x)={px=11px=0
and gX() can be defined like this:
gX(x)={0.5x=10.5x=0
Therefore, you can simulate X via distribution gX(), but all the observations with X=1 will have the weight p0.5=2p and all the observations with X=0 will have the weight 1p0.5=2(1p). Simulation of the ARIMA process will not be affected.
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.