Next Article in Journal
Rising Tides: Election Cycles, Economic Uncertainty, Equity and Commodity Markets Fluctuations
Previous Article in Journal
Does “Paper Oil” Matter? Energy Markets’ Financialization and Co-Movements with Equity Markets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

What Insights Do Short-Maturity (7DTE) Return Predictive Regressions Offer about Risk Preferences in the Oil Market?

Fox School of Business, Temple University, Philadelphia, PA 19122, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Commodities 2024, 3(2), 225-247; https://doi.org/10.3390/commodities3020014
Submission received: 27 March 2024 / Revised: 14 May 2024 / Accepted: 15 May 2024 / Published: 28 May 2024
(This article belongs to the Topic Energy Market and Energy Finance)

Abstract

:
In this study, we investigate the ability of three higher-order risk-neutral return cumulants to predict short maturity (weekly) returns of oil futures. Our data includes weekly West Texas Crude Oil futures options that expire in 7 days (7DTE). Using a model-free approach, we estimate these risk-neutral return cumulants at the beginning of each options expiration cycle. Our results suggest that the third risk-neutral return cumulant consistently predicts the returns of various oil futures (including WTI, Brent, Dubai, Heating Oil, and RBOB Gasoline). We compare our findings with 14 other predictors and offer a theoretical explanation for the negative coefficient observed for the 7DTE third risk-neutral return cumulant. Our theory connects higher-order risk-neutral return cumulants with the risk premiums of oil futures. Furthermore, our quantitative investment strategy favors the predictability of oil futures returns.

1. Introduction

The ever-changing oil market can be influenced by a range of variables, such as economic data, OPEC’s production decisions, geopolitical events, and oil discoveries. To stay ahead, investors closely monitor commodities contracts for oil, particularly the West Texas Intermediate (WTI) crude oil futures and options, which allow them to anticipate these developments and realign their exposures. The recent introduction of weekly options contracts with seven days to expiration (7DTE) for WTI crude oil futures offers new opportunities to manage short maturity price fluctuations.
Our study analyzes the predictability of short maturity returns in these WTI oil futures contracts with options-based information, which are primarily influenced by sudden price changes (jumps) in oil. We demonstrate that our findings are also applicable to futures contracts based on related crude oil products, such as Brent, Dubai, Heating Oil, and RBOB Gasoline.
Empirical measure of 7DTE return asymmetries. Options market-inferred asymmetries of oil futures returns can be a relevant attribute in predicting oil returns, making it a significant area of study for financial economists, investors, and policymakers. In contrast to metrics such as volatility, short maturity return asymmetry provides information about the likely direction of crude oil futures price changes due to its measurement of upside versus downside risk imbalance.
This paper considers the relevance of 7DTE conditional asymmetry of oil futures returns in understanding oil price dynamics and oil futures risk premiums. Using model-free measures of asymmetry, particularly forward-looking asymmetry—the third risk-neutral return cumulant—implied by prices of options, we find that conditional oil asymmetry displays significant variation. This variable contains predictive power for short maturity 7DTE returns of oil futures.
Motivating the role of oil futures return asymmetries and jump risks. There is a scarcity of studies examining the asymmetry of small maturity return distributions in oil markets. Models of oil futures curve and options are often ambivalent about return asymmetries. However, when studying the empirical conditional distributions of small maturity oil futures returns, notable asymmetries emerged, with the sample third risk-neutral return cumulant—inferred from 7DTE options—being stochastic and switches sign.
There are two methods we use to calculate conditional asymmetry, using prices of oil futures and 7DTE options on futures. Realized asymmetry is determined through the actual moments of intraday observed oil futures returns. We utilize returns sampled at five-minute intervals. Conversely, implied asymmetry is based on risk-neutral return cumulants derived from options on oil futures. We concentrate on option-implied third risk-neutral return cumulant for our analysis as it is indicative of subsequent developments, can be updated from the crude oil options market, and is founded in our theory.
Shifts in asymmetry are consistent over periods of time, indicating that the level of downside versus upside risks—as reflected in the relative pricing of 7DTE puts and calls—in the oil futures market is significantly different during these periods. This suggests that investors become more aware about potential market changes and downside versus upside risks to oil prices.
Changes in conditional asymmetry reflect shifts in how oil is perceived and often anticipate changes in supply and demand conditions. This measure contains information for predicting oil price changes. It is a predictor of oil futures returns, even when controlling for other candidate predictors. The COVID-19 pandemic is an example of this predictive power, as implied asymmetry indicated the heightened downside risk to oil and later predicted the rise in oil prices.
Our empirical contribution from 7DTE options and futures returns. Our key finding is that the (model-free) third risk-neutral return cumulant can predict variations in weekly oil futures returns. This finding is new and the third risk-neutral return cumulant survives the inclusion of other considered predictors from oil and equity markets. We provide empirical support based on futures written on WTI, Brent, Dubai, Heating Oil, and RBOB Gasoline.
The combination of the third risk-neutral return cumulant with additional variables, such as realized variance of oil futures returns (calculated from five-minute intraday intervals), leads to an improved ability to predict returns and increases the statistical significance of the predictors. The skewness derived from weekly options on the S&P 500 equity index is another variable that has predictive value, suggesting that both volatility and asymmetric returns in both oil and equity markets play a role in determining subsequent oil futures returns. Additionally, the weekly percentage changes in crude oil stock (available from the U.S. Energy Information Administration) also contributes to predicting short maturity oil futures returns. We also develop a quantitative trading strategy that takes advantage of past options momentum—based on 7DTE third risk-neutral return cumulant—to potentially generate profits in subsequent periods. This strategy is tested in real-world out-of-sample conditions to evaluate its effectiveness.
The remainder of this paper is organized as follows. Section 2 provides a review of the literature, followed, in Section 3, by a presentation of key features of 7DTE options data on WTI oil futures and the S&P 500 index. The question of oil futures return predictability is then examined through various angles (including R os 2 statistic and trading strategy) in Section 4. Next, a framework is presented in Section 5 that connects risk-neutral return cumulants to oil futures risk premiums.

2. Forecasting Oil Futures Returns and Relation to the Literature

A common question that has been addressed in various studies is the prediction of oil prices. For instance, [1,2] examine different models and data related to oil prices. Their aim is to determine whether there is predictability in oil price movements and how effective oil futures markets and survey forecasts are in predicting oil prices. Our approach differs in that we focus on using options data to better understand the predictability of weekly (7DTE) futures returns, taking into account the impact of jumps.
Our empirical study is motivated by the insights from a tractable model of oil futures risk premiums without imposing any restrictive assumptions about the stochastic discount factor (SDF) projection. Several recent papers, including [3,4,5,6,7], have examined the impact of different predictors, focusing on both their ability to accurately forecast and their performance on data that was not used for forecasting.
Ref. [7] use an approach commonly seen in empirical studies to investigate the predictability of oil prices. They construct forecast combinations from several individual predictor forecasts and find that these combined forecasts are more accurate than assuming no change in prices. In particular, they observe significant reductions in mean square forecast errors and improved directional accuracy using combination forecasts compared to the no-change forecast. However, when examining forecast accuracy for oil prices, the authors find no significant improvements in forecast errors or directional accuracy. They also note that using these forecasts to inform investment and hedging decisions does not result in superior economic value for investors. This suggests that the statistical and economic significance of oil price forecasts may be affected by how the underlying price series is constructed. In contrast, our study does not involve predicting oil prices but instead centers on forecasting oil futures returns.
This study is distinct from previous research on the impact of oil price shocks on the economy (such as [8,9]. Specifically, we examine short maturity futures return predictive regressions and examine their significance in understanding risk preferences embedded in the properties of risk-neutral return cumulants of oil futures returns.
Our approach involves incorporating information from weekly oil options, while also taking into consideration the impact of jump risks on futures returns and their expected behavior according to risk-neutral cumulants. This differs from other studies (such as [10,11,12,13,14,15]), which do not take into account the significance of option market variables in predicting oil futures returns.
While [16] examine energy risk premiums, they do not specifically address time-series predictability. Our study connects developments in 7DTE option prices to the predictability of subsequent short maturity oil futures returns.
Ref. [17] investigate the relationship between the skewness of commodity futures returns and expected returns. They find that a strategy focusing on buying futures with highly negative skewness and selling those with highly positive skewness can generate significant excess returns, even after controlling for risk factors. This suggests that the presence of a tradeable skewness factor plays a role in explaining variations in commodity returns, beyond traditional risk premiums. These findings align with investors’ preferences for skewness based on cumulative prospect theory, as well as their selective hedging practices.
Ref. [18] conduct a study on the behavior of WTI and Brent crude oil prices, specifically focusing on the debate surrounding the regionalization and globalization of the two oil markets. They utilize a Dynamic Time Warping technique to identify instances of decoupling between the two oil price series. The results of their study suggest that WTI and Brent prices can decouple and recouple based on local market conditions, indicating that the WTI/Brent market is not consistently integrated.
The study undertaken by [19] investigates the influence of various fundamental factors and financial indicators on the realized volatility of WTI. A time-varying parameter vector autoregression model with stochastic volatility is used to analyze weekly data from January 2008 to October 2021. The results revealed that the volatility of WTI oil price is affected by shocks in oil production, oil inventories, the U.S. dollar index, and VIX. On the other hand, it decreases in response to shocks in the U.S. economic activity. Of all the factors, the VIX index has the most significant impact on the volatility of WTI oil price. These findings shed light on the evolving nature and determinants of WTI oil price volatility.
The idea of characterizing risk-neutral moments—as presented by [20,21]—has been applied in a number of recent studies. Examples of such research include [22,23,24]. Our study places emphasis on distinguishing between downside and upside risks extracted from 7DTE put and call options for WTI crude oil. We show that this distinction has implications for investment management in the 7DTE energy commodities markets.

3. Data on Short Maturity (7DTE) Options and Oil Futures

There has been limited investigation into the information contained in crude oil options with short maturity. To investigate further, let F t represent the date-t price of the oil futures contract at the beginning of the options expiration cycle. Additionally, F t + Δ represents the futures price at the expiration date of options at t + Δ , where Δ corresponds to a week (seven days). The (excess) return of a fully collateralized long position in oil futures contract over seven days is F t + Δ F t 1 .
WTI crude oil serves as the underlier for the oil futures contract, which is listed on the New York Mercantile Exchange. Price data for WTI crude oil is provided by the Chicago Mercantile Exchange and includes spot prices, futures prices, and options on oil futures. Our dataset combines information on futures and options prices, specifically excluding in-the-money options.
The dataset on oil futures prices, including those for the S&P 500 equity index, is constructed over the weekly options expiration cycles (Friday to Friday). Our data covers the period from 12 August 2016, to 24 February 2023, which includes a total of 341 option expiration cycles. This period encompasses the time-series of model-free risk-neutral return cumulants and captures the COVID-19 crisis period, which poses some challenges due to unprecedented volatility levels, but also presents informative insights on uncertainty and investor preferences. To address these challenges, we analyze empirical models that predict the level of the oil futures returns.
The options available have a maturity of seven days, with an average of 45 puts and 46 calls on WTI futures at the beginning of each option expiration cycle, making it suitable for accurately estimating the higher-order risk-neutral return cumulants. The data on S&P 500 equity index options is from the Chicago Board of Trade. See Table 1 for the summary statistics on the number of put and call options.
Table 2 displays the characteristics of five futures contracts related to oil, namely (i) WTI, (ii) Brent, (iii) Dubai, (iv) Heating Oil, and (v) RBOB Gasoline. In Panel A, we present the summary of returns of a long position in these contracts, while in Panel B, we show the annualized volatility over weekly intervals, computed from intraday returns measured at five-minute intervals. As we analyze the correlation patterns of futures returns in Panels C and D, we observe that the first extracted principal component, which represents the average return factor, explains 78.5% of the common variation and has a positive relationship with the returns of all five oil futures contracts. The second factor, responsible for 10.5% of the common variation, is most influenced by Dubai, while the third principal component, accounting for 4.5% of the variation, is more affected by Brent in a positive direction and by Heating Oil in a negative direction. Figure 1 provides a visual representation of crude oil futures prices at the beginning of oil option expiration cycles.

4. Short Maturity (7DTE) Oil Futures Return Predictability

We define the time-varying 7DTE risk-neutral return cumulants as follows:
(1) κ 1 , t Q , oil E t Q F t + Δ F t 1 = 0 , (2) κ 2 , t Q , oil E t Q { F t + Δ F t 1 κ 1 , t Q } 2 , (3) κ 3 , t Q , oil E t Q { F t + Δ F t 1 κ 1 , t Q } 3 , and (4) κ 4 , t Q , oil E t Q { F t + Δ F t 1 κ 1 , t Q } 4 3 ( κ 2 , t Q ) 2 ,
where E t Q ( · ) denotes conditional expectation under the risk-neutral probability measure Q . Since E t Q ( F t + Δ ) = F t (the martingale condition for futures), the first Q cumulant is zero. Hence, κ 1 , t Q , oil contains no information about forecasting short maturity oil futures returns F t + Δ F t 1 .
In our predictive regressions, we focus on measuring risk-neutral return cumulants empirically by utilizing model-free analogs. These measures do not rely on a specific model and are derived from current prices of oil options with one-week maturity and varying strike prices. This approach is based on the notion that weekly option prices reflect the market’s anticipation of return payoffs and, specifically, the jump risks associated with oil, which drive higher-order return cumulants.
Following [20], we have, for n = 2 , 3 , 4 , as follows:
(5) E t Q { F t + Δ F t 1 } n = K < F t w [ K ] put t [ K ] d K + K > F t w [ K ] call t [ K ] d K , (6)     with w [ K ] = R t rf ( n 1 ) n K F t 1 n 2 F t 2 .
Here, R t rf is the gross risk-free return over one week, and put t [ K ] ( call t [ K ] ) are the date-t price of the 7DTE put (call) option, with strike price K, that expires at date t + Δ . These higher-order risk-neutral return cumulants are computed at the beginning of the weekly option expiration cycle.
In economic terms, the second risk-neutral return cumulant represents the conditional return variance when evaluated under the risk-neutral measure, while the third risk-neutral return cumulant reflects conditional return asymmetries. Unlike the S&P 500 index, which demonstrates negative return asymmetries, the oil futures market displays both negative and positive asymmetries. From Table 3, we observe that the 7DTE third risk-neutral return cumulant is negative in 68% of the data points and positive in 32%.
The econometric model that we consider is
F t + Δ F t 1 weekly   oil   futures   returns = constant + b s q u a r e d κ 2 , t Q , o i l + b c u b i c κ 3 , t Q , o i l + b q u a r t i c κ 4 , t Q , o i l + ϵ t + Δ .
We use data from options prices at the beginning of a weekly option expiration period to make predictions about the returns of oil futures. The values and characteristics of these futures prices and return cumulants are influenced by the potential for large oil price fluctuations.
Our predictive exercises are based on linear regressions that use lagged predictor variables to model the (excess) returns of oil futures. The summary statistics and empirical properties for these predictors, which come from both the oil and equity markets, are shown in Table 3. The first-order autocorrelation coefficients—displayed in the “Acf” column—for the predictors are not high. This matters because typically, predictors that have near-unit-root dynamics can lead to inference problems in predictive regressions.
We examine the estimated slope coefficients in predictive regressions and their statistical significance using the p-values of the HAC estimator of Newey and West, with the lag selected automatically. These p-values, designated as NW[p], allow us to determine the level of significance of the predictive coefficients. Additionally, we assess the adequacy of the forecasts by calculating the goodness-of-fit adjusted R 2 ( R ¯ 2 ) value.
The field of commodities has researched methods for forecasting futures returns, but evidence on forecasting 7DTE oil futures returns using 7DTE risk-neutral return cumulants is lacking. This predictive association in Section 5. As a result, our focus is on models in this class that rely on theory. Developing our estimates, our model specifications also include nonstandard models that utilize high-frequency realized returns. These models emphasize the importance of persistence as well as the added predictive power of prior variances and quantities that reflect the perception of jump risks. This can help inform the development of oil models.
Table 4 presents the outcomes of univariate regressions (aligning with multivariate regressions) of the type F t + Δ F t 1 = constant + b × Predictor t + ϵ t + Δ . Figure 2 plots the time-variation in F t + Δ F t 1 of WTI futures whereas Figure 3 displays the 7DTE Skewness t Q , o i l . We use seven predictors from the oil market, seven predictors from the equity market, and three predictors from U.S. Energy Information Administration (EIA). We construct the following predictors:
  • Q skewness from weekly oil options ( Skewness t Q , o i l ). This weekly variable is κ 3 , t Q , o i l { κ 2 , t Q , o i l } 3 / 2 .
  • Q excess kurtosis from weekly oil options ( Kurtosis t Q , o i l ). This weekly variable is κ 4 , t Q , o i l { κ 2 , t Q , o i l } 2 .
  • Realized variance of oil futures returns ( Variance t P , o i l , 5 m i n ). This weekly variable is based on oil futures returns sampled at five-minute intervals.
  • Realized skewness from oil futures returns ( Skewness t P , o i l , 5 m i n ). This weekly variable is based on oil futures returns sampled at five-minute intervals.
  • Q second equity return cumulant ( κ 2 , t Q , e q u i t y ). This is based on weekly S&P 500 equity index options prices. We use E t Q { S t + Δ S t R t r f } n = K < S t R t r f w [ K ] put t [ K ] d K + K > S t R t r f w [ K ] call t [ K ] d K , with w [ K ] = R t r f ( n 1 ) n K S t R t r f n 2 S t 2 .
  • Q third equity return cumulant ( κ 3 , t Q , e q u i t y ). This is based on weekly S&P 500 equity index options prices.
  • Q fourth equity return cumulant ( κ 4 , t Q , e q u i t y ). This is based on weekly S&P 500 equity index options prices.
  • Q skewness from weekly equity options ( Skewness t Q , e q u i t y ). This weekly variable is κ 3 , t Q , e q u i t y { κ 2 , t Q , e q u i t y } 3 / 2 . We use S&P 500 equity index options prices.
  • Q excess kurtosis from weekly equity options ( Kurtosis t Q , e q u i t y ). This weekly variable is κ 4 , t Q , e q u i t y { κ 2 , t Q , e q u i t y } 2 . We use S&P 500 equity index options prices.
  • Realized variance of equity futures returns ( Variance t P , e q u i t y , 5 m i n ). This weekly variable is based on S&P 500 E-mini equity futures returns sampled at five-minute intervals.
  • Realized skewness from equity futures returns ( Skewness t P , e q u i t y , 5 m i n ). This weekly variable is based on S&P 500 E-mini equity futures returns sampled at five-minute intervals.
  • Growth rate of crude oil stock ( EIA t s t o c k g r o w t h ). This is constructed based on EIA releases of petroleum status reports (https://www.eia.gov/petroleum/supply/weekly/, accessed on 14 May 2024). The underlying quantity is crude oil stock.
  • Growth rate of crude oil production ( EIA t p r o d u c t i o n g r o w t h ). The underlying variable is domestic crude oil production (EIA estimates).
  • Growth rate of crude oil imports ( EIA t i m p o r t g r o w t h ). The underlying variable is crude oil imports (EIA estimates).
Adhering to our approach, we analyze the influence of various predictive factors such as oil-specific variables and equity market-specific variables on the forecasting performance of our model. We find that the coefficient estimates for the third risk-neutral return cumulant, κ 3 , t Q , o i l , and the fourth risk-neutral return cumulant, κ 4 , t Q , o i l , are significant with negative and positive values, respectively. The estimate for the second risk-neutral return cumulant, κ 2 , t Q , o i l , is close to zero and is not significant. Additionally, many other potential predictors are individually insignificant.
Despite this, the NW[p]-value for the slope coefficient related to the third risk-neutral return cumulant is below the two-sided 5% significance level, indicating significance. Similarly, the predictors Skewness t Q , e q u i t y and EIA t s t o c k g r o w t h also have NW[p] values below 0.05 and are considered significant predictors. However, at the short maturity weekly forecasting horizon, these univariate predictors explain less than 2% of the variability in crude oil futures returns.
To gain a better understanding of the results in a wider empirical context, we investigate whether the third risk-neutral return cumulant ( κ 3 , t Q , o i l ) maintains its predictive power when other variables are included in bivariate predictive regressions. Specifically, we examine short maturity regressions of the following form:
F t + Δ F t 1 = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ .
The results in Table 5 show a consistent pattern that κ 3 , t Q , o i l remains a predictor of futures price changes, with b c u b i c negative and statistically significant. This aligns with our theoretical framework presented in Section 5, which suggests that the third risk-neutral return cumulant captures important aspects of asymmetry as well as downside and upside jump risks. We find that κ 3 , t Q , oil is not highly correlated with other predictors.
Additionally, the NW[p] values for four other variables are below 0.05 , making them stand out as predictors: (i) κ 4 , t Q , o i l , (ii) Variance t P , o i l , 5 m i n , (iii) Skewness t Q , e q u i t y , and (iv) EIA t s t o c k g r o w t h . The combination of κ 3 , t Q , o i l and Variance t P , o i l , 5 m i n has the highest significance and an R ¯ 2 of 3.47%. Comparing these results with other predictors, the slope coefficients ( b a l t e r ) are not statistically significant and the R ¯ 2 values are lower. However, incorporating S&P 500 equity index based predictors does increase the R ¯ 2 values, suggesting that additional variables from the equity market may be helpful in explaining the variations in oil futures returns.
We further investigate the relationship between oil prices and the third risk-neutral return cumulant by using various approaches to construct the returns of a crude oil basket, including the use of several individual crude oil futures, equal weight basket of five crude oil futures, and regression analysis. Specifically, we construct an equal weight basket of five crude oil futures, including WTI, Brent, Dubai, Heating Oil, and RBOB Gasoline, as follows:
z t + Δ b a s k e t , o i l = i = 1 5 1 5 ( F t + Δ i F t i 1 ) .
We use this basket, labeled z t + Δ b a s k e t , o i l , to smooth out idiosyncratic variations in the five individual oil futures prices. The bivariate regressions in Table 6 shows that the coefficients for 7DTE κ 3 , t Q , o i l are negative, suggesting a significant relationship between z t + Δ b a s k e t , o i l and crude oil options-based κ 3 , t Q , o i l . Furthermore, the R ¯ 2 values tend to improve on the bivariate counterparts in Table 5.
Table 7, Table 8, Table 9 and Table 10 analyze the impact of individual crude oil futures, respectively, for Brent, Dubai, Heating Oil, and RBOB Gasoline, on predicting weekly oil returns and finds that 7DTE κ 3 , t Q , o i l remains a relevant predictor. All in all, this suggests that the third risk-neutral return cumulant has predictive power for oil futures returns.
Next, we assess the out-of-sample performance using the [25] statistic, which is calculated as R o s 2 = 1 t = 1 T z t z ^ t , f u l l 2 t = 1 T z t z ^ t , n e s t e d 2 . We use the portion of the data up until the first weekly cycle of May 2020 (i.e., 194 (out of 341) cycles) as the initial in-sample estimates for coefficients. Our implementation relies on the expanding scheme updated weekly.
  • When the nested model is the historical average and κ 3 , t Q , o i l is the full model, the R o s 2 values are 1.13% and 0.97% for WTI and the oil basket, respectively. These values suggest that using κ 3 , t Q , o i l as a predictor adds to the predictive ability beyond using the historical average.
  • When the alternative predictor is the nested model and the bivariate predictor constitutes the full model, the resulting R o s 2 values reported in Table 5 and Table 6 are all positive and range from 0.25% to 4.1%. These R o s 2 values align with the notion that κ 3 , t Q , o i l has additional predicting power over the considered alternative predictor.
The conclusion to draw is that, in an out-of-sample setting, κ 3 , t Q , o i l is an informative predictor.
Table 11 builds on the empirical results from Table 5 and Table 6 and provides further evidence supporting the relevance of 7DTE κ 3 , t Q , o i l as a predictor for oil futures returns. First, to evaluate the statistical significance of our results, we include the effect of all three higher-order risk-neutral return cumulants (as in Equation (7)). Second, we analyze whether there is a significant distinction in oil futures returns between three sample periods: pre-COVID, COVID, and post-COVID. This was done by conducting a regression on weekly oil futures returns using dummy variables to indicate different sample periods, as follows:
F t + Δ F t 1 W T I o r b a s k e t = constant + ( b s q u a r e d , p r e - c o v i d 1 p r e - c o v i d + b s q u a r e d , c o v i d 1 c o v i d + b s q u a r e d , p o s t - c o v i d 1 p o s t - c o v i d ) κ 2 , t Q , o i l     + ( b c u b i c , p r e - c o v i d 1 p r e - c o v i d + b c u b i c , c o v i d 1 c o v i d + b c u b i c , p o s t - c o v i d 1 p o s t - c o v i d ) κ 3 , t Q , o i l (10)     + ( b q u a r t i c , p r e - c o v i d 1 p r e - c o v i d + b q u a r t i c , c o v i d 1 c o v i d + b q u a r t i c , p o s t - c o v i d 1 p o s t - c o v i d ) κ 4 , t Q , o i l + ϵ t + Δ .
The results from the multiple regressions in Table 11 (Panel A) reveal that b c u b i c displays a high level of significance, while the coefficients for the risk-neutral second and fourth return cumulants become insignificant. When κ 2 , t Q , o i l , κ 3 , t Q , o i l , and κ 4 , t Q , o i l are combined, the overall explanatory power, measured by R ¯ 2 , is 1.31%, compared to 3.47% using just κ 3 , t Q , o i l with variance t P , o i l , 5 m i n (see Table 5). This suggests that other variables—related to realized oil futures return variance, the equity market, and EIA data—contribute to explaining the variation in the oil futures returns, compared to using information only from the oil options market.
The results, presented in Table 11 (Panel B), show that the returns of oil futures are affected differently by risk-neutral return cumulants. The dummy variables 1 p r e - c o v i d , 1 covid , and 1 p o s t - c o v i d correspond to the time period before the COVID-19 pandemic (12 August 2016, to 31 January 2020), during the pandemic (7 February 2020, to 5 June 2020), and after the pandemic (12 June 2020, to 23 February 2023). We find that b c u b i c ,   c o v i d and b c u b i c ,   p o s t - c o v i d are highly significant, indicating that the changes in the third risk-neutral return cumulant exert an effect on the oil futures returns. The second risk-neutral return cumulant does not show a significant impact on the futures returns.
Overall, we present an approach for predicting futures returns in the oil market using the 7DTE third risk-neutral return cumulant, which is an indicator that varies with economic conditions. This variable is derived from the crude oil options market and, if reliable, can provide insights for predicting oil price movements. However, we must consider that this finding may potentially be a chance discovery and, therefore, it is essential to consider the model’s performance on out-of-sample data and trading rule.
Mindful of such possibilities, we implement a trading strategy in which we take a long or short position in the basket of crude oil commodities. The rule that we adhere to is as follows:
Trading position t = Long the oil futures basket i f j = 1 J κ 3 , t j Δ Q , o i l < 0 , Short the oil futures basket i f j = 1 J κ 3 , t j Δ Q , o i l > 0 .
In words, we take a long position in each of the basket components when the prior momentum in third risk-neutral return cumulant is negative (i.e., puts remain more expensive than calls). The return on this signal over the subsequent week is z t + Δ b a s k e t , o i l . To the contrary, we take a short position in the basket of oil futures when the prior momentum is positive. The return to this signal is z t + Δ b a s k e t , o i l . It remains to be seen if this empirical model can translate into a viable investment strategy, as it will require consideration of the signal-to-noise ratio and our ability to diversify and rebalance when the signal changes direction.
The results in Table 12 show the mean returns of our strategy, along with the corresponding 90% block bootstrap confidence intervals. While this method shows a positive mean return, it also faces a high level of return variability over the entire sample. However, in the post COVID-19 sample, the strategy produces significantly positive mean returns. For instance, using the prior six week momentum in κ 3 , t Q , o i l results in an average weekly return of 0.77% with a 90% confidence interval of 0.16 % 1.42 % . This strategy yields reliable results with a frequency of 62% and is supported by the reported sizable Sharpe Ratio’s. Though the simple structure of the model is appealing, we acknowledge that the predictability of oil prices is influenced by multiple factors beyond just economic cycles related to the 7DTE third risk-neutral return cumulant.
In summary, market expectations of asymmetric behavior in oil futures returns is a consistent predictor of oil returns. This is due to the fact that the third risk-neutral return cumulant, which measures asymmetry and risk imbalances, provides insight into the likely direction of oil changes. In contrast, metrics like volatility offer information only about the magnitude of changes. This paper highlights the significance of conditional asymmetry in explaining oil futures dynamics and risk premiums. By using model-free measures of the third risk-neutral return cumulant, specifically the forward-looking asymmetry derived from option prices, our study reveals that this measure exhibits cyclical variation. Furthermore, this asymmetry measure is a predictor of subsequent oil futures returns and is not subsumed by information obtained from other predictors.
Other papers connect to our analysis, including, for instance, [26,27,28,29,30,31,32,33,34,35,36,37].

5. Risk Preferences and the Third Risk-Neutral Return Cumulant

The third risk-neutral return cumulant κ 3 , t Q , o i l is linked to the risk premium of oil futures, and the sign of b c u b i c in the regressions can provide insights into the characteristics of the risk-neutral oil density and the corresponding futures risk premium. The estimated predictive coefficient, b c u b i c , has shown a negative pattern in Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11. A negative (positive) κ 3 , t Q , o i l indicates a left-skewed (right-skewed) futures return distribution, where extreme negative (positive) events are more likely to occur, leading to a higher (lower) oil futures risk premium.
The negative sign of b c u b i c poses a potential challenge for models explaining the behavior of oil futures returns. This section proposes a framework to reconcile the negative value of b c u b i c by exploring restrictions on the properties of the SDF of the oil market. Define
gross return of oil futures Z t + Δ F t + Δ F t and net ( excess ) return z t + Δ = F t + Δ F t 1 .
Additionally, let P represent the real-world probability measure and p [ z t + Δ ] the density of oil futures returns. Define Ω = { z t + Δ > 1 } and note that p [ z t + Δ ] satisfies 1 = Ω p [ z t + Δ ] d z t + Δ .

5.1. Oil Futures Risk Premiums

Refs. [11,14] have investigated the properties of the SDFs for the crude oil market. However, their studies did not explore the characteristics of 7DTE SDFs. Additionally, the content of weekly oil options and jump behaviors have yet to be fully examined in the existing literature. Furthermore, the potential implications of the negative predictive coefficient on κ 3 , t Q , o i l has not been analyzed.
Our substantive result is a theoretical link between the risk premium of oil futures, namely, E t P ( F t + Δ F t 1 ) , and the higher-order risk-neutral return cumulants. We present a theory that matches the sign of the estimated coefficient b c u b i c in Equation (7). Furthermore, we show that the constructed 7DTE κ 3 , t Q , o i l can switch direction in the oil market.
We assume that there exists a SDF projection of the oil market, m [ z t + Δ ] , satisfying
m [ z t + Δ ] > 0 , E t P ( m [ z t + Δ ] ) = 1 R t r f < , and E t P ( { m [ z t + Δ ] } 2 ) < .
In line with [38,39], we only need to model m [ z t + Δ ] without incorporating state-dependent specifications. We write Q P if the measure Q , corresponding to the density q [ z t + Δ ] , is absolutely continuous with respect to measure P . We assume both Q P and P Q hold, and, thus, Q and P are equivalent.
The n-th order raw conditional return moments of the distribution, under P and Q , are
μ n , t P Ω { z t + Δ } n p [ z t + Δ ] d z t + Δ and μ n , t Q Ω { z t + Δ } n q [ z t + Δ ] d z t + Δ , respectively .
We denote the moment-generating functions under P and Q , as mgf t P [ λ ] and mgf t Q [ λ ] , respectively, for some constant parameter λ . Assume mgf t P [ λ ] < and mgf t Q [ λ ] < . For brevity, we omit the oil superscript on both κ n , t P , o i l and κ n , t Q , o i l . To develop expressions that rely on variables constructed from option prices, we assume that m [ z t + Δ ] is continuous and infinitely differentiable in z t + Δ ; that is, m [ z t + Δ ] C . We define
G [ z t + Δ ] 1 m [ z t + Δ ] .
Pertinent to our formalizations, the properties of G [ z t + Δ ] relate to the shape of m [ z t + Δ ] . The class G [ z t + Δ ] C encompasses empirically plausible specifications of m [ z t + Δ ] . We denote
G [ z t + Δ ] d G [ z t + Δ ] d z t + Δ , G [ z t + Δ ] d 2 G [ z t + Δ ] d z t + Δ 2 , and G [ z t + Δ ] d 3 G [ z t + Δ ] d z t + Δ 3 .
Additionally,
(17) G [ 0 ] = G [ z t + Δ ] | z t + Δ = 0 , G [ 0 ] = G [ z t + Δ ] | z t + Δ = 0 , and (18) G [ 0 ] = G [ z t + Δ ] | z t + Δ = 0 , G [ 0 ] = G [ z t + Δ ] | z t + Δ = 0 .
The remaining step is to formulate the expression for the real-world density of oil futures returns in terms of the risk-neutral return density, as follows:
p [ z t + Δ ] = G [ z t + Δ ] q [ z t + Δ ] Ω G [ z t + Δ ] q [ z t + Δ ] d z t + Δ .
Then, Equation (20) is the oil futures risk premium counterpart of that in [39] for equity market index (proof available from authors):
E t P F t + Δ F t 1 o i l f u t u r e s r i s k p r e m i u m = G [ 0 ] G [ 0 ] κ 2 , t Q + 1 2 G [ 0 ] G [ 0 ] κ 3 , t Q + 1 6 G [ 0 ] G [ 0 ] κ 4 , t Q + .
The expression for E t P ( F t + Δ F t 1 ) is an infinite series and is determined by its sensitivity to the risk-neutral return cumulants ( κ n , t Q ). This representation is based on the existence and well-definedness of G [ 0 ] G [ 0 ] , G [ 0 ] G [ 0 ] , and G [ 0 ] G [ 0 ] . Our analysis does not require a parametric assumption about the real-world oil futures return distributions or the risk-neutral oil futures return distributions.

5.2. Sign of G [ 0 ] G [ 0 ] (Theoretical Counterpart to b c u b i c ) in Theoretical Economies

To clarify, our focus is on the m [ z t + Δ ] and its reciprocal, the G [ z t + Δ ] function, in reproducing the negative value of predictive coefficient b c u b i c (which coincides with the sign of G [ 0 ] G [ 0 ] ). In order to ascertain the restrictions needed for this result, our study examines the example economies for the sign of G [ 0 ] G [ 0 ] .
Example 1
(Monotonic SDFs). Suppose that the SDF is monotonic in oil futures returns. That is, for constants m 0 > 0 and ϕ,
m [ z t + Δ ] = m 0 exp ϕ ( 1 + z t + Δ ) .
Here G [ z t + Δ ] = 1 m [ z t + Δ ] = 1 m 0 exp ϕ ( 1 + z t + Δ ) . Hence, G [ z t + Δ ] = ϕ m 0 exp ϕ ( 1 + z t + Δ ) , G [ z t + Δ ] = ϕ 2 m 0 exp ϕ ( 1 + z t + Δ ) , and G [ z t + Δ ] = ϕ 3 m 0 exp ϕ ( 1 + z t + Δ ) . In this illustration, it follows that
G [ 0 ] G [ 0 ] = ϕ , 1 2 G [ 0 ] G [ 0 ] = 1 2 ϕ 2 > 0 , and 1 6 G [ 0 ] G [ 0 ] = 1 6 ϕ 3 .
Monotonic SDFs ( ϕ > 0 or ϕ < 0 ) cannot support the estimated negative b c u b i c .
Example 2
(Variance-dependent SDFs). Suppose that for constants m 0 > 0 and η > 0 ,
m [ z t + Δ ] = m 0 exp η { log ( 1 + z t + Δ ) } 2 .
The essence of this model is that
G [ 0 ] G [ 0 ] = 0 , 1 2 G [ 0 ] G [ 0 ] = η < 0 , and 1 6 G [ 0 ] G [ 0 ] = η > 0 .
This form of variance-dependent SDF can support the negative b c u b i c . Hence, variance-dependent SDFs—that symmetrically weight down and up returns—can be consistent with our findings.
Example 3
(SDFs that incorporate asymmetry in down and up oil futures return movements). Suppose that for constants m 0 > 0 and η > 0 ,
m [ z t + Δ ] = m 0 exp η 1 + z t + Δ + η 2 ( 1 + z t + Δ ) 2 .
It follows that
G [ 0 ] G [ 0 ] = 0 , 1 2 G [ 0 ] G [ 0 ] = 3 2 η < 0 , and 1 6 G [ 0 ] G [ 0 ] = η > 0 .
This form of the SDF gives rise to a negative b c u b i c and a positive b q u a r t i c . Unlike Example 2, this model does not impose equal and opposite signed effects on b c u b i c and b q u a r t i c .
Example 4
(SDFs sensitive to the cubic oil futures return payoff). Suppose that for constant η > 0 ,
m [ z t + Δ ] = 1 4 3 η exp η 1 + z t + Δ + η 3 ( 1 + z t + Δ ) 3 .
For this model, which implies risk imbalance to the downside versus the upside of returns, it follows that
G [ 0 ] G [ 0 ] = 0 , 1 2 G [ 0 ] G [ 0 ] = 2 η < 0 , and 1 6 G [ 0 ] G [ 0 ] = 2 3 η > 0 .
This form of the SDF gives rise to a negative b c u b i c and a positive b q u a r t i c . The SDF m [ z t + Δ ] is an asymmetric U-shaped function, with greater weights to the downside and is centered at one.
Example 5
(SDFs exponential quadratic in log futures returns). Suppose that
m [ z t + Δ ] = exp a × log ( 1 + z t + Δ ) + b × { log ( 1 + z t + Δ ) } 2 .
For constant parameters a > 0 and b > 0 , it follows that
G [ 0 ] G [ 0 ] = a , 1 2 G [ 0 ] G [ 0 ] = 1 2 { a + a 2 2 b } , 1 6 G [ 0 ] G [ 0 ] = 1 6 { 2 a + a 3 3 a ( a + 2 b ) + 6 b } .
For a = 1.5 and b = 9.5 , this SDF implied b c u b i c = 9.12 , b q u a r t i c = 4.81 , and b s q u a r e d = 1.5 . The SDF m [ z t + Δ ] is an asymmetric U-shaped function with greater weighting to the downside of short maturity oil futures returns. This parameterization mimics features of Table 11 (Panel A).
In summary, the theory proposed suggests that incorporating higher-order return dependence in the SDFs can result in negative predictive coefficients for the third risk-neutral return cumulant. This perspective takes into account the impact of down and up jumps and has not been previously investigated in the analysis of short maturity SDFs.

6. Conclusions

This study reveals that the third risk-neutral return cumulant—inferred from short maturity (weekly expiring (7DTE)) crude oil option prices—is a consistent predictor of short maturity oil futures returns. This finding is supported by evidence across different model specifications, even when other predictors are taken into account. The return predictability is strongest when the 7DTE third risk-neutral return cumulant is combined with measures of (i) realized variance of intraday oil futures returns (based on five-minute returns), (ii) skewness-derived variable from S&P 500 equity options, and (iii) growth rate of crude oil stock. The adjusted R 2 value for weekly horizon oil returns achieves a maximum of 3.47%. We also implement a quantitative trading strategy that reflect the measurement—its prior momentum—of the 7DTE third risk-neutral return cumulant.
These empirical observations from our return predictive regressions support a theoretical framework that includes an SDF projection whose dispersion reflects the effects of varying levels of oil futures return volatility. The estimated negative predictive coefficient for the third risk-neutral return cumulant is also in line with this theory, as it considers the inclusion of squared and cubic return dependence in the SDF projection.

Author Contributions

Conceptualization, G.B., X.G. and Z.Z.; Methodology, G.B.; Formal analysis, G.B., X.G. and Z.Z.; Writing—original draft, G.B., X.G. and Z.Z.; Writing—review & editing, G.B., X.G. and Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Available upon request.

Acknowledgments

First Draft: January 2024. We are immensely grateful to two anonymous referees and the Editor for guiding the work done in this paper and for suggesting numerous improvements.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Alquist, R.; Kilian, L. What do we learn from the price of crude oil futures? J. Appl. Econom. 2010, 25, 539–573. [Google Scholar] [CrossRef]
  2. Alquist, R.; Kilian, L.; Vigfusson, R. Forecasting the price of oil. In Handbook of Economic Forecasting; Elsevier: Amsterdam, The Netherlands, 2013; Volume 2, pp. 427–507. [Google Scholar]
  3. Miao, H.; Ramchander, S.; Wang, T.; Yang, D. Influential factors in crude oil price forecasting. Energy Econ. 2017, 68, 77–88. [Google Scholar] [CrossRef]
  4. Funk, C. Forecasting the real price of oil—Time-variation and forecast combination. Energy Econ. 2018, 76, 288–302. [Google Scholar] [CrossRef]
  5. Cotter, J.; Eyiah-Donkor, E.; Poti, V. Commodity Futures Return Predictability and Intertemporal Asset Pricing; Working Paper; University College Dublin: Dublin, Ireland, 2020. [Google Scholar]
  6. Ellwanger, R.; Snudden, S. Forecasts of the real price of oil revisited: Do they beat the random walk? J. Bank. Financ. 2023, 154, 106962. [Google Scholar] [CrossRef]
  7. Conlon, T.; Cotter, J.; Eyiah-Donkor, E. Forecasting the price of oil: A cautionary note. J. Commod. Mark. 2024, 33, 100378. [Google Scholar] [CrossRef]
  8. Killian, L. Not all oil price shocks are alike: Disentangling demand and supply shocks in the crude oil market. Am. Econ. Rev. 2009, 99, 1053–1069. [Google Scholar] [CrossRef]
  9. Baumeister, C.; Hamilton, J. Structural interpretation of vector autoregressions with incomplete identification: Revisiting the role of oil supply and demand shocks. Am. Econ. Rev. 2019, 109, 1873–1910. [Google Scholar] [CrossRef]
  10. Routledge, B.; Seppi, D.; Spatt, C. Equilibrium forward curves for commodities. J. Financ. 2000, 55, 1297–1338. [Google Scholar] [CrossRef]
  11. Trolle, A.; Schwartz, E. Unspanned stochastic volatility and the pricing of commodity derivatives. Rev. Financ. Stud. 2009, 22, 4423–4461. [Google Scholar] [CrossRef]
  12. Chiang, I.; Hughen, K.; Sagi, J. Estimating Oil Risk Factors Using Information from Equity and Derivatives Markets. J. Financ. 2015, 70, 769–804. [Google Scholar] [CrossRef]
  13. Ready, R. Oil consumption, economic growth, and oil futures: The impact of long-run oil supply uncertainty on asset prices. J. Monet. Econ. 2018, 94, 1–26. [Google Scholar] [CrossRef]
  14. Christoffersen, P.; Jacobs, K.; Pan, X. The State Price Density Implied by Crude Oil Futures and Option Prices. Rev. Financ. Stud. 2022, 35, 1064–1103. [Google Scholar] [CrossRef]
  15. Gao, L.; Hitzemann, S.; Shaliastovich, I.; Xu, L. Oil Volatility Risk. J. Financ. Econ. 2022, 144, 456–494. [Google Scholar] [CrossRef]
  16. Kang, W.; Rouwenhorst, G.; Tang, K. A tale of two premiums: The role of hedgers and speculators in commodity futures markets. J. Financ. 2014, 75, 377–417. [Google Scholar] [CrossRef]
  17. Fernandez-Perez, A.; Frijns, B.; Fuertes, A.M.; Miffre, J. The skewness of commodity futures returns. J. Bank. Financ. 2018, 86, 143–158. [Google Scholar] [CrossRef]
  18. Mastroeni, L.; Mazzoccoli, A.; Quaresima, G.; Vellucci, P. Decoupling and recoupling in the crude oil price benchmarks: An investigation of similarity patterns. Energy Econ. 2021, 94, 105036. [Google Scholar] [CrossRef]
  19. Le, T.H.; Boubaker, S.; Bui, M.T.; Park, D. On the volatility of WTI crude oil prices: A time-varying approach with stochastic volatility. Energy Econ. 2023, 117, 106474. [Google Scholar] [CrossRef]
  20. Bakshi, G.; Madan, D. Spanning and Derivative-Security Valuation. J. Financ. Econ. 2000, 55, 205–238. [Google Scholar] [CrossRef]
  21. Bakshi, G.; Kapadia, N.; Madan, D. Stock Return Characteristics, Skew Laws, and the Differential Pricing of Individual Equity Options. Rev. Financ. Stud. 2003, 16, 101–143. [Google Scholar] [CrossRef]
  22. Prokopczuk, M.; Symeonidis, L.; Simen, C. Variance risk in commodity markets. J. Bank. Financ. 2017, 81, 136–149. [Google Scholar] [CrossRef]
  23. Ruan, X.; Zhang, J. Risk-neutral moments in the crude oil market. Energy Econ. 2018, 72, 583–600. [Google Scholar] [CrossRef]
  24. Wong, P. Predicting intraday crude oil returns with higher order risk-neutral moments. J. Commod. Mark. 2023, 100331. [Google Scholar] [CrossRef]
  25. Campbell, J.; Thompson, S. Predicting Excess Stock Returns Out of Sample: Can Anything Beat the Historical Average? Rev. Financ. Stud. 2008, 21, 1509–1531. [Google Scholar] [CrossRef]
  26. Fuertes, A.M.; Miffre, J.; Rallis, G. Tactical allocation in commodity futures markets: Combining momentum and term structure signals. J. Bank. Financ. 2010, 34, 2530–2548. [Google Scholar] [CrossRef]
  27. Casassus, J.; Liu, P.; Tang, K. Economic Linkages, Relative Scarcity, and Commodity Futures Returns. Rev. Financ. Stud. 2013, 26, 1324–1362. [Google Scholar] [CrossRef]
  28. Ahmed, S.; Tsvetanov, D. The predictive performance of commodity futures risk factors. J. Bank. Financ. 2016, 71, 20–36. [Google Scholar] [CrossRef]
  29. Gao, X.; Li, B.; Liu, R. The relative pricing of WTI and Brent crude oil futures: Expectations or risk premia? J. Commod. Mark. 2023, 30, 100274. [Google Scholar] [CrossRef]
  30. Yang, Y.; Zhang, J.; Chen, S. Information effects of monetary policy announcements on oil price. J. Commod. Mark. 2023, 30, 100268. [Google Scholar] [CrossRef]
  31. Chang, C.L. Extreme events, economic uncertainty and speculation on occurrences of price bubbles in crude oil futures. Energy Econ. 2024, 130, 107318. [Google Scholar] [CrossRef]
  32. Cui, J.; Maghyereh, A. Unveiling interconnectedness: Exploring higher-order moments among energy, precious metals, industrial metals, and agricultural commodities in the context of geopolitical risks and systemic stress. J. Commod. Mark. 2024, 33, 100380. [Google Scholar] [CrossRef]
  33. Dai, Z.; Zhang, X.; Liang, C. Efficient predictability of oil price: The role of VIX-based panic index shadow line difference. Energy Econ. 2024, 129, 107234. [Google Scholar] [CrossRef]
  34. Gronwald, M.; Jin, X. Measuring world oil market integration with a Thick Pen. Energy Econ. 2024, 130, 107315. [Google Scholar] [CrossRef]
  35. Gupta, R.; Nielsen, J.; Pierdzioch, C. Stock market bubbles and the realized volatility of oil price returns. Energy Econ. 2024, 132, 107432. [Google Scholar] [CrossRef]
  36. Luo, Q.; Ma, F.; Wang, J.; Wu, Y. Changing determinant driver and oil volatility forecasting: A comprehensive analysis. Energy Econ. 2024, 129, 107187. [Google Scholar] [CrossRef]
  37. Zhang, L.; Chen, Y.; Bouri, E. Time-varying jump intensity and volatility forecasting of crude oil returns. Energy Econ. 2024, 129, 107236. [Google Scholar] [CrossRef]
  38. Bakshi, G.; Gao, X.; Xue, J. Recovery with Applications to Forecasting Equity Disaster Probability and Testing the Spanning Hypothesis in the Treasury Market. J. Financ. Quant. Anal. 2023, 58, 1808–1842. [Google Scholar] [CrossRef]
  39. Bakshi, G.; Crosby, J.; Gao, X.; Zhou, W. A New Formula for the Expected Excess Return of the Market; Working Paper; Temple University: Philadelphia, PA, USA, 2024. [Google Scholar]
Figure 1. Crude oil prices at the beginning of each option expiration cycle on Friday. We graph the value of oil futures at the start of each weekly options expiration cycle, specifically at the close of trading on Fridays. We have data for 341 cycles, starting from 12 August 2016 and ending on 24 February 2023.
Figure 1. Crude oil prices at the beginning of each option expiration cycle on Friday. We graph the value of oil futures at the start of each weekly options expiration cycle, specifically at the close of trading on Fridays. We have data for 341 cycles, starting from 12 August 2016 and ending on 24 February 2023.
Commodities 03 00014 g001aCommodities 03 00014 g001b
Figure 2. Weekly (Friday to Friday) WTI oil futures returns. We plot the change in the value of oil futures over a weekly options expiration cycle (Friday to Friday), calculated as F t + Δ F t 1 . The market for these futures is based on West Texas Intermediate (WTI) crude oil. Our dataset tracks the prices for 341 weekly option expiration cycles, spanning from 12 August 2016, to 24 February 2023.
Figure 2. Weekly (Friday to Friday) WTI oil futures returns. We plot the change in the value of oil futures over a weekly options expiration cycle (Friday to Friday), calculated as F t + Δ F t 1 . The market for these futures is based on West Texas Intermediate (WTI) crude oil. Our dataset tracks the prices for 341 weekly option expiration cycles, spanning from 12 August 2016, to 24 February 2023.
Commodities 03 00014 g002
Figure 3. Measure of asymmetry of the risk-neutral return distribution of oil futures. We display the risk-neutral skewness of oil futures over weekly options expiration cycles, with WTI crude oil serving as the underlying asset for the futures. The skewness is calculated for each weekly options expiration cycle from 12 August 2016, to 24 February 2023, which covers a total of 341 cycles. Our formula for skewness is Skewness t Q , o i l = E t Q ( { F t + Δ F t 1 } 3 ) / { E t Q ( { F t + Δ F t 1 } 2 ) } 3 / 2 . Specifically, E t Q ( { F t + Δ F t 1 } n ) = K < F t w [ K ] put t [ K ] d K + K > F t w [ K ] call t [ K ] d K , with w [ K ] = R t r f ( n 1 ) n ( K F t 1 ) n 2 F t 2 .
Figure 3. Measure of asymmetry of the risk-neutral return distribution of oil futures. We display the risk-neutral skewness of oil futures over weekly options expiration cycles, with WTI crude oil serving as the underlying asset for the futures. The skewness is calculated for each weekly options expiration cycle from 12 August 2016, to 24 February 2023, which covers a total of 341 cycles. Our formula for skewness is Skewness t Q , o i l = E t Q ( { F t + Δ F t 1 } 3 ) / { E t Q ( { F t + Δ F t 1 } 2 ) } 3 / 2 . Specifically, E t Q ( { F t + Δ F t 1 } n ) = K < F t w [ K ] put t [ K ] d K + K > F t w [ K ] call t [ K ] d K , with w [ K ] = R t r f ( n 1 ) n ( K F t 1 ) n 2 F t 2 .
Commodities 03 00014 g003
Table 1. Number of OTM puts and calls employed in the construction of 7DTE risk-neutral return cumulants.
Table 1. Number of OTM puts and calls employed in the construction of 7DTE risk-neutral return cumulants.
MeanSDMin.5th25th50th75th95thMax.
WTI crude oil futures: Number of OTM puts452981219466595182
WTI crude oil futures: Number of OTM calls46321113204462101203
S&P 500 equity index: Number of OTM puts985314365783136193322
S&P 500 equity index: Number of OTM calls352481419264487158
The time frame of our study is confined to 12 August 2016, to 24 February 2023, as we use the weekly (7DTE) options on crude oil futures and the S&P 500 equity index. This period was chosen due to the need for consistent and accurate data at regular intervals, as well as reliable options-derived quantities. Tabulated are the summary statistics on OTM puts and calls at the beginning of the options expiration cycles (on Friday). The options data for WTI crude oil futures are obtained from the Chicago Mercantile Exchange, while the data for matching S&P 500 equity options are obtained from the Chicago Board of Trade. The standard deviation is shown in the column “SD”. Our analysis covers a total of 341 weekly option expiration cycles.
Table 2. Summary statistics of weekly oil futures returns and volatility of futures returns.
Table 2. Summary statistics of weekly oil futures returns and volatility of futures returns.
Mean
(%)
SD
(%)
Block
Bootstrap
NW[ p ]Min.Max.AcfSkewnessKurtosis 1 { z > 0 }
(%)
LowerUpper
Panel A: Oil futures returns, F t + Δ F t 1 (weekly, %)
WTI futures0.245.86⌊−0.320.78⌋0.46−32.324.70.11−0.75.955
Brent futures0.305.08⌊−0.180.74⌋0.29−20.923.30.000.02.956
Dubai futures0.314.89⌊−0.260.85⌋0.32−20.935.60.111.213.155
Heating Oil futures0.345.31⌊−0.110.80⌋0.23−24.133.3−0.010.56.353
RBOB Gasoline futures0.406.13⌊−0.220.97⌋0.27−33.024.40.08−0.45.253
Equal weight basket0.324.84⌊−0.180.81⌋0.28−22.425.10.090.014.155
Panel B: Volatility of oil futures returns (from intraday returns, annualized (%))
WTI futures29.032.8⌊24.635.0⌋0.009.9539.20.4711.9175.4
Brent futures27.015.4⌊24.330.3⌋0.0010.9163.60.734.731.4
Dubai futures27.734.3⌊23.133.5⌋0.001.2349.60.364.428.4
Heating Oil futures24.613.1⌊22.027.3⌋0.0011.7121.40.793.315.8
RBOB Gasoline futures27.617.2⌊24.531.5⌋0.0013.9155.80.884.828.4
Panel C: Return correlationsPanel D: Correlation between return volatilities
WTIBrentDubaiHeating Oil WTIBrentDubaiHeating Oil
Brent futures0.86 0.81
Dubai futures0.610.58 0.480.70
Heating Oil futures0.850.770.62 0.710.910.70
RBOB Gasoline futures0.810.810.570.79 0.760.920.650.87
The time frame of our study is confined to 12 August 2016, to 24 February 2023, as we use the weekly (7DTE) options on crude oil futures and weekly options on the S&P 500 equity index. All reported returns data is expressed in weekly units. The volatility is annualized and computed from intraday returns, as follows: Weekly ( intraday ) variance = i = 1 N { z t , i } 2 , where z t , i = log ( F t , i N ) log ( F t , i 1 N ) , where N is the number of return observations in a trading week (i.e., N = 12 × 6.5 × 5 ). We compute the equal weighted return across the five crude oil futures contracts (i.e., WTI, Brent, Dubai, Heating Oil, and RBOB Gasoline), as z t + Δ basket , oil i = 1 5 1 5 ( F t + Δ i F t i 1 ) . The realized futures returns volatility are based on return observations over five-minute intervals (except for Dubai for which intraday data is unavailable). The source of the matched futures and options data is Chicago Mercantile Exchange. The standard deviation is shown in the column “SD”, and Acf is the first-order autocorrelation. Our analysis covers a total of 341 weekly option expiration cycles.
Table 3. Summary statistics of weekly predictors constructed from options and intraday data.
Table 3. Summary statistics of weekly predictors constructed from options and intraday data.
MeanSDBootstrap
Block
NW[ p ]Min.Max.AcfSkewnessKurtosis 1 { a > 0 }
(%)
LowerUpper
κ 2 , t Q , oil 0.00420.0105⌊0.00260.0065⌋0.000.00030.13280.668.1280.3
κ 3 , t Q , oil −0.000120.00185⌊−0.000330.00001⌋0.24−0.032110.005810.17−15.10263.632
κ 4 , t Q , oil 0.000320.00360⌊0.000030.00083⌋0.20−0.001210.064420.2116.92299.9
Skewness t Q , oil −0.17730.47⌊−0.26−0.09⌋0.00−1.371.870.620.431.2732
Kurtosis t Q , oil 6.22.4⌊5.86.5⌋0.002.826.20.433.5621.15
Variance t P , oil , 5 min 0.00370.0306⌊0.00130.0077⌋0.0590.00020.55910.1217.76322.81
Skewness t P , oil , 5 min −0.12091.25⌊−0.25−0.02⌋0.08−9.326.63−0.02−0.6012.5144
κ 2 , t Q , equity 0.000650.00116⌊0.000470.00089⌋0.000.000060.013590.786.8761.56
κ 3 , t Q , equity −0.0000270.00010⌊−0.00005−0.00001⌋0.00−0.001070.000000.71−8.5479.770
κ 4 , t Q , equity 0.00000240.00001⌊0.000000.00000⌋0.020.000000.000210.5113.20200.50
Skewness t Q , equity −1.30.5⌊−1.3−1.2⌋0.00−2.7−0.10.61−0.660.220
Kurtosis t Q , equity 7.44.8⌊6.68.2⌋0.002.148.30.513.3518.49
Variance t P , equity , 5 min 0.000480.00116⌊0.000320.00071⌋0.000.000020.014430.798.3384.51
Skewness t P , equity , 5 min −0.170.66⌊−0.24−0.11⌋0.00−4.781.710.05−2.4114.5540
EIA t stock growth × 100 −0.100.56⌊−0.18−0.02⌋0.03−1.722.260.370.611.4441
EIA t production growth × 100 0.112.23⌊−0.020.23⌋0.19−13.9812.31−0.20−0.8516.0245
EIA t import growth × 100 −0.109.64⌊−0.420.28⌋0.70−28.3832.60−0.530.220.5247
The time frame of our study is confined to 12 August 2016, to 24 February 2023, as we use the weekly (7DTE) options on crude oil futures and the S&P 500 equity index. The source of the matched futures and options data is Chicago Mercantile Exchange and Chicago Board of Trade. All data is expressed in weekly units. The 7DTE higher-order cumulants, κ n , t Q , are constructed based on both WTI crude oil options and the S&P 500 equity options. The realized variances are based on return observations over five-minute intervals. Our estimates imply an annualized P ( Q ) measure oil volatility of 39.2% (28.0%). Acf is the first-order autocorrelation. Our analysis covers a total of 341 weekly option expiration cycles.
Table 4. Short-maturity (weekly) predictive regression results for WTI crude oil futures returns.
Table 4. Short-maturity (weekly) predictive regression results for WTI crude oil futures returns.
Predictor
(Univariate)
ConstantNW[p] b NW[p] R ¯ 2
(%)
CORRPredict
(Yes or No)
Panel A: Three 7DTE higher-order risk-neutral return cumulants from oil market
κ 2 , t Q , o i l 0.000.450.070.88−0.280.01No
κ 3 , t Q , o i l 0.000.55−3.570.010.98−0.11Yes
κ 4 , t Q , o i l 0.000.520.950.010.050.06Yes
Panel B: Other predictors from oil markets
Skewness t Q , o i l 0.000.64−0.010.40−0.11−0.04No
Kurtosis t Q , o i l −0.010.480.000.210.040.06No
Variance t P , o i l , 5 m i n 0.000.560.080.07−0.130.04Yes
Skewness t P , o i l , 5 m i n 0.000.57−0.000.33−0.10−0.05No
Panel C: Predictors from equity markets
κ 2 , t Q , e q u i t y 0.010.11−4.750.330.59−0.09No
κ 3 , t Q , e q u i t y 0.000.1188.920.141.880.15No
κ 4 , t Q , e q u i t y 0.000.32−295.60.300.13−0.07No
Skewness t Q , e q u i t y −0.010.16−0.010.020.40−0.08Yes
Kurtosis t Q , e q u i t y −0.000.580.000.140.140.07No
Variance t P , e q u i t y , 5 m i n 0.010.09−6.20.231.22−0.12No
Skewness t P , e q u i t y , 5 m i n 0.000.490.000.81−0.290.01No
Panel D: Predictors from the Energy Information Administration (EIA)
EIA t s t o c k g r o w t h 0.000.79−1.390.031.51−0.15Yes
EIA t p r o d u c t i o n g r o w t h 0.000.48−0.060.69−0.25−0.02No
EIA t i m p o r t g r o w t h 0.000.50−0.030.290.01−0.06No
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on WTI crude oil. The results are presented in the form of univariate predictive regressions expressed as follows: F t + Δ F t 1 w e e k l y W T I f u t u r e s r e t u r n = constant + b × Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The correlation between the predictor and F t + Δ F t 1 is reported as CORR. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday. The data is collected from 12 August 2016, to 24 February 2023, covering a total of 341 option expiration cycles.
Table 5. Evidence from the bivariate predictive regressions (WTI crude oil futures).
Table 5. Evidence from the bivariate predictive regressions (WTI crude oil futures).
Alternative PredictorConstantNW[ p ]Predictor Is
κ 3 , t Q , oil
Alternative
Predictor
R ¯ 2
(%)
Joint p R os 2
(%)
b cubic NW[ p ] b alter NW[ p ]
κ 2 , t Q , o i l 0.00 0.15−5.410.01−0.520.281.230.010.39
κ 4 , t Q , o i l 0.00 0.48−9.460.00−3.400.011.580.002.78
Skewness t Q , o i l 0.000.64−3.430.010.000.650.740.020.52
Kurtosis t Q , o i l −0.010.46−3.560.010.000.211.020.011.25
Variance t P , o i l , 5 m i n 0.000.32−16.10.00−0.820.003.470.004.10
Skewness t P , o i l , 5 m i n 0.000.64−3.450.010.000.530.790.031.08
κ 2 , t Q , e q u i t y 0.010.11−3.970.01−5.520.261.880.040.90
κ 3 , t Q , e q u i t y 0.000.13−4.080.0297.10.103.280.041.23
κ 4 , t Q , e q u i t y 0.000.38−3.840.02−357.90.211.320.041.21
Skewness t Q , e q u i t y −0.010.12−3.650.01−0.010.021.460.000.96
Kurtosis t Q , e q u i t y −0.000.50−3.630.010.000.121.190.011.00
Variance t P , e q u i t y , 5 m i n 0.010.10−3.980.02−6.880.192.530.051.08
Skewness t P , e q u i t y , 5 m i n 0.000.56−3.570.010.000.740.720.031.18
EIA t s t o c k g r o w t h −0.000.90−3.970.00−0.780.013.370.000.95
EIA t p r o d u c t i o n g r o w t h 0.000.59−3.550.01−0.060.680.750.031.12
EIA t i m p o r t g r o w t h 0.000.62−3.640.01−0.040.231.080.031.25
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on WTI crude oil. The results are presented in the form of bivariate predictive regressions expressed as follows: F t + Δ F t 1 w e e k l y W T I f u t u r e s r e t u r n = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday. The out-of-sample performance is reported as R o s 2 = 1 t = 1 T z t z ^ t , f u l l 2 t = 1 T z t z ^ t , n e s t e d 2 .
Table 6. Evidence from the bivariate predictive regressions (equal weight basket across five oil futures contracts).
Table 6. Evidence from the bivariate predictive regressions (equal weight basket across five oil futures contracts).
Alternative PredictorConstantPredictor Is
κ 3 , t Q , oil
Alternative
Predictor
NW[ p ] R ¯ 2
(%)
Joint p R os 2
NW[ p ] b cubic NW[ p ] b alter
κ 2 , t Q , o i l 0.000.15−4.920.01−0.270.521.960.000.51
κ 4 , t Q , o i l 0.000.33−6.520.00−1.460.212.000.001.56
Skewness t Q , o i l 0.000.40−3.910.000.000.751.790.000.25
Kurtosis t Q , o i l −0.010.47−3.990.000.000.172.200.001.13
Variance t P , o i l , 5 m i n 0.000.25−11.100.00−0.470.003.100.002.21
Skewness t P , o i l , 5 m i n 0.000.44−3.860.000.000.391.940.000.88
κ 2 , t Q , e q u i t y 0.010.03−4.350.00−4.970.193.230.000.68
κ 3 , t Q , e q u i t y 0.000.06−4.420.0080.790.064.450.000.97
κ 4 , t Q , e q u i t y 0.000.24−4.200.00−268.280.182.330.001.06
Skewness t Q , e q u i t y −0.010.11−4.070.00−0.010.012.810.000.60
Kurtosis t Q , e q u i t y 0.000.48−4.060.000.000.092.620.000.95
Variance t P , e q u i t y , 5 m i n 0.010.04−4.340.00−5.840.123.770.000.95
Skewness t P , e q u i t y , 5 m i n 0.000.39−4.010.000.000.651.860.000.76
EIA t s t o c k g r o w t h 0.000.60−4.190.00−0.880.092.850.000.81
EIA t p r o d u c t i o n g r o w t h 0.000.38−3.990.00−0.030.771.800.000.82
EIA t i m p o r t g r o w t h 0.000.40−4.080.00−0.040.102.490.001.04
We use regression analysis based on an equal weight basket of five crude oil futures: (i) WTI, (ii) Brent, (iii) Dubai, (iv) Heating Oil, and (v) RBOB Gasoline, as follows: z t + Δ b a s k e t , o i l = i = 1 5 1 5 ( F t + Δ i F t i 1 ) . The return of the crude oil basket serves to smooth idiosyncratic variations in the oil prices. The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on WTI crude oil. The results are presented in the form of bivariate predictive regressions expressed as i = 1 5 1 5 ( F t + Δ i F t i 1 ) e q u a l w e i g h t b a s k e t = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ i F t i 1 represents the return of the aforementioned individual futures position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The out-of-sample performance is reported as R o s 2 = 1 t = 1 T z t z ^ t , f u l l 2 t = 1 T z t z ^ t , n e s t e d 2 .
Table 7. Evidence from the bivariate predictive regressions (Brent oil futures).
Table 7. Evidence from the bivariate predictive regressions (Brent oil futures).
Alternative PredictorConstantPredictor Is
κ 3 , t Q , oil
Alternative
Predictor
NW[ p ] R ¯ 2
(%)
Joint p
NW[ p ] b cubic NW[ p ] b alter
κ 2 , t Q , o i l 0.00 0.11−3.590.08−0.420.180.460.21
κ 4 , t Q , o i l 0.00 0.28−9.970.00−4.530.002.110.00
Skewness t Q , o i l 0.000.51−1.900.20−0.010.380.230.22
Kurtosis t Q , o i l 0.000.69−2.120.170.000.340.200.28
Variance t P , o i l , 5 m i n 0.000.16−15.590.00−0.880.004.300.00
Skewness t P , o i l , 5 m i n 0.000.41−2.020.180.000.440.120.32
κ 2 , t Q , e q u i t y 0.000.17−2.300.17−2.370.610.310.37
κ 3 , t Q , e q u i t y 0.000.11−2.400.1852.120.341.010.28
κ 4 , t Q , e q u i t y 0.000.26−2.260.17−172.50.490.210.32
Skewness t Q , e q u i t y −0.010.09−2.210.15−0.010.011.000.01
Kurtosis t Q , e q u i t y 0.000.54−2.190.150.000.130.590.15
Variance t P , e q u i t y , 5 m i n 0.000.09−2.360.17−3.940.370.820.28
Skewness t P , e q u i t y , 5 m i n 0.000.32−2.160.160.000.360.190.26
EIA t s t o c k g r o w t h 0.000.55−2.300.12−0.750.140.720.12
EIA t p r o d u c t i o n g r o w t h 0.000.36−2.130.17−0.020.870.020.37
EIA t i m p o r t g r o w t h 0.000.38−2.220.15−0.050.050.870.06
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on Brent oil futures. The results are presented in the form of bivariate predictive regressions expressed as follows: F t + Δ F t 1 B r e n t = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday.
Table 8. Evidence from the bivariate predictive regressions (Dubai oil futures).
Table 8. Evidence from the bivariate predictive regressions (Dubai oil futures).
Alternative PredictorConstantPredictor Is
κ 3 , t Q , oil
Alternative
Predictor
NW[ p ] R ¯ 2
(%)
Joint p
NW[ p ] b cubic NW[ p ] b alter
κ 2 , t Q , o i l 0.010.04−10.450.00−0.820.129.560.00
κ 4 , t Q , o i l 0.000.45−9.430.01−1.070.567.810.00
Skewness t Q , o i l 0.000.36−7.740.000.000.487.800.00
Kurtosis t Q , o i l 0.000.59−7.570.000.000.247.950.00
Variance t P , o i l , 5 m i n 0.000.40−11.480.00−0.260.258.100.00
Skewness t P , o i l , 5 m i n 0.000.58−7.380.000.000.318.060.00
κ 2 , t Q , e q u i t y 0.010.00−8.350.00−10.460.0013.790.00
κ 3 , t Q , e q u i t y 0.010.05−8.250.00123.50.0013.690.00
κ 4 , t Q , e q u i t y 0.000.24−8.040.00−595.80.0010.220.00
Skewness t Q , e q u i t y −0.010.08−7.680.00−0.010.019.160.00
Kurtosis t Q , e q u i t y −0.010.28−7.680.000.000.049.050.00
Variance t P , e q u i t y , 5 m i n 0.010.02−8.120.00−8.790.0012.030.00
Skewness t P , e q u i t y , 5 m i n 0.000.68−7.550.000.000.677.860.00
EIA t s t o c k g r o w t h 0.000.71−7.760.00−0.810.258.590.00
EIA t p r o d u c t i o n g r o w t h 0.000.51−7.580.000.010.917.710.00
EIA t i m p o r t g r o w t h 0.000.51−7.650.00−0.040.248.240.00
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on Dubai oil. The results are presented in the form of bivariate predictive regressions expressed as follows: F t + Δ F t 1 Dubai = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday.
Table 9. Evidence from the bivariate predictive regressions (Heating Oil futures).
Table 9. Evidence from the bivariate predictive regressions (Heating Oil futures).
Alternative PredictorConstantPredictor Is
κ 3 , t Q , oil
Alternative
Predictor
NW[ p ] R ¯ 2
(%)
Joint p
NW[ p ] b cubic NW[ p ] b alter
κ 2 , t Q , o i l 0.000.20−3.910.01−0.150.630.850.00
κ 2 , t Q , o i l 0.000.26−6.410.03−1.760.201.090.00
Skewness t Q , o i l 0.000.34−3.320.000.000.880.800.00
Kurtosis t Q , o i l 0.000.54−3.360.000.000.211.100.00
Variance t P , o i l , 5 m i n 0.000.23−9.510.00−0.400.021.630.00
Skewness t P , o i l , 5 m i n 0.000.38−3.160.000.000.171.130.00
κ 2 , t Q , e q u i t y 0.000.08−3.500.00−2.000.461.010.00
κ 3 , t Q , e q u i t y 0.000.11−3.560.0038.390.211.310.00
κ 4 , t Q , e q u i t y 0.000.25−3.340.0024.170.870.820.00
Skewness t Q , e q u i t y −0.010.24−3.430.00−0.010.041.520.00
Kurtosis t Q , e q u i t y 0.000.72−3.420.000.000.161.230.00
Variance t P , e q u i t y , 5 m i n 0.000.08−3.520.00−2.770.281.180.00
Skewness t P , e q u i t y , 5 m i n 0.000.26−3.380.000.000.580.890.00
EIA t s t o c k g r o w t h 0.000.59−3.640.00−1.240.032.550.00
EIA t p r o d u c t i o n g r o w t h 0.000.32−3.360.00−0.030.820.820.00
EIA t i m p o r t g r o w t h 0.000.33−3.450.00−0.040.171.430.00
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on Heating Oil. The results are presented in the form of bivariate predictive regressions expressed as follows: F t + Δ F t 1 H e a t i n g O i l = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday.
Table 10. Evidence from the bivariate predictive regressions (RBOB Gasoline futures).
Table 10. Evidence from the bivariate predictive regressions (RBOB Gasoline futures).
Alternative PredictorConstantPredictor Is
κ 3 , t Q , oil
Alternative
Predictor
NW[ p ] R ¯ 2
(%)
Joint p
NW[ p ] b cubic NW[ p ] b alter
κ 2 , t Q , o i l 0.000.68−1.260.630.580.401.030.00
κ 4 , t Q , o i l 0.000.352.670.473.460.051.260.00
Skewness t Q , o i l 0.000.38−3.170.000.000.620.480.00
Kurtosis t Q , o i l −0.010.32−3.310.000.000.081.100.00
Variance t P , o i l , 5 m i n 0.000.35−2.860.550.030.910.440.00
Skewness t P , o i l , 5 m i n 0.000.37−3.300.000.000.910.430.00
κ 2 , t Q , e q u i t y 0.010.08−3.650.00−4.480.471.180.00
κ 3 , t Q , e q u i t y 0.010.05−3.820.0092.80.182.660.00
κ 4 , t Q , e q u i t y 0.000.24−3.500.00−239.40.440.710.00
Skewness t Q , e q u i t y −0.010.43−3.380.00−0.010.140.830.00
Skewness t Q , e q u i t y 0.000.57−3.400.000.000.161.080.00
Variance t P , e q u i t y , 5 m i n 0.010.04−3.730.00−6.830.272.150.00
Skewness t P , e q u i t y , 5 m i n 0.000.33−3.370.000.000.340.660.00
EIA t s t o c k g r o w t h 0.000.40−3.350.00−0.130.840.440.00
EIA t p r o d u c t i o n g r o w t h 0.000.34−3.310.00−0.080.570.510.00
EIA t i m p o r t g r o w t h 0.000.36−3.420.00−0.050.110.980.00
The data analyzed involves a predictor based on oil and equity markets to forecast the returns of a fully collateralized long position in weekly futures contract on RBOB Gasoline. The results are presented in the form of bivariate predictive regressions expressed as follows: F t + Δ F t 1 R B O B G a s o l i n e = constant + b c u b i c κ 3 , t Q , o i l + b a l t e r Predictor t + ϵ t + Δ , where F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday.
Table 11. Predictive regression results allowing for COVID-19 dummy variables.
Table 11. Predictive regression results allowing for COVID-19 dummy variables.
Panel A: F t + Δ F t 1 = constant + b s q u a r e d κ 2 , t Q , o i l + b c u b i c κ 3 , t Q , o i l + b q u a r t i c κ 4 , t Q , o i l + ϵ t + Δ
constantNW[p] b s q u a r e d NW[p] b c u b i c NW[p] b q u a r t i c NW[p] R ¯ 2 (%)
WTI0.000.33−0.150.87−9.130.01−2.910.421.31
Oil basket0.000.20−0.130.84−6.230.01−1.020.651.74
Panel B: Regression with dummy variables
WTIconstantNW[p] b s q u a r e d ,   p r e - c o v i d NW[p] b s q u a r e d ,   c o v i d NW[p] b s q u a r e d ,   p o s t - c o v i d NW[p]
−0.010.127.310.07−0.210.826.150.01
b c u b i c ,   p r e - c o v i d NW[p] b c u b i c ,   c o v i d NW[p] b c u b i c ,   p o s t - c o v i d NW[p]
53.070.62−7.960.04−27.940.00
b q u a r t i c ,   p r e - c o v i d NW[p] b q u a r t i c ,   c o v i d NW[p] b q u a r t i c ,   p o s t - c o v i d NW[p] R ¯ 2 (%)
−188.340.48−1.960.59−109.460.052.22
Oil basketconstantNW[p] b s q u a r e d ,   p r e - c o v i d NW[p] b s q u a r e d ,   c o v i d NW[p] b s q u a r e d ,   p o s t - c o v i d NW[p]
−0.010.305.660.18−0.160.804.510.04
b c u b i c ,   p r e - c o v i d NW[p] b c u b i c ,   c o v i d NW[p] b c u b i c ,   p o s t - c o v i d NW[p]
83.030.39−4.950.05−24.380.01
b q u a r t i c ,   p r e - c o v i d NW[p] b q u a r t i c ,   c o v i d NW[p] b q u a r t i c ,   p o s t - c o v i d NW[p] R ¯ 2 (%)
−152.420.59−0.170.94−81.530.102.33
The dummy variables 1 p r e - c o v i d , 1 c o v i d , and 1 p o s t - c o v i d correspond to the time period before the COVID-19 pandemic (12 August 2016, to January 31, 2020), during the pandemic (7 February 2020, to 5 June 2020), and after the pandemic (12 June 2020, to 23 February 2023). The data analyzed involves predictors based on the oil options market to forecast the returns of a fully collateralized long position in weekly futures contract on WTI crude oil or the basket of oil futures. The results are presented based on the predictive regressions expressed, as follows:
F t + Δ F t 1 W T I o r b a s k e t = constant + ( b s q u a r e d , p r e - c o v i d 1 p r e - c o v i d + b s q u a r e d , c o v i d 1 c o v i d + b s q u a r e d , p o s t - c o v i d 1 p o s t - c o v i d ) κ 2 , t Q , o i l     + ( b c u b i c ,   p r e - c o v i d 1 p r e - c o v i d + b c u b i c ,   c o v i d 1 c o v i d + b c u b i c ,   p o s t - c o v i d 1 p o s t - c o v i d ) κ 3 , t Q , o i l + ( b q u a r t i c ,   p r e - c o v i d 1 p r e - c o v i d + b q u a r t i c ,   c o v i d 1 c o v i d + b q u a r t i c ,   p o s t - c o v i d 1 p o s t - c o v i d ) κ 4 , t Q , o i l + ϵ t + Δ .
Here F t + Δ F t 1 represents the return of the aforementioned position. Reported are the b estimates, as well as two-sided p-values for the null hypothesis b = 0 , denoted by NW[p], based on the procedure in Newey and West with the lag automatically selected. The adjusted R 2 , denoted by R ¯ 2 , is also presented as a measure of goodness-of-fit, expressed as a percentage. The predictive regressions are conducted over weekly options expiration cycles, spanning from Friday to Friday.
Table 12. Investible trading strategy.
Table 12. Investible trading strategy.
J
Weeks
Full Sample
12 August 2016, to 23 February 2023
Weekly Returns (%)
Subsample
12 August 2016, to 5 June 2020
Weekly Returns (%)
Subsample
5 June 2020, to 23 February 2023
Weekly Returns (%)
Mean
Return
Block
Bootstrap
NW[ p ] 1 z > 0
(%)
Mean
Return
Block
Bootstrap
NW[ p ] 1 z > 0
(%)
Mean
Return
Block
Bootstrap
NW[ p ] 1 z > 0
(%)
Sharpe
Ratio
80.13⌊−0.390.67⌋0.6854−0.24⌊−1.130.34⌋0.58480.63⌊0.041.28⌋0.10630.90
70.09⌊−0.420.61⌋0.7753−0.37⌊−1.240.18⌋0.39470.71⌊0.141.34⌋0.06631.01
60.33⌊−0.170.85⌋0.28550.01⌊−0.830.54⌋0.98500.77⌊0.161.42⌋0.05621.09
50.22⌊−0.310.75⌋0.4855−0.01⌊−0.860.53⌋0.97500.54⌊−0.161.28⌋0.17610.76
40.04⌊−0.470.56⌋0.8952−0.19⌊−1.030.35⌋0.66470.36⌊−0.291.02⌋0.33600.51
30.14⌊−0.380.65⌋0.65530.05⌊−0.780.58⌋0.90510.25⌊−0.400.91⌋0.51560.35
20.07⌊−0.440.60⌋0.8354−0.13⌊−0.980.43⌋0.76510.34⌊−0.301.01⌋0.36580.48
10.09⌊−0.410.64⌋0.7754−0.13⌊−0.960.46⌋0.76510.40⌊−0.261.09⌋0.29580.57
Our approach follows a quantitative strategy where we take a long position in the five oil futures contracts if the momentum in third risk-neutral return cumulant is negative. Otherwise, we adopt a short position in the five oil commodities. The momentum is calculated over various time intervals. The profits or losses are measured over the next seven days, using an equal weighting method. We present the average weekly return (%), along with the 90% confidence intervals obtained using the block bootstrap technique. Additionally, we report the two-sided p-values, denoted by NW[p], which reflect the hypothesis of a zero mean return, based on the Newey-West procedure with an automatic lag selection. Our strategy is adjusted weekly, in line with options expiration cycles, from Friday to Friday. We collect data from 12 August 2016, to 24 February 2023, capturing a total of 341 options expiration cycles. The entries under the column “Sharpe Ratio” reflect the computed annualized Sharpe Ratio of the strategy.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bakshi, G.; Gao, X.; Zhang, Z. What Insights Do Short-Maturity (7DTE) Return Predictive Regressions Offer about Risk Preferences in the Oil Market? Commodities 2024, 3, 225-247. https://doi.org/10.3390/commodities3020014

AMA Style

Bakshi G, Gao X, Zhang Z. What Insights Do Short-Maturity (7DTE) Return Predictive Regressions Offer about Risk Preferences in the Oil Market? Commodities. 2024; 3(2):225-247. https://doi.org/10.3390/commodities3020014

Chicago/Turabian Style

Bakshi, Gurdip, Xiaohui Gao, and Zhaowei Zhang. 2024. "What Insights Do Short-Maturity (7DTE) Return Predictive Regressions Offer about Risk Preferences in the Oil Market?" Commodities 3, no. 2: 225-247. https://doi.org/10.3390/commodities3020014

APA Style

Bakshi, G., Gao, X., & Zhang, Z. (2024). What Insights Do Short-Maturity (7DTE) Return Predictive Regressions Offer about Risk Preferences in the Oil Market? Commodities, 3(2), 225-247. https://doi.org/10.3390/commodities3020014

Article Metrics

Back to TopTop