আমি মনে করি আপনি এটি সঠিকভাবে পেয়েছেন তবে ভিএআর মডেলটি তৈরি করার সময় আমি সাধারণত নিশ্চিত হয়েছি যে আমি এই পদক্ষেপগুলি অনুসরণ করেছি:
1. ভেরিয়েবল নির্বাচন করুন
এটি আপনার মডেলটি তৈরির সবচেয়ে গুরুত্বপূর্ণ অংশ। আপনি যদি কোনও সম্পদের দামের পূর্বাভাস দিতে চান তবে আপনাকে দামের গঠনের পদ্ধতির সাথে সম্পর্কিত ভেরিয়েবলগুলি অন্তর্ভুক্ত করতে হবে। এটি করার সর্বোত্তম উপায়টি একটি তাত্ত্বিক মডেলের মাধ্যমে। যেহেতু আপনি সম্পদটি কী এবং আপনি আপনার মডেলটিতে অন্তর্ভুক্ত অন্যান্য পরিবর্তনগুলি কী তা উল্লেখ করেননি আমি এই আইটেমটি সম্পর্কে সত্যই বেশি কিছু বলতে পারি না তবে আপনি এখানে সম্পদ মূল্যের মডেলগুলির একটি সংক্ষিপ্তসার জানতে পারেন ।
2. ডেটা পরীক্ষা করুন এবং সঠিক সামঞ্জস্য করুন
একবার আপনি ভেরিয়েবলগুলি নির্বাচন করেন, আপনি ডেটাতে কিছু সামঞ্জস্য করতে পারেন যা মডেলটির অনুমান এবং ব্যাখ্যা উন্নত করবে। সংক্ষিপ্ত পরিসংখ্যান ব্যবহার করা এবং বিদেশী, নিখোঁজ ডেটা এবং অন্যান্য অদ্ভুত আচরণগুলি সনাক্ত করতে সিরিজের একটি প্লট দেখতে দরকারী। দামের ডেটা নিয়ে কাজ করার সময়, লোকেরা সাধারণত প্রাকৃতিক লগগুলি গ্রহণ করে, যা একটি বৈকল্পিক-স্থিতিশীল রূপান্তর এবং এটির একটি ভাল ব্যাখ্যাও রয়েছে (লগের দামের পার্থক্যটি ক্রমাগত যৌগিক রিটার্নে পরিণত হয়)। মডেলটি অনুমান করার আগে আপনি লগগুলি গ্রহণ করেছেন কিনা তা সম্পর্কে আমি নিশ্চিত নই, তবে আপনি যদি সম্পদের দাম নিয়ে কাজ করছেন তবে তা করা ভাল ধারণা।
৩. ডেটাতে স্থিতিশীল উপাদান নেই কিনা তা পরীক্ষা করুন
আপনার সিরিজটি স্থির কিনা তা পরীক্ষা করতে এখন আপনি ইউনিট রুট পরীক্ষাগুলি ব্যবহার করতে পারেন। আপনি যদি কেবলমাত্র পূর্বাভাস দিতে আগ্রহী হন, যেমন @ জ্যাকোবিএইচ-র দ্বারা উল্লিখিত হয়েছে, আপনার সিরিজটি স্থির নয় এমনকী আপনি স্তরগুলিতে ভিএআর চালাতে পারেন, তবে তারপরেও আপনার স্ট্যান্ডার্ড ত্রুটিগুলি বিশ্বাস করা যায় না, যার অর্থ আপনি এর মান সম্পর্কে অনুমান করতে পারবেন না সহগ। আপনি এডিএফ পরীক্ষাটি ব্যবহার করে স্টেশনারি পরীক্ষা করেছেন, যা এই অ্যাপ্লিকেশনগুলিতে খুব বেশি ব্যবহৃত হয়, তবে নোট করুন যে আপনি যদি পরীক্ষা দিয়ে চালাতে চান তবে নির্দিষ্ট করা উচিত i) কোন ধ্রুবক এবং প্রবণতা নেই; ii) একটি ধ্রুবক এবং কোন প্রবণতা; এবং iii) একটি ধ্রুবক এবং একটি প্রবণতা। সাধারণত মূল্য সিরিজের স্টোকাস্টিক প্রবণতা থাকে, সুতরাং লিনিয়ার প্রবণতাটি সঠিক হবে না। এই ক্ষেত্রে আপনি স্পেসিফিকেশন চয়ন করতে পারেন ii। আপনার কোড আপনি ব্যবহার করেছেনndiffs
পূর্বাভাস প্যাকেজ ফাংশন। আমি নিশ্চিত নই যে এই তিনটি বিকল্পের মধ্যে কোনটি এই ফাংশনটি প্রয়োগ করে তারতম্যের সংখ্যা গণনা করার জন্য (আমি এটি ডকুমেন্টেশনে খুঁজে পাইনি)। আপনার ফলাফল যাচাই করতে আপনি ur.df
"urca" প্যাকেজে ফাংশনটি ব্যবহার করতে পারেন :
adf <- ur.df(x[, "VAR1"], type = "drift", lags = 10, selectlags = "AIC")
লক্ষ্য করুন এই কমান্ড আপনি ফলাফল ব্যাখ্যা শুধু তাকান সমস্যা থেকে থাকে 10. সর্বোচ্চ ল্যাগ সঙ্গে একটি ধ্রুবক এবং এআইসি কমান্ড দ্বারা নির্বাচিত lags সঙ্গে ADF পরীক্ষা চালানো হবে, এই প্রশ্নের । সিরিজটি যদি আমি (1) হয় তবে কেবলমাত্র পার্থক্যটি ব্যবহার করুন, যা ক্রমাগত চক্রবৃদ্ধি ফেরতের সমান হবে equal যদি পরীক্ষাটি নির্দেশ করে যে সিরিজটি আমি (2) এবং আপনি সন্দেহের মধ্যে থাকেন যে আপনি অন্যান্য পরীক্ষাও ব্যবহার করতে পারেন, যেমন ফিলিপস-পেরোন পরীক্ষা (PP.test
আর এ ফাংশন)। যদি সমস্ত পরীক্ষাগুলি নিশ্চিত করে যে আপনার সিরিজটি আমি (২) (পরীক্ষা চালানোর আগে সিরিজের লগটি ব্যবহার করতে ভুলবেন না) তবে দ্বিতীয় পার্থক্যটি গ্রহণ করুন, তবে নোট করুন যে ফলাফলগুলির আপনার ব্যাখ্যাটি পরিবর্তিত হবে, যেহেতু এখন থেকে আপনি কাজ করছেন ধারাবাহিকভাবে যৌগিক রিটার্নের পার্থক্য। সম্পদের দামগুলি সাধারণত আমি (1) হয় যেহেতু তারা এলোমেলো হাঁটার কাছাকাছি হয়, যা প্রথম পার্থক্যটি প্রয়োগ করার সময় একটি সাদা শব্দ।
4. মডেল ক্রম নির্বাচন করুন
এটি আকাইকে, শোয়ার্জ (বিআইসি) এবং হান্নান-কুইন এর মতো সাধারণভাবে ব্যবহৃত মানদণ্ডগুলির সাহায্যে করা যেতে পারে। আপনি VARselect
ফাংশনটি দিয়ে এটি করেছেন এবং এটি সঠিক, তবে মনে রাখবেন যে আপনি নিজের সিদ্ধান্ত নেওয়ার জন্য যে মানদণ্ডটি ব্যবহার করেছিলেন। সাধারণত বিভিন্ন মানদণ্ড ভিএআর জন্য বিভিন্ন আদেশ নির্দেশ করে।
৫) কন্টিগ্রেটিং সম্পর্ক রয়েছে কিনা তা পরীক্ষা করুন
যদি আপনার সমস্ত সিরিজ আমি (1) বা আমি (2) হয় তবে ভিএআর মডেল চালানোর আগে সাধারণত সিরিজের মধ্যে কোনও সমন্বয় সম্পর্ক নেই কিনা তা পরীক্ষা করা ভাল ধারণা, বিশেষত যদি আপনি এর সাথে আবেগ প্রতিক্রিয়া বিশ্লেষণ করতে চান অবশিষ্টাংশ। আপনি জোহেনসেন পরীক্ষা বা এনগেল-গ্রেঞ্জার ব্যবহার করে এটি করতে পারেন (কেবলমাত্র বাইভারিয়েট মডেলের জন্য)। আর-তে আপনি ca.jo
"urca" প্যাকেজটির কার্যকারিতা সহ জোহেনসেন পরীক্ষা চালাতে পারেন । দ্রষ্টব্য যে এই পরীক্ষারও আলাদা আলাদা স্পেসিফিকেশন রয়েছে। দামের সিরিজের জন্য আমি সাধারণত নিম্নলিখিত কোডটি ব্যবহার করি ( p
আইটেম 4 এর ল্যাগ দৈর্ঘ্যটি, স্তরটিতে সিরিজের সাথে সম্পাদিত হয়):
jo_eigen <- ca.jo(x, type = "eigen", ecdet = "const", K = p)
jo_trace <- ca.jo(x, type = "trace", ecdet = "const", K = p)
The. মডেলটি অনুমান করুন
যদি আপনার সিরিজটি সমন্বিত না হয়, আপনি VAR
আপনার কোড অনুসারে কমান্ডের সাহায্যে মডেলটি সহজেই অনুমান করতে পারেন । যদি সিরিজটি সমন্বিত হয় তবে নিম্নলিখিত কোডের সাথে ভেক্টর ত্রুটি সংশোধন মডেলটি (যেখানে k
সংযুক্তির ক্রমটি রয়েছে) অনুমান করে আপনাকে দীর্ঘকালীন সম্পর্কের বিষয়টি বিবেচনা করতে হবে :
vecm <- cajorls(joeigen, r = k)
Diagn. ডায়াগনস্টিক পরীক্ষা চালান
আপনার মডেলটি সুনির্দিষ্টভাবে নির্দিষ্ট করা আছে কিনা তা পরীক্ষা করতে আপনি অবশিষ্টাংশগুলিতে সিরিয়াল পারস্পরিক সম্পর্কের একটি পরীক্ষা চালাতে পারেন। আপনার কোডে আপনি ফাংশনটির সাথে একটি পোর্টম্যান্টিউ পরীক্ষা ব্যবহার করেছেন serial.test
। আমি এই ফাংশনটি কখনও ব্যবহার করি নি তবে আমি মনে করি এটি ঠিক আছে। এছাড়াও প্যাকেজ এমটিএসে প্রয়োগ করা লজং-বক্স পরীক্ষার একটি মাল্টিভিয়ারেট সংস্করণ রয়েছে যা আপনি ফাংশনটি দিয়ে চালাতে পারেন mq
।
8. ভবিষ্যদ্বাণী করুন
আপনার মডেলটি সুনির্দিষ্টভাবে নির্দিষ্ট হওয়ার পরে আপনি আপনার কোডটিতে যেমন predict
ফাংশনটি ব্যবহার করতে পারেন তা নিশ্চিত হওয়ার পরে । এমনকি ভেরিয়েবলগুলি কীভাবে কোনও বিশেষ ধাক্কায় irf
ফাংশনটি ব্যবহার করে তা প্রতিক্রিয়া জানায় তা পরীক্ষা করতে আপনি ইমপ্লাস প্রতিক্রিয়া ফাংশন প্লট করতে পারেন ।
9. পূর্বাভাস মূল্যায়ন
একবার আপনি আপনার পূর্বাভাস দেওয়ার পরে আপনাকে সেগুলি মূল্যায়ন করতে হবে এবং অন্যান্য মডেলের সাথে তুলনা করতে হবে। পূর্বাভাসের নির্ভুলতার মূল্যায়নের জন্য কিছু পদ্ধতি এখানে পাওয়া যাবে , তবে এটি করার জন্য যে লিঙ্কে ব্যাখ্যা করা হয়েছে, আপনি নিজের সিরিজটিকে একটি প্রশিক্ষণ এবং একটি পরীক্ষার সেটে ভাগ করে নেওয়া অত্যন্ত গুরুত্বপূর্ণ।