আমি দুটি পর্যায়ক্রমিক টাইম-সিরিজের মধ্যে পর্বের পার্থক্যটি কীভাবে অনুমান করতে পারি?

Question

আমার 2 টি দৈনিক সময়-সিরিজ রয়েছে, প্রতি 6 বছর দীর্ঘ। কোলাহল করার সময় এগুলি উভয়ই স্পষ্টভাবে পর্যায়ক্রমিক (a 1 বছরের ফ্রিকোয়েন্সি সহ), তবে পর্যায়ক্রমের বাইরে রয়েছে বলে মনে হয়। আমি এই সময়-সিরিজের মধ্যে পর্বের পার্থক্যটি অনুমান করতে চাই।

আমি প্রতিটি টাইম-সিরিজে ফর্মের ফিটিং বক্ররেখাকে এবং খ এর জন্য দুটি পৃথক মানের তুলনা করছি, তবে আমি সন্দেহ করি যে এখানে আরও মার্জিত ( এবং কঠোর!) পদ্ধতিগুলি করার জন্য (সম্ভবত ফুরিয়ার ট্রান্সফর্ম ব্যবহার করছেন?)। আমি যদি সম্ভব হয় তবে আমার পর্বের পার্থক্যের প্রাক্কলনটিতে অনিশ্চয়তা সম্পর্কে এক ধরণের ধারণা রাখতেও পছন্দ করব।$a\sin(\frac{2\pi}{365}t - b)$

আপডেট :

ছায়াযুক্ত অঞ্চলগুলি 95% সিআই হয়।

দুটি সময়-সিরিজের মধ্যে নমুনা ক্রসক্রোলেশন:

Comments

একটি প্লট আকর্ষণীয় এবং সম্ভাব্য, সহায়ক হবে। দুটি সিরিজ সাইনোসয়েডালের কতটা কাছাকাছি?

— কার্ডিনাল

হাই কার্ডিনাল। আমার একটি প্লট আছে, তবে এটি আপলোড করার জন্য 10 খ্যাতি পয়েন্ট দরকার, দুর্ভাগ্যক্রমে! যা হওয়ার সাথে সাথে আমি তা করবো। সময়-সিরিজটি তাপমাত্রা এবং উদ্ভিদের বৃদ্ধির সাথে সম্পর্কিত, সুতরাং এটি বেশ স্বচ্ছ মরসুমী আচরণ অনুসরণ করুন, যদিও এটি শোরগোল। প্লট না দেখে, কীভাবে এই সমস্যাটির কাছে যাবেন সে সম্পর্কে আপনার কোনও চিন্তা আছে?

— পল কেটিং

হ্যাঁ. আমার বেশ কয়েকটি ধারণা রয়েছে তবে আপনি কীভাবে কাজ করছেন সে সম্পর্কে আরও ভাল ধারণা পেতে আমি প্রথমে একটি প্লট দেখতে চাই। যদি আপনি লিগটি সরবরাহ করার জন্য পোস্টটি ইগুর করতে এবং সম্পাদনা করতে আপলোড করেন তবে আমি ইমেজটিকে ইনলাইন রাখতে আবারও এটি সম্পাদনা করতে পারি। প্রতিনিধি শীঘ্রই যথেষ্ট আসবে। সাইটে স্বাগতম।

— কার্ডিনাল

কোনও কারণে, সম্ভবত আমার বর্তমানে পঙ্গু মেশিন, আমি কোনও লিঙ্কে ক্লিক করতে বা প্লট দেখতে পাচ্ছি না।

— jboman

খুব প্রথম তাত্পর্যপূর্ণ ও নোংরা চেক হিসাবে, আপনি কি দুটি সিরিজের মধ্যে নমুনা ক্রসকোরেলেশন প্লট করার এবং শিখরটি আবিষ্কার করার চেষ্টা করেছেন? কিছু কৌতূহল বিচ্ছিন্নতা রয়েছে, উদাহরণস্বরূপ, প্রথম সিরিজে কিছুটা ফেব্রুয়ারী 2005 বা এর কাছাকাছি সময়ে।

— কার্ডিনাল

Answer 1

ক্রস বর্ণালী বিশ্লেষণের জন্য এটি খুব সমস্যা problem এর পরে আপনার কাছে ভোক্তার মূল্য (তফাতের মধ্যে) এবং তেলের দাম ব্যবহার করে কোডের একটি সহনশীলতা (মোটামুটিভাবে, ফ্রিকোয়েন্সি ব্যান্ড দ্বারা ভাঙা একটি স্কোয়ার পারস্পরিক সম্পর্ক সহগ) এবং ফেজ (রেডিয়ানে পিছিয়ে, আবার ফ্রিকোয়েন্সি ব্যান্ডের মাধ্যমে) কোডের উদাহরণ রয়েছে।

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")

এগুলি শেষ নির্দেশাবলী দ্বারা উত্পাদিত গ্রাফগুলি s আপনি সম্ভবত এটি আপনার সেটআপের সাথে মানিয়ে নিতে পারেন।