R- तुच्छ कोड उदाहरण में कण फ़िल्टर


9

मैं एक साधारण कोड उदाहरण की तलाश कर रहा हूं कि आर में एक कण फ़िल्टर कैसे चलाया जाए। धूमधाम पैकेज राज्य स्पेस मैथ को थोड़ा सपोर्ट करता प्रतीत होता है, लेकिन उदाहरणों के लिए एक साधारण OO डेवलपर के लिए प्रोग्रामेटिक रूप से फॉलो करना थोड़ा मुश्किल होता है, विशेष रूप से अवलोकन किए गए डेटा को एक धूमधाम वस्तु में कैसे लोड किया जाए।

कहते हैं कि मेरे पास इनपुट के रूप में शोर डेटा के 1 कॉलम के साथ एक सीएसवी फ़ाइल है, और मैं इसे एक और फिल्टर करने के लिए एक कण फ़िल्टर के माध्यम से चलाना चाहता हूं, ताकि इसे दूसरी सीएसवी फ़ाइल में अनुमानों के आउटपुट के साथ इसे साफ किया जा सके।

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

उदाहरणों के साथ मुख्य कठिनाई सीएसवी डेटा को एक धूमधाम वस्तु में लोड करना है।

एक बहुत ही सरल स्टेट स्पेस मॉडल अभी के लिए पर्याप्त होना चाहिए।

आर-उत्सुक के लिए कोई विचार?


यह उन सभी के लिए उपयोगी हो सकता है जो जोड़े के व्यापार, या एल्गोरिथम के व्यापार को सामान्य रूप से देखते हैं, जहां ये शोर से घिरा एक आर्थिक संबंध है।

IMHO आप फ़िल्टर को स्वयं कोड करने से बेहतर हैं ...
डॉ। जी

बहुत बढ़िया! इसे देखने वाले सभी लोगों के नमूने / सुझाव / संकेत? एक वैकल्पिक समाधान बिना समाधान से बेहतर है।
14:11 पर user1180428

@ user1180428: मैंने अपना उत्तर संपादित कर लिया है, जो अब आपके लिए एक संभावित विकल्प प्रदान कर सकता है।
वेन

जवाबों:


1

संपादित करें: ऐसा लगता है कि अधिकांश कण फिल्टर पैकेज अब चले गए हैं। हालाँकि, मैं लाप्लासडेमोन (एक बायेसियन एमसीएमसी पैकेज) के साथ खेल रहा हूं और इसमें पीएमसी (जनसंख्या मोंटे कार्लो) फ़ंक्शन है जो पीएमसी को लागू करता है, जो एक प्रकार का कण फिल्टर है। शायद एक त्वरित कण फिल्टर तरह के लिए बहुत अधिक मशीनरी, लेकिन सीखने के लायक पैकेज।

आप CRAN पर पैकेज और ट्यूटोरियल पा सकते हैं ।

मूल: ईमानदार होने के लिए, सरलतम मामले में, pompउपयोग करना कठिन है। यह कुछ भी आप कर सकते हैं के लिए बहुत लचीला है, लेकिन यह किराने की दुकान पर जाने के लिए एक अंतरिक्ष जहाज का उपयोग करने जैसा है।

क्या आपने कलमन फ़िल्टर (यदि आपका डेटा कलमन फ़िल्टर की मान्यताओं को पूरा कर सकता है) को देखने की कोशिश की है, जिसमें आधार फ़ंक्शंस tsSmoothऔर StructTSकेवल (यूनीवेरिएट), और पैकेज शामिल हैं dlm? मैं भी loessऔर अन्य धूम्रपान करने वालों पर एक नज़र डालूँगा।

मुझे आशा है कि मैं गलत हूं और कोई व्यक्ति जल्दी से यहां आकर खड़ा हो जाता है, "यहां यह सरल सरल डेटा के लिए कैसे करना है जैसे कि आपके पास कुछ मामूली धारणाएं हैं।" मुझे पैकेज का उपयोग करने में सक्षम होना पसंद है।


1
वहाँ हैक किया गया। दुर्भाग्य से, एक सरल चलती औसत इस में कलमन फ़िल्टर की तुलना में एक प्रयोग करने योग्य संकेत को बेहतर ढंग से कील करने के लिए लगता है, और कई अन्य उदाहरण - कलमन: लिंक , एसएमए: लिंक डेटा स्थिर है, अब तक एक डिक्की फुलर पी <0.01 है। शायद मैं गलत कर रहा हूं। इस डेटा पर एक कण फ़िल्टर चलाने के लिए एक स्क्रिप्ट, और अन्य जोड़े व्यापारिक उम्मीदवार भयानक होंगे (मुझे लगता है)।
user1180428
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.