شناسنامه علمی شماره

نویسندگان

1 عضو هیات علمی مؤسسه عالی آموزش و پژوهش مدیریت و برنامه‌ریزی

2 دانشجوی دکتری اقتصاد نفت و گاز دانشگاه علامه طباطبایی

چکیده

تلاطم قیمت نفت در بازارهای بین‌المللی، بازیگران این بازار را در معرض ریسک‌های بالقوه زیادی قرار می‌دهد. روش ارزش در معرض ریسک (VaR) از مهم‌ترین شیوه‌های اندازه‌گیری ریسک بازار می‌باشد که در سال‌های اخیر توسعه زیادی پیدا کرده است. در این تحقیق ریسک‌های فراسوی و فروسوی بازدهی قیمت نفت خام WTI با استفاده از مدل GED- GARCH که برای توزیع‌های کشیده با دنباله چاق مناسب است، برآورده شده است. داده‌های روزانه قیمت نقدی و آتی‌های WTI درون نمونه از ژانویه 1986 تا دسامبر 2010 و داده‌های خارج از نمونه نیز از ژانویه 2011 تا جولای 2012 در نظر گرفته شده است. برای بررسی صحت مدل VaR برآورد شده از آزمون کوپیک استفاده شده و همچنین با استفاده از روش علیت گرنجری در ریسک به بررسی اثر سرریز ریسک میان بازدهی قیمت آتی‌ها و نقدی برای شاخص نفت خام WTI در بازار نایمکس پرداخته شده است. یافته‌های تحقیق نشان می‌دهد که توزیع آماری بازدهی قیمت در بازار نقدی و آتی‌های نفت WTI، توزیعی با دنباله چاق می‌باشد. بررسی‌ها همچنین اثر سرریز ریسک فراسوی معنی‌داری در سطح اطمینان 99 درصد از آتی‌ها به نقدی را تأیید می‌کند که نشان می‌دهد در زمان‌های افزایش قیمت نفت مانند دهه 2000، ریسک بازار آتی‌ها به نقدی منتقل ‌شده است.

کلیدواژه‌ها

عنوان مقاله [English]

Risk Spillover Effect between Oil Spot and Futures Price Returns

نویسندگان [English]

  • Ahmadreza Jalali Naiini 1
  • Vahid Ghorbani Pashakolae 2
  • Mohamad Sayadi 2

1

2

چکیده [English]

Due to price volatility in the oil market, market players are exposed to large risks. Value at Risk (VaR) is one of the main methods to measure market risk in various asset markets including commodities.,. In this study, Upside and Downside Risks are estimated by using the GED-GARCH method that is appropriate for leptokurtic distributions with fat tail. The daily spot and Futures oil prices data from January 1986 to December 2010 data for "in sample" and from January 2011 to July 2012 for "out of sample" are our data sample. To test the reliability of estimated VaR, the Kupiec test is used. Also by using Granger Causality analysis, the spillover effect risk between spot and futures oil price returns are investigated. Results show that spot and futures returns have leptokurtic distribution with fat tails. There is also a significant upside spillover effect risk from futures to spot price returns at 99% confidence level as for oil price increases during 2000s.

کلیدواژه‌ها [English]

  • Value at Risk
  • GED-GARCH Model
  • Upside and Downside Risk Spillover Effect
  • Kupiec Test del

اثر سرریز ریسک بین بازدهی قیمت در بازارهای نقدی و آتی‌های نفت خام

سیداحمدرضا جلالی نائینی*، وحید قربانی پاشاکلایی** و محمد صیادی***

 

تاریخ دریافت: 8 آبان 1392                       تاریخ پذیرش: 21 اسفند 1392

 

چکیده

تلاطم قیمت نفت در بازارهای بین‌المللی، بازیگران این بازار را در معرض ریسک‌های بالقوه زیادی قرار می‌دهد. روش ارزش در معرض ریسک (VaR) از مهم‌ترین شیوه‌های اندازه‌گیری ریسک بازار می‌باشد که در سال‌های اخیر توسعه زیادی پیدا کرده است. در این تحقیق ریسک‌های فراسوی و فروسوی بازدهی قیمت نفت خام WTI با استفاده از مدل GED- GARCH که برای توزیع‌های کشیده با دنباله چاق مناسب است، برآورده شده است. داده‌های روزانه قیمت نقدی و آتی‌های WTI درون نمونه از ژانویه 1986 تا دسامبر 2010 و داده‌های خارج از نمونه نیز از ژانویه 2011 تا جولای 2012 در نظر گرفته شده است. برای بررسی صحت مدل VaR برآورد شده از آزمون کوپیک استفاده شده و همچنین با استفاده از روش علیت گرنجری در ریسک به بررسی اثر سرریز ریسک میان بازدهی قیمت آتی‌ها و نقدی برای شاخص نفت خام WTI در بازار نایمکس پرداخته شده است. یافته‌های تحقیق نشان می‌دهد که توزیع آماری بازدهی قیمت در بازار نقدی و آتی‌های نفت WTI، توزیعی با دنباله چاق می‌باشد. بررسی‌ها همچنین اثر سرریز ریسک فراسوی معنی‌داری در سطح اطمینان 99 درصد از آتی‌ها به نقدی را تأیید می‌کند که نشان می‌دهد در زمان‌های افزایش قیمت نفت مانند دهه 2000، ریسک بازار آتی‌ها به نقدی منتقل ‌شده است.

 

واژه‌های کلیدی: ارزش در معرض ریسک، مدل GED-GARCH، اثر سرریز ریسک فراسوی و فروسوی، آزمون کوپیک.

طبقه‌بندی JEL: Q40، G15، C13.

 

1. مقدمه

پیش از دهه 80 میلادی، ساختار ساماندهی و شرایط حاکم بر بازار بین‌المللی نفت به گونه‌ای بود که قیمت‌های نفت خام نوسانات چندانی نداشت و بنابراین ریسک چندانی را متوجه فعالین در بازیگران بازار نفت نمی‌کرد. پس از آن که به علت افزایش نوسانات قیمت در دهه 1980 میلادی قیمت‌گذاری نفت به بازار فیزیکی واگذار شد زمینه برای تشکیل بازار کاغذی[1] نفت ایجاد گردید.

از اواخر دهة80 میلادی، پس از راه‌اندازی بازار بورس نایمکس[2] در نیویورک و تعمیق تدریجی این بازار در دهه بعدی کشف قیمت‌ نفت‌های شاخص جهان (مانند (Brent , WTI متأثر از شرایط عرضه و تقاضا در بازارهای کاغذی آتی‌ها و نیز نقدی[3]  گردید. با توسعه و گسترش بازار مذکور و مالی شدن این بازار[4] اثرگذاری عوامل غیربنیادین دیگر، قیمت جهانی نفت، با گستره زیادی متلاطم گردید.[5] قیمت‌های جهانی نفت تغییرپذیری[6] زیادی را طی دوره‌های مختلف تجربه کرده که به عنوان نمونه می‌توان به تغییرپذیری شدید قیمت‌ها در سال 2008 اشاره نمود که طی آن، قیمت نفت پس از ثبت رکورد 148 دلار در هر بشکه در ماه جولای، به حدود 40 دلار در هر بشکه در اواخر دسامبر تنزل پیدا کرد. ریسک‌های بازاری فراوان ایجاد شده بر اثر نوسانات شدید قیمتی‌، فعالین در صنعت و بازار بین‌المللی نفت را در معرض زیان‌های بالقوه سنگین قرار داده است.

ارزش در معرض ریسک (VaR) یکی از معیارهای معتبر و رایج اندازه‌گیری ریسک بازار است. ارزش در معرض ریسک به عنوان «برآوردی از بیشترین زیانی که احتمالاً یک سبد دارایی در طول کل دوره متحمل می‌شود، یا به صورت دقیق‌تر، حداکثر زیانی است که یک بنگاه می‌تواند مطمئن باشد در کسر معینی از زمان طی یک دوره مشخص خواهد داشت» تعریف می‌شود[7]. استفاده از روش ارزش در معرض ریسک برای اندازه‌گیری ریسک بازار، مزایای متعددی به همراه دارد. این روش، ریسک بازار را با استفاده از توزیع احتمال یک متغیر تصادفی و ارزیابی ریسک را با یک عدد حقیقی منفرد، توصیف می‌کند. لذا ارزش در معرض ریسک تبدیل به یک ابزار ضروری در بازارهای مالی شده است. به طور مثال، اگر نوسانات قیمت نفت در بازار آتی‌ها بالا باشد و این نوسانات به بازارهای نقدی  منتقل شود تولید کننده و یا تاجر نفت خام می‌تواند از این شیوه برای سنجش ریسک موجودی نفت خام خود استفاده کند. استفاده از این روش برای سنجش ریسک سرمایه در بازارهای نفت کاربرد دارد و مقالات متعددی از این شیوه برای سنجش ریسک استفاده کرده‌اند.

کبدو و مویا[8] (2003) از سه روش برای محاسبه ارزش در معرض ریسک برای قیمت نفت استفاده نمودند و نشان دادند که روش HSAF[9] بهترین کارایی را دارد. صادقی و شوالپور (2006) نیز ارزش در معرض ریسک را در بازار نفت برای سبد نفتی اوپک بررسی و نشان داده‌اند که روش HSAF روشی کاراتر نسبت به روش گارچ معمولی می‌باشد ولی توزیع چاق در بازارهای مالی را در نظر نگرفته‌اند در حالیکه شواهد زیادی دلالت بر این نکته دارند که توزیع بازده دارایی‌های مالی دارای کشیدگی و توزیع با دنباله چاق هستند و عدم توجه به این نکته باعث می‌شود که برآوردهای ارزش در معرض ریسک، بیشتر یا کمتر از مقدار واقعی برآورد شوند. در برخی مطالعات استفاده از توزیع t برای داده‌های با دنباله توزیع چاق پیشنهاد شده است مانند انجلیدیس، بنوس و دجیناکی[10] (2004) و هانگ و همکاران[11] (2007) ؛ مطالعات بعدی روش [12]GED-GARCH را مبنای کار برآورد ارزش در معرض ریسک دانسته‌اند.

پیگیلدین[13] (2009) ارزش در معرض ریسک را برای نفت خام WTI و گاز طبیعی هنری هاب محاسبه و نشان داد که با توجه به دنباله چاق بودن توزیع، روش GED-GARCH بهترین نتایج را دارد. وی همچنین برای بررسی اثر سرریز ریسک از آزمون علیت گرنجر استفاده و نشان داده که اثر سرریز ریسک از بازار نفت خام به گاز وجود دارد ولی اثر سرریز بالعکسی وجود ندارد. فان و همکاران (2008) برای محاسبه ارزش در معرض ریسک بازدهی نفت خام نقدی  WTI و برنت از رویکرد  GED-GARCHو برای اثر سرریز ریسک از رویکرد آزمون علیت گرنجر بر مبنای توزیع‌های کرنل[14] استفاده نمودند. نتایج مطالعه آن‌ها اثر سرریز ریسک دو طرفه میان بازارهای فوق را تأیید می‌نماید. بررسی اثرسرریز ریسک با مطالعات هونگ و همکاران (2001،2003، 2007) آغاز گردید. هونگ[15] (2001) علیت گرنجری ریسک بین دو سری زمانی را معرفی، هونگ و همکاران[16] (2003) اثرسرریز ریسک را بین بازار سهام چین و بازارهای سهام بین‌المللی با استفاده از علیت گرنجری و براساس توابع دنیل کرنل بررسی کردند. هونگ و همکاران[17] (2007) در مطالعه دیگری اثر سرریز ریسک را بین بازارهای مالی با استفاده از سه روش آزمون علیت گرنجری، گرنجر براساس توابع کرنل و دنیل کرنل بررسی نمودند.

در مطالعات داخلی، کاظمی (1385) ارزش در معرض ریسک را برای ریسک بانکی، پیکارجو و حسین‌پور (1389) نیز برای شرکت‌های بیمه آنرا محاسبه نمودند. خیابانی و ساروقی (1390) با استفاده از شبیه‌سازی مونت‌کارلو برای برآورد ارزش در معرض ریسک در بازه خارج نمونه، نشان دادند بازده روزانه در بورس اوراق بهادار تهران دارای توزیع با دنباله چاق می‌باشد که در نظر داشتن توزیع t در برآورد ارزش در معرض ریسک، کاراترین نتیجه را در پی دارد. شاهمرادی و زنگنه (1386) نیز در برآورد ارزش در معرض ریسک در بورس اوراق بهادار تهران با استفاده از مدل‌های گارچ، توزیع با دنباله چاق را در نظر داشته و توزیع t را برای محاسبات انتخاب نمودند.

اغلب مطالعات صورت گرفته در ادبیات موضوع، مبتنی بر این فرض می‌باشد که سری‌های اجزاء اخلال دارای توزیع نرمال و یا توزیع t است، در حالی که بازدهی‌های قیمت نفت اغلب دارای توزیع کشیده با دنباله چاق می‌باشد که از فرض مذکور کاملاً متفاوت بوده و در نتیجه اغلب مدل‌های ارزش در معرض خطر برآورد شده، نامناسب و دارای اریب بوده و قضاوت‌های نادرستی را در رابطه با وضعیت ریسک بازار نفت موجب شده است.[18]

در این مطالعه، ابتدا مدل نوسانات بازدهی بازار نقدی  و آتی‌ها با استفاده از مدل‌های خانواده گارچ تعیین گردید و سپس ریسک‌های فراسوی[19] و فروسوی[20] بازدهی قیمت نفت خام WTI با استفاده از مدل گارچ بر مبنای توزیع خطای تعمیم‌یافته (GED) که برای توزیع‌های کشیده با دنباله چاق مدل مناسبی است، برآورد شده است در ادامه برای ارزیابی صحت مدل‌های ارزش در معرض ریسک فراسوی و فروسوی بکار گرفته شده، از آزمون کوپیک[21] (1995) استفاده شده است. در نهایت نیز با استفاده از روش آزمون علیت گرنجر (هونگ و همکاران، 2003) به بررسی اثر سرریز ریسک میان بازدهی قیمت آتی‌ها و نقدی  بازار نفت WTI پرداخته شده است.

 

2. متدولوژی

برای محاسبه ارزش در معرض ریسک به طور عمده از سه روش تئوری ارزش حدی، شبیه‌سازی تصادفی و روش واریانس- کوواریانس استفاده می‌شود. در این مقاله ، روش واریانس- کوواریانس و به طور خاص، روش GED-GARCH از خانواده مدل‌های GARCH مورد استفاده قرار گرفته است. توزیع جزء اخلال تعمیم‌یافته (GED) در این روش، امکان توصیف بهتر توزیع‌های کشیده با دنباله چاق را که در اغلب بازارهای مالی مصداق دارد، فراهم می‌کند. در شکل 1 انواع روش‌های محاسبه ارزش در معرض ریسکنشان داده شده است (پیگیلدین، 2009).

 

 

شکل 1. روش‌های محاسبه ارزش در معرض ریسک

2-1. مدل VaR بر مبنای GED-GARCH

مطابق با کارهای هندریک (1996)[22] و هیلتون (2003)[23]، ارزش در معرض ریسک بیانگر مقدار حداکثر پولی است که در یک پروتفولیو طی دوره زمانی معین و سطح اطمینان مشخص با توجه به تابع توزیع از دست می‌رود. به لحاظ آماری ارزش در معرض ریسک بیانگر چارک تابع توزیع است. در این مقاله، چارک چپ بازدهی‌های بین‌المللی قیمت نفت برای اندازه‌گیری ریسک فروسوی که نشانگر زیان فروش تولیدکنندگان نفت خام به دلیل کاهش بازدهی‌های نفت و چارک راست برای توصیف ریسک فراسوی که نشانگر زیان خریداران نفت به دلیل افزایش سریع در قیمت نفت استفاده شده است.

از آن‌جا که بازده قیمت نفت اغلب دارای توزیع کشیده[24] بوده که توزیع نرمال استاندارد ندارد (فرض توزیع نرمال استاندارد منجر به برآورد کمتر از حد ریسک حدی می‌شود)، توزیع خطای تعمیم یافته (GED) برای محاسبه سری‌های باقیمانده مدل‌های GARCH مورد استفاده قرار می‌گیرد.

مدل‌های استاندارد ARCH توسط انگل[25] (1982) توسعه داده است. همچنین هنگامی وقفه مدل‌های ARCH بسیار بزرگ می‌گردد، بولرسلیو[26] (1986) مدل ARCH تعمیم‌یافته را که به مدل GARCH شهرت دارد پیشنهاد می‌کند. مدل GARCH بصورت زیر نمایش داده می‌شود:

(1)                                                               

(1.1)                                      

در رابطه (1)  بیانگر بازدهی قیمت نفت،  برداری ستونی شامل متغیرهای مستقل توضیحی، و β بردار ستونی ضرایب است. همچنین‌ محدودیت‌های همانند ، ، ، ، ،  و  و  که منعکس‌کننده مدت تلاطم بازده است.

غالباً اثر اهرمی قابل توجهی در تلاطم بازدهی‌های قیمت نفت مشاهده می‌شود که با استفاده از رابطه (1.1) نمی‌توان آن را استخراج نمود. به عبارت دیگر تلاطم بازدهی‌های قیمت نفت می‌تواند نامتقارن باشد، یعنی اخبار خوب و اخبار بد بر نوسانات بازار می‌تواند اثر نامتقارن داشته باشد که در اینصورت مدل TGARCH برای بررسی این موضوع بکار می‌رود و واریانس شرطی آن می‌تواند به صورت زیر نشان داده شود:

(2)                             

چنانچه  باشد  و در غیر اینصورت  خواهد بود. با توجه به استفاده از  در رابطه فوق، اثر افزایش بازده () و کاهش () در واریانس شرطی () قابل تفکیک شدن است.

همانطور که بیان شد، توزیع بازدهی در بازار نفت، توزیعی دارای کشیدگی و دنباله‌های چاق می‌باشد که توزیع نرمال استاندارد قادر به توضیح رفتار آن نمی‌باشد. برای رفع این مشکل از توزیع GED در مدل‌های GARCH استفاده می‌گردد. تابع چگالی احتمال GED بصورت زیر می‌باشد:

(3)                                

که در آن  و  تابع گاما است و k پارامتر GED و همچنین درجه آزادی آن نیز نامیده می‌شود. به طور خاص، زمانی که  است، GED توزیع نرمال استاندارد را نشان می‌دهد؛ برای  دنباله نازک‌تر از توزیع نرمال استاندارد و برای  دنباله پهن‌تر و چاق‌تر از توزیع نرمال استاندارد خواهد بود. براساس تعریف ارزش در معرض ریسک، ریسک فراسوی و فروسوی بصورت زیر تعریف می‌گردد:

(4)                                                  

(5)                                           

که در آن  بازده انتظاری شرطی بازار m؛  بیانگر چارک αام چپ توزیع GED که در سری‌های باقیمانده مدل GARCH در بازار m استفاده می‌گردد و  نیز سری‌های واریانس شرطی در بازار m است.

پیش از محاسبه ارزش در معرض ریسک، ضروری است تا آزمون اطمینان[27] برای این‌که آیا مدل VaR بکار گرفته شده به میزان کافی ریسک واقعی حدی را نشان می‌دهد یا خیر، انجام شود. برای این منظور، آزمون اطمینان که توسط کوپیک[28] (1995) معرفی شده، بکار گرفته می‌شود. ایده اصلی بدین صورت است که چنانچه فرض کنیم سطح اطمینان α-1، اندازه نمونه T، تعداد شکست‌ها N و تناوب شکست به صورت  در نظر گرفته شود، آنگاه کوپیک آزمون معروف نسبت راستنمایی[29] با فرض صفر  را با آماره زیر معرفی می‌کند:

(6)                    

تحت فرض صفر،  بوده و مقادیر بحرانی سطوح اطمینان 95 درصد و 99 درصد آن برابر با 84/3 و 64/6 است. بر اساس تعریف توزیع ، چنانچه مقدار LR بزرگ‌تر از مقدار بحرانی باشد، آنگاه فرض صفر می‌بایست رد شود. به عبارت دیگر می‌توان گفت که مدل VaR ناکافی است.

 

2-2. آزمون اثر سرریز ریسک[30]

سرریز ریسک به معنای انتقال واریانس‌ها و کوواریانس‌های بازدهی از یک دارایی معین به دارایی معین دیگر تعریف می‌شود (بوناتو و همکاران، 2011). متدلوژی استفاده از علیت گرنجر در ریسک بوسیله هونگ[31] (2001) و هونگ و همکاران[32] (2003) ارائه گردید که می‌توان برای آزمون اثر سرریز ریسک میان بازارهای مختلف نفت از آن بهره برد که این مهم نیازمند محاسبه تغییر زمانی VaR برای هر بازده است. تعیین این‌که اطلاعات تاریخی درباره ریسک در یک بازار توانایی پیش‌بینی وقوع آن را در بازار دیگر افزایش می‌دهد از طریق مفهوم علیت گرنجر فراهم شده است.

 

2-2-1. آزمون علیت گرنجری ریسک فروسوی

برای آزمون اثر یک طرفه ریسک فروسوی از بازار 2 به بازار 1، فرض صفر و یک می‌تواند به صورت زیر بیان شود:

(7)                     

تابع مشخصه ریسک فروسوی بصورت مقابل تعریف می‌گردد:

(8)                                             

در این تابع،  تابع مشخصه‌ای است که اگر زیان واقعی، بیشتر از ارزش در معرض ریسک باشد تابع  مقدار یک را خواهد گرفت و در غیر اینصورت مقدار صفر را خواهد گرفت و  نیز در مطالعه فوق بازار نقدی  و آتی‌های WTI را نشان می‌دهد.

حال با توجه به تعریف تابع فوق، می‌توان فرض صفر و مقابل را بصورت زیر و براساس تابع مشخصه تعریف نمود:

(9.1)                                       

(9.2)                                       

که در آن  بیانگر مجموعه اطلاعات موجود در زمان  و  مجموعه اطلاعات بازار 1 است. چنانچه فرض صفر تأیید شود، هیچ علیت گرنجری یک طرفه در ریسک از بازار 2 به بازار 1 وجود ندارد، آن‌گاه هنگامی که ریسک حدی در بازار 2 مشاهده شد، نمی‌تواند برای پیش‌بینی ریسک محتمل در بازار 1 در آینده مورد استفاده قرار گیرد. به طریق مشابه آزمون اثر سرریز ریسک فراسوی نیز تعریف می‌گردد.

 

2-2-2. آزمون علیت گرنجری ریسک فراسوی

برای آزمون اثر یک طرفه ریسک فراسوی از بازار 2 به بازار 1، فرض صفر و یک می‌تواند به صورت زیر بیان شود:

(10.1)                        

(10.2)                        

تابع مشخصه برای ریسک فراسوی بصورت مقابل می‌باشد:

(11)                                             

در این حالت اگر بازدهی واقعی بیشتر از VaR باشد، تابع  مقدار یک و در غیر اینصورت مقدار صفر را خواهد گرفت. همانند قبل، می‌توان فرض صفر و مقابل را براساس تابع مشخصه تعریف نمود: (هونگ و همکاران (2003)، لیو و همکاران (2008))

(12.1)                                     

(12.2)                                     

فرض می‌شود که  سری‌های VaR بازار m در سطح معنی‌داری α از معادلات (4) و (5) بدست می‌آید. چنانچه ، آن‌گاه تابع همبستگی متقاطع نمونه[33] (CCF) میان  و  می‌تواند از روابط زیر حاصل شود:

(13)

که در آن ،  وT اندازه نمونه سری‌های بازده  و همبستگی متقاطع نمونه به صورت زیر است:

(14)                                  

که  واریانس نمونه  است.

هونگ و دیگران (2003)، دسته‌ای از آزمون‌های مبتنی بر آزمون کرنل را برای بررسی علیت گرنجر در ریسک پیشنهاد دادند. جهت آزمون نمودن وجود علیت یک طرفه در ریسک از بازار 2 به بازار 1، از آماره زیر می‌توان استفاده نمود.

(15)                        

 یک تابع کرنل است که وزن‌ها را به وقفه‌های مختلف با توجه به این‌که M درجه وقفه و فاکتور مرکزی و فاکتور مقیاس به صورت زیر است، اختصاص می‌دهد:

 

(16)                      

 

3. برآورد مدل

3-1. محاسبه ارزش در معرض ریسک

جهت بررسی اثر سرریز ریسک بین بازدهی در بازارهای نقدی  و آتی‌های نفت خام، از داده‌های روزانه نقدی  و آتی‌های یک ماهه بازار نفت نایمکس برای نفت شاخص  WTIاستفاده شده است. داده‌های مورد نظر از اداره اطلاعات انرژی[34] استخراج و برای تبدیل سری قیمت‌ها به بازدهی نیز از رابطه (1) استفاده شده است.

(17)                                                                            

در رابطه فوق،  و  قیمت نفت‌خام در دوره t و 1-t و  نیز بازدهی قیمت را در دوره t نشان می‌دهد. بازه زمانی دادها از 2 ژانویه 1986 تا 24 جولای 2012 می‌باشد. داده‌های درون نمونه[35] از 2 ژانویه 1986 تا 31 دسامبر 2010 و داده‌های خارج از نمونه[36] نیز از 3 ژانویه 2011 تا 24 جولای 2012 انتخاب شده است.

ویژگی‌های توصیفی مربوط به سری بازدهی قیمت آتی‌ها و نقدی  نفت‌خام در جدول 1 نشان داده شده است. با توجه به داده‌های جدول، میانگین بازدهی‌ دو بازار بسیار به هم نزدیک بوده، به عبارتی بازدهی دو بازار تقریباً یکسان است. انحراف معیار نشان می‌دهد که نوسانات بازدهی در بازار نقدی  اندکی بیشتر از قرار داد آتی‌های  WTIاست. برای توزیع نرمال چولگی[37] و کشیدگی[38] به ترتیب صفر و 3 خواهد بود و اگر کشیدگی بیشتر از 3 باشد توزیع کشیده با دنباله چاق[39] خواهد بود (فان و همکاران، 2008). همانطور که اطلاعات جدول فوق نشان می‌دهد، عدد کشیدگی هم برای بازار نقدی و هم برای بازار آتی‌ها بیشتر از 3 بوده که دنباله چاق بودن توزیع بازدهی بازار نقدی و آتی‌ها را تأیید می‌نماید. آماره جارک برا[40] که برای بررسی نرمال بودن توزیع‌ها بکار می‌رود معنی‌دار شده است که نرمال استاندارد نبودن توزیع را تأیید می‌نماید. آزمون دیکی فولر تعمیم یافته نیز مانایی متغیرها را تأیید می‌نماید.

 

جدول 1. آماره‌های توصیفی و آزمون‌های مربوطه برای بازدهی نفت خام نقدی و آتی‌های WTI

بازدهی

میانگین

انحراف معیار

چولگی

کشیدگی

آزمون جارگ برا (JB)

آزمون ADF

نقدی

000185/0

0258/0

756/0-

462/17

(000/0) 65/59033

(00/0) 7742/60-

آتی‌ها

000186/0

0252/0

797/0-

236/17

(000/0) 12/56963

(00/0) 3349/61-

مأخذ: یافته‌های تحقیق

 

جهت استفاده از مدل‌های خانواده گارچ، آزمون  ARCH-LMانجام شده است. آماره آزمون فوق برای بازدهی در بازار نقدی  و آتی‌های نفت خام WTI به ترتیب (00/0) 74/78 و (00/0) 39/116 بوده که وجود ناهمسانی شرطی در باقیمانده‌ها را تأیید می‌نماید. با در نظر داشتن معیار آکائیک (AIC) که باید حداقل گردد، برای بازدهی در بازار نقدی  فرایند ARMA(1,1)-GARCH(1,1) و برای بازدهی در بازار آتی‌ها  AR(2)-GARCH(1,1)انتخاب و نتایج آن در زیر گزارش گردیده است.

 

 

 

 

جدول 2. مدل گارچ برای بازدهی در بازار نقدی نفت خام WTI - ARMA(1,1)-GARCH(1,1)

پارامتر

ضریب

انحراف معیار

آماره z

معادله میانگین :

 

 

 

C

00057/0

00019/0

(003/0)977/2

AR(1)

939/0

0241/0

(00/0)91/38

MA(1)

947/0-

0223/0

(00/0)47/42-

معادله واریانس :

 

 

 

C

000006/0

000001/0

(00/0)71/4

 

077/0

0065/0

(00/0)82/11

 

914/0

007/0

(00/0) 130

GED PARAMETER

337/1

024/0

(00/0)63/55

AIC: -4.8091, Log likelihood: 15170.1

 

جدول 3. مدل گارچ برای بازدهی در بازار آتی‌های نفت خام WTI - AR(2)-GARCH(1,1)

پارامتر

ضریب

انحراف معیار

آماره z

معادله میانگین :

 

 

 

C

00052/0

00021/0

(015/0)415/2

AR(2)

028/0-

0123/0

(02/0)273/2-

معادله واریانس :

 

 

 

C

000005/0

0000012/0

(00/0)666/4

 

0723/0

0064/0

(00/0)239/11

 

9194/0

007/0

(00/0)46/129

GED PARAMETER

37/1

027/0

(00/0)96/49

AIC: -4.8462, Log likelihood: 15191.6

 

همانطور که از نتایج مدل نوسانات بازدهی قیمت نقدی  و آتی‌ها (جداول 2 و 3) مشخص است پارامتر GED برای بازدهی بازار نقدی  و آتی‌ها به ترتیب 33/1 و 37/1 بدست آمده است. اگر پارامتر GED کمتر از 2 باشد توزیع با دنباله چاق خواهد بود (فان و همکاران، 2008) لذا معنی‌دار بودن و کوچکتر از 2 بودن پارامترهای محاسبه شده نشان می‌دهد توزیع بازده قیمت نفت خام نقدی و قرارداد آتی‌ها در بازار نایمکس دارای توزیع با دنباله‌های چاق می‌باشد که نتایج حاصل از پارامترهای چولگی و کشیدگی را نیز تأیید می‌نماید.

با در نظر داشتن دنباله چاق بودن توزیع بازدهی در بازار نقدی و آتی‌های نفت خام، ارزش در معرض ریسک فراسوی[41] و فروسوی[42] به واسطه استفاده از داده های درون نمونه برای پیش بینی داده های خارج از نمونه محاسبه گردیده است. ریسک فراسوی و فروسوی در واقع نتیجه ای است که از پیش بینی عددی داده های درون نمونه ای بدست آمده است. برای این منظور از رابطه (4) و (5) استفاده شده است. به عبارتی پس از محاسبه انحراف معیارهای پیش بینی و با استفاده از داده های خارج از نمونه، همانطور که در رابطه (4) و (5) توضیح داده شد، ریسک فراسوی و فروسوی محاسبه گردید. شکل (2) و (3) ارزش در معرض ریسک فراسوی و فروسوی را برای بازدهی بازار نقدی و شکل (4) و (5) نیز بازدهی قرارداد آتی‌نفت خام در سطح اطمینان 95 و 99 درصد نشان می‌دهد.

 

 

شکل 2. ارزش در معرض ریسک (VaR) برای بازدهی بازار نقدی نفت خام WTI در سطح 95 درصد

 

شکل 3. ارزش در معرض ریسک (VaR) برای بازدهی بازار نقدی نفت خام WTI در سطح 99 درصد

 

شکل 4. ارزش درمعرض ریسک (VaR) برای بازدهی قرارداد آتی‌های‌ یک ماهه در سطح 95 درصد

 

شکل 5. ارزش در معرض ریسک (VaR) برای بازدهی بازار قرارداد آتی یک ماهه در سطح 99 درصد

 

اگر بازه سری‌های استخراج شده ارزش در معرض ریسک، کوچک باشد، در اینصورت تعداد خطا (تعداد دفعاتی که قیمت واقعی از بازه فراسوی و فروسوی مربوطه خارج شده) افزایش می‌یابد و اگر بازه سری‌های فوق بزرگ باشد، در این صورت محاسبات، کارایی لازم را نخواهد داشت. لذا سری‌های ارزش در معرض ریسک فراسوی و فروسوی استخراج شده مناسب، بازه‌ای است که هم تعداد خطای پایینی داشته باشد و هم کارایی بالایی داشته باشد. لذا نیاز است تا آزمون اطمینان کوپیک برای بررسی صحت سری‌های ارزش در معرض ریسک استخراج شده صورت پذیرد.

پس از استخراج سری‌های ارزش در معرض ریسک فراسوی و فروسوی در سطح 95 و 99 درصد برای بازدهی نقدی  و آتی‌ها که در شکل (3) تا (6) صورت پذیرفت، تعداد خطا، نسبت خطا ( نسبت تعداد خطا به کل نمونه) و نهایتاً آماره LR (آماره مربوطه در معادله (6)) محاسبه شد که نتایج آن درجدول (4) و (5) بیان گردیده است.

سطح معنی‌داری برای آماره LR در سطح 95 درصد و 99 درصد به ترتیب 841/3 و 635/6 می‌باشد. فرض صفر در این آزمون مناسب بودن سری‌های ارزش در معرض ریسک استخراج شده است و فرض مقابل نامناسب بودن آن است. نتایج جدول نشان می‌دهد که آماره‌های LR محاسبه شده از مقادیر بحرانی پایین‌تر بوده، لذا فرض صفر برای تمامی VaRهای محاسبه شده پذیرفته می‌باشد. تنها برای ارزش در معرض ریسک فراسوی برای بازار نقدی  و آتی‌ها و در سطح اطمینان 95 درصد، مقدار آماره فوق بالاتر از مقدار بحرانی بوده که این امر در مطالعه فان و همکاران (2008) که بر روی نفت‌های شاخص WTI و برنت انجام شده نیز دیده می‌شود که دلیل آن افزایش قیمت نفت از سال 2000 تا 2008 بوده که سری ارزش در معرض ریسک فراسوی فوق بدلیل جهش بسیار زیاد قیمت واقعی در این سال‌ها نتوانسته آنرا شامل شود.

 

جدول 4. آزمون کوپیک برای بررسی صحت مدل‌ VaR محاسبه شده در بازار نقدی  WTI

 

سطح اطمینان

میانگین

انحراف معیار

تعداد خطا

درصد خطا

آمارهLR

فراسوی

95%

034/0

0079/0

10

025/0

986/5

99%

0525/0

0123/0

3

007/0

237/0

فروسوی

95%

034/0-

0079/0

21

053/0

102/0

99%

0525/0-

0123/0

6

015/0

959/0

مأخذ: یافته‌های تحقیق

 

جدول 5. آزمون کوپیک برای بررسی صحت مدل VaR محاسبه شده در بازار آتی‌های WTI

 

سطح اطمینان

میانگین

انحراف معیار

تعداد خطا

درصد خطا

آمارهLR

فراسوی

95%

033/0

008/0

11

028/0

689/4

99%

052/0

012/0

7

017/0

981/1

فروسوی

95%

033/0-

008/0

20

051/0

008/0

99%

052/0-

012/0

3

007/0

237/0

مأخذ: یافته‌های تحقیق

 

 

3-2. اثرسرریز ریسک بین بازار نقدی  و قراردادهای آتی برای نفت WTI

پس ازحصول اطمینان از صحت ارزش‌های در معرض ریسک فراسوی و فروسوی محاسبه شده، می‌توان اثر سرریز را در بین سری‌ها با استفاده از آزمون گرنجر مورد بررسی قرار داد لذا آزمون علیت گرنجری بر روی سری‌ Z های محاسبه شده (معادلات (8) و (11)) صورت پذیرفته است که نتایج آن در جدول (4) نشان داده شده است. (هونگ و همکاران (2003))

نتایج جدول (6) نشان می‌دهد که اثر سرریز ریسک فراسوی در سطح اطمینان 95 درصد از آتی‌ها به نقدی  وجود ندارد به عبارتی هیچگونه ریسک فراسویی از بازار آتی‌ها با نقدی  در سطح اطمینان 95 درصد منتقل نمی‌شود. اثر سرریز ریسک فراسوی از نقدی  به آتی‌ها در تمامی وقفه‌ها جز در وقفه 20 بی‌معنی بوده است و به عبارتی ریسک فراسویی از بازار نقدی  به آتی‌ها در سطح اطمینان 95 درصد انتقال نمی‌یابد. معنی‌داری تنها در وقفه بیستم وجود دارد که نشان می‌دهد پس از زمانی بسیار طولانی ریسک از بازار نقدی  به آتی‌ها منتقل خواهد شد که با توجه به عدم صحت ارزش در معرض ریسک محاسبه شده برای آن می‌توان از آن چشم پوشی نمود.

نتایج همچنین حاکی از آن است که اثر سرریز ریسک فراسوی شدیدی در سطح اطمینان 99 درصد از آتی‌ها به نقدی  وجود دارد. آماره مورد نظر در تمامی وقفه‌ها برای اثر سرریز ریسک از بازار آتی‌ها به نقدی  قویاً معنی‌دار بوده است و این امر نشان می‌دهد در زمانی‌که قیمت‌ها برای قراردادهای آتی‌های WTI افزایشی است، ریسکی برای معامله‌گران در هر دو بازار نقدی  و آتی‌ها ایجاد می‌شود و بازار نقدی  علاوه بر ریسکی که در خود دارد تأثیرپذیری شدیدی نیز از ریسک بازار آتی‌ها خواهد داشت. نتایج همچنین نشان داده که در سطح اطمینان 99 درصد اثر سرریز ریسک فراسویی از بازار نقدی  به آتی‌ها وجود ندارد.

نتایج همچنین نشان داده که تمامی ضرایب مربوط به اثر سرریز ریسک فروسوی از بازار نقدی  به آتی‌ها و از بازار آتی‌ها به نقدی  در سطوح 95 و 99 درصد بی‌معنی بوده است لذا هیچگونه ریسک فروسویی بین بازارهای نقدی  و آتی‌ها وجود در سطوح اطمینان 95 و 99 درصد وجود ندارد. این امر نشان می‌دهد زمانی‌که قیمت‌ها در بازار نقدی  و آتی‌های WTI نزولی است، ریسکی را برای معامله‌گران در هر دو بازار ایجاد خواهد نمود ولی ریسک از بازارهای به یکدیگر منتقل نخواهد شد.

 

 

جدول 6. اثر سرریز ریسک بین بازار نقدی و آتی‌های نفت WTI

اثر سرریز

5=M

10=M

15=M

20=M

فراسوی 95 درصد :

 

 

 

 

F→S

NS(99/0) 035/0

NS (606/0) 823/0

NS (314/0) 144/1

NS (346/0) 101/1

S→F

NS (99/0) 098/0

NS (789/0) 628/0

NS (08/0) 568/1

**(027/0) 727/1

فراسوی 99 درصد :

 

 

 

 

F→S

*** (00/0) 08/31

*** (00/0) 269/16

*** (00/0) 429/11

*** (00/0) 23/9

S→F

NS (00/1) 001/0

NS (00/1) 00/0

NS (00/1) 00/0

NS (00/1) 00/0

فروسوی 95 درصد :

 

 

 

 

F→S

NS (99/0) 024/0

NS (00/1) 03/0

NS (00/1) 015/0

NS (874/0) 648/0

S→F

NS (99/0) 022/0

NS (00/1) 02/0

NS (00/1) 045/0

NS (00/1) 103/0

فروسوی 99 درصد :

 

 

 

 

F→S

NS (91/0) 293/0

NS (00/1) 015/0

NS (00/1) 017/0

NS (00/1) 06/0

S→F

NS (99/0) 013/0

NS (99/0) 141/0

NS (00/1) 091/0

NS (00/1) 02/0

M: تعداد وقفه     **: معنی‌داری در سطح 95 درصد         ***: معنی‌داری در سطح 99 درصد         NS: بی‌معنی

مأخذ: یافته‌های تحقیق

 

4. نتیجه‌گیری

یافته‌های مقاله که مربوط به محاسبات ارزش در معرض ریسک و اثر سرریز ریسک است را می‌تواند بصورت زیر خلاصه نمود:

-      نتایج نشان می‌دهد که دو بازار نقدی  و آتی‌های نفت خام WTI، ویژگی‌های نزدیک به هم از نظر بازدهی و ریسک بازدهی دارند. بازدهی در بازار نقدی  و آتی‌ها بسیار نزدیک به یک دیگر است و نوسانات بازدهی در بازار نقدی  اندکی بیشتر از نوسانات در بازار آتی‌های WTI است.

-      توزیع آماری بازده قیمت نفت خام نقدی  و قرار داد آتی یک ماهه‌ در بازار WTI دارای توزیع با دنباله‌های چاق می‌باشد که نتایج حاصل از پارامترهای GED، چولگی و کشیدگی آنرا تأیید می‌نماید. توزیع چاق بازدهی در بازار نقدی  و آتی‌ها نشان می‌دهد که احتمال تغییرات و جهش‌های بزرگ قیمتی در این بازارها وجود دارد.

-      نتایج بررسی اثر سر ریز ریسک فراسوی نشان داده که در سطح اطمینان 99 درصد از بازار آتی‌ها به نقدی  اثر سرریز معنی‌داری وجود دارد. آماره مورد نظر در تمامی وقفه‌ها برای اثر سرریز ریسک از بازار آتی‌ها به نقدی  معنی‌دار بوده است. در سطح 95 درصد نیز از بازار نقدی  به آتی‌ها و تنها در وقفه بیستم اثر سرریز مشاهده شده است که با توجه به آزمون کوپیک می‌توان از این نتیجه چشم پوشی نمود.

-      وجود اثر سرریز ریسک فراسوی از بازار آتی‌ها به نقدی  نشان می‌دهد که در بازار نفت خام WTI، در زمان افزایش قیمت، ریسک از بازار آتی‌ها به نقدی  هدایت می‌شود. به این ترتیب می‌توان از ریسک در بازار آتی‌ها برای پیش‌بینی ریسک در بازار نقدی  در زمان افزایش قیمت‌ها استفاده نمود.

-      نتایج بررسی اثر سر ریز ریسک فروسوی نشان داده که در سطح اطمینان 95 و 99 درصد اثر سرریز ریسک وجود ندارد. لذا نمی‌توان از اطلاعات بازار آتی‌ها برای پیش‌بینی ریسک در بازار نقدی  در زمان کاهش قیمت‌ها استفاده نمود و بالعکس.

-      صحت ارزش در معرض ریسک فراسوی و فروسوی با استفاده از آزمون کوپیک مورد تأیید بوده است و تنها برای ریسک فراسوی در سطح اطمینان 95 درصد برای بازار نقدی  و آتی‌ها صحت آن تأیید نشده که دلیل آن نیز جهش‌های بزرگ قیمتی خصوصاً در دهه 2000 میلادی بوده است و در نتیجه این عدم صحت، نتایج اثر سرریز ریسک فراسوی در سطح 95 درصد نادیده گرفته شده است لذا پیشنهاد می‌شود تا پژوهش‌های آینده سری ارزش در معرض ریسکی را استخراج نمایند که مشکل فوق را نداشته و بتوان به نتایج حاصل از اثر سرریز آن نیز اعتماد نمود.

 

منابع

الف- فارسی

پیکارجو، پ و حسین‌پور، ب (1389). اندازه‌گیری ارزش در معرض ریسک در شرکت‌های بیمه با استفاده از مدل  GARCH. پژوهشنامه بیمه (صنعت بیمه)، شماره 25.

خیابانی، ن و ساروقی، م.، (1390). «ارزش‌گذاری برآورد VaR براساس مدل‌های خانواده ARCH (مطالعه موضوعی برای بازار اوراق بهادار تهران)». فصلنامه پژوهش‌های اقتصادی ایران، سال شانزدهم، شماره 47.

شاهمرادی، الف و زنگنه، م.، (1387). «محاسبه ارزش در معرض ریسک برای شاخص‌های عمده بورس اوراق بهادار تهران با استفاده از روش پارامتریک». مجله تحقیقات اقتصادی، دروه 42، شماره 79.

کاظمی، ر.، (1385). «معرفی ارزش در معرض ریسک و کاربرد آن در اندازه‌گیری ریسک بانک». مجموعه مقالات هفدهمین همایش بانکداری اسلامی، مؤسسه عالی بانکداری ایران.

مشیری، س و مروت، ح.، (1385). «پیش‌بینی شاخص کل بازدهی سهام ایران با استفاده از مدل‌های خطی و غیرخطی». فصلنامه پژوهشنامه بازرگانی، شماره 41.

 

ب- انگلیسی

Angelidis, T., Benos,A., Degiannakis,S., (2004). “The Use of GARCH Models in VaR Estimation”, working paper

Bollerslev, T., (1986). “Generalized Autoregressive Conditional         Heteroskedasticity”. Journal of Econometrics, No 31. 

Bonato, Matteo, Caporin Massimiliano, Ranaldo, Angelo., (2011). “Hedging risk spillovers in international equity Portfolios”, working paper

Cabedo, J.D., Moya, I., (2003). “Estimating Oil Price ‘Value at Risk’ Using the Historical Simulation Approach”. Energy Economics, No 25.

Energy Information Administration, www.eia.gov.

Engle, R.F., (1982). “Autoregressive Conditional Heteroskedasticity with Estimates of the Variance of United Kingdom Inflation”.Econometrica, No 50.

Fan, Y., (2000).”VaR Methodology and its Application in Stock Market Risk Analysis”.Chinese Journal of Management Science, No 8.

Fan, Ying, Yue-Jun Zhang, Hsien- Tang Tsai, and Yi- Ming Wei., (2008). “Estimating Value at Risk of Crude Oil Price and its Spillover Effect Using the GED-GARCH Approach”.Energy Economics, No 30.

Hendricks, D., (1996). “Evaluation of Value at Risk Modeling Using Historical Data”.Economic Policy Review.Federal Reserve Bank of New York.

Hilton, G.A., (2003). Value-at-Risk, Theory and Practice. New York.

Hong, Y.M., (2001). “Granger Causality in Risk and Detection of Risk Transmission between Financial Markets”.Working paper.Department of Economics and Department of Statistical Science, Cornell University.

Hong, Y.M., Cheng, S.W., Liu, Y.H., Wang, S.Y., (2003). “Extreme Risk Spillover between Chinese Stock Markets and International Stock Markets”.Working paper.Department of Economics and Department of Statistical Science, Cornell University.

Hopper, G.P., (1996). “Value at Risk: a New Methodology for Measuring Portfolio Risk”. Federal Reserve Bank of Philadelphia Business Review, July August.

Jalali-Naini, Ahmad R., (2009). The Impact of Financial Markets on the Price of Oil and Volatility: Since 2007. OPEC Bulletin.

Kupiec, P.H., (1995). “Techniques for Verifying the Accuracy of Risk Measurement Models”.Journal of Derivatives, No 3.

Pigildin, D.(2009). Value-at-Risk estimation and extreme risk spillover between oil and natural gas markets. Master thesis, department of economics, Budapest, Hungry

Sadeghi, M., Shavvalpour, S., (2006). “Energy Risk Management and Value at Risk Modeling”.Energy Policy, No 34.

Timotheos, A., Benos, A., and Degiannakis, S. (2004). “The Use of GARCH Models in VaR Estimation”.Statistical Methodology, vol. 1, No 2.

 



* عضو هیات علمی مؤسسه عالی آموزش و پژوهش مدیریت و برنامه‌ریزی                                  ahmad_jalali@hotmail.com

** دانشجوی دکتری اقتصاد نفت و گاز دانشگاه علامه طباطبایی                                                             v.ghorbani@st.atu.ac.ir

*** دانشجوی دکتری اقتصاد نفت و گاز دانشگاه علامه طباطبایی                                                               m.sayadi@st.atu.ac.ir

[1]. Paper Market

[2] . Nymex

[3] . Spot

[4]. Financialization

[5] . برای توضیحات بیشتر رجوع کنید به .Jalali-Naini (2009)

[6]. Volatility

[7]. Hopper (1996)

[8]. Cabedo, J. and Moya, I.

[9]. Historical Simulation with ARMA Forecasts

[10]. Angelidis, Benos and Degiannakis

[11]. Hung,J., Lee, M., Liu, H.

[12] . Generalized Error Distribution GARCH

[13] . Pigildin, D.

[14] . kernel Distributions

[15]. Hong,Y.

[16]. Hong, Y., Cheng, S. Liu, Y.,Wang, S.

[17]. Hong, Y., Liu, Y. ,Wang, S

[18]. Fan, et al (2008)

[19]. Up Side Risk

[20]. Down Side Risk

[21]. Kupiec, (1995)

[22]. Hendricks (1996)

[23]. Hilton (2003)

[24]. Leptokurtic and Fat Tail Distribution

[25]. Engle

[26]. Bollerslev

[27]. Backtest

[28]. Kupiec

[29]. Likelihood Ratio

[30]. Risk Spillover Effect

[31]. Hong

[32]. Hong and et al

[33]. Sample Cross-Correlation Function

[34]. Energy Information Administration, www.eia.gov

[35]. in-Sample Period

[36]. Out-of-Sample Period

[37]. Skewness

[38]. kurtosis

[39]. Fat tail

[40]. JarqueBera (JB)

[41]. Upside Value at Risk

[42]. Downside Value at Risk

الف- فارسی
پیکارجو، پ و حسین‌پور، ب (1389). اندازه‌گیری ارزش در معرض ریسک در شرکت‌های بیمه با استفاده از مدل  GARCH. پژوهشنامه بیمه (صنعت بیمه)، شماره 25.
خیابانی، ن و ساروقی، م.، (1390). «ارزش‌گذاری برآورد VaR براساس مدل‌های خانواده ARCH (مطالعه موضوعی برای بازار اوراق بهادار تهران)». فصلنامه پژوهش‌های اقتصادی ایران، سال شانزدهم، شماره 47.
شاهمرادی، الف و زنگنه، م.، (1387). «محاسبه ارزش در معرض ریسک برای شاخص‌های عمده بورس اوراق بهادار تهران با استفاده از روش پارامتریک». مجله تحقیقات اقتصادی، دروه 42، شماره 79.
کاظمی، ر.، (1385). «معرفی ارزش در معرض ریسک و کاربرد آن در اندازه‌گیری ریسک بانک». مجموعه مقالات هفدهمین همایش بانکداری اسلامی، مؤسسه عالی بانکداری ایران.
مشیری، س و مروت، ح.، (1385). «پیش‌بینی شاخص کل بازدهی سهام ایران با استفاده از مدل‌های خطی و غیرخطی». فصلنامه پژوهشنامه بازرگانی، شماره 41.
 
ب- انگلیسی
Angelidis, T., Benos,A., Degiannakis,S., (2004). “The Use of GARCH Models in VaR Estimation”, working paper
Bollerslev, T., (1986). “Generalized Autoregressive Conditional         Heteroskedasticity”. Journal of Econometrics, No 31. 
Bonato, Matteo, Caporin Massimiliano, Ranaldo, Angelo., (2011). “Hedging risk spillovers in international equity Portfolios”, working paper
Cabedo, J.D., Moya, I., (2003). “Estimating Oil Price ‘Value at Risk’ Using the Historical Simulation Approach”. Energy Economics, No 25.
Energy Information Administration, www.eia.gov.
Engle, R.F., (1982). “Autoregressive Conditional Heteroskedasticity with Estimates of the Variance of United Kingdom Inflation”.Econometrica, No 50.
Fan, Y., (2000).”VaR Methodology and its Application in Stock Market Risk Analysis”.Chinese Journal of Management Science, No 8.
Fan, Ying, Yue-Jun Zhang, Hsien- Tang Tsai, and Yi- Ming Wei., (2008). “Estimating Value at Risk of Crude Oil Price and its Spillover Effect Using the GED-GARCH Approach”.Energy Economics, No 30.
Hendricks, D., (1996). “Evaluation of Value at Risk Modeling Using Historical Data”.Economic Policy Review.Federal Reserve Bank of New York.
Hilton, G.A., (2003). Value-at-Risk, Theory and Practice. New York.
Hong, Y.M., (2001). “Granger Causality in Risk and Detection of Risk Transmission between Financial Markets”.Working paper.Department of Economics and Department of Statistical Science, Cornell University.
Hong, Y.M., Cheng, S.W., Liu, Y.H., Wang, S.Y., (2003). “Extreme Risk Spillover between Chinese Stock Markets and International Stock Markets”.Working paper.Department of Economics and Department of Statistical Science, Cornell University.
Hopper, G.P., (1996). “Value at Risk: a New Methodology for Measuring Portfolio Risk”. Federal Reserve Bank of Philadelphia Business Review, July August.
Jalali-Naini, Ahmad R., (2009). The Impact of Financial Markets on the Price of Oil and Volatility: Since 2007. OPEC Bulletin.
Kupiec, P.H., (1995). “Techniques for Verifying the Accuracy of Risk Measurement Models”.Journal of Derivatives, No 3.
Pigildin, D.(2009). Value-at-Risk estimation and extreme risk spillover between oil and natural gas markets. Master thesis, department of economics, Budapest, Hungry
Sadeghi, M., Shavvalpour, S., (2006). “Energy Risk Management and Value at Risk Modeling”.Energy Policy, No 34.
Timotheos, A., Benos, A., and Degiannakis, S. (2004). “The Use of GARCH Models in VaR Estimation”.Statistical Methodology, vol. 1, No 2.