मेरे पास एक समय श्रृंखला है जिसमें दोहरे मौसमी घटक शामिल हैं और मैं निम्नलिखित समय श्रृंखला घटकों (प्रवृत्ति, मौसमी घटक 1, मौसमी घटक 2 और अनियमित घटक) में श्रृंखला को विघटित करना चाहूंगा। जहाँ तक मुझे पता है, R में एक श्रृंखला को विघटित करने की STL प्रक्रिया केवल एक मौसमी घटक की अनुमति देती है, इसलिए मैंने श्रृंखला को दो बार विघटित करने की कोशिश की है। पहले, निम्न कोड का उपयोग करके पहले मौसमी घटक होने की आवृत्ति सेट करके:
ser = ts(data, freq=48)
dec_1 = stl(ser, s.window="per")
फिर, मैंने विघटित श्रृंखला के अनियमित घटक ( dec_1) को दूसरे मौसमी घटक होने की आवृत्ति सेट करके, जैसे:
ser2 = ts(dec_1$time.series[,3], freq=336)
dec_2 = stl(ser2, s.window="per")
मैं इस दृष्टिकोण के साथ बहुत आश्वस्त नहीं हूं। और मैं जानना चाहूंगा कि क्या किसी सीरीज़ को विघटित करने के कोई अन्य तरीके हैं, जिनमें कई मौसमी हैं। इसके अलावा, मैंने देखा है कि tbats()आर फोरकास्ट पैकेज में फंक्शन एक मॉडल को कई मौसमी श्रृंखलाओं के साथ फिट करने की अनुमति देता है, हालांकि, यह नहीं कहता कि इसके साथ श्रृंखला को कैसे विघटित किया जाए।