ตัวกรองอนุภาคในตัวอย่างโค้ดเล็ก ๆ น้อย ๆ

ฉันกำลังมองหาตัวอย่างรหัสง่ายๆของวิธีการเรียกใช้ตัวกรองอนุภาคในอาร์แพคเกจ pomp ปรากฏขึ้นเพื่อสนับสนุนบิตสเปซทางคณิตศาสตร์ของรัฐ แต่ตัวอย่างนั้นค่อนข้างยุ่งยากที่จะติดตามโดยทางโปรแกรม วิธีการโหลดข้อมูลที่สังเกตลงในวัตถุเอิกเกริก

• ตัวอย่างที่นี่: http://cran.r-project.org/web/packages/pomp/vignettes/intro_to_pomp.pdf

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

 y <- read.csv("C:/Dev/VeryCleverStatArb/inputData.csv", header=FALSE)
 #CSV to Pomp object ???
 #Run Particle Filter
 #Write estimates to csv.

ปัญหาหลักของตัวอย่างคือการโหลดข้อมูล csv ลงในวัตถุเอิกเกริก

แบบจำลองอวกาศของรัฐที่ง่ายมากน่าจะเพียงพอสำหรับตอนนี้

ความคิดใด ๆ สำหรับ R- อยากรู้อยากเห็น?

แก้ไข: ดูเหมือนว่าแพคเกจตัวกรองอนุภาคส่วนใหญ่หายไปแล้ว อย่างไรก็ตามฉันได้เล่นกับ LaplacesDemon (แพ็คเกจ Bayesian MCMC) และมันมีฟังก์ชั่น PMC (Population Monte Carlo) ซึ่งใช้ PMC ซึ่งเป็นตัวกรองอนุภาคชนิดหนึ่ง อาจจะเป็นเครื่องจักรมากเกินไปสำหรับสิ่งที่กรองอนุภาคอย่างรวดเร็ว แต่แพคเกจก็คุ้มค่าที่จะเรียนรู้

คุณสามารถค้นหาแพคเกจและแบบฝึกหัดที่CRAN

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

คุณได้ลองดูตัวกรองคาลมานแล้วหรือยัง (ถ้าข้อมูลของคุณอาจเป็นไปตามสมมติฐานของตัวกรองคาลมาน) รวมถึงฟังก์ชั่นพื้นฐานtsSmoothและStructTS(ไม่รวมตัวแปรเท่านั้น) และแพ็คเกจdlm? ฉันจะดูloessและเครื่องปั่นอื่น ๆ

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