একাধিক মৌসুমী উপাদানগুলির সাথে একটি সময় সিরিজ কীভাবে পচে যায়?


22

আমার একটি টাইম সিরিজ রয়েছে যাতে দ্বৈত মৌসুমী উপাদান থাকে এবং আমি সিরিজটি নিম্নলিখিত সময়ের সিরিজ উপাদানগুলিতে (প্রবণতা, মৌসুমী উপাদান 1, মৌসুমী উপাদান 2 এবং অনিয়মিত উপাদান) বিভক্ত করতে চাই। যতদূর আমি জানি, আরে সিরিজটি পচানোর জন্য এসটিএল পদ্ধতিটি কেবল একটি মৌসুমী উপাদানকে অনুমতি দেয়, তাই আমি সিরিজটি দু'বার পচানোর চেষ্টা করেছি। প্রথমত, নিম্নলিখিত কোডটি ব্যবহার করে প্রথম মৌসুমী উপাদান হিসাবে ফ্রিকোয়েন্সি সেট করে:

ser = ts(data, freq=48)
dec_1 = stl(ser, s.window="per")

তারপরে, আমি পচনশীল সিরিজের অনিয়মিত উপাদানটি ( dec_1) theতুকে দ্বিতীয় মৌসুমী উপাদান হিসাবে ফ্রিকোয়েন্সিটি সেট করে পচন করেছি:

ser2 = ts(dec_1$time.series[,3], freq=336)
dec_2 = stl(ser2, s.window="per")

আমি এই পদ্ধতির সাথে খুব আত্মবিশ্বাসী নই। এবং আমি জানতে চাই যে কোনও একাধিক মৌসুমী সিরিজটি পচে যাওয়ার অন্য কোনও উপায় আছে কিনা। এছাড়াও, আমি লক্ষ করেছি যে tbats()আর পূর্বাভাস প্যাকেজটিতে ফাংশনটি একাধিক মৌসুমীর সাথে একটি সিরিজে একটি মডেল ফিট করতে দেয় তবে এটি কীভাবে এটির সাথে কোনও সিরিজ পচিয়ে ফেলবে তা বলে না।


হাই এবং সাইটটিতে আপনাকে স্বাগতম। আপনার দুটি মৌসুমী উপাদানগুলির জন্য, তাদের কি বিভিন্ন সময়সীমার রয়েছে, যেমন একটি সাপ্তাহিক এবং অন্য মাসিক?
মিশেল

1
রব হ্যান্ডম্যান, কোহলার, অর্ড ও স্নাইডার "অধ্যায় 14 এর পূর্বাভাস" তাত্পর্যপূর্ণ স্মুথিংয়ের সাথে এটি অন্তর্ভুক্ত। হেন্ডম্যানের আর-তে একটি পূর্বাভাস প্যাকেজও রয়েছে I আমি মনে করি হিনডম্যান এই সাইটে এই সাইটে পোস্ট করেছিলেন তবে এটি সম্ভবত তার ব্লগে।
zbicyclist

উত্তরটি দেওয়ার জন্য @ মিশেল হাই ধন্যবাদ thanks হ্যাঁ দুটি মৌসুমী উপাদানগুলির বিভিন্ন পর্যায়ক্রমিকতা রয়েছে। প্রথমটির পর্যায়কালটি 48 (দৈনিক seasonতু), এবং দ্বিতীয়টির পর্যায়কাল 336 (সাপ্তাহিক seasonতু) থাকে। এটি একটি আধ ঘন্টা সময় সিরিজ।
এসেই

@ জাইসাইক্লিস্ট আমি বিশ্বাস করি যে আপনি যে পূর্বাভাস প্যাকেজ সম্পর্কে রয়েছেন তা হ'ল 'পূর্বাভাস' প্যাকেজ যা আমি মূল পোস্টে উল্লেখ করেছি। আমি এই প্যাকেজের টিবিটস ফাংশনটি দেখেছি তবে এটি কীভাবে পচে যাওয়ার জন্য এটি ব্যবহার করবেন তা বলে না। আমি আরও কোনও চিত্র খুঁজে পেতে পারি কিনা তা দেখার জন্য আমি বইটি পরীক্ষা করে দেখব।
এসেস

2
আমি যা ভাবছিলাম তা এখানে। এটি হেন্ডম্যানের ব্লগে ছিল। robjhyndman.com
কমপ্লেক্স-

উত্তর:


13

forecastbats()tbats()এক্স(টি)টি

দেখুন http://robjhyndman.com/papers/complex-seasonality/ সূত্র এবং Hyndman এট (2008) ETS মডেলের একটি উন্নততর বর্ণনার জন্য জন্য। বিএটিএস এবং টিবিএটিএস হ'ল ইটিএসের একটি বর্ধিতাংশ।

উদাহরণ স্বরূপ:

fit <- bats(myTimeseries)
fit$x

এই ক্ষেত্রে, প্রতিটি সারি xফুরিয়ার-জাতীয় সুরেলাতে থাকবে।

এছাড়াও আছে plot.tbats()এবং plot.bats()স্বয়ংক্রিয়ভাবে পচা ফাংশন এবং উপাদান দেখতে।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.