SharpeR: statistics concerning Sharpe ratio and Markowitz portfolio

Description Sharpe Ratio Optimal Sharpe Ratio Spanning and Hedging Legal Mumbo Jumbo Note Author(s) References

Description

Inference on Sharpe ratio and Markowitz portfolio.

Sharpe Ratio

Suppose xi are n independent draws of a normal random variable with mean mu and variance sigma^2. Let xbar be the sample mean, and s be the sample standard deviation (using Bessel's correction). Let c0 be the 'risk free' or 'disastrous rate' of return. Then

z = (xbar - c0)/s

is the (sample) Sharpe ratio.

The units of z are per root time. Typically the Sharpe ratio is annualized by multiplying by sqrt(d), where d is the number of observations per year (or whatever the target annualization epoch.) It is not common practice to include units when quoting Sharpe ratio, though doing so could avoid confusion.

The Sharpe ratio follows a rescaled non-central t distribution. That is, z/K follows a non-central t-distribution with m degrees of freedom and non-centrality parameter ζ / K, for some K, m and zeta.

We can generalize Sharpe's model to APT, wherein we write

x_i = alpha + sum_j beta_j F_j,i + epsilon_i,

where the F_{j,i} are observed 'factor returns', and the variance of the noise term is sigma^2. Via linear regression, one can compute estimates alpha, and sigma, and then let the 'Sharpe ratio' be

z = (alpha - c0)/sigma.

As above, this Sharpe ratio follows a rescaled t-distribution under normality, etc.

The parameters are encoded as follows:

Optimal Sharpe Ratio

Suppose xi are n independent draws of a q-variate normal random variable with mean mu and covariance matrix Sigma. Let xbar be the (vector) sample mean, and S be the sample covariance matrix (using Bessel's correction). Let

Z(w) = (w'xbar - c0)/sqrt(w'Sw)

be the (sample) Sharpe ratio of the portfolio w, subject to risk free rate c0.

Let w* be the solution to the portfolio optimization problem:

max {Z(w) | 0 < w'Sw <= R^2},

with maximum value z* = Z(w*). Then

w* = R S^-1 xbar / sqrt(xbar' S^-1 xbar)

and

z* = sqrt(xbar' S^-1 xbar) - c0/R

The variable z* follows an Optimal Sharpe ratio distribution. For convenience, we may assume that the sample statistic has been annualized in the same manner as the Sharpe ratio, that is by multiplying by d, the number of observations per epoch.

The Optimal Sharpe Ratio distribution is parametrized by the number of assets, q, the number of independent observations, n, the noncentrality parameter,

zeta* = sqrt(mu' Sigma^-1 mu),

the 'drag' term, c0/R, and the annualization factor, d. The drag term makes this a location family of distributions, and by default we assume it is zero.

The parameters are encoded as follows:

Spanning and Hedging

As above, let

Z(w) = (w'xbar - c0)/sqrt(w'Sw)

be the (sample) Sharpe ratio of the portfolio w, subject to risk free rate c0.

Let G be a g x q matrix of 'hedge constraints'. Let w* be the solution to the portfolio optimization problem:

max {Z(w) | 0 < w'Sw <= R^2, G S w = 0},

with maximum value z* = Z(w*). Then z*^2 can be expressed as the difference of two squared optimal Sharpe ratio random variables. A monotonic transform takes this difference to the LRT statistic for portfolio spanning, first described by Rao, and refined by Giri.

Legal Mumbo Jumbo

SharpeR is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.

Note

The following are still in the works:

  1. Corrections for standard error based on skew, kurtosis and autocorrelation.

  2. Tests on Sharpe under positivity constraint. (c.f. Silvapulle)

  3. Portfolio spanning tests.

  4. Tests on portfolio weights.

This package is maintained as a hobby.

Author(s)

Steven E. Pav shabbychef@gmail.com

References

Sharpe, William F. "Mutual fund performance." Journal of business (1966): 119-138. https://ideas.repec.org/a/ucp/jnlbus/v39y1965p119.html

Johnson, N. L., and Welch, B. L. "Applications of the non-central t-distribution." Biometrika 31, no. 3-4 (1940): 362-389. doi: 10.1093/biomet/31.3-4.362

Lo, Andrew W. "The statistics of Sharpe ratios." Financial Analysts Journal 58, no. 4 (2002): 36-52. https://www.ssrn.com/paper=377260

Opdyke, J. D. "Comparing Sharpe Ratios: So Where are the p-values?" Journal of Asset Management 8, no. 5 (2006): 308-336. https://www.ssrn.com/paper=886728

Ledoit, O., and Wolf, M. "Robust performance hypothesis testing with the Sharpe ratio." Journal of Empirical Finance 15, no. 5 (2008): 850-859. doi: 10.1016/j.jempfin.2008.03.002

Giri, N. "On the likelihood ratio test of a normal multivariate testing problem." Annals of Mathematical Statistics 35, no. 1 (1964): 181-189. doi: 10.1214/aoms/1177703740

Rao, C. R. "Advanced Statistical Methods in Biometric Research." Wiley (1952).

Rao, C. R. "On Some Problems Arising out of Discrimination with Multiple Characters." Sankhya, 9, no. 4 (1949): 343-366. https://www.jstor.org/stable/25047988

Kan, Raymond and Smith, Daniel R. "The Distribution of the Sample Minimum-Variance Frontier." Journal of Management Science 54, no. 7 (2008): 1364–1380. doi: 10.1287/mnsc.1070.0852

Kan, Raymond and Zhou, GuoFu. "Tests of Mean-Variance Spanning." Annals of Economics and Finance 13, no. 1 (2012) https://econpapers.repec.org/article/cufjournl/y_3a2012_3av_3a13_3ai_3a1_3akanzhou.htm

Britten-Jones, Mark. "The Sampling Error in Estimates of Mean-Variance Efficient Portfolio Weights." The Journal of Finance 54, no. 2 (1999): 655–671. https://www.jstor.org/stable/2697722

Silvapulle, Mervyn. J. "A Hotelling's T2-type statistic for testing against one-sided hypotheses." Journal of Multivariate Analysis 55, no. 2 (1995): 312–319. doi: 10.1006/jmva.1995.1081

Bodnar, Taras and Okhrin, Yarema. "On the Product of Inverse Wishart and Normal Distributions with Applications to Discriminant Analysis and Portfolio Theory." Scandinavian Journal of Statistics 38, no. 2 (2011): 311–331. doi: 10.1111/j.1467-9469.2011.00729.x


shabbychef/SharpeR documentation built on Aug. 21, 2021, 8:50 a.m.