मैं दो आवधिक समय-श्रृंखलाओं के बीच चरण अंतर का अनुमान कैसे लगा सकता हूं?

Question

मेरे पास 2 दैनिक समय-श्रृंखलाएं हैं, प्रत्येक 6 साल लंबी हैं। शोर करते हुए, वे दोनों स्पष्ट रूप से आवधिक हैं (~ 1 वर्ष की आवृत्ति के साथ), लेकिन चरण से बाहर दिखाई देते हैं। मैं इन समय-श्रृंखलाओं के बीच के चरण अंतर का अनुमान लगाना चाहूंगा।

मैंने प्रत्येक टाइम-सीरीज़ के लिए फॉर्म के फिटिंग कर्व्स पर विचार किया है और बस b के लिए दो अलग-अलग मूल्यों की तुलना कर रहा हूं, लेकिन मुझे संदेह है कि अधिक सुरुचिपूर्ण हैं ( और यह करने के लिए कठोर!) तरीके (शायद फूरियर रूपांतरण का उपयोग?)। यदि संभव हो तो मैं अपने चरण अंतर अनुमान में अनिश्चितता के बारे में किसी तरह का विचार रखना पसंद करूंगा।$a\sin(\frac{2\pi}{365}t - b)$

अपडेट :

छायांकित क्षेत्र 95% सीआई हैं।

दो समय-श्रृंखला के बीच नमूना पारगमन:

Comments

एक भूखंड दिलचस्प और संभवतः, सहायक होगा। साइनसॉइडल के कितने करीब हैं दो श्रृंखला?

— कार्डिनल

हाय कार्डिनल। मेरे पास एक भूखंड है, लेकिन दुर्भाग्य से इसे अपलोड करने के लिए 10 प्रतिष्ठा बिंदुओं की आवश्यकता है! ऐसा होते ही मैं करूंगा। समय-श्रृंखला तापमान और वनस्पति विकास से संबंधित है, इसलिए बहुत लगातार मौसमी व्यवहार का पालन करें, हालांकि यह शोर है। भूखंडों को देखे बिना, क्या आपके पास इस समस्या से संपर्क करने के बारे में कोई विचार है?

— पॉल कीटिंग

हाँ। मेरे पास कुछ विचार हैं, लेकिन मैं सबसे पहले एक कथानक देखना चाहता हूं जो आप के साथ काम कर रहे हैं। यदि आप लिंक देने के लिए अपना प्लॉट imgur और पोस्ट को संपादित करते हैं, तो मैं छवि इनलाइन डालने के लिए इसे फिर से संपादित कर सकता हूं। जल्द ही रिपीट आएगा। साइट पर आपका स्वागत है।

— कार्डिनल

किसी कारण से, शायद मेरी वर्तमान में अपंग मशीन, मैं एक लिंक पर क्लिक नहीं कर सकता या एक प्लॉट नहीं देख सकता।

— 1

बहुत पहले त्वरित और गंदे चेक के रूप में, क्या आपने दो श्रृंखलाओं के बीच नमूना क्रॉसकोरेल्टेशन की साजिश रचने और चोटी खोजने की कोशिश की है? कुछ उत्सुक असंतोष हैं, उदाहरण के लिए, पहली श्रृंखला में कुछ समय फरवरी 2005 या उसके आसपास।

— कार्डिनल

Answer 1

यह बहुत ही समस्या है क्रॉस-स्पेक्ट्रल विश्लेषण के लिए अच्छा है। आगे आपके पास उपभोक्ता मूल्यों (तेल के अंतर) और तेल की कीमत का उपयोग करते हुए कोड का एक उदाहरण है, और सुसंगतता का अनुमान है (मोटे तौर पर, आवृत्ति बैंड द्वारा टूटा हुआ एक सहसंबंध गुणांक) और चरण (रेडियन में अंतराल, फिर से आवृत्ति बैंड द्वारा)।

Crudo <- dget(file="Crudo.dge")
IPC <- dget(file="ipc2001.dge")[,1]
dIPC <- diff(IPC)
datos <- ts.union(dIPC,
           Crudo)
datos <- window(datos,
           start=c(1979,1),
           end=c(2002,1))
sp <- spectrum(datos,
           main="Petróleo e IPC",
           spans=rep(3,5))
par(mfrow=c(2,1))
plot(sp,plot.type="coh")
plot(sp,plot.type="phase")

ये अंतिम निर्देशों द्वारा निर्मित रेखांकन हैं। आप शायद इसे अपने सेटअप में ढाल सकते हैं।