Editorial
Author
Abstract
One of the most important roles of a futures market is to provide the means of risk reduction. Optimal hedge strategy is determined via calculation of the hedge ratio. Estimation of hedge ratio and hedging effectiveness depend on correct specification of relation between spot and futures prices. Thus in this paper hedge ratio is estimated for the natural gas futures market by different methods e.g. OLS, VAR, VECM, and GARCH and their effectiveness is compared. In GARCH method, hedge ratio is time-varying so the time series of hedge ratio are estimated while the in other methods a fixed hedge ratio is estimated. Results show that GARCH hedge ratio has higher effectiveness rather than other methods and the effectiveness of other methods are ranked as: VECM, OLS and VAR respectively.
Keywords
برآورد نسبتهای بهینه پوشش ریسک ایستا و پویا و مقایسه میزان اثربخشی آنها در بازار آتیهای گاز طبیعی
محمد علیمرادی*
تاریخ دریافت: 7 آذر 1390 تاریخ پذیرش: 21 اسفند 1392
چکیده
یکی از مهمترین نقشهای بازار آتیها، فراهم کردن ابزاری برای پوشش ریسک است. استراتژی بهینه برای پوشش ریسک از طریق تخمین نسبت پوشش ریسک، مشخص میشود. محاسبه نسبت پوشش ریسک و همچنین میزان اثربخشی پوشش به تصریح درست رابطه بین قیمت آتیها و قیمت نقطهای بستگی دارد. از این رو در این مقاله نسبت بهینه پوشش ریسک با روشهای مختلف OLS، VAR، VECM و GARCH چندمتغیره، برای دادههای ماهانه بازار آتیهای گاز طبیعی در دوره زمانی 2011-2000 برآورد شده و سپس میزان اثربخشی آنها مورد مقایسه قرار میگیرد. در روش GARCH چندمتغیره، یک سری زمانی از نسبتهای بهینه پوشش به دست میآید در حالی که در سایر روشها نسبت بهینه پوشش منفردی برای کل دوره به دست میآید. از این رو روش GARCH یک روش پویا و سایر روشها ایستا هستند. براساس نتایج تحقیق، نسبت بهینه پوشش به دست آمده از روش GARCH چندمتغیره نسبت به سایر روشها از اثربخشی بالاتری برخوردار است.
واژههای کلیدی: بازار آتیها، نسبت پوشش ریسک، اثربخشی پوشش، گاز طبیعی، GARCH چندمتغیره
طبقهبندی JEL: G13، C13.
1. مقدمه
یکی از مهمترین کاربردهای مشتقات مالی از قبیل قرارداد آتیها، پوشش ریسک است. در سالهای گذشته اهالی دانشگاه و اهالی کسب و کار، علاقه بسیار زیادی به موضوع پوشش ریسک با آتیها نشان دادهاند. این امر از مقالات بسیار زیادی که در این حوزه نوشته شده است، کاملاً مشهود است. مفهوم پایهای پوشش ریسک، ترکیب سرمایهگذاریها در بازار نقطهای و آتیها برای شکل دادن سبدی از داراییها است که نوسانات در ارزش آن را کم کرده و یا حذف میکند. قیمت جاری یک قرارداد آتیها، تخمینی از قیمت نقطهای در سررسید قرارداد براساس اطلاعات قابل دسترس در بازار ارائه میدهد. در بلندمدت، قیمت آتیها با قیمت نقطهای به اضافه هزینه حمل در طی زمان شامل هزینه انبارداری، حق بیمه ریسک و ثمرات رفاهی برابر میشود از این رو مانع از شکلگیری آربیتراژ میشود.
بازار آتیها دو نقش مهم بر عهده دارد یکی انتقال ریسک و دیگری کشف قیمت. میتوان برای کاهش ریسک ناشی از شوکهای غیرقابل انتظار در قیمت نقطهای به طور همزمان در بازار آتیها شرکت کرد. اگر شرکتکنندگان در یک بازار آتیها، از بازار آتیها برای کاهش ریسک قرار گرفتن در موضعی[1] خاص در بازار نقطهای استفاده کنند تحت عنوان پوششدهنده ریسک[2] شناخته میشوند. پس از تعیین موضع آتیها یکی از مهمترین موضوعاتی که مطرح میشود این است که چه تعداد قرارداد آتیها برای پوشش ریسک موضع نقطهای، کفایت میکند. تعداد بهینه در موضع آتیها از طریق محاسبه نسبت پوشش ریسک (کواریانس شرطی بین قیمت آتیها و نقطهای تقسیم بر واریانس شرطی قیمت آتیها) تعیین میشود. پس از تعیین نسبت پوشش ریسک، استراتژی پوشش بهینه مشخص میشود.
اگر پوششدهنده ریسک به ازای نگهداری هر موضع نقطهای، به اندازه نسبت پوشش ریسک، در بازار آتیها در موضع فروش قرار گیرد میتواند ریسک خود را حداقل کند. محاسبه نسبت پوشش ریسک و همچنین میزان اثربخشی پوشش، بستگی به تصریح درست رابطه بین قیمت آتیها و قیمت نقطهای دارد.
برخی از شرکتکنندگان در بازار آتیها تلاش میکنند تا از فرصتهای آربیتراژ ناشی از انحراف کوتاهمدت از شرایط تعادلی بلندمدت به سود خود بهرهبرداری کنند که این گروه تحت عنوان سفتهبازان[3] شناخته میشوند. سفتهبازان در بازار تنها زمانی به انحرافات کوتاهمدت پاسخ میدهند که معتقد باشند این انحرافات هزینههای حمل در طی زمان را به درستی نشان نمیدهند. بنابراین میتوان رابطه پویای بین قیمت آتیها و نقطهای را به وسیله تفاوت بین قیمتهای نقطهای و آتیها (مبنا)[4] تخمین زد.
یکی از مهمترین موضوعات تئوریکی در پوشش ریسک، تعیین نسبت پوشش بهینه است. نسبت پوشش بهینه بستگی زیادی به نوع تابع هدفی که بهینهیابی میشود، دارد. اخیراً در همین راستا توابع هدف متعددی مورد استفاده قرار گرفتهاند. به عنوان مثال یکی از پرکاربردترین استراتژیهای پوشش، ریسک مبتنی بر حداقلسازی واریانس سبد پوشش داده شده است[5]. این روش تحت عنوان نسبت پوشش حداقل واریانس[6] شناخته میشود.
تخمین نسبت پوشش ریسک در وهله اول به نوع تابع هدف مورد استفاده بستگی دارد و در وهله دوم به نحوه تصریح مدل استفاده شده برای تخمین تجربی رابطه بین قیمت نقطهای و آتیها بستگی خواهد داشت. به عنوان مثال در روش حداقل واریانس یک روش تخمین ساده و معمول از نسبت پوشش ریسک بهینه، روش حداقل مربعات معمولی (OLS)[7] است. ضریب قیمت آتیها در تخمین OLS معادلهای که در آن قیمت نقطهای به صورت تابعی از قیمت آتیها تعریف شده است، اغلب به عنوان نسبت پوشش ریسک بهکار گرفته میشود. این روش چندان صحیح به نظر نمیرسد زیرا در این روش روند تاریخی قیمت به عنوان یکی از تعیینکنندههای قیمت جاری نادیده گرفته میشود. بنابراین مدل خودرگرسیون برداری (VAR) نمایش بهتری از این رابطه ارائه میدهد.
معمولاً سطح قیمتهای آتیها و نقطهای نامانا و تفاضل مرتبه اول آنها ماناست. بنابراین تخمین مدل VAR باید با استفاده از دیفرانسیل مرتبه اول قیمتها صورت بگیرد. استفاده از دیفرانسیل مرتبه اول در تخمین مدل VAR باعث میشود که این مدل نتواند همانباشتگی بین قیمتهای آتیها و نقطهای را به حساب آورد.
انگل و گرنجر (1987) نشان دادند که اگر همانباشتگی بین متغیرها وجود داشته باشد، نمایش VAR استاندارد، صحیح نیست. اضافه کردن یک جزء تصحیح خطا میتواند این مشکل را حل نماید از این رو، مدل تصحیح خطای برداری (VECM) برای تصریح رابطه بین قیمت نقطهای و آتیها بکار گرفته میشود که در آن از مبنا (تفاوت بین قیمت نقطهای و آتیها) به عنوان جزء تصحیح خطا استفاده میشود. قوش[8] (1993) عملکرد نسبتهای پوشش محاسبه شده از طریق مدل VAR استاندارد و مدل VECM را با هم مقایسه کرده است. ما نیز در این مقاله با بهرهگیری از روش حداقل واریانس با استفاده از روشهای مختلف، نسبت بهینه پوشش را برآورد کرده و سپس به مقایسه میزان اثربخشی نسبتهای پوشش ریسک محاسبه شده برای بازار آتیهای گاز طبیعی میپردازیم.
2. نسبت بهینه پوشش ریسک
در یک تعریف ساده، نسبت بهینه پوشش ریسک، نشاندهنده تعداد قرارداد آتیهای مورد نیاز برای جبران تغییر در ارزش یک قرارداد دارایی پایه[9] به دلیل تغییرات قیمت آن دارایی میباشد. این نسبت به این علت ساخته میشود که این اطمینان را در زمانی که قیمت نقدی دارایی کاهش یا افزایش مییابد، به معاملهگر بدهد که به اندازه کافی قرارداد آتیها برای حمایت مالی دارد. به عنوان مثال شرکتی که قصد دارد در زمانی مشخص در آینده میزان مشخصی نفت خام خریداری نماید این شرکت نگران افزایش قیمت نفت خام است، از این رو تصمیم میگیرد با خرید آتیهای نفت خام، افزایش احتمالی قیمت را به نحو مناسبی پوشش دهد. خطری که این شرکت را تهدید میکند آن است که هماکنون به درستی نمیداند قیمت جاری و قیمت آتیهای نفت خام، در روزی که قرار است نفت خام را از بازار نقطهای خریداری نماید، چیست. تفاوت قیمت نقطهای و قیمت آتیها را اصطلاحاً «ریسک مبنا»[10] میگویند.
اکنون اگر فرض شود که شرکتی قرار است Na واحد از یک دارایی پایه را در زمان 2t بخرد، برای پوشش دادن ریسک افزایش قیمت، این شرکت تصمیم میگیرد که Nf واحد قرارداد آتیها خریداری کند. نسبت پوشش ریسک را که با h نشان داده میشود برابر است با[11]:
اگر فرض کنیم که فردی که یک موضع غیرقابل معامله نقطهای در زمان t در اختیار دارد بخواهد برای یک دوره زمانی، استراتژی پوشش ریسک با حداقل واریانس اتخاذ کند که ریسکش را حداقل کند همانطور که اشاره شد باید موضع بهینه در بازار آتیها در پیش بگیرد که این موضع به وسیله محاسبه نسبت پوشش ریسک یعنی نسبتی از قرارداد آتیها که باید توسط این فرد در موضع فروش قرار بگیرد، مشخص میشود. به عنوان مثال اگر فرض شود St و 1+St نشاندهنده قیمت نقطهای در زمانهای t و 1+t باشد و Ft و 1+Ftهم قیمت قرارداد آتیها در دو زمان t و 1+t باشد و همچنین فرض شود که پوششدهنده ریسک قرارداد آتیها میفروشد. دارایی این فرد در پایان دوره به صورت زیر خواهد بود:
(1)
در این معادله نشاندهنده نسبت پوشش ریسک ()، تغییر در قیمت نقطهای و تغییر در قیمت آتیها است. نسبت پوشش ریسک بهینه، نسبتی است که واریانس را که مبتنی بر مجموعه اطلاعات در دسترس در زمان t (It) است و به صورت زیر قابل نمایش است، حداقل میکند:
(2)
نسبت بهینه پوشش ریسک که واریانس را حداقل میکند به وسیله نشان داده میشود که برابر است با[12]:
(3)
نسبت بهینه پوشش ریسک است زیرا در واقع آن مقداری از است که معادله (2) را حداقل میسازد. از آنجا که پوششدهنده ریسک ممکن است علاوه بر حداقل کردن واریانس بازدهیها به دنبال حداقل کردن تغییرپذیری و نوسانات عایدی سبد خود نیز باشد از این رو به جای سطح متغیر به لگاریتم متغیر توجه خواهد کرد. بنابراین اگر فرض کنیم که و باشند، آنگاه پوششدهنده ریسک واریانس را در معادله حداقل خواهد کرد. به عبارت بهتر به جای تفاضل متغیرها از تفاضل لگاریتم متغیرها استفاده میکند. از آنجایی که هم قیمت نقطهای و هم قیمت آتیها از ویژگی ناهمسانی واریانس[13] برخوردار هستند از این رو روش دوم مناسبتر خواهد بود. همانند قبل هدف حداقل کردن واریانس rm میباشد و مشابه حالت بالا میتوان نسبت بهینه پوشش ریسک را برآورد کرد که در این حالت نسبت پوشش بهینه به صورت زیر خواهد بود:
(4)
3. روشهای تخمین نسبت بهینه پوشش ریسک
نسبت بهینه پوشش ریسک به لحاظ تجربی اغلب از طریق رگرسیون خطی ساده تخمین زده میشود. بر این اساس اگر را متغیر مستقل و را متغیر وابسته در نظر بگیریم نرخ بهینه پوشش ریسک با تخمین β در معادله رگرسیونی زیر به دست خواهد آمد:
(5)
در این معادله جزء اخلال است. به طور مشخص این معادله را میتوان از طریق روش حداقل مربعات معمولی (OLS) تخمین زد. به راحتی ثابت میشود که تخمین β در این معادله برابر با در معادله (3) میباشد. اگر فرض شود که و آنگاه تخمین β در معادله (4) به صورت زیر خواهد بود:
که در آن و مقادیر انحراف از میانگین هستند. با توجه به این که و ، داریم:
همانطور که مشخص است تخمین از معادله (4) برابر با نسبت بهینه پوشش است. اگر چه از لحاظ تجربی تخمین معادله (4) بسیار آسان است اما روش حداقل مربعات معمولی از یک طرف از مجموعه اطلاعات موجود به طور کامل استفاده نمیکند و از طرف دیگر از گشتاورهای دوم غیرشرطی استفاده میکند که در نتیجه آن، نسبت پوشش بهینه ثابت برای دوره مورد مطالعه به دست میآورد. بنابراین تصریح معادله به شکل معادله (5) تصریح دقیقی از رابطه دو متغیر نیست.
برای استفاده کامل از تمام اطلاعات موجود در تخمین رابطه بین دو متغیر، روش خودرگرسیون برداری در مقایسه با روش OLS، روش مناسبتری است. این مدل بر پایه روند تاریخی قیمتها بنا میشود از این رو اطلاعات بیشتری در مدل وارد میشود. مدل خودرگرسیون برداری به صورت زیر قابل نمایش است:
(6)
(7)
در این دو معادله و اجزای اخلال بوده که مستقل از هم هستند. اجزای اخلال در معادلات (6) و (7) بیانگر میزان تغییراتی از قیمتهای نقطهای و آتیها هستند که توضیح داده نشده باقیماندهاند. بنابراین میتوان از ماتریس واریانس- کواریانس اجزای اخلال استفاده کرد و براساس معادله (4) نسبت بهینه پوشش ریسک را محاسبه کرد. اگر واریانس اجزای اخلال را به ترتیب با و و کواریانس بین این دو را با نشان دهیم نسبت بهینه پوشش ریسک خواهد شد.
معمولاً قیمت آتیها و نقطهای نامانا و انباشته از مرتبه یک هستند. در معادلات بالا تفاضل مرتبه اول سریهای زمانی استفاده شده بنابراین هر دو مانا بوده و برای تحلیل سری زمانی باکس-جنکینز[14] مناسب هستند. بنابراین معادلات بالا تصریح درستی دارند اما اگر قیمت نقطهای و آتیها همانباشتگی داشته باشند آنگاه اجزای اخلال در معادلات VAR کارایی خود را از دست خواهند داد چرا که در این حالت قسمتی از تغییرات قیمتهای نقطهای و آتیها از طریق یک رابطه خطی بین قیمت نقطهای و قیمت آتیها قابل توضیح است. از این رو این جزء را باید به معادلات اضافه کرد تا در اجزای اخلال اطلاعات اضافی باقی نمانده باشد. بنابراین مدل صحیحتر، مدل تصحیح خطا است که به صورت زیر خواهد بود:
(8)
(9)
در این معادلات، متغیر متغیر همانباشتگی است[15] و حضور آن، مدل تصحیح خطا را از مدل VAR متمایز میکند. وقتی که هر دو مدل دارای وقفههای یکسانی باشند، مدل تصحیح خطا، مدل VAR را هم دربر میگیرد. مشابه با حالت و مدل VAR میتوان نسبت پوشش ریسک را از مدل تصحیح خطا استخراج کرد. بنابراین به نظر میرسد که با توجه به تصریح دقیقتر در روش VECM باید اثربخشی نسبت بهینه پوشش در این روش نسبت به روش VAR بیشتر باشد.
یکی دیگر از مشکلاتی که در تخمین مدل با روش OLS وجود دارد و این مشکل در روشهای VAR و VECM همچنان باقی ماند، استفاده از گشتاورهای غیرشرطی مرتبه دوم در این مدلها است که به همین دلیل یک نسبت پوشش بهینه منفرد برای کل دوره به دست آمد.
مدلهایی که بتوانند هم از تمام اطلاعات موجود استفاده کنند و همچنین از گشتاورهای شرطی مرتبه دوم استفاده کنند مدلهای GARCH و ARCH سیستمی و یا چندمتغیره میباشد. خصوصاً این که در دادههای مربوط به بازدهی در بازارهای مالی مشاهده میشود که معمولاً تغییرات بزرگ در بازدهی، تغییرات بزرگتری را در بازدهی به دنبال دارد. به عبارت دیگر زمانی که بازدهی دارایی در بازار مالی تغییرات غیرمنتظره زیادی در جهت کاهش یا افزایش داشته باشد، آنگاه معاملهگر واریانس انتظاری بازدهی دارایی را برای دوره بعد افزایش خواهد داد. این ویژگیها در مدلهای ARCH و GARCH دو متغیره به خوبی لحاظ شده بنابراین تصریح درستی برای محاسبه نسبت پوشش بهینه خواهد بود.
از زمان توسعه مدلهای ARCH و GARACH، روش تخمین نسبت پوشش به گونهای تعدیل شده است که بتواند طبیعت ناهمسانی واریانس جزء خطا در معادله (4) را به حساب آورد. در این حالت به جای استفاده ساده از واریانس و کواریانس غیرشرطی در تخمین نسبت پوشش، از واریانس و کواریانس شرطی مدل GARCH استفاده میشود. این روش باعث پویایی نسبت پوشش در طی زمان و تغییر آن در دوره پوشش میشود.
مدل GARCH دو متغیره زیر مدلی است که میتوان از آن برای محاسبه نسبت پوشش ریسک استفاده کرد:
(10)
(11)
نسبت پوشش حداقل واریانس شرطی در زمان t را میتوان از معادله به دست آورد. این مدل این امکان را فراهم میآورد که نسبت پوشش در طی زمان تغییر کرده و در نتیجه به جای یک نسبت پوشش واحد به یک سری نسبت پوشش دست پیدا کنیم. پس از انجام مراحل فوق نسبتهای پوشش ریسک برای هر چهار روش OLS، VAR، تصحیح خطا و GARCH محاسبه شده و مقایسه خواهد شد.
4. پیشینه تحقیق
در داخل کشور مطالعه چندانی در این زمینه صورت نگرفته است. نائینی و کاظمیمنش (1383)، در مقالهای نسبت بهینه پوشش ریسک را در دوره پنجساله، 1999 تا 2003 با استفاده از سری زمانی هفتگی قیمت نفت خام WTI و قراردادهای آتیهای NYMEX مورد مطالعه قرار دادند. بدین منظور، از مدلهای ARCH و GARCH استفاده شده است و نتایج تخمین نشان میدهد که با افزایش دوره قرارداد آتیها، نسبت بهینه پوشش ریسک بزرگتر میشود. ابراهیمی و قنبری هم در مطالعهای با استفاده از قراردادهای آتیهای یک تا چهار ماهه بورس نایمکس و روشهای OLS، VAR و ECM نسبت پوشش بهینه را محاسبه کرده و حالات مختلفی برای استراتژی پوشش ریسک به دست آوردهاند. در پایان نتیجه میگیرند که قراردادهای آتیهای چهارماهه با نسبت پوشش مدل ECM برای پوشش ریسک و ثبات درآمدهای نفتی ایران مناسبترین گزینه است.
در سطح بینالمللی هم مطالعات مرتبطی صورت گرفته است که در ادامه به بررسی برخی از مهمترین این مطالعات پرداخته خواهد شد. در یکی از مهمترین مطالعاتی که در این زمینه صورت گرفته است، شنگ- سیان چن، چنگ- فو لی و کشاب سرسدا (2003) در مقالهای به بازبینی روشهای مختلف نظری درباره نسبتهای بهینه پوشش ریسک آتیها میپردازد. در این مقاله روشهای OLS، ARCH/GARCH، ضریب تصادفی و همانباشتگی و تصحیح خطا مورد بررسی قرار میگیرد و در نهایت نتیجه میگیرد که نسبتهای بهینه پوشش ریسک براساس روشهای مختلف متفاوت بوده و هیچ نسبت پوشش بهینه منفردی وجود ندارد که به طور مشخص بر سایر نسبتها تفوق داشته باشد.
در مطالعه دیگری توفیق چودری[16] (2003) به بررسی و آزمون آثار رابطه بلندمدت بین شاخص نقدی و شاخص آتیهای سهام بر اثربخشی پوشش شش بازار آتیهای سهام میپردازد. اثربخشی پنج نسبت پوشش مختلف مورد مقایسه قرار میگیرد و در نهایت به این نتیجه میرسد که در میان نسبتهای پوشش، نسبتهای پوشش متغیر که از روش GARCH به دست میآیند نسبت به نسبتهای پوشش ثابت عملکرد بهتری دارند.
توفیق چودری (2004) در مطالعه دیگری به بررسی میزان اثربخشی پوشش در بازار آتیهای ژاپن، استرالیا و هنگکنگ میپردازد و در پایان باز هم نتیجه میگیرد که برای اغلب موارد عملکرد نسبت پوشش متغیر به دست آمده از روش GARCH نسبت به نسبت پوشش ثابت عملکرد بهتری دارد.
چونرانگ آی، آرجون چاترا و فرانک سونگ (2007) در مقاله خود با تکیه بر روش جدید نیمه پارامتریک و با استفاده از آمار بازارهای غلات، پنبه و سویای آمریکا به تخمین نسبت پوشش بهینه ریسک میپردازد و نتایج را با تخمینهای حاصل از روشهای پارامتریک اعم از OLS، VAR و GARCH مورد مقایسه قرار میدهد و نتیجه میگیرد که برای افزایش اثربخشی پوشش ریسک در بازار آتیها باید بر روی مدلهای ناپارامتریک تمرکز کرد هر چند مدلهای نیمهپارامتریک ممکن است با کنترل برخی شرایط عملکرد بهتری نسبت به سایر مدلها داشته باشند.
در همین ارتباط مقاله دیگری توسط دونالد لین (2004) در فصلنامه بررسی اقتصاد و مالی به چاپ رسیده است. در این مقاله، نویسنده به ارزیابی آثار نادیده گرفتن رابطه همانباشتگی بین قیمتهای نقطهای و آتیها بر نسبت پوشش ریسک بهینه و اثربخشی پوشش پرداخته است. وی در پایان نتیجه میگیرد که نادیده گرفتن همانباشتگی بین این دو متغیر باعث میشود که نسبت پوشش بهینه کوچکتری تخمین زده شود.
گپرت (1995) نیز در مقاله خود به تخمین نسبت بهینه پوشش ریسک حداقل واریانس با روش OLS و روشهای همانباشتگی برای طولهای متفاوت از افق برنامهریزی میپردازد. برای نمونه مورد نظر نتیجه میگیرد که اثربخشی در هر دو روش با افزایش طول افق برنامهریزی بهبود مییابد.
کرونر و سلتان (1993) در مقالهای با استفاده از مدل تصحیح خطا و GARCH به تخمین نسبت پوشش ریسک حداقل واریانس برای پنج ارز میپردازند. در این مقاله با استفاده از مدلهای شرطی، سریهای زمانی از نسبتهای پوشش ریسک تخمین زده میشود. شواهد حاکی از آن است که استراتژی پوشش ریسک توصیه شده در این مقاله نسبت به استراتژیهای مرسوم عملکرد بهتری دارند.
5. دادهها و برآورد نسبت پوشش بهینه ریسک
در این مقاله از آمار و اطلاعات مربوط به قیمتهای نقطه ای و آتیهای بازار گاز طبیعی استفاده شده است که از دفتر اطلاعات انرژی آمریکا (EIA)[17] استخراج شده است. منبع این قیمتها در این سازمان آماری برای قیمتهای نقطهای، تامسون رویترز[18] و برای قیمت آتیها، بورس نایمکس[19] اعلام شده است. همچنین قابل ذکر است که قیمت آتیها در این مقاله مربوط به چهار نوع قرارداد آتیها است که با قرارداد یک، قرارداد دو، قرارداد سه و قرارداد چهار نامگذاری شدهاند. هر قرارداد آتیها یک روز تحویل مشخص دارد که با آن شناخته میشود. قرارداد آتیهای گاز طبیعی سه روز کاری قبل از اولین روز تقویمی ماه تحویل منقضی میشوند. از این رو قرارداد یک، قراردادی است که ماه تحویل در آنها اولین ماه تقویمی بعد از روز معامله است و به همین ترتیب، قراردادهای دو تا چهار، قراردادهایی هستند که ماه تحویل در آنها به ترتیب دومین، سومین و چهارمین ماه تقویمی بعد از روز معامله است.
در این مقاله برای مقایسه میزان اثربخشی نسبتهای پوشش بهینه از روشی ساده استفاده میشود. برای نیل به این هدف از معادله (1) استفاده میشود. با در اختیار داشتن برآوردهای مختلف از نسبت بهینه پوشش ریسک و همچنین داشتن سری زمانی قیمت نقطهای و آتیهای گاز طبیعی، تنها مجهول معادله تعداد قراردادهای نقطهای میباشد از این رو بازدهی دارایی مورد نظر برحسب تعداد قراردادهای نقطهای حاصل خواهد شد. با انجام این فرآیند میتوان به یک سری زمانی از بازدهی با تناوب ماهانه دست پیدا کرد. در ادامه، انحراف معیار این سری زمانی مورد محاسبه قرار میگیرد. انحراف معیار کمتر بازدهی، حاکی از اثربخشی بیشتر نسبت پوشش ریسک خواهد بود از این رو اثربخشی براساس نسبت انحراف معیارها قابل محاسبه است.
در ادامه این بخش به بررسی مدلهای استفاده شده برای برآورد نسبت بهینه پوشش حداقل واریانس و نتایج حاصل از آنها پرداخته میشود:
1. روش OLS: در سادهترین روش برای تخمین نسبت پوشش بهینه، معادله (5) با استفاده از روش OLS تخمین زده شده است. همانطور که پیشتر توضیح داده شده و اثبات شد ضریب در این معادله، نسبت بهینه پوشش ریسک است. این تخمین برای قیمت آتیهای چهار نوع قرارداد در چهار مدل متفاوت صورت گرفته است که نتایج حاصل از این تخمین، برای در جدول 1 آمده است. همچنین در جدول 1 میزان اثربخشی نسبت بهینه پوشش برحسب انواع قراردادها آورده شده است. در میان قراردادهای مختلف و روش OLS، قرارداد چهار اثربخشی بالاتری نسبت به سایر انواع قراردادها دارد.
جدول 1. برآوردهای نسبت بهینه پوشش از روش OLS
|
قرارداد یک |
قرارداد دو |
قرارداد سه |
قرارداد چهار |
برآورد |
008491/1 |
011888/1 |
005470/1 |
048077/1 |
اثربخشی |
9871/0 |
9785/0 |
9826/0 |
9926/0 |
2. روش VAR و VECM: در این دو روش معادلات 6 تا 9 تخمین زده میشوند. در بدو امر ضروری است مانا یا نامانا بودن متغیرها مورد بررسی قرار بگیرد. از این رو آزمون ریشه واحد براساس آزمون دیکی- فولر در مورد متغیرها صورت گرفته است که براساس نتایج آن که در جدول 2 آمده است، قیمت نقطهای و قیمت آتیهای گاز طبیعی در سطح اطمینان 95 درصد هر دو نامانا و از مرتبه یک هستند از این رو تفاضل مرتبه اول آنها مانا خواهد بود. آنچه که ما در تخمین معادلات VAR به دنبال آن هستیم ماتریس واریانس- کواریانس باقیماندهها برای محاسبه نسبت بهینه پوشش ریسک میباشد. بنابراین پس از تخمین معادلات و استخراج باقیماندههای معادلات (6) و (7)، واریانس و کواریانس باقیماندهها به دست میآیند که این ماتریس در جدول 3 نمایش داده شده است.
با توجه به اینکه متغیرهای موجود در معادله VAR به صورت تفاضلی هستند از این رو اگر بین متغیرها همانباشتگی وجود داشته باشد این نوع تصریح معادله صحیح نخواهد بود و باید جزء تصحیح خطا به معادله افزوده شود. از این رو در ادامه به بررسی و آزمون همانباشتگی دو متغیر قیمت آتیها و قیمت نقطهای پرداخته شده است. آزمون مورد نظر با استفاده از آزمون جوهانسون بررسی شده است که نتایج آن در جدول 2 آمده است که براساس آن حداقل یک رابطه همانباشتگی بین دو متغیر قابل تصور است. با رد فرض عدم وجود رابطه همانباشتگی بین دو متغیر مورد نظر ضروری است که جزء تصحیح خطا به معادله افزوده شود. از این رو در ادامه تخمین VECM صورت گرفته است و در جدول 2 ماتریس واریانس- کواریانس باقیماندههای تخمین معادلات (8) و (9) در این روش آمده است.
همانطور که پیشتر توضیح داده شد، اگر ماتریس واریانس- کواریانس باقیماندههای مدل VAR یا VECM با نشان داده شود آنگاه نسبت بهینه پوشش ریسک خواهد شد. بنابراین با توجه به ماتریسهای واریانس– کواریانس حاصل از تخمین این دو مدل که در جداول 3 و 4 نشان داده شده است نسبت پوشش بهینه برای قراردادهای نوع اول به صورت زیر خواهد بود:
در مقایسه بین نسبتهای بهینه پوشش محاسبه شده برای قراردادهای نوع اول در سه روش فوق که یک نسبت پوشش منفرد را میدهند میتوان مشاهده کرد که رابطه بین آنها برقرار است. برای سایر انواع قراردادها نسبت بهینه پوشش ریسک به همین ترتیب محاسبه میشود که نتایج آن در جداول 3 و 4 آمده است. در این جداول میزان اثربخشی نسبتهای بهینه پوشش هم آمده است. همانطور که مشخص است اثربخشی نسبت پوشش ریسک در روش VECM و قراردادهای یک ماهه بیشترین حالت ممکن است.
جدول 2. نتایج آزمون ریشه واحد
|
قیمت آتیهای قرارداد یک |
قیمت آتیهای قرارداد دو |
قیمت آتیهای قرارداد سه |
قیمت آتیهای قرارداد چهار |
قیمت نقطهای |
|
ریشه واحد |
آماره t |
7073/1- |
8040/2- |
6065/2- |
5932/2- |
5647/2- |
احتمال |
4252/0 |
0603/0 |
0941/0 |
0968/0 |
1028/0 |
|
نتیجه آزمون |
نامانا |
نامانا |
نامانا |
نامانا |
نامانا |
|
همانباشتگی با قیمت نقطهای |
آماره |
129/6 |
621/5 |
762/5 |
399/5 |
|
احتمال |
0133/0 |
0177/0 |
01644/0 |
0201/0 |
|
|
نتیجه آزمون |
همانباشته |
همانباشته |
همانباشته |
همانباشته |
|
منبع: یافتههای تحقیق
جدول 3. ماتریس واریانس- کواریانس باقیماندهها در روش VAR
واریانس- کواریانس |
|||||||||
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
||
1 |
01918/0 |
|
|
|
01728/0 |
|
|
|
|
2 |
|
01949/0 |
|
|
|
01369/0 |
|
|
|
3 |
|
|
02071/0 |
|
|
|
01278/0 |
|
|
4 |
|
|
|
0210/0 |
|
|
|
0096/0 |
|
1 |
01728/0 |
|
|
|
01690/0 |
|
|
|
|
2 |
|
01369/0 |
|
|
|
01291/0 |
|
|
|
3 |
|
|
01278/0 |
|
|
|
01210/0 |
|
|
4 |
|
|
|
0096/0 |
|
|
|
0084/0 |
|
نسبت بهینه پوشش ریسک |
0225/1 |
0604/1 |
0562/1 |
1429/1 |
|||||
اثربخشی |
9859/0 |
9743/0 |
9788/0 |
9848/0 |
منبع: یافتههای تحقیق
جدول 4. ماتریس واریانس- کواریانس باقیماندهها در روش VECM
واریانس - کواریانس |
|||||||||
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
||
1 |
01875/0 |
|
|
|
01695/0 |
|
|
|
|
2 |
|
01794/0 |
|
|
|
01341/0 |
|
|
|
3 |
|
|
01948/0 |
|
|
|
01284/0 |
|
|
4 |
|
|
|
02/0 |
|
|
|
0101/0 |
|
1 |
01695/0 |
|
|
|
01685/0 |
|
|
|
|
2 |
|
01341/0 |
|
|
|
01292/0 |
|
|
|
3 |
|
|
01284/0 |
|
|
|
01219/0 |
|
|
4 |
|
|
|
0101/0 |
|
|
|
0085/0 |
|
نسبت بهینه پوشش ریسک |
0059/1 |
0379/1 |
0533/1 |
1906/1 |
|||||
اثربخشی |
9873/0 |
9763/0 |
979/0 |
9835/0 |
منبع: یافتههای تحقیق
3. روش GARCH: معادله (10) و (11) تصریحی از یک مدل GARCH دو متغیره را نشان میدهد. پس از تخمین این دو معادله ماتریس واریانس کواریانس قابل استخراج خواهد بود. همانطور که قبلاً نیز توضیح داده شد، این روش یک سری زمانی از واریانسها و کواریانسها را خواهد داد. مشابه قبل با دستیابی به ماتریس واریانس کواریانس میتوان نسبت بهینه پوشش ریسک را محاسبه کرد که خلاصه نتایج آن تنها برای ماههای ژانویه و برای آتیهای نوع اول در جدول 5 آمده است. همانطور که از جدول مشخص است در این حالت نسبتهای بهینه پوشش مختلفی به دست میآید که این امر مزیت محاسبه از روش GARCH را نشان میدهد. هر چند باید میزان اثربخشی این روش را هم سنجید تا معلوم گردد که آیا واقعاً نسبت پوشش متغیر در طی زمان به دست آمده از روش GARCH، عملکرد بهتری نسبت به سایر روشها دارد یا خیر که این موضوع در ادامه بررسی خواهد شد.
جدول 5. بخشی از سری زمانی نسبت بهینه پوشش ریسک در روش GARCH
|
ژانویه 2000 |
ژانویه 2001 |
ژانویه 2002 |
ژانویه 2003 |
ژانویه 2004 |
ژانویه 2005 |
ژانویه 2006 |
ژانویه 2007 |
ژانویه 2008 |
ژانویه 2009 |
ژانویه 2010 |
ژانویه 2011 |
0154/0 |
0233/0 |
0153/0 |
0175/0 |
0182/0 |
0161/0 |
0391/0 |
0164/0 |
0164/0 |
0183/0 |
0167/0 |
0155/0 |
|
0144/0 |
0263/0 |
0224/0 |
0176/0 |
0197/0 |
0185/0 |
397/0 |
0179/0 |
0165/0 |
0196/0 |
0174/0 |
0151/0 |
|
065/1 |
887/0 |
682/0 |
996/0 |
926/0 |
874/0 |
986/0 |
918/0 |
995/0 |
934/0 |
038/1 |
023/1 |
6. اثربخشی نسبت بهینه پوشش
همانطور که پیشتر نیز اشاره شد، میتوان براساس نسبتهای پوشش مختلفی که از روشهای OLS، VAR، VECM و GARCH چندمتغیره، به دست میآید سری زمانی بازدهیها و انحراف معیار مربوط به آنها را به دست آورد. برای مقایسه میزان اثربخشی نسبتهای پوشش مختلف، از حاصل تقسیم انحراف معیارها استفاده میشود. به این صورت که هر یک از انحراف معیارها بر سایر انحراف معیارها تقسیم میشود. از این رو نسبت کمتر نشاندهنده اثربخشی بیشتر خواهد بود.
نتایج حاصل از محاسبه انحراف معیارها و ماتریس نسبت اثربخشی آنها در جدول 6 آمده است. در این ماتریس هر یک از درایهها نشاندهنده نسبت اثربخشی عناصر سطر به عناصر ستون است. همانطور که از عناصر این جدول مشخص است اعداد ستون نهایی همگی بزرگتر از یک هستند که این امر نشان میدهد که اثربخشی نسبت بهینه پوشش ریسک با استفاده از قراردادهای چهارماهه و روش GARCH چندمتغیره در نسبت به سایر روشها بیشتر است.
جدول 6. ماتریس اثربخشی نسبتهای بهینه پوشش
|
OLS |
VAR |
VECM |
GARCH |
|||||||||||||
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
||
OLS
|
1 |
1/0000 |
0/9913 |
0/9955 |
1/0056 |
0/9988 |
0/9871 |
0/9916 |
0/9977 |
1/0002 |
0/9891 |
0/9918 |
0/9964 |
1/0026 |
0/9917 |
0/9987 |
1/0131 |
2 |
1/0087 |
1/0000 |
1/0042 |
1/0144 |
1/0075 |
0/9957 |
1/0003 |
1/0065 |
1/0090 |
0/9977 |
1/0005 |
1/0051 |
1/0114 |
1/0004 |
1/0075 |
1/0220 |
|
3 |
1/0045 |
0/9958 |
1/0000 |
1/0101 |
1/0033 |
0/9915 |
0/9961 |
1/0022 |
1/0047 |
0/9935 |
0/9963 |
1/0008 |
1/0071 |
0/9962 |
1/0032 |
1/0177 |
|
4 |
0/9945 |
0/9858 |
0/9900 |
1/0000 |
0/9933 |
0/9816 |
0/9861 |
0/9922 |
0/9947 |
0/9836 |
0/9863 |
0/9908 |
0/9970 |
0/9862 |
0/9932 |
1/0075 |
|
VAR
|
1 |
1/0012 |
0/9925 |
0/9967 |
1/0068 |
1/0000 |
0/9882 |
0/9928 |
0/9989 |
1/0014 |
0/9903 |
0/9930 |
0/9975 |
1/0038 |
0/9929 |
0/9999 |
1/0143 |
2 |
1/0131 |
1/0043 |
1/0086 |
1/0188 |
1/0119 |
1/0000 |
1/0046 |
1/0108 |
1/0133 |
1/0020 |
1/0048 |
1/0094 |
1/0157 |
1/0047 |
1/0118 |
1/0264 |
|
3 |
1/0085 |
0/9997 |
1/0040 |
1/0141 |
1/0073 |
0/9954 |
1/0000 |
1/0062 |
1/0087 |
0/9975 |
1/0002 |
1/0048 |
1/0111 |
1/0001 |
1/0072 |
1/0217 |
|
4 |
1/0023 |
0/9936 |
0/9978 |
1/0079 |
1/0011 |
0/9893 |
0/9939 |
1/0000 |
1/0025 |
0/9913 |
0/9941 |
0/9986 |
1/0049 |
0/9940 |
1/0010 |
1/0154 |
|
VECM |
1 |
0/9998 |
0/9911 |
0/9953 |
1/0053 |
0/9986 |
0/9868 |
0/9914 |
0/9975 |
1/0000 |
0/9888 |
0/9916 |
0/9961 |
1/0023 |
0/9915 |
0/9985 |
1/0129 |
2 |
1/0110 |
1/0023 |
1/0065 |
1/0167 |
1/0098 |
0/9980 |
1/0025 |
1/0087 |
1/0113 |
1/0000 |
1/0028 |
1/0074 |
1/0137 |
1/0026 |
1/0098 |
1/0243 |
|
3 |
1/0037 |
0/9950 |
0/9992 |
1/0093 |
1/0025 |
0/9907 |
0/9952 |
1/0014 |
1/0039 |
0/9927 |
0/9955 |
1/0000 |
1/0062 |
0/9953 |
1/0024 |
1/0168 |
|
4 |
1/0037 |
0/9950 |
0/9992 |
1/0093 |
1/0025 |
0/9907 |
0/9952 |
1/0014 |
1/0039 |
0/9927 |
0/9955 |
1/0000 |
1/0062 |
0/9953 |
1/0024 |
1/0168 |
|
GARCH |
1 |
0/9974 |
0/9888 |
0/9930 |
1/0030 |
0/9962 |
0/9845 |
0/9890 |
0/9952 |
0/9977 |
0/9865 |
0/9893 |
0/9938 |
1/0000 |
0.9891 |
0/9962 |
1/0105 |
2 |
1/0084 |
0/9996 |
1/0039 |
1/0140 |
1/0072 |
0/9953 |
0/9999 |
1/0061 |
1/0086 |
0/9974 |
1/0001 |
1/0047 |
1/0110 |
1/0000 |
1.0071 |
1/0216 |
|
3 |
1/0013 |
0/9926 |
0/9968 |
1/0068 |
1/0001 |
0/9883 |
0/9928 |
0/9990 |
1/0015 |
0/9903 |
0/9931 |
0/9976 |
1/0038 |
0/9929 |
1/0000 |
1/0144 |
|
4 |
0/9871 |
0/9785 |
0/9826 |
0/9926 |
0/9859 |
0/9743 |
0/9788 |
0/9848 |
0/9873 |
0/9763 |
0/9790 |
0/9835 |
0/9896 |
0/9789 |
0/9858 |
1/0000 |
براساس سایر ارقام موجود میتوان نتیجه گرفت که نسبت بهینه پوشش حاصل از روش OLS نسبت به روشهای VAR و VECM و نسبت بهینه پوشش حاصل از روش VECM نسبت به روش VAR اثربخشی بیشتری دارد. در مقایسه بین نوع قراردادها هم کاملاً مشخص است که قراردادهای نوع چهارم نسبت به سایر روشها برای پوشش ریسک مناسبتر هستند که در روش GARCH بیشترین اثربخشی و در روش OLS در رتبه دوم اثربخشی قرار دارد.
7. نتیجهگیری
در این مقاله ما به بررسی و مقایسه میزان اثربخشی نسبتهای مختلف پوشش ریسک (به دست آمده از روش حداقل واریانس) در بازار گاز طبیعی و برای دادههای ماهانه در طی دوره ژانویه 2000- اکتبر 2011 پرداختیم. در راستای تخمین نسبت پوشش در روش حداقل واریانس از چهار مدل OLS، VAR، VECM و GARCH استفاده شد. در ابتدا از مدل OLS استفاده شد اما این روش به دلیل استفاده نکردن از تمام اطلاعات موجود ناقص بود از این رو در ادامه از مدلهای VAR و VECM استفاده شد. این روشها اگر چه نسبت به روش اول کارایی بیشتری در استفاده از تمام اطلاعات داشتند اما وجود واریانس ناهمسانی در متغیرها باعث میشود که از کارایی این مدلها کاسته شود. بنابراین برای رفع این نقصانها از یک سیستم GARCH دومتغیره برای مدلسازی این رابطه استفاده شد.
در تعداد زیادی از مقالات مرتبط با بازار آتیها ادعا شده است که باید نسبت بهینه پوشش در طی زمان متغیر باشد نه ثابت. از این رو استفاده از مدل GARCH در این مقاله باعث شد که علاوه بر مرتفع شدن نقصان روشهای فوق به جای دستیابی به یک نسبت بهینه پوشش منفرد به یک سری از نسبتهای پوشش دست پیدا کنیم.
در این مقاله مقایسه میزان اثربخشی از طریق محاسبه میزان انحراف معیار بازدهی دارایی پس از پوشش ریسک با استفاده از نسبت پوشش به دست آمده از هر روش صورت گرفت و نتایج نشان داد که مدلسازی رابطه بین قیمتهای نقطهای و آتیها در قالب یک سیستم GARCH دومتغیره، اثربخشی نسبت پوشش بهینه را در بازار آتیهای گاز طبیعی بهبود بخشید. به عبارت دیگر نتایج نشان میدهد که نسبتهای پوشش متغیر در طی زمان عملکرد بهتری نسبت به نسبتهای پوشش ثابت دارند. از این رو پوششدهندگان ریسک که میخواهند در بازار ریسک نوسانات قیمت گاز طبیعی در بازار را پوشش دهند میتوانند بهترین نسبت پوشش ریسک را در قالب سیستم GARCH دومتغیره پیشبینی کنند.
منابع
الف- فارسی
ابراهیمی، محسن و علیرضا قنبری (1388)، «پوشش ریسک نوسانات درآمدهای نفتی با استفاده از قراردادهای آتی در ایران»، تهران، پژوهشنامه اقتصادی، سال نهم، شماره سوم
جلالی نائینی، سید احمدرضا و مریم کاظمیمنش (1383)، «بررسی تغییرات نرخ بهینه پوشش ریسک در بازار نفت»، تهران، فصلنامه مطالعات اقتصاد انرژی، سال اول، شماره یک
درخشان، مسعود (1390)، مشتقات و مدیریت ریسک در بازارهای نفت، تهران، مؤسسه مطالعات بینالمللی انرژی، چاپ دوم 1390
ب- انگلیسی
Ai, C., Chatrath, A. & F. Song (2007), “A Semiparametric Estimation of the Optimal Hedge Ratio”, The Quarterly Review of Economics and Finance, No. 47, PP. 366-381.
Chen, S. S., Lee, C. F. & K. Shrestha (2003), “Futures Hedge Ratios: A Review”, The Quarterly Review of Economics and Finance, No. 43, PP. 433-465.
Choudhry, T. (2003), “Short-run Deviations and Optimal Hedge Ratio: Evidence from Stock Futures”, Journal Multinational Financial Management, No. 13, PP. 171-192.
Choudhry, T. (2004), “The Hedging Effectiveness of Constant and Time-varying Hedge Ratios Using Three Pacific Basin Stock Futures”, International Review of Economics & Finance, No. 13, PP. 371-385.
Geppert, J. M. (1995), “A Statistical Model for the Relationship between Futures Contract Hedging Effectiveness and Investment Horizon Length”, Journal of Futures Markets, No. 15, PP. 507-536.
Ghosh, A. (1993), “Hedging with Stock Index Futures: Estimation and Forecasting with Error Correction Model”, Journal of Futures Markets, No. 13, PP. 295-305.
Kroner, K. F. & J. Sultan (1993), “Time-varying Distributions and Dynamic Hedging with Foreign Currency Futures”, Journal of Financial and Quantitive Analysis, No. 28, PP. 535-551.
Lein, D. (2004), “Cointegration and the Optimal Hedge Ratio: the General Case”, The Quarterly Review of Economics and Finance, No. 44, PP. 654-658.
* دانشجوی دکتری اقتصاد نفت و گاز دانشگاه علامه طباطبائی mdalimoradi@yahoo.com
[1]. Position
[2]. Hedger
[3]. Speculator
[4]. Basis
[5]. Ederington (1979)
[6]. Minimum Variance (MV)
[7]. Ordinary Least Square
[8]. Ghosh
[9]. منظور از قرارداد، قرارداد استانداردی است که بسته به نوع کالا متفاوت بوده و در بورس معامله میشود.
[10]. Basis Risk
[11]. درخشان (1390)
[12]. برای توضیحات بیشتر و اثبات معادله مراجعه شود به: درخشان (1390)
[13]. Heteroscedastic
[14]. Box-Jenkins
[15]. این متغیر همانباشتگی، یک ترکیب خطی از دو متغیر قیمت نقطهای و قیمت آتیها (به صورت لگاریتمی) است و از آنجایی که رابطه بین این دو متغیر یک رابطه یک به یک است بنابراین میتوان این نتیجه را گرفت که این رابطه خطی، به طور معناداری نزدیک به تفاضل این دو متغیر و به بیان دیگر ریسک مبنا باشد.
[16]. Taufiq Choudhry
[17]. Energy Information Administration
[19]. New York Mercantile Exchange (NYMEX)