Deriving the Gibbons-Ross-Shanken CAPM Test from a Market Model
Summary
The document sets up a multivariate excess-return market model for several assets, with normally distributed, independent errors and a shared market excess return. It writes the likelihood and differentiates with respect to the intercept vector, yielding an intercept estimate based on sample mean asset returns, estimated market exposure, and the sample mean market return. The question concerns why the estimated intercept variance includes an adjustment involving the market mean and variance, and how that variance leads to a Wald test of the joint hypothesis that all asset alphas are zero.
The supplied response does not complete the derivation. It points to the general Wald statistic for linear restrictions in a regression and suggests using the regressor cross-product matrix, but explicitly leaves unresolved how its inverse produces the expression in the question. Thus the document offers a model setup and a useful connection to regression-based Wald testing, rather than a full derivation or empirical result. Its stated normality and temporal independence assumptions delimit the setup, and readers would need another source to verify the missing algebra and finite-sample test details.
Key ideas
- The setup models multiple asset excess returns with intercepts, market exposures, and jointly distributed errors.
- Maximum likelihood differentiation gives an intercept estimate formed from sample means and estimated market exposure.
- The question asks how uncertainty in the intercept estimate produces the Gibbons-Ross-Shanken Wald statistic.
- The provided response relates the problem to the regression Wald test but does not finish the variance derivation.
Tags
Full text
# Gibbons, Ross, Shanken Test derivation by MLE
# Gibbons, Ross, Shanken Test derivation by MLE
I Am trying to derive the expression for the GRS test of the CAPM. I am following the book: The Econometrics Of Financial Markets by Campbell, Lo, McKinley (1997).
Define $Z_t$ as an $N×1$ vector of excess returns for N assets. We assume that the excess returns can be described by the following excess-return market model:
$$Z_t = \alpha + \beta Z_{mt} + \epsilon_t$$ We assume that excess returns are jointly normal, with: $$E[\epsilon_t]=0 $$ N×1 vector $$E[\epsilon_t \epsilon_t']=\Sigma$$
Accordingly, because excess returns are normally distributed conditionally on the excess return of the market and assuming they are temporally IID, given T observations, we get the following log-likelihood function:
$$L(\alpha,\beta,\Sigma)=-NTlog(2\pi)-T/2log(det(\Sigma))-1/2 \sum_{t=1}^{T} (Z_t-\alpha-\beta Z_{mt})'\Sigma^{-1}(Z_t-\alpha-\beta Z_{mt})$$
The partial first derivative w.r.t. alpha is: (1) $$\partial L/\partial \alpha=\Sigma^{-1}\sum_{t=1}^{T}(Z_t-\alpha-\beta Z_{mt}) $$
From which, by setting it equal to 0, we get the MLE of alpha:
$$\hat{\alpha}=\hat{\mu}-\hat{\beta}\hat{\mu_{m}}$$
Where $\hat{\mu}=1/T\sum_{t=1}^{T} Z_t$ and $\hat{\mu_m}=1/T\sum_{t=1}^{T} Z_{mt}$
The authors claim that the variance of the MLE estimator of alpha is $$Var[\hat{\alpha}]=1/T[1+\hat{\mu_m}^2/\hat{\sigma_m}^2]\Sigma$$ Where $\hat{\sigma_m}^2=1/T\sum_{t=1}^{T} (Z_{mt}-\hat{\mu_m})^2 $
So that the GRS test is simply the Wald statistics:
$$J= \hat{\alpha}'[var[\hat{\mu}]]^{-1}\hat{\alpha}=T[1+\hat{\mu_m}^2/\hat{\sigma_m}^2]^{-1}\hat{\alpha}'\Sigma^{-1}\hat{\alpha}$$
Of the null hypothesis that the alphas are jointly zero.
I know that the variance of the estimates can be derived using the inverse of the Fisher information matrix. However, if I compute the derivative of (1), namely the second derivative of the LogLik w.r.t. alpha, change sign and then take its expectation, I can not obtain the expression of the variance claimed by the authors. Can you help me with this last step , please?
## Answer by mark leeds (score 1)
https://quant.stackexchange.com/a/47146
Hi: This is an incomplete answer but I needed room. The Wald statistic for testing a linear constraint , $Rb = r$ is ,
$(Rb - r )^{\prime}[R(X^{\prime}X)^{-1} R^{\prime}]^{-1}(Rb - r)/s^2$
$X^{\prime}X$ can be obtained from P4 and, in your case, $R = 1$ and $b = \alpha$. But I still don't see how the expression for $(X^\prime X)^{-1}$ results in what you wrote. Hopefully someone else can help here because I don't see it. Note that $\Sigma$ is just a scale factor so don't worry about that.Shown in full with attribution under the source's licence. Licence: CC BY-SA 4.0 (Stack Exchange)
This summary was written by Stratmill's research agent from the original; it is not a copy of the source.