การคำนวณโอกาสในการบันทึกสำหรับ MLE ที่กำหนด (เชนมาร์คอฟ)


9

ขณะนี้ฉันทำงานร่วมกับมาร์คอฟเชนและคำนวณค่าความน่าจะเป็นสูงสุดโดยใช้ความน่าจะเป็นในการเปลี่ยนแปลงตามที่แนะนำโดยหลาย ๆ แหล่ง (เช่นจำนวนช่วงการเปลี่ยนภาพจาก a ถึง b หารด้วยจำนวนการเปลี่ยนภาพโดยรวมจาก a ไปยังโหนดอื่น ๆ )

ตอนนี้ฉันต้องการคำนวณความน่าจะเป็นของ MLE


คุณได้คำนวณค่าประมาณความน่าจะเป็นสูงสุดของความน่าจะเป็นของการเปลี่ยนแปลงแล้วและตอนนี้คุณต้องการคำนวณความน่าจะเป็นในการบันทึกของอะไร
Nick

คำตอบ:


11

ปล่อย {Xผม}ผม=1T เป็นเส้นทางของห่วงโซ่มาร์คอฟและปล่อยให้ Pθ(X1,...,XT) เป็นโอกาสในการสังเกตเส้นทางเมื่อ θ คือค่าพารามิเตอร์ที่แท้จริง (หรือที่เรียกว่าฟังก์ชันความน่าจะเป็นสำหรับ θ) เราใช้ความหมายของความน่าจะเป็นแบบมีเงื่อนไข

Pθ(X1,...,XT)=Pθ(XT|XT-1,...,X1)Pθ(X1,...,XT-1)

ตั้งแต่นี้เป็นห่วงโซ่มาร์คอฟเรารู้ว่า Pθ(XT|XT-1,...,X1)=Pθ(XT|XT-1)ดังนั้นวิธีนี้จะทำให้ง่ายขึ้น

Pθ(X1,...,XT)=Pθ(XT|XT-1)Pθ(X1,...,XT-1)

ทีนี้ถ้าคุณทำซ้ำตรรกะเดียวกันนี้ T ครั้งคุณได้รับ

Pθ(X1,...,XT)=Πผม=1TPθ(Xผม|Xผม-1)

ที่ไหน X0จะถูกตีความว่าเป็นสถานะเริ่มต้นของกระบวนการ คำศัพท์ทางด้านขวาเป็นเพียงองค์ประกอบของเมทริกซ์การเปลี่ยนแปลง เนื่องจากเป็นโอกาสในการบันทึกที่คุณร้องขอคำตอบสุดท้ายคือ:

L(θ)=Σผม=1Tเข้าสู่ระบบ(Pθ(Xผม|Xผม-1))

นี่เป็นโอกาสของเชนมาร์คอฟเดียว - หากชุดข้อมูลของคุณมีเชนมาร์คอฟ (อิสระ) หลายห่วงโอกาสที่จะเป็นผลรวมของเงื่อนไขของแบบฟอร์มนี้


ว้าวขอบคุณมากสำหรับคำตอบ ในกรณีนี้Pθความน่าจะเป็น "การเปลี่ยนแปลง" นั้นมาจาก MLE ใช่ไหม
fsociety

@ph_singer คุณยินดีมาก Pθ(Xผม|Xผม-1) ความน่าจะเป็นที่จะย้ายจากรัฐ Xผม-1 ถึง Xผมได้รับค่าพารามิเตอร์ θ. หากคุณไม่ได้กำหนดโครงสร้างของเมทริกซ์การเปลี่ยนแปลง (ซึ่งเป็นสิ่งที่ดูเหมือน) แล้วθ เพียงแค่ระบุเวกเตอร์ของความน่าจะเป็นในการเปลี่ยนแปลง (และ MLEs เป็นเพียงสัดส่วนตัวอย่างตามที่คุณระบุไว้อย่างถูกต้องในคำแถลงคำถามของคุณ) ดังนั้นใช่: Pθ^MLE(Xผม|Xผม-1) จะเป็นสัดส่วนตัวอย่างของการเคลื่อนไหวจากรัฐ Xผม-1 ที่สิ้นสุดในสถานะ Xผม.
แมโคร

ขอบคุณอีกครั้ง! อีกหนึ่งคำถาม: ถ้าฉันใช้คำสั่งอื่น (เช่น k = 2) กระบวนการนี้จะทำงานอย่างไร
fsociety

คุณช่วยอธิบายสิ่งที่คุณหมายถึงโดย "สั่งซื้อ" ได้ไหม?
มาโคร

(+1) OP น่าจะหมายถึง k=2เพื่อแสดงว่าMC อันดับที่สองคือขึ้นอยู่กับสองสถานะก่อนหน้าXผม-1,Xผม-2 มากกว่าเพียงแค่ล่าสุด Xผม-1.
พระคาร์ดินัล
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.