Inference for the Sharpe Ratio using a Likelihood-Based Approach
|
|
- Mildred Scott
- 6 years ago
- Views:
Transcription
1 Inference for the Sharpe Ratio using a Likelihood-Based Approach Ying Liu Marie Rekkas Augustine Wong Abstract The Sharpe ratio is the prominent risk-adjusted performance measure used by practitioners. Statistical testing of this ratio using its asymptotic distribution has lagged behind its use. In this paper, highly accurate likelihood analysis is applied for inference on the Sharpe ratio. Both the one and two-sample problems are considered. The methodology has O(n 3/ ) distributional accuracy and can be implemented using any parametric return distribution structure. Simulations are provided to demonstrate the method s superior accuracy over existing methods used for testing in the literature. Keywords: Canonical Parameter; Lagrange Multiplier; Likelihood Analysis; p-value; Sharpe Ratio; Tilted Log-Likelihood Department of Economics, York University, 4700 Keele Street, Toronto, Ontario M3J 1P3 Corresponding Author. Department of Economics, Simon Fraser University, 8888 University Drive, Burnaby, British Columbia V5A 1S6, mrekkas@sfu.ca, phone: (778) , fax: (778) Department of Mathematics and Statistics, York University, 4700 Keele Street, Toronto, Ontario M3J 1P3. It is noted that part of this paper originates from a chapter of Liu s doctoral dissertation. 1
2 1 Introduction The measurement of fund performance is an integral part of investment analysis. Investments are often ranked and evaluated on the basis of their risk-adjusted returns. Several risk-adjusted performance measures are available to money managers of which the Sharpe ratio is the most popular. Introduced by William Sharpe in 1966 (Sharpe (1966)), this ratio provides a measure of a fund s excess returns relative to its volatility. Expressed in its usual form, the Sharpe ratio for an asset with an expected return given by µ and standard deviation given by σ is given by SR = µ R f σ where R f is the risk-free rate of return. From this expression, it is clear to see how this ratio provides a measure of a fund s excess return per unit of risk. The Sharpe ratio has been extensively studied in the literature., The main criticism leveled against this measure concerns its reliance on only the first two moments of the returns distribution. If investment returns are normally distributed then the Sharpe ratio can be justified. On the other hand, if returns are asymmetric then it can be argued that the measure may not accurately describe the fund s performance as moments reflecting skewness and kurtosis are not captured by the ratio. To address this issue, several measures exist in the literature which integrate higher moments into the performance measure. The Omega measure is one such measure that uses all the available information in the returns distribution. Keating and Shadwick (00) provide an introduction to this measure. While various methods are available, they are also more complex and often very difficult to implement in practice. To gauge the trade-off between the attractiveness of such measures and their cost, Eling and Schuhmacher (007) compared the Sharpe ratio with 1 other approaches to performance measurement. Eling and Schuhmacher (007) focussed on the returns of,763 hedge funds. Hedge funds are known to have return distributions which differ significantly from the normal distribution and, as such, provide a rich and relevant environment for such a comparison. Their study concluded that the Sharpe ratio produced rankings that were largely identical to those obtained from the 1 other performance measures. In other words, the choice of performance measure did not matter much in terms of ranking funds. From a practical perspective, this imparts credibility to the use of the Sharpe ratio for ranking funds with highly non-normal returns. Our focus in this paper is on the statistical properties of the Sharpe ratio. While this ratio may
3 be the most widely known and used measure of risk-adjusted performance for an investment fund, fund managers rarely (if ever) indicate any measure of associated statistical significance of their rankings. This point was made in Opdyke (007). We use a likelihood-based statistical method to obtain highly accurate p-values and confidence intervals for evaluating the significance of the Sharpe ratio. While the methodology we use is applicable under any parametric distributional assumption, for expositional clarity we demonstrate the use of the method under the assumption of normal returns. We investigate inference for both the one sample problem as well as the two-sample problem. We compare our method to currently existing methods in the literature. Simulations indicate the superior performance of the method compared to existing methods in the literature. The paper is organized as follows. Section provides an overview of the asymptotic distribution of the Sharpe ratio under various assumptions. Section 3 lays out the framework for the likelihoodbased approach we will use for testing. Section 4 applies this methodology for inference on the Sharpe ratio. Examples and simulations for this one sample case are then presented in Section 5. The two-sample problem is given in Section 6, with example and simulations presented in Section 7. Section 8 concludes. Distribution of the Sharpe Ratio Jobson and Korkie (1981) and Lo (00) derived the asymptotic distribution of the Sharpe ratio assuming identically and independently normally distributed returns. Consider a fund with a return at time t given by R t, t = 1,, T. Further assume that the returns are identically and independently distributed (IID) as N(µ, σ ). If we let the vector θ = (µ, σ ), then as given in Lo (00) we have the following asymptotically normal distribution for ˆθ: T (ˆθ θ) N(0, V ), (1) where V σ 0 0 σ 4 and ˆθ = (ˆµ, ˆσ ) is the maximum likelihood estimator (MLE) with ˆµ and ˆσ given by the formulas ˆµ = 1 T ΣT t=1r t ˆσ = 1 T ΣT t=1(r t ˆµ). 3
4 Lo s result is straightforward to prove. For large samples, by the Central Limit Theorem, we have (ˆθ θ) [var(ˆθ)] 1 (ˆθ θ) χ p, where k = dim(θ) and var(ˆθ) = var(ˆµ) cov(ˆµ, ˆσ ) cov(ˆµ, ˆσ ) var(ˆσ ). See for instance, Cox and Hinkley (1974) for this result. Given µ = (1, 0)θ and σ = (0, 1)θ, we have (ˆµ µ) N(0, 1) and (ˆσ σ ) N(0, 1). var(ˆµ) var(ˆσ ) Given the Sharpe ratio is a function of θ, the Delta method can be used to derive its asymptotic distribution. Let SR = g(θ) be a continuous function with Jacobian denoted by g(θ)/ θ and let the Sharpe ratio estimate be given by around the true parameter value θ, we have ˆ SR = g(ˆθ). Applying a first-order Taylor expansion on g(ˆθ) g(ˆθ) =. [ ] g(θ) g(θ) + (ˆθ θ). θ Rearranging terms and multiplying by T yields [ ]. g(θ) T (g(ˆθ) g(θ)) = T (ˆθ θ). θ Using the above result and applying the Delta method, we have the asymptotic distribution of g(ˆθ): [ ]. g(θ) ( [ g(θ) T (g(ˆθ) g(θ)) = T (ˆθ θ) N 0, θ θ ] V g(θ) θ ). Note that the Jacobian term is readily computed as g(θ) = g(θ)/ µ = θ g(θ)/ σ 1/σ (µ R f )/σ 3. Combining terms yields the variance of the asymptotic distribution: [ ] g(θ) V g(θ) θ θ = 1 + (µ R f ) σ = SR. 4
5 By the Central Limit Theorem, the distribution of the Sharpe ratio assuming identically and independently distributed normal returns is therefore T ( SR ˆ SR) N(0, 1) () ˆ SR From this distribution, 100(1 α)% confidence intervals for the Sharpe ratio can be constructed in the usual fashion: ˆ SR ± z α/ ( ˆ SR ) /T, (3) where z α/ is the 100 (1 α/)-th percentile of the standard normal distribution. The above derivation only holds under the IID normal assumption. Under the IID assumption but without assuming normally distributed returns, Mertens (00) derived the following asymptotic distribution of the Sharpe ratio: T ( SR ˆ SR) ˆ SR SR ˆ ˆγ 3 + SR ˆ ˆγ where γ 3 and γ 4 are defined as N(0, 1) (4) γ 3 = E(R t µ) 3 σ 3, γ 4 = E(R t µ) 4 σ 4 with σ = SR. (5) Christie (005) derived the asymptotic distribution of the Sharpe ratio under the more relaxed assumption of stationarity and ergodicity. Opdyke (007) interestingly showed that the derivation provided by Christie (005) under the non-iid returns condition was in fact identical to the one provided by Mertens (00). For a complete discussion and proof see Opdyke (007). For our purposes, we will use asymptotic distributions () and (4) to compare to our proposed likelihoodbased approach. In the next section below we review the likelihood-based approach to inference and in the following section we apply the methodology to the Sharpe ratio. 3 Likelihood-Based Approach to Inference In this section we review standard first-order likelihood methods and then turn our attention to higher-order methods. For this section, the underlying model under consideration is as follows: Y = (Y 1,..., Y n ) is a vector of IID random variables with density given by f( ; θ) where θ is a k-dimensional vector of parameters. Statistical inference concerns drawing conclusions about θ or 5
6 a function of θ, namely ψ = ψ(θ), given an observed sample y = (y 1,..., y n ) and in this paper we assume that the interest parameter is ψ(θ). We further assume that this interest parameter is scalar. Let λ = λ(θ) be a (k 1)-dimensional nuisance parameter vector. The log-likelihood function of θ is ( n ) l(θ) = l(θ; y) = log f(y i ; θ). i=1 Likelihood analysis typically involves the maximum likelihood estimator (MLE) ˆθ = ( ˆψ, ˆλ) = arg max θ l(θ; y) and the constrained maximum likelihood estimator θ = (ψ, λ) = arg max λ l(θ; y) where the maximum is taken for fixed values of ψ. The constrained MLE can be solved by maximizing l(θ) subject to ψ(θ) = ψ. This is often done using the Lagrange multiplier method where the Lagrangean is given by H(θ, α) = l(θ; y) + α[ψ(θ) ψ]. (6) The Lagrange multiplier is given by ˆα. The tilted log-likelihood is defined as l(θ) = l(θ) + ˆα[ψ(θ) ψ]. (7) Using the above information, standard first-order departure measures can be obtained. These measures include the Wald departure (the standardized maximum likelihood estimator) and the signed square root log-likelihood ratio statistic: q q(ψ) = ( ˆψ ψ)[var( ˆψ)] 1/ (8) R R(ψ) = sgn( ˆψ ψ){[l(ˆθ; y) l( θ; y)]} 1/. (9) Note that var( ˆψ) can be estimated by using the Delta method var( ˆ ˆψ) ψ θ (ˆθ)[var(ˆθ)]ψ θ (ˆθ), where the variance of ˆθ, var(ˆθ), can be estimated by using either the Fisher expected information matrix [ ] E l(θ) θ θ 6
7 or the observed information matrix evaluated at the MLE j θθ (ˆθ) = l(θ) θ θ. θ=ˆθ The latter is preferred however because of the simplicity in calculation. The p-value function for these methods are p(ψ) = Φ(q) and p(ψ) = Φ(R), where Φ represents the standard normal cumulative distribution function. These p-value functions have order of convergence O(n 1/ ) and are correspondingly referred to as first-order methods. 1 A 100(1 α)% confidence interval is then given by: { [ min p 1 ( α ) (, p 1 1 α )] [, max p 1 ( α ) (, p 1 1 α )]}. (10) Approximations aimed at improving the accuracy of the first-order methods have been worked on during the past three decades. Higher-order asymptotic approximations for the distribution of an estimate can be constructed using expansion methods of which Edgeworth expansions and saddlepoint approximations are most common. Edgeworth approximations tend to produce good approximations near the center of a distribution, with relatively poorer tail approximations, while saddlepoint approximations produce remarkably accurate tail approximations. Given we are interested in tail probabilities, our focus will thus be on approximation methods that use the saddlepoint method. We will not review the saddlepoint approximation literature here but the interested reader is directed to Daniels (1954, 1987) and Barndorff-Nielsen and Cox (1994). For a canonical exponential family model with canonical parameter vector θ = (ψ, λ ) and density f(y; θ) = exp{ψu(y) + λ v(y) c(θ) + d(y)}, where (U(Y ), V (Y )) is a minimal sufficient statistic for θ, the saddlepoint-based approximation for the density of ˆθ is given as f(ˆθ; θ). = c j θθ (ˆθ) 1/ exp[l(θ; y) l(ˆθ; y)], (11) where c is a normalizing constant. This extremely important approximation to the density of the maximum likelihood estimator is referred to as Barndorff-Nielsen s p formula and can be found in 1 Note that (9) is invariant to reparameterization whereas (8) is not. However, in practice, (8) is preferred because of its simplicity of use. Doganaksoy and Schmee (1993) illustrate that (9) has better coverage than (8) in the cases they examined. 7
8 Barndorff-Nielsen (1980). Thus the marginal density for ˆψ can be obtained. In general however, the closed form of the marginal density is not available. Fraser, Reid and Wong (1991) showed that the result by Lugannani and Rice (1980) can be applied to approximate the p-value function and it takes the form: The statistic Q is given by p(ψ) =. ( 1 Φ(R) + φ(r) R 1 ). (1) Q Q Q(ψ) = ( ˆψ ψ) j θθ (ˆθ) 1/ j λλ ( θ) 1/, (13) where j θθ (ˆθ) is the observed information matrix and j λλ ( θ) is the observed nuisance information matrix defined as follows: j θθ (ˆθ) = l(θ; y) θ θ j λλ ( θ) = l(θ; y) λ λ θ=ˆθ θ= θ (14). (15) The term j θθ (ˆθ) j λλ ( θ) 1 can be viewed as an estimate of the variance of ˆψ where the elimination of the nuisance parameter λ has been taken into consideration. Note that Q is a standardized maximum likelihood departure in the canonical parameter scale. The statistic R is the log-likelihood ratio given by (9). An alternate approximation to (1) is given by Barndorff-Nielsen (1986): p(ψ) =. ( Φ R 1 ( R.= R Q)) log Φ(R ). (16) These approximations are asymptotically equivalent (see Jensen (199)) and have order of convergence O(n 3/ ). They are correspondingly referred to as third-order methods. Hence 100(1 α)% confidence interval for ψ can be obtained from (10). Consider now a general full rank exponential family model with parameter vector θ and density given by f(y; θ) = exp{ϕ (θ)t(y) c(θ) + h(t(y))}, where ϕ(θ) is the canonical parameter and t(y) is the canonical variable. In order to use either the Lugannani and Rice approximation given in (1) or the Barndorff-Nielsen approximation given in (16) for inference on our interest parameter ψ(θ), we need to calculate both R and Q in ϕ(θ) This approximation is valid for more general models but requires the existence of an ancillary statistic. 8
9 scale. Note that R calculated in the original θ scale is equivalent to R calculated in ϕ(θ) scale. Moreover, Fraser and Reid (1995) derive a Q in the ϕ(θ) scale. Their methodology involves replacing the parameter of interest ψ by a linear function of the ϕ(θ) coordinates. 3 parameter is given by This newly calibrated χ(θ) = ψ θ ( θ)ϕ 1 θ ( θ)ϕ(θ). (17) Then χ(ˆθ) χ( θ) measures the departure of ˆψ ψ in the ϕ(θ) scale. The function Q is then constructed as [ ( 1/ Q = sgn(ψ(ˆθ) ψ) χ(ˆθ) χ( θ) var ˆ χ(ˆθ) χ( θ))]. (18) Fraser and Reid (1995) show that ( ) var ˆ χ(ˆθ) χ( θ) ψ θ ( θ) j 1 θθ ( θ)ψ θ ( θ) j θθ ( θ) ϕ θ ( θ) j θθ (ˆθ) ϕ θ (ˆθ) where j θθ (ˆθ) is defined in (14) and j θθ ( θ) is defined as j θθ ( θ) = l(θ) θ θ where l(θ) is the tilted log-likelihood function defined in (7). The function Q in (18) along with R is used in the Lugannani and Rice approximation in (1) or the Barndorff-Nielsen expression in (16) which implicitly yields a new R. In the next section we apply the above methodology and use Barndorff-Nielsen s approximation to obtain highly accurate p-values to test particular hypothesized values of the Sharpe ratio. We further note that for all our examples and simulations in this paper, the Lugannani and Rice expression given in (1) was also used but we do not report results from this method as they are θ= θ virtually identical to those produced by the Barndorff-Nielsen expression. 4 Methodology for the Sharpe Ratio, We emphasize that while the third-order methodology discussed in the previous section is applicable under any parametric distributional assumptions, for expositional clarity we demonstrate the use of the method under the assumption of IID normal returns. 4 Consider a fund with return at time 3 For a detailed discussion and derivation see Fraser and Reid (1995). 4 See Fraser, Reid and Wu (1999) for the general model set-up. 9
10 t given by R t, (t = 1,, T ). And let R f be the mean return for the risk-free asset. We assume that R 1, R T IID N(µ, σ ). The interest parameter is the Sharpe Ratio: SR = µ R f σ = ψ (19) and θ = (µ, σ ). The log-likelihood and related first and second derivatives, denoted by l (θ) and l (θ) respectively, are given as: l(θ) = l(µ, σ ) = T log σ 1 σ Σ(R t µ) (0) l µ (θ) = 1 σ Σ(R t µ) l σ (θ) = T σ + 1 σ 4 Σ(R t µ) l µµ (θ) = T σ l µσ (θ) = 1 σ 4 Σ(R t µ) l σ σ (θ) = T σ 4 1 σ 6 Σ(R t µ). To obtain the maximum likelihood estimators, the first order conditions l µ (θ) = 0 l σ (θ) = 0 are solved, and we have The estimated Sharpe ratio is then ˆµ = R (1) ˆσ = Σ(R t R). () T SR ˆ = ˆµ R f. ˆσ The observed information matrix is given as j θθ (θ) = l θθ (θ) which can be expressed as follows: j θθ (θ) = l µµ(θ) l σ µ(θ) l µσ (θ) l σ σ (θ) = T σ 1 σ 4 Σ(R t µ) 1 σ 4 Σ(R t µ) T σ σ 6 Σ(R t µ). (3) 10
11 The observed information matrix is then evaluated at the MLE: j θθ (ˆθ) = Ṱ 0 σ 0 T ˆσ (4) with determinant equal to j θθ (ˆθ) = ( ) ( ) T T ˆσ ˆσ 4 = T ˆσ 6. (5) To solve for the constrained MLE, the Lagrange multiplier method from (6) is used [ µ Rf H(θ, α) = l(θ) + α σ ] ψ (6) with first derivatives equal to H µ (θ, α) = l µ (θ) + α σ H σ (θ, α) = l σ (θ) α(µ R f ) σ 3 H α (θ, α) = µ R f σ Setting these first derivatives equal to zero, we have: ψ. ˆα = σl µ ( θ) (7) µ = R f + ψ σ (8) Working on (9) we have which further simplifies to l σ ( θ) + σl µ ( θ) µ R f σ 3 = 0. (9) T σ + 1 σ 4 Σ(R t µ) + σ 1 σ Σ(R t µ) µ R f σ 3 = 0 T σ ψσ(r t R f ) σ + Σ(R t R f ) = 0. 11
12 If we let r t = R t R f and r = Σr t /T, then we have with σ + ψ r σ Σr t T = 0 σ = ψ r + (ψ r) + 4(Σrt /T ), since σ > 0. The constrained MLE are therefore given by: σ = ψ r + (ψ r) + 4(Σr t /T ) (30) µ = R f + ψ σ (31) ˆα = ( T ψ r σ ) (3) so that θ = ( µ, σ). Using these constrained maximum likelihood estimators we are able to obtain the titled log-likelihood as defined in (7), and the corresponding first and second derivatives are: l(θ) = l(θ) + ˆα ( µ Rf σ lµ (θ) = l µ (θ) + ˆα σ lσ (θ) = l σ (θ) ˆα ( µ Rf σ 3 lµµ (θ) = l µµ (θ) ) ψ lµσ (θ) = l µσ (θ) ˆα σ 3 lσ σ (θ) = l σ σ (θ) + 3 ˆα(µ R f ) 4 σ 5. ) We can use these second derivatives to calculate j θθ ( θ) as j θθ ( θ) = l µµ ( θ) l µσ ( θ) l σ µ( θ) l σ σ ( θ) For the normal model, the canonical parameter is given by: ϕ(θ) = ( µ σ, 1 ) σ. (33) 1
13 Two related matrices derived from the canonical parameter are given below ϕ θ (θ) = 1 σ µ σ σ 4 and With our parameter of interest ϕ 1 θ (θ) = σ µσ 0 σ 4. (34) we have ψ(θ) = µ R f σ ψ θ (θ) = (ψ µ (θ), ψ σ (θ)) = Quantities (33), (34) and (35) can be used to calculate (17). The quantity ( 1 σ, µ R ) f σ 3. (35) ( ) var ˆ χ(ˆθ) χ( θ) can further be calculated. We now have all the ingredients to calculate R (ψ) given in (16). 5 Examples and Simulations In this section, we provide examples and simulations for inference on the Sharpe ratio. We compute confidence intervals and p-values using our proposed third-order method given in (16). We label this method Proposed in the tables below. We compare our results to those obtained by using the existing methods, in particular, using both the Jobson and Korkie (1981) asymptotic distribution in () and the Mertens (00) distribution given in (4). These methods are labeled as Jobson and Korkie and Mertens, respectively. Results from the signed square root log-likelihood ratio statistic in (9) are additionally provided and labeled Likelihood Ratio. 5.1 Examples The dataset for our examples consists of monthly return prices for three time series. The first series represents return prices for a large-cap mutual fund (Fund), the second for a market index (Market), and the third for 90-day Treasury bills (Cash). This data spans a period of one year. The data is listed below in Table 1 and originates from the Matlab Financial Toolbox TM User s Guide. 5 5 Data in the User s Guide is provided for a 5-year period. We focus on the most recent year of data. See Section 4: Investment Performance Metrics of the manual for further information. 13
14 Table 1: Monthly Return Prices Month Fund Market Cash Table below reports 95% confidence intervals for the Sharpe ratio separately for the large-cap mutual fund and the market index for the four methods discussed above. Monthly returns are calculated from the return prices listed in Table 1 and are used to calculate the Sharpe ratio. The mean of the 90-day Treasury returns is used as the risk-free rate. Table shows that the confidence intervals obtained from the four methods produce rather different results. Theoretically, the proposed method has third-order accuracy whereas the remaining three methods do not. Although the confidence intervals produced by the Jobson and Korkie (1981) approximation have the best concordance with those produced by the proposed method, the difference is still noticeable. This result will be borne out in the simulations as well. Table : 95% Confidence Intervals for the Sharpe Ratio Method 95% CI for SR (Fund) 95% CI for SR (Market) Jobson and Korkie (-0.38, ) ( , 1.099) Mertens (-0.838, 1.004) ( , 1.066) Likelihood Ratio (-0.087, 0.964) ( , 1.161) Proposed (-0.319, ) (-0.11, 1.094) The p-value functions calculated from the methods discussed in this paper for the Sharpe Ratio for the large-cap mutual fund and the Sharpe Ratio for the market index are plotted in Figures 1 and, respectively. These significance functions can be used to obtain p-values for specific hypothesized values of the Sharpe ratio. As we are typically interested in tail probabilities which tend to be small, it is important to estimate such probabilities with precision. A few values with 14
15 their corresponding p-values are provided in Tables 3 and 4 below for the market fund and market index, respectively. From these tables we can see that the p-values vary across the methods. If for instance, interest is on testing whether SR = 1 for either fund, then the corresponding p-values for such a test are given in the last column of Tables 3 and 4. Focussing on the market index and using a 5% level of significance, the Sharpe ratio may or may not be statistically significant depending on the method chosen for the hypothesis test. Figure 1: p-value Function for Fund Probability Jobson and Korkie Mertens Likelihood Ratio Proposed SR 15
16 Figure : p-value Function for Market Probability Jobson and Korkie Mertens Likelihood Ratio Proposed SR Table 3: Fund: p-values for SR SR Method Jobson and Korkie Mertens Likelihood Ratio Proposed
17 Table 4: Market: p-values for SR SR Method Jobson and Korkie Mertens Likelihood Ratio Proposed
18 5. Simulations Two simulation studies of size 10,000 were performed to compare the three existing methods to the proposed third-order method. The first simulation was constructed to mimic the mutual fund return data from the above example and the second to mimic the market index returns: X N( , ), Y N( , ) and R f = %, 95% and 99% confidence intervals for the Sharpe ratio were obtained for each sample and the lower probability error, the upper probability error, and the central coverage were recorded. Lower (upper) probability error refers to the proportion of times the true parameter value falls below (above) the confidence interval. Central coverage is defined as the proportion of times the true parameter values falls within the confidence intervals. Tables 5 and 6 report the results from these simulations for the fund and market returns, respectively. For reference, the nominal values and the corresponding standard errors have been included for each the lower probability error, the upper probability error, and the central coverage. Table 5: Results for Simulation Study for Fund for n=1 CI Method Lower Error Upper Error Central Coverage 90% Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error % Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error % Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error From these results tables it is clear that the third-order method outperforms the other three methods based on the criteria we examined. It must be noted however, that the Jobson and Korkie (1981) method provides surprisingly good results. We note that other values for the parameters of 18
19 Table 6: Results for Simulation Study for Market for n=1 CI Method Lower Error Upper Error Central Coverage 90% Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error % Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error % Jobson and Korkie Mertens Likelihood Ratio Proposed Nominal Standard Error the simulation were chosen (but not reported) and the results were consistent with the reported results. Standard errors for the simulations can easily be computed. This information has been included in the tables and labeled Standard Error. From these standard errors, it can be seen that the proposed method produces results that are uniformly within three standard deviations of the nominal value. The other three methods produce less satisfactory results. 6 The Two-Sample Case Suppose one is interested in testing hypotheses concerning the Sharpe ratios of two funds X and Y. For instance, one may be interested in testing the null hypothesis SR X SR Y or SR X = SR Y against the alternative hypothesis SR X < SR Y or SR X SR Y, respectively. In this section we apply the third-order methodology described in Section 3 to test the difference between the Sharpe ratios of two funds. Consider two funds with returns given by (x 1,..., x n ) and (y 1,..., y m ). Further assume for expositional clarity that these returns are identically and independently distributed as N(µ X, σx ) and N(µ Y, σy ), respectively. The mean return for the risk-free asset is given by R f. Our interest is in testing the difference in Sharpe ratios of the two funds. The interest parameter 19
20 ψ is then defined as ψ = µ X R f σ X µ Y R f σ Y. (36) The log-likelihood function for this problem is written as: l(θ) = l(µ X, σx, µ Y, σy ) = n log σ X 1 σx Σ(x i µ X ) m log σ Y 1 σy Σ(y i µ Y ). (37) The corresponding first derivatives are: l µx = 1 σx Σ(x i µ X ) l σ X = n σ X l µy = 1 σy Σ(y j µ Y ) l σ Y = m σ Y + 1 σx 4 Σ(x i µ X ) + 1 σy 4 Σ(y j µ Y ). These first derivatives are set equal to zero and simultaneously solved for the overall MLE: j θθ (θ) = ˆµ X = x ˆσ X = Σ(x i x) n ˆµ Y = ȳ ˆσ Y = Σ(y j ȳ). m The second derivatives are required to obtain the observed information matrix: n 1 Σ(x σx σx 4 i µ X ) Σ(x σx 4 i µ X ) n σ 4 X + 1 Σ(x σx 6 i µ X ) 0 0 m 0 0 σy Σ(y σy 4 j µ Y ) m σ 4 Y 1 Σ(y σy 4 j µ Y ) + 1 Σ(y σy 6 j µ Y ) Evaluating this observed information matrix at the MLE produces j θθ (ˆθ) with corresponding determinant denoted by j θθ (ˆθ). To solve for the constrained MLE, the Lagrange multiplier method is used H(θ) = l(θ) + α µ X R f σ X µ Y R f σy. ψ (38) 0
21 with the following first derivatives: H µx = 1 σx Σ(x i µ X ) + α σx H σ X = n σ X H µy = 1 σy Σ(y j µ Y ) H σ X = m σ Y H α = µ X R f σx + 1 σx 4 Σ(x i µ X ) α(µ X R f ) (σx )3/ α σ Y + 1 σy 4 Σ(y j µ Y ) + α(µ Y R f ) (σy )3/ µ Y R f σy Setting these derivatives equal to zero and solving the resulting system produces the constrained ψ. MLE. The constrained MLE can be obtained by solving the following iteratively: and Lagrangian σ X = Σ(x i µ X ) + ( µ X R f )Σ(x i µ X ) n σ Y = Σ(y j µ Y ) + ( µ Y R f )Σ(y j µ Y ) m µ X = m σ Y R f + n x σ Y + mȳ σ X m σ X R f + mψ σ X σ Y (n + m) σ X µ Y = n x σ Y n σ Y R f + n σ X R f nψ σ X σ Y + mȳ σ X (n + m) σ X ˆα = Σ(y j µ Y ) σ Y. So we have ˆθ ψ = ( µ X, σ X, µ Y, σ Y ). We are now able to define the tilted log-likelihood: l(θ) = l(θ) + ˆα [ µx R f σ X µ Y R f σ Y ] ψ 1
22 with first derivatives: lµx = l µx + ˆα σ X [ µx R f lσ = l X σ ˆα X σx 3 lµy = l µy ˆα σ Y lσ Y = l σ Y + ˆα [ µy R f σ 3 Y lµx µ X = l µx µ X lµx σ X = l µ X σ X ˆα σ 3 Y lµx µ Y = 0. ] ] The second derivatives are needed to calculate j θθ (θ): j θθ (θ) = j θθ (θ) + 0 ˆα σ 3 X ˆα σ 3 X 3ˆα(µ X R f ) 4σ 5 X ˆα 0 0 ˆα σ 3 Y σy 3 3ˆα(µ Y R f ) 4σY 5. The canonical parameter for this model is given by ϕ θ (θ) = ( µx σx, 1 σ X, µ Y σy, The matrices associated with this vector can now be calculated: 1 µ σ X 0 0 X σx σ ϕ θ (θ) = X µ σ Y Y σy σy 4 1 σ Y ). (39) and ϕ 1 θ (θ) = σx µ XσX σx σy µ Y σy. (40) σ 4 Y
23 Given our parameter of interest ψ(θ) = µ X R f σx µ Y R f σy we have ψ θ (θ) = 1, (µ X R f ) σx (σx, )3/ 1 σ Y, (µ Y R f ) (σ Y )3/. (41) Quantities ϕ θ (θ), ϕ 1 θ (θ) and ψ θ(θ) can be used to calculate the new parameter χ(θ) given in (17). The additional ingredients necessary for the construction of the standardized maximum likelihood departure Q given in (18) are readily available from the above information. 7 Example and Simulations In this section we provide an example and a simulation study for the two-sample Sharpe ratio case. We compare the third-order likelihood-based inference method to the classical methods used for testing, namely the maximum likelihood statistic and the likelihood ratio statistic. These latter two statistics are analogous to expressions (8) and (9). 7.1 Example For our example, we use the data presented in Table 1 to compare Sharpe ratios. We may for instance, be interested in whether the mutual fund s risk-adjusted return as captured by the Sharpe ratio is significantly better than the market s return. In Table 7 we present the 95% confidence interval for the difference between the Sharpe ratios for the mutual fund and market index. We can see that while the intervals produced using the MLE and likelihood ratio methods are similar, the interval obtained from the third-order method differs from these methods. As this is an example, we cannot comment on which interval is more accurate, but it is relevant to note the differences between the intervals which may be important in real world settings. The p-values for testing a null hypothesis of a zero difference between the Sharpe ratios of the mutual fund and market index are given by: , , and for the MLE, Likelihood Ratio and Proposed, respectively. As tail probabilities tend to be small probabilities, it is important to approximate these as accurately as we can. 3
24 Table 7: 95% Confidence Intervals for Sharpe Ratio Difference Method 95% CI MLE ( , ) Likelihood Ratio (-0.986, ) Proposed ( , ) 7. Simulations In this section we provide a simulation study to assess the performance of the third-order method relative the MLE and likelihood ratio. The size of each simulation is 10,000 and again the parameter values were chosen to mimic the example data: X N( , ), Y N( , ) and R f = Confidence intervals for the difference between the Sharpe ratios for funds X and Y were obtained for various sample sizes of each fund. Table 8 records the results from this simulation. Lower and upper probability errors are recorded as is central coverage. The nominal values and standard errors of the simulation are also reported for reference. As in the one sample case, these simulation results generally indicate that the propsed method outperforms the other methods based on the criteria we examined. 8 Conclusion A higher order likelihood method was applied for inference on the Sharpe ratio. The two-sample problem for comparing Sharpe ratios was also considered. This statistical method is known to be extremely accurate and has O(n 3/ ) distributional accuracy. The methodology was demonstrated and worked out explicitly for independently and identically normally distributed returns. Simulations were provided to show the exceptional accuracy of the method even for very small sample sizes. While our assumption of independently and identically normally distributed returns may seem restrictive, we stress that the methodology can be applied more generally for any parametric return distribution structure and our choice of normality was made to clearly illustrate the methodology s merits to the applied practitioner. A natural next step would be to consider inference when returns are not independently and identically distributed. For the two-sample problem, it would be fruitful to compare the approach taken by Ledoit and Wolf (008) which involves the studentized bootstrap with the proposed methodology. 4
25 Table 8: Results for Simulation Study CI n, m Method Lower Error Upper Error Central Coverage 90% n = m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error n = 1, m = 0 MLE Likelihood Ratio Proposed Nominal Standard Error n = 0, m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error % n = m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error n = 1, m = 0 MLE Likelihood Ratio Proposed Nominal Standard Error n = 0, m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error % n = m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error n = 1, m = 0 MLE Likelihood Ratio Proposed Nominal Standard Error n = 0, m = 1 MLE Likelihood Ratio Proposed Nominal Standard Error
26 References [1] Barndorff-Nielsen, O.E., 1980, Conditionality Resolutions, Biometrika 67, [] Barndorff-Nielsen, O.E., 1983, On a Formula for the Distribution of the Maximum Likelihood Estimator, Biometrika 70, [3] Barndorff-Nielsen, O.E., 1986, Inference on Full or Partial Parameters based on the Standardized Signed Log Likelihood Ratio, Biometrika 73, [4] Barndorff-Nielsen, O.E., Cox, D.R., 1994, Inference and Asymptotics, Chapman and Hall, UK. [5] Christie, S., 005, Is the Sharpe Ratio Useful in Asset Allocation? MAFC Research Papers No. 31, Applied Finance Centre, Macquarie University. [6] Cox, D.R., Hinkley, D.V., 1974, Theoretical Statistics, CRC Press, USA. [7] Daniels, H.E., 1954, Saddlepoint Approximations in Statistics, Annals of Mathematical Statistics 5, [8] Daniels, H.E., 1987, Tail Probability Approximations, International Statistical Review 55, [9] Doganaksoy, N., Schmee, J., 1993, Comparisons of Approximate Confidence Intervals for Distributions used in Life-Data Analysis, Technometrics 35, [10] Eling, M., Schuhmacher, F., 007, Does the Choice of Performance Measure Influence the Evaluation of Hedge Funds?, Journal of Banking & Finance 31(9), [11] Fraser, D., Reid, N., 1995, Ancillaries and Third Order Significance, Utilitas Mathematica 47, [1] Fraser, D., Reid, N., Wong, A., 1991, Exponential Linear Model: a Two-Pass Procedure for Saddlepoint Approximation, Journal of the Royal Statistical Society Series B 53(), [13] Fraser, D., Reid, N., Wu, J., 1999, A Simple General Formula for Tail Probabilities for Frequentist and Bayesian Inference, Biometrika 86(), [14] Jensen, J.L., 199, The Modified Signed Likelihood Statistic and Saddlepoint Approximations, Biometrika 79(4),
27 [15] Jobson, J.D., Korkie, B.M., 1981, Performance Hypothesis Testing with the Sharpe and Treynor Measures, The Journal of Finance 36(4), [16] Keating C., Shadwick, W.F., 00, An Introduction to Omega, The Finance Development Centre Limited. [17] Ledoit, O., Wolf, M., 008, Robust Performance Hypothesis Testing with the Sharpe Ratio, Journal of Empirical Finance 15, [18] Lo, A.W., 00, The Statistics of Sharpe Ratios, Financial Analysts Journal 58(4), [19] Lugannani, R., Rice, S.O., 1980, Saddlepoint Approximation for the Distribution of the Sum of Independent Random Variables, Advances in Applied Probability 1, [0] Mertens, E., 00, Comments on Variance of the IID estimator in Lo (00), Research Note, [1] Opdyke, J.D., 007, Comparing Sharpe ratios: So where are the p-values?, Journal of Asset Management 8(5), [] Sharpe, W.F., 1966, Mutual Fund Performance, The Journal of Business 39(1),
Comparing the Means of. Two Log-Normal Distributions: A Likelihood Approach
Journal of Statistical and Econometric Methods, vol.3, no.1, 014, 137-15 ISSN: 179-660 (print), 179-6939 (online) Scienpress Ltd, 014 Comparing the Means of Two Log-Normal Distributions: A Likelihood Approach
More informationSharpe Ratio Practice Note
Sharpe Ratio Practice Note Geng Deng, PhD, FRM Tim Dulaney, PhD Craig McCann, PhD, CFA Securities Litigation and Consulting Group, Inc. 3998 Fair Ridge Drive, Suite 250, Fairfax, VA 22033 June 26, 2012
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
More informationLecture 10: Point Estimation
Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,
More informationChapter 4: Asymptotic Properties of MLE (Part 3)
Chapter 4: Asymptotic Properties of MLE (Part 3) Daniel O. Scharfstein 09/30/13 1 / 1 Breakdown of Assumptions Non-Existence of the MLE Multiple Solutions to Maximization Problem Multiple Solutions to
More informationLikelihood Methods of Inference. Toss coin 6 times and get Heads twice.
Methods of Inference Toss coin 6 times and get Heads twice. p is probability of getting H. Probability of getting exactly 2 heads is 15p 2 (1 p) 4 This function of p, is likelihood function. Definition:
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More informationApplied Statistics I
Applied Statistics I Liang Zhang Department of Mathematics, University of Utah July 14, 2008 Liang Zhang (UofU) Applied Statistics I July 14, 2008 1 / 18 Point Estimation Liang Zhang (UofU) Applied Statistics
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More informationInterval estimation. September 29, Outline Basic ideas Sampling variation and CLT Interval estimation using X More general problems
Interval estimation September 29, 2017 STAT 151 Class 7 Slide 1 Outline of Topics 1 Basic ideas 2 Sampling variation and CLT 3 Interval estimation using X 4 More general problems STAT 151 Class 7 Slide
More informationChapter 8: Sampling distributions of estimators Sections
Chapter 8 continued Chapter 8: Sampling distributions of estimators Sections 8.1 Sampling distribution of a statistic 8.2 The Chi-square distributions 8.3 Joint Distribution of the sample mean and sample
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationA New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationCS340 Machine learning Bayesian model selection
CS340 Machine learning Bayesian model selection Bayesian model selection Suppose we have several models, each with potentially different numbers of parameters. Example: M0 = constant, M1 = straight line,
More information1 Residual life for gamma and Weibull distributions
Supplement to Tail Estimation for Window Censored Processes Residual life for gamma and Weibull distributions. Gamma distribution Let Γ(k, x = x yk e y dy be the upper incomplete gamma function, and let
More informationABSTRACT. involved therein. This paper has been motivated by the desire to meet the challenge of statistical estimation. A new estimator for
A Shorter-Length Confidence-Interval Estimator (CIE) for Sharpe-Ratio Using a Multiplier k* to the Usual Bootstrap-Resample CIE and Computational Intelligence Chandra Shekhar Bhatnagar 1, Chandrashekhar.Bhatnagar@sta.uwi.edu
More informationA Bayesian Control Chart for the Coecient of Variation in the Case of Pooled Samples
A Bayesian Control Chart for the Coecient of Variation in the Case of Pooled Samples R van Zyl a,, AJ van der Merwe b a PAREXEL International, Bloemfontein, South Africa b University of the Free State,
More informationAn Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process
Computational Statistics 17 (March 2002), 17 28. An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Gordon K. Smyth and Heather M. Podlich Department
More informationPoint Estimators. STATISTICS Lecture no. 10. Department of Econometrics FEM UO Brno office 69a, tel
STATISTICS Lecture no. 10 Department of Econometrics FEM UO Brno office 69a, tel. 973 442029 email:jiri.neubauer@unob.cz 8. 12. 2009 Introduction Suppose that we manufacture lightbulbs and we want to state
More informationOn Performance of Confidence Interval Estimate of Mean for Skewed Populations: Evidence from Examples and Simulations
On Performance of Confidence Interval Estimate of Mean for Skewed Populations: Evidence from Examples and Simulations Khairul Islam 1 * and Tanweer J Shapla 2 1,2 Department of Mathematics and Statistics
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationSharpe Ratio over investment Horizon
Sharpe Ratio over investment Horizon Ziemowit Bednarek, Pratish Patel and Cyrus Ramezani December 8, 2014 ABSTRACT Both building blocks of the Sharpe ratio the expected return and the expected volatility
More information**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More informationConfidence Intervals Introduction
Confidence Intervals Introduction A point estimate provides no information about the precision and reliability of estimation. For example, the sample mean X is a point estimate of the population mean μ
More informationChapter 7 - Lecture 1 General concepts and criteria
Chapter 7 - Lecture 1 General concepts and criteria January 29th, 2010 Best estimator Mean Square error Unbiased estimators Example Unbiased estimators not unique Special case MVUE Bootstrap General Question
More informationIntroduction to the Maximum Likelihood Estimation Technique. September 24, 2015
Introduction to the Maximum Likelihood Estimation Technique September 24, 2015 So far our Dependent Variable is Continuous That is, our outcome variable Y is assumed to follow a normal distribution having
More informationSTAT 830 Convergence in Distribution
STAT 830 Convergence in Distribution Richard Lockhart Simon Fraser University STAT 830 Fall 2013 Richard Lockhart (Simon Fraser University) STAT 830 Convergence in Distribution STAT 830 Fall 2013 1 / 31
More informationTesting Out-of-Sample Portfolio Performance
Testing Out-of-Sample Portfolio Performance Ekaterina Kazak 1 Winfried Pohlmeier 2 1 University of Konstanz, GSDS 2 University of Konstanz, CoFE, RCEA Econometric Research in Finance Workshop 2017 SGH
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.
More informationPoint Estimation. Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage
6 Point Estimation Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage Point Estimation Statistical inference: directed toward conclusions about one or more parameters. We will use the generic
More informationHomework Problems Stat 479
Chapter 10 91. * A random sample, X1, X2,, Xn, is drawn from a distribution with a mean of 2/3 and a variance of 1/18. ˆ = (X1 + X2 + + Xn)/(n-1) is the estimator of the distribution mean θ. Find MSE(
More informationOn Some Test Statistics for Testing the Population Skewness and Kurtosis: An Empirical Study
Florida International University FIU Digital Commons FIU Electronic Theses and Dissertations University Graduate School 8-26-2016 On Some Test Statistics for Testing the Population Skewness and Kurtosis:
More information6. Genetics examples: Hardy-Weinberg Equilibrium
PBCB 206 (Fall 2006) Instructor: Fei Zou email: fzou@bios.unc.edu office: 3107D McGavran-Greenberg Hall Lecture 4 Topics for Lecture 4 1. Parametric models and estimating parameters from data 2. Method
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationEstimating the Distribution of Sharpe Ratios
Estimating the Distribution of Sharpe Ratios Nadeem Aftab, Ingo Jungwirth, Tomás Sedliacik, Nadir Virk June 30, 2008 Abstract The Sharpe ratio is the most common measure for risk adjusted return of a nancial
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationTwo-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models
Economic Review (Otaru University of Commerce), Vo.59, No.4, 4-48, March, 009 Two-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models Haruhiko
More informationStatistical estimation
Statistical estimation Statistical modelling: theory and practice Gilles Guillot gigu@dtu.dk September 3, 2013 Gilles Guillot (gigu@dtu.dk) Estimation September 3, 2013 1 / 27 1 Introductory example 2
More informationPoint Estimation. Principle of Unbiased Estimation. When choosing among several different estimators of θ, select one that is unbiased.
Point Estimation Point Estimation Definition A point estimate of a parameter θ is a single number that can be regarded as a sensible value for θ. A point estimate is obtained by selecting a suitable statistic
More informationEstimation after Model Selection
Estimation after Model Selection Vanja M. Dukić Department of Health Studies University of Chicago E-Mail: vanja@uchicago.edu Edsel A. Peña* Department of Statistics University of South Carolina E-Mail:
More informationA Saddlepoint Approximation to Left-Tailed Hypothesis Tests of Variance for Non-normal Populations
UNF Digital Commons UNF Theses and Dissertations Student Scholarship 2016 A Saddlepoint Approximation to Left-Tailed Hypothesis Tests of Variance for Non-normal Populations Tyler L. Grimes University of
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Further Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Outline
More informationMonitoring Processes with Highly Censored Data
Monitoring Processes with Highly Censored Data Stefan H. Steiner and R. Jock MacKay Dept. of Statistics and Actuarial Sciences University of Waterloo Waterloo, N2L 3G1 Canada The need for process monitoring
More informationExercise. Show the corrected sample variance is an unbiased estimator of population variance. S 2 = n i=1 (X i X ) 2 n 1. Exercise Estimation
Exercise Show the corrected sample variance is an unbiased estimator of population variance. S 2 = n i=1 (X i X ) 2 n 1 Exercise S 2 = = = = n i=1 (X i x) 2 n i=1 = (X i µ + µ X ) 2 = n 1 n 1 n i=1 ((X
More informationInference of Several Log-normal Distributions
Inference of Several Log-normal Distributions Guoyi Zhang 1 and Bose Falk 2 Abstract This research considers several log-normal distributions when variances are heteroscedastic and group sizes are unequal.
More informationBayesian Linear Model: Gory Details
Bayesian Linear Model: Gory Details Pubh7440 Notes By Sudipto Banerjee Let y y i ] n i be an n vector of independent observations on a dependent variable (or response) from n experimental units. Associated
More informationImproved Inference for Signal Discovery Under Exceptionally Low False Positive Error Rates
Improved Inference for Signal Discovery Under Exceptionally Low False Positive Error Rates (to appear in Journal of Instrumentation) Igor Volobouev & Alex Trindade Dept. of Physics & Astronomy, Texas Tech
More informationOn Some Statistics for Testing the Skewness in a Population: An. Empirical Study
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017), pp. 726-752 Applications and Applied Mathematics: An International Journal (AAM) On Some Statistics
More informationmay be of interest. That is, the average difference between the estimator and the truth. Estimators with Bias(ˆθ) = 0 are called unbiased.
1 Evaluating estimators Suppose you observe data X 1,..., X n that are iid observations with distribution F θ indexed by some parameter θ. When trying to estimate θ, one may be interested in determining
More informationMeasuring Financial Risk using Extreme Value Theory: evidence from Pakistan
Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan Dr. Abdul Qayyum and Faisal Nawaz Abstract The purpose of the paper is to show some methods of extreme value theory through analysis
More informationIntroduction Dickey-Fuller Test Option Pricing Bootstrapping. Simulation Methods. Chapter 13 of Chris Brook s Book.
Simulation Methods Chapter 13 of Chris Brook s Book Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 April 26, 2017 Christopher
More informationJackknife Empirical Likelihood Inferences for the Skewness and Kurtosis
Georgia State University ScholarWorks @ Georgia State University Mathematics Theses Department of Mathematics and Statistics 5-10-2014 Jackknife Empirical Likelihood Inferences for the Skewness and Kurtosis
More informationSample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method
Meng-Jie Lu 1 / Wei-Hua Zhong 1 / Yu-Xiu Liu 1 / Hua-Zhang Miao 1 / Yong-Chang Li 1 / Mu-Huo Ji 2 Sample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method Abstract:
More informationDealing with forecast uncertainty in inventory models
Dealing with forecast uncertainty in inventory models 19th IIF workshop on Supply Chain Forecasting for Operations Lancaster University Dennis Prak Supervisor: Prof. R.H. Teunter June 29, 2016 Dennis Prak
More informationA New Multivariate Kurtosis and Its Asymptotic Distribution
A ew Multivariate Kurtosis and Its Asymptotic Distribution Chiaki Miyagawa 1 and Takashi Seo 1 Department of Mathematical Information Science, Graduate School of Science, Tokyo University of Science, Tokyo,
More informationWeek 7 Quantitative Analysis of Financial Markets Simulation Methods
Week 7 Quantitative Analysis of Financial Markets Simulation Methods Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 November
More informationMuch of what appears here comes from ideas presented in the book:
Chapter 11 Robust statistical methods Much of what appears here comes from ideas presented in the book: Huber, Peter J. (1981), Robust statistics, John Wiley & Sons (New York; Chichester). There are many
More informationINSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 SPECIMEN EXAMINATION
INSTITUTE AND FACULTY OF ACTUARIES Curriculum 2019 SPECIMEN EXAMINATION Subject CS1A Actuarial Statistics Time allowed: Three hours and fifteen minutes INSTRUCTIONS TO THE CANDIDATE 1. Enter all the candidate
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationLecture 10: Performance measures
Lecture 10: Performance measures Prof. Dr. Svetlozar Rachev Institute for Statistics and Mathematical Economics University of Karlsruhe Portfolio and Asset Liability Management Summer Semester 2008 Prof.
More informationFinancial Econometrics (FinMetrics04) Time-series Statistics Concepts Exploratory Data Analysis Testing for Normality Empirical VaR
Financial Econometrics (FinMetrics04) Time-series Statistics Concepts Exploratory Data Analysis Testing for Normality Empirical VaR Nelson Mark University of Notre Dame Fall 2017 September 11, 2017 Introduction
More informationExperience with the Weighted Bootstrap in Testing for Unobserved Heterogeneity in Exponential and Weibull Duration Models
Experience with the Weighted Bootstrap in Testing for Unobserved Heterogeneity in Exponential and Weibull Duration Models Jin Seo Cho, Ta Ul Cheong, Halbert White Abstract We study the properties of the
More informationModelling Returns: the CER and the CAPM
Modelling Returns: the CER and the CAPM Carlo Favero Favero () Modelling Returns: the CER and the CAPM 1 / 20 Econometric Modelling of Financial Returns Financial data are mostly observational data: they
More informationApplications of Good s Generalized Diversity Index. A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK
Applications of Good s Generalized Diversity Index A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK Internal Report STAT 98/11 September 1998 Applications of Good s Generalized
More informationChapter 7: Point Estimation and Sampling Distributions
Chapter 7: Point Estimation and Sampling Distributions Seungchul Baek Department of Statistics, University of South Carolina STAT 509: Statistics for Engineers 1 / 20 Motivation In chapter 3, we learned
More information12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.
12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance
More informationGPD-POT and GEV block maxima
Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,
More informationQuantitative Risk Management
Quantitative Risk Management Asset Allocation and Risk Management Martin B. Haugh Department of Industrial Engineering and Operations Research Columbia University Outline Review of Mean-Variance Analysis
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationDoes Naive Not Mean Optimal? The Case for the 1/N Strategy in Brazilian Equities
Does Naive Not Mean Optimal? GV INVEST 05 The Case for the 1/N Strategy in Brazilian Equities December, 2016 Vinicius Esposito i The development of optimal approaches to portfolio construction has rendered
More informationAnalysis of truncated data with application to the operational risk estimation
Analysis of truncated data with application to the operational risk estimation Petr Volf 1 Abstract. Researchers interested in the estimation of operational risk often face problems arising from the structure
More informationMulti-armed bandit problems
Multi-armed bandit problems Stochastic Decision Theory (2WB12) Arnoud den Boer 13 March 2013 Set-up 13 and 14 March: Lectures. 20 and 21 March: Paper presentations (Four groups, 45 min per group). Before
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationThe Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp
The Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp. 351-359 351 Bootstrapping the Small Sample Critical Values of the Rescaled Range Statistic* MARWAN IZZELDIN
More informationAn Information Based Methodology for the Change Point Problem Under the Non-central Skew t Distribution with Applications.
An Information Based Methodology for the Change Point Problem Under the Non-central Skew t Distribution with Applications. Joint with Prof. W. Ning & Prof. A. K. Gupta. Department of Mathematics and Statistics
More informationEE641 Digital Image Processing II: Purdue University VISE - October 29,
EE64 Digital Image Processing II: Purdue University VISE - October 9, 004 The EM Algorithm. Suffient Statistics and Exponential Distributions Let p(y θ) be a family of density functions parameterized by
More informationOptimizing Portfolios
Optimizing Portfolios An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Introduction Investors may wish to adjust the allocation of financial resources including a mixture
More informationRisk Measurement in Credit Portfolio Models
9 th DGVFM Scientific Day 30 April 2010 1 Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 9 th DGVFM Scientific Day 30 April 2010 2 Quantitative Risk Management Profit
More informationDecision theoretic estimation of the ratio of variances in a bivariate normal distribution 1
Decision theoretic estimation of the ratio of variances in a bivariate normal distribution 1 George Iliopoulos Department of Mathematics University of Patras 26500 Rio, Patras, Greece Abstract In this
More informationBack to estimators...
Back to estimators... So far, we have: Identified estimators for common parameters Discussed the sampling distributions of estimators Introduced ways to judge the goodness of an estimator (bias, MSE, etc.)
More informationMTH6154 Financial Mathematics I Stochastic Interest Rates
MTH6154 Financial Mathematics I Stochastic Interest Rates Contents 4 Stochastic Interest Rates 45 4.1 Fixed Interest Rate Model............................ 45 4.2 Varying Interest Rate Model...........................
More informationFinancial Econometrics
Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value
More informationModelling financial data with stochastic processes
Modelling financial data with stochastic processes Vlad Ardelean, Fabian Tinkl 01.08.2012 Chair of statistics and econometrics FAU Erlangen-Nuremberg Outline Introduction Stochastic processes Volatility
More informationTesting for the martingale hypothesis in Asian stock prices: a wild bootstrap approach
Testing for the martingale hypothesis in Asian stock prices: a wild bootstrap approach Jae H. Kim Department of Econometrics and Business Statistics Monash University, Caulfield East, VIC 3145, Australia
More informationContinuous random variables
Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),
More informationUnit 5: Sampling Distributions of Statistics
Unit 5: Sampling Distributions of Statistics Statistics 571: Statistical Methods Ramón V. León 6/12/2004 Unit 5 - Stat 571 - Ramon V. Leon 1 Definitions and Key Concepts A sample statistic used to estimate
More informationBIO5312 Biostatistics Lecture 5: Estimations
BIO5312 Biostatistics Lecture 5: Estimations Yujin Chung September 27th, 2016 Fall 2016 Yujin Chung Lec5: Estimations Fall 2016 1/34 Recap Yujin Chung Lec5: Estimations Fall 2016 2/34 Today s lecture and
More informationUnit 5: Sampling Distributions of Statistics
Unit 5: Sampling Distributions of Statistics Statistics 571: Statistical Methods Ramón V. León 6/12/2004 Unit 5 - Stat 571 - Ramon V. Leon 1 Definitions and Key Concepts A sample statistic used to estimate
More informationMVE051/MSG Lecture 7
MVE051/MSG810 2017 Lecture 7 Petter Mostad Chalmers November 20, 2017 The purpose of collecting and analyzing data Purpose: To build and select models for parts of the real world (which can be used for
More informationChapter 6 Forecasting Volatility using Stochastic Volatility Model
Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using SV Model In this chapter, the empirical performance of GARCH(1,1), GARCH-KF and SV models from
More informationChapter 7. Inferences about Population Variances
Chapter 7. Inferences about Population Variances Introduction () The variability of a population s values is as important as the population mean. Hypothetical distribution of E. coli concentrations from
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationA Two-Step Estimator for Missing Values in Probit Model Covariates
WORKING PAPER 3/2015 A Two-Step Estimator for Missing Values in Probit Model Covariates Lisha Wang and Thomas Laitila Statistics ISSN 1403-0586 http://www.oru.se/institutioner/handelshogskolan-vid-orebro-universitet/forskning/publikationer/working-papers/
More informationWeek 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals
Week 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg :
More informationPractice Exam 1. Loss Amount Number of Losses
Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000
More informationEVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS. Rick Katz
1 EVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu
More information