Uncovering the Distribution of Option Implied Risk Aversion

This paper explores the dynamics of risk aversion of a representative agent with an iso-elastic utility function. In contrast to most of the existing literature, we estimate the coefficient of relative risk aversion from option prices. To do this, we transform the risk-neutral density function obtained from a cross-section of option prices to an objective distribution function that reflects individuals’ risk aversion through a CRRA utility function. The dynamics of the relative risk-aversion coefficient are obtained by repeating the same estimation procedure over rolling windows. This procedure uncovers strong variation in risk aversion over time. We also propose a simulation procedure to construct confidence intervals for the risk-aversion coefficient in each period. We assess the robustness of these confidence intervals under different assumptions on the data generating process of stock prices. The results imply a strong influence of volatility on the variation of risk aversion. In an empirical application, we compare the forecasting performance of our approach based on our risk-aversion estimates against the method proposed in [1]. Overall, we find that our simulation based approach obtains better forecasting results than bootstrap methods.


Introduction
Stock options are priced using risk-neutral expectations of the payoff of the underlying asset. These risk-neutral expectations usually incorporate in the option pricing formula a flexible specification of the density function for describing the price of the underlying asset. Stock prices are priced, instead, by discounting the expected value of the future cash flows of the asset under an objective distribution function. Both approaches incorporate individuals' attitude towards risk. The risk-neutral approach embeds individuals' risk aversion in the risk-neutral probability distribution of the risky payoffs, that are discounted using the risk-free rate. In contrast, the objective approach uses objective probabilities, inferred from historical prices or alternative parametric methods, and discount future payoffs using a stochastic discount factor that incorporates a risk premium over the risk-free rate.
There are different approaches to connect the risk-neutral valuation with the objective valuation. A convenient formulation that relies on an utilitarian approach is to assume a parametric relationship between the distribution functions through a utility function. The analytical tractability of power utility functions offers a convenient modelling device for describing individuals' preferences. In this paper we exploit this parametric relationship between the risk-neutral and objective distribution functions and the power utility function to extract consistent estimates of the coefficient of relative risk aversion of a representative agent. There are two reasons for extracting risk aversion from risk-neutral valuations of option prices. First, risk-aversion estimates are economically important and provide relevant information for understanding the performance and dynamics of financial markets. Second, the empirical finance literature confirms the significant improvement in the ability of parametric models for forecasting stock prices that incorporate reliable estimates of risk aversion.
The link between expected returns and investors' risk-preferences in financial markets is the foundation of research concerning the risk-return relationship.
Several methods have been developed over time in order to determine the risk-premia underlying asset markets. Popular approaches are, for example, the standard and consumption-based capital asset pricing model. These pioneering asset pricing methods are usually estimated based on asset price returns, treasury yields and consumption data. Although these models are still the building blocks of current asset pricing models and the foundation of today's financial economics, their success in estimating the underlying risk premia in financial markets has been limited. A prominent example is the seminal article by [2] on the equity premium puzzle. These authors retrieve extremely high values of the risk aversion coefficient in order for the asset pricing model to be compatible with observed asset prices. These values are not economically justified. Other studies have tried to correct these estimates by imposing different assumptions on individuals' preferences. For example, [3] makes strong assumptions on the flexibility of the underlying utility function. This author finds estimates of the risk aversion coefficient equal to one and consistent with the presence of a log-utility function for modelling the preferences of a representative agent. [4] find values of the risk aversion coefficient in the range of two. [5] in their M. Kyriacou [2], discussed above, find values of the risk aversion coefficient around 55 that are difficult to reconcile with practitioners' and financial economists' expectations on the underlying risk aversoin in financial markets. [6] propose a recursive utility function for modelling individuals' preferences and find under this novel approach a risk aversion coefficient that oscillates between 0.4 and 1.4. It is important to remark that the Epstein-Zin utility function is able to disentangle individuals' risk aversion from the rate of intertemporal substitution. [7] propose habit formation models in which individuals' preferences depend on a reference point under which perceived utility of the asset is zero. This approach provides more variability in the estimates of the risk aversion coefficient. Thus, these authors find a risk aversion coefficient between 0 and 12. [8] n an asset pricing context find values in the range of 40 -50 consistent with those found in [2]. [9] and [10] find values between 5.4 and 11.9 and smaller than 3, respectively. More recently, [11] derives a coefficient of risk aversion of 12.7 and [12] using option prices close to their asset price returns and consumption data find a value of the relative risk aversion coefficient equal to 3.52. These authors exploit option prices to extract implied volatility measures for estimating the risk aversion coefficient. Our aim is to follow the third approach, and more specifically, [1]. One of our objectives is to gain insight on the variation of risk-aversion estimates over time and, hence, through the business cycle. [1] provide quantitative methods to evaluate risk-aversion estimates. These authors introduce two different approaches. First, a Monte-Carlo simulation method to retrieve the standard error Journal of Mathematical Finance of risk-aversion estimates. Second, a Bootstrapping method to capture the influence of actual data and potential model misspecification on the reliability of risk-aversion estimates. Despite the appealing of both approaches, these methods face challenges due to the limited availability of asset prices and option prices.
To overcome this issue, we explore different parametric data generating processes for constructing the risk-neutral distribution function for describing the dynamics of option and asset prices. Variation in the estimates of risk aversion is obtained by constructing rolling windows of 120 monthly observations and obtaining estimates of risk aversion for each window. This number is sufficiently large to avoid too much variation in risk-aversion estimates at the same time small enough to capture changes in market conditions that trigger shifts in risk aversion. Confidence intervals of the risk-aversion estimates can be constructed by simulating the prices from a parametric data generating process.
The simulated distribution of risk-aversion estimates provides valuable insights on the variation of the risk-aversion coefficient. Furthermore, in the proposed set-up the underlying asset price process allows for different characteristics starting with the base case of a Geometric Brownian Motion and its extension to accommodate stochastic volatility processes. The simulation results are compared with the classic approach followed by [1] nd [27] These authors derive their risk-aversion estimates based on past series of risk-neutral densities and realised asset prices.
Our extensive simulation exercise provides the following insights. First, our results suggest that there is more variation in the risk-aversion coefficient than implied by [1]. Not surprisingly, the volatility component has a major impact on the coefficient of relative risk aversion. The choice of parametric model for describing asset prices also has a major effect on the estimates of risk aversion for a representative agent. Thus, we note that imposing a lognormal density on asset prices produces greater variation in the coefficient of risk aversion compared to the mixture-lognormal and generalised beta distribution of the second kind.
Lastly, we also perform a forecasting exercise to assess the relative out-of-sample ability of the different data generating processes for estimating the risk-aversion coefficient for forecasting stock prices. The forecasting performance is evaluated using the subjective density derived from the average simulated risk-aversion estimate. The findings are compared against the classic approach obtained from [1]. Overall, we find that the simulation based approach obtains better forecasting results than bootstrap methods. This fact is particularly apparent in the period after the sub-prime crisis.
The paper is structured in six sections. Section 2 introduces the considered density types and transformation from risk-neutral to subjective density. Section

Option Implied Densities
The risk-aversion estimates obtained in this section are retrieved from option prices. This approach goes back to the density forecasting literature for option pricing as, for example, [27]. The foundations to this literature on option pricing is due to [28]. These authors show that when there are no-arbitrage opportunities and options are available for all strikes K, a call-option c can be evaluated by: where S is the stock price, r the risk-free rate, T the time to maturity and Q f a unique risk-neutral density. All parameters can be observed in the market except for the risk-neutral density. Assuming a parametric form of Q f , its parameters To reflect the dynamics of the underlying asset, the risk-neutral density needs to be transformed into a subjective counterpart density function that incorporates individuals' risk-aversion profile. We follow the transformation in [1], who relate the risk-neutral Q f to a subjective density function P f by making an assumption about the preferences of a representative investor. This parametric relationship goes back to [11]. In our study, we assume a power utility function The constant relative risk-aversion parameter γ , inherited from the power utility and responsible for the transformation into the subjective density, is the central aspect of the subsequent analysis. As in [27], its value is determined by a series of M densities, maximising the log-likelihood of the realization of asset  [29] introduces a mixture of lognormal densities, MLN herafter, to evaluate option prices. Following [27], this distribution is defined as a weighted combination of two lognormal densities. A European call price is evaluated as the sum of two [30]models:

Mixtures of Lognormal Densities
The risk-neutrality constraint for the MLN reduces the amount of free parameters to four and by stating that: with F the forward price. Transforming the density into its subjective counterpart as in (3) leads to another mixture-lognormal density. The subjective parameters become: The lognormal density (LN) is a special case of the above for 1 w = . As a result, we do not describe this distribution separately.

Generalized Beta Densities
The generalized beta density of the second kind, GB2 hereafter, was introducted by [31] and has been applied to option prices by [27]. They show that a European call option price can be evaluated by:  and F β the distribution function of a random variable following a beta distribution.
The risk-neutrality constraint reduces the amount of free parameters from four to three: The transformation of the GB2 risk-neutral density as in (3) leads a different type of GB2 density function. The subjective parameters are then

Risk Aversion Evaluation
One of the biggest challenges in the evaluation of risk-aversion estimates is the  (4). This approach requires, however, access to a sufficiently long dataset. In situations where the time horizon is too short, estimates might be misleading and will not accurately reflect the risk preferences of the representative agent.
[1] introduce a Monte-Carlo simulation method and a Bootstrapping procedure to gain better insight into the variation of risk-aversion estimates. However, their approach does not overcome the problem of single price realisations. More importantly, risk aversion is only estimated from past data, which may accurately reflect individuals' past preferences but not necessarily those relevant for the forthcoming period.
We propose an alternative Monte-Carlo simulation exercise that addresses all these problems. Our simulation setting is structured in a way that only considers the most recent risk-neutral dynamics and simulates potential asset price outcomes based on different stochastic processes. This allows us to estimate risk aversion by taking into account current market dynamics. Importantly, it is possible to simulate several asset price outcomes for the same time period providing a more detailed outlook of the anticipated variation of the underlying risk aversion.

Simulation Framework
In our simulation framework current market dynamics are captured by a set of assumptions on risk-neutral and subjective processes. Given these, we can simulate a whole range of potential asset price outcomes for a specific period and use these to estimate risk aversion. Simulation-based risk-aversion estimates provide better insight into the full spectrum of potential risk-aversion estimates. Accordingly, the utility function can be validated. The Monte-Carlo simulation framework consists of the following five steps 1 : 1) Make assumptions about the subjective process that includes the utility function of the representative agent and the risk-neutral process.
2) Obtain options from the risk-neutral process and fit a risk-neutral density into option prices. Since, the risk-neutral process does not change, estimated risk-neutral parameter values do not change for the time horizon of M periods.
3) Generate M asset price outcomes based on the subjective process for M risk-neutral densities.
4) Estimate the risk-aversion parameter γ over the estimation horizon of M periods based on (4). 1 Note that for the empirical application, Step 2 does not draw option prices from the risk-neutral process since options can be observed directly in the market.

Likelihood Criteria
Log-likelihood based tests on the asset price outcomes can be useful in providing insights on the quality of the retrieved risk-aversion values, as used in the work of [32], [27] and [24] to compare the quality of their density methods on equity indices. The log-likelihood function is defined as: where f is the density function, m the density type and i the number of densities.
The log-likelihood function serves two different testing purposes. First, we construct a log-likelihood ratio test to assess statistically whether the additional risk-aversion parameter γ significantly improves the forecasting performance of the densities.
The corresponding test statistic is defined as: Second, we test whether the two competing methods, denoted as m and n, provide equal expected log-likelihoods as described: This test is proposed by [33], AG hereafter, and has been applied to option based densities by [24]. The AG test is an out-of sample test, which uses a rolling window estimation scheme to rank density forecasts. Its loss function is based on a weighted Likelihood ratio test and offers the advantage that it can be robust to heterogeneous data sets, while it can be applied to densities obtained from both parametric and non-parametric models.
Under the null hypothesis 0 H , it follows that

Market Assumptions and Estimation
We impose two different parametric forms to model market dynamics. First, market dynamics are assumed to follow a Geometric Brownian Motion as in [34]. Second, we incorporate stochastic volatility into the model by implementing  [35]. In both cases, we assume that the volatility parameters are derived from option prices including also the subjective dynamics. Only the drift parameter is changed between specifications in order to account for the riskpremium of the subjective approach. This transformation is similar to the risk-premium transformation in [24]. The drift term of the subjective process is estimated using historical asset price returns under the assumption of lognormality of stock prices. Finally, we illustrate how expected risk aversion can be obtained by taking into account the risk-neutral and subjective dynamics.

Geometric Brownian Motion
In the first specification, we assume that markets follow a Geometric Brownian Motion as in [34]: with µ the drift term and σ the constant volatility. Option prices can be evaluated by the well-known Black and Scholes formula as follows.
The parameter estimation is performed as in (2) by only considering a single at-the-money option.

Stochastic Volatility
An alternative to the geometric Brownian motion is the stochastic volatility model of [35]. This specification allows for more flexibility and captures more accurately observed market dynamics. As in the previous model, the change in asset prices follows a Geometric Brownian Motion as in (17), whereas the volatility parameter σ follows an Ornstein-Uhlenbeck process.
where µ denotes the drift term, σ the volatility of stock price, κ the mean reversion speed, θ the long-run variance, t v the current variance and ξ the volatility of volatility parameter.
In contrast to the estimation of the risk-neutral density, the risk aversion parameters are not obtained directly from option prices as in (2). The Heston model is estimated by minimising the squared error between market implied volatility and model implied volatility 2 2 The required calculations for this method are computationally extensive and the weighting approach by [36] helps to reduce the amount of computations.

Expected Risk-Aversion
In line with [11] [13] and [17], the risk-aversion coefficient can be determined by the relation of the utility function to the risk-neutral and objective density function. The absolute risk-aversion, RA, can be then defined as: Using this result, the relative risk-aversion, γ , can be retrieved by: Since it is assumed that the transformation from risk-neutral to subjective dynamics only impacts the mean return, the risk-aversion estimate is the same at every point x.

Simulations
This simulation exercise aims to unveil the influence of different characteristics of the underlying asset price process on risk aversion. As a starting point, we use the conventional Geometric Brownian Motion to model asset price dynamics and we then extend this to accommodate stochastic volatility, as discussed in Section 3. The simulations provide useful insights about how and to what extent risk-aversion estimates are influenced by different factors within the considered models. The main point of interest is the influence of stochastic volatility on the risk-aversion estimates. Therefore, our simulations below include six different scenarios and examine the impact on risk aversion, as depicted in Panels A-F in the Table 1, Table 2, and Figures 1-3. Simulation A evaluates the change in the mean return, while Simulation B evaluates the impact of volatility, under a Geometric Brownian Motion. The remainder panels focus on the presence of stochastic volatility. The volatility process is modelled by an Ornstein-Uhlenbeck process. Simulation C investigates on the effect of the mean reversion κ of volatility. Simulation D evaluates the impact of fat-tailed distributions via the volatility of volatility ξ . Lastly, Simulation E and F evaluate the impact of skewness in conjunction with a fat-tailed distribution. The latter might be the empirically most relevant case which can provide guidance for practical purposes.
Each simulation scheme derives risk aversion over an estimation horizon of 120 H = monthly periods. As described in Section 3, the risk-neutral and subjective parameters are kept constant over the entire estimation horizon to ensure only the current market dynamics are considered. The simulated option prices are obtained one month prior to expiry. One month consists of 21 trading days, which results in 252 trading days per year. Each day is discretised by ten sub-steps to assure a sufficient precision of the simulated asset prices. The entire simulation is repeated 5000 N = times.

Geometric Brownian Motion
Within the Geometric Brownian Motion, both the mean, µ , and the variance, respectively. Figures 1-3 A displays the case of an stepwise increase of µ from −0.5% to 4%. It is apparent that risk-aversion increases when the expected return rises since the compensation for risk becomes higher. Consequently, the mean of the subjective density shifts to reflect the premium in the underlying asset price process.  However, this relation is inverse, meaning when variance is high (low) the variation of risk aversion is low (high). This feature can be confirmed by (7), which clearly illustrates that the risk-premium is linked to the product of

Stochastic Volatility
Extending the Geometric Brownian motion to accommodate stochastic volatility allows to further analyse the impact of volatility on the underlying asset. We focus on the effect of the following parameters: the speed of mean reversion to the long run volatility, κ , the volatility of volatility, ξ , and the correlation between the two Wiener processes of the stock price and volatility, ρ .  In contrast, the parameter ξ has an opposite effect on risk aversion, as shown in Figure 1-3 D. Increasing ξ results in more variation of risk-aversion. This can be attributed to the fact that ξ models the tail thickness of the asset price distribution. The simulation unveils that a higher frequency of extreme outcomes outweighs the impact of higher volatility, resulting in higher variation of risk-aversion.
Simulations E and F elaborate on the remaining parameter ρ , which captures asymmetry in the asset price distribution. Both plots show that with an increase in ρ the asymmetry of the relative risk-aversion coefficient increases. This can be seen by the diverging mean and median parameters but also in the shrinking upper percentiles. The higher ξ is the greater the absolute difference between median and average is. The influence of ρ on the mean squared error is only minor as this is primarily driven by ξ .
In addition to Mean Square Error consideration, we also provide simulation results on the maximum likelihood estimates from the transformation from risk-neutral to subjective density and compare them to the risk-neutral likelihoods.
The results are tabulated in Table 2 and draw a more complete picture along with results on mean squared error. In simulation A and B values are fairly close to each other. However, this changes from simulation C and ongoing. It shows that the GB2 density clearly fits the data the best, enforcing that a flexible density type is preferable when the underlying asset price distribution moves further away from lognormal. Furthermore, when the subjective density is compared to the risk-neutral density the likelihood is in all cases higher suggesting an improved fit.
Collectively, simulation results indicate that the main influence on the certainty of risk-aversion may be attributed to the volatility and the volatility of volatility. The mean mainly shifts the risk-aversion estimate in one or the other direction. In addition, flexible class of densities enable to obtain more accurate estimates for risk-aversion when the asset price distribution is skewed and has excess kurtosis. The fact that more flexible density types are favourable could also be shown by the likelihood comparison of the simulation. Finally, variation of risk-aversion strongly depends on the volatility component. Beyond that it could be shown that asymmetry in the asset price distribution also skews risk-aversion estimates. All these results provide useful insights to practitioners and highlight the importance of accounting for volatility when evaluating risk aversion estimates.

Data and Estimation
This study considers European call and put options, which can only be exercised at the maturity date. We obtain option prices from Optionmetrics for the period  Table 3).
The risk-neutral density parameters are estimated minimising (2) at each observation date. Multiple start values were supplied into the minimisation problem to ensure the minimum is robust. Table 4 provides the summary statistics of the option pricing errors. The reported values are in a similar range as, for example, [24]. There are remarkable differences between the more flexible densities and the lognormal approach. The pricing errors can attributed to the greater flexibility of the MLN and GB2, which better replicate the non-normal distribution contained in option prices.

Results
The empirical results, applying the simulation method to observable S&P 500 index options, are presented in Figure 4 and Figure 5. The graphs display the mean and median simulated risk-aversion estimate enclosed by confidence intervals. Furthermore, the estimation using a sequence of risk-neutral densities and asset price outcomes as in [1] is displayed. Their procedure is denoted as "historical approach''. For the case of the Geometric Brownian Motion, the results of the different density types are similar with no major discrepancies. Overall, the simulation method suggests a symmetric distribution of risk aversion since the mean and median are close to each other and the percentiles are evenly distributed. Furthermore, during times of high volatility, as for example at the subprime crisis, the variation in risk aversion lowers. These observations confirm the results from the simulation.
More complex are the findings assuming stochastic volatility for the price process. The greater flexibility of the underlying process allows for an asymmetric shape. While in the previous case the mean and median for all density types are close together, they diverge for the lognormal case. The median tends to be in excess towards the median, imposing a negative skew. This observation is a result of the negative skew induced by the stochastic volatility process. The lognormal distribution is not able to accommodate the more complex shape of the underlying probability distribution. Therefore, the estimation of the risk-aversion parameter compensates for this inflexibility, resulting in a skewed distribution of risk aversion. This observation highlights the importance of flexible density types when trying to estimate risk aversion. This finding lines up with the characteristics of the simulation, that more flexible density types yield a lower mean square error. Furthermore, the confidence  intervals have a stronger variation compared to the Geometric Brownian Motion. This is not a surprise due to the greater complexity of stochastic volatility. However, remarkable are the differences in the lower confidence intervals. For the MLN and GB2, they swing far less into negative values of risk aversion, which further highlights the importance of flexible density types.
Overall, it is remarkable how uncertain the estimates are even within basic assumptions as displayed in Figure 4 and Figure 5. Generally, the simulationbased methods tend to exceed the risk-aversion estimates of the historical approach. Nevertheless, the confidence intervals mostly encompass the historical estimates, particularly when assuming stochastic volatility.
To gain further insight, we use the obtained mean risk-aversion estimate and transform the risk-neutral density into its subjective counterpart. Afterwards, we evaluate their predictability based on the log-likelihood of the final asset price outcomes. This results in 122 evaluation periods since the historical method retrieves the estimates using a 10 year rolling window. The log-likelihood is reported in Table 5 in excess to the risk-neutral LN. The reported values can clearly show that the more flexible density types are preferable towards the LN.
Furthermore, the GB2 is clearly in excess towards the MLN. The additional parameter γ is further tested using the log-likelihood ratio test. All methods are significant at the 20% level. The reason for the low significance level could be rooted in the short sample size. [27] encounter similar difficulties and argue the insignificance is a type II error. These authors base this claim on the findings of [37] and others since the risk-premium is small relative to volatility and, therefore, difficult to estimate.
Finally, we test two competing forecasts for equal expected log-likelihoods following [33]. The results are reported in Table 6. Comparing the different density types with each other clearly shows that the GB2 is the preferable method. The test rejects the null hypothesis at the 1% level in all cases. Within each density type findings are less clear. Only for the historical risk aversion, the subjective density is significantly different from the risk-neutral density. In none Figure 5. Risk-aversion simulation results and historical estimates. The simulation assumes the asset price follows stochastic volatility. The dark blue line denotes the average estimate of risk-aversion whereas the green dashed-dotted line denotes the median. The inner condence interval (dark blue) de ned by the 25 th and 75 th percentile. The outer con dence interval (light blue) is de ned by the 5 th and 95 th percentile. The red line de nes the risk-aversion obtained by the classic method.  of the cases, the proposed simulation method is significantly different to the risk-neutral density. A possible explanation for this observation could be that the volatility dynamics are assumed to be the same in the risk-neutral and subjective model dynamics. It is well known that option prices contain a volatility premium and the volatility in the subjective process tends to be lower. Therefore, when assuming the risk-neutral and subjective volatility dynamics to be the same, we implicitly assume similar expected log-likelihoods (Table 6).

Conclusions
This paper provides a novel methodology for uncovering the presence of risk aversion in asset prices. Risk aversion is interpreted as the coefficient of relative risk aversion characterizing a power utility function. In contrast to most of the literature, we extract risk aversion from the risk-neutral distribution underlying the risk-neutral valuation of option prices. More specifically, we exploit a parametric expression relating the risk-neutral distribution function with an objective distribution describing the true probability distribution of asset prices and a power utility function describing individuals' risk attitude. This expression provides a closed-form solution that links risk-neutral expectations with a subjective valuation that incorporates the preferences of a representative agent.
This approach allows us to estimate consistently the coefficient of relative risk aversion as the sample size increases.
We have also proposed a simulation scheme to introduce uncertainty into the pricing models. By doing so, we have constructed confidence intervals that Journal of Mathematical Finance reflect the uncertainty around the parameter estimates of relative risk aversion.
We have also accommodated the presence of dynamics in risk aversion by proposing a parametric model for describing the dynamics of asset prices.
Finally, as an additional contribution, we have explored different types of density functions for describing the probability law of asset prices. By doing so, we have explored the predictive ability of the model and assessed, out of sample, the reliability of the different estimates of the relative risk-aversion coefficient for prediction purposes.
Our results, obtained from an extensive simulation exercise, suggest that there is more variation in the risk-aversion coefficient than implied by [1]. Not surprisingly, the volatility component has a major impact on the coefficient of relative risk aversion. This is due to the direct link between volatility, risk aversion and the risk premium. A way to make the risk-aversion estimate less varying and more stable across time is to assume more flexible density types.
Especially, when stochastic volatility is present more flexible density types obtain more stable values. The choice of a parametric model for describing asset prices also has a major effect on the estimates of risk aversion. Thus, we note that imposing a lognormal density on asset prices produces greater variation in the coefficient of risk aversion compared to the mixture-lognormal and generalised beta distribution of the second kind. Lastly, we also perform a forecasting exercise to assess the relative predictive ability of each estimation method out of sample. For all density types, there is no statistically significant difference between the risk-aversion coefficient based on the Monte-Carlo simulation and the classic approach developed by [1]. In contrast, comparing both approaches against the risk-neutral density reveals statistically significant differences.
Within this paper we centred our attention on the classic case of a power utility function. Therefore, it can be extended in different directions to obtain further insights into the degree of risk aversion in financial markets and its variation over time. First, an interesting avenue for further research is to evaluate the influence of the volatility premium on risk-aversion estimates.
Second, exploring alternative utility functions that provide a different characterization of individuals' risk aversion can be a fruitful strategy to uncover different patterns of risk aversion under different assumptions on individuals' risk profile. Finally, our simulation exercise clearly shows that individuals' risk aversion changes over time. It might be of much interest to applied economists and financial practitioners to explore suitable time series models for forecasting asset prices constructed from reliable estimates of the risk-aversion coefficient.

Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.