David VEREDASOptimal portfolios with end-of-period target

We study the estimation of optimal portfolios for a Reserve Fund with an end-of-period target and when the returns of the assets that constitute the Reserve Fund portfolio follow two specifications. In the first one, assets are split into short memory (bonds) and long memory (equity), and the optimality of the portfolio is based on maximizing the Sharpe ratio. In the second, returns follow a conditional heteroskedasticity autoregressive nonlinear model, and we study when the distribution of the innovation vector is heavy-tailed stable. For this specification, we consider appropriate estimation methods, which include bootstrap and empirical likelihood. © 2012 Hiroshi Shiraishi et al.


Introduction
The Government Pension Investment Fund (GPIF) of Japan was established in April 1st 2006 as an independent administrative institution with the mission of managing and investing the Reserve Fund of the employees' pension insurance and the national pension (www.gpif.go.jp for more information).It is the world's largest pension fund ($1.4 trillions in assets under management as of December 2009) and it has mission of managing and investing the Reserve Funds in safe and efficient investment with a long-term perspective.Business management targets to be achieved by GPIF are set by the Minister of health, labour and welfare based on the Law on the general rules of independent administrative agencies.In the Actuarial science, "required resurve fund" for pension insurance has been investigated for a long time.The traditional approach focuses on the expected value of future obligations and interest rate.Then, the investment strategy is determined for exceeding the expected value of interest rate.Recently, solvency for the insurer is defined in terms of random values of future obligations (e.g.Olivieri and Pitacco (2000)).In this article, we assume that the reserve fund is defined in terms of the random interest rate and the expected future obligations.Then, we propose optimal portfolios by optimaizing the randomized reserve fund.
The GPIF invest in a portfolio of domestic and international stocks and bonds.In this article we consider the optimal portfolio problem of the Reserve Fund under two econometric specifications for the asset's returns.
First, we select the optimal portfolio weights based on the maximization of the Sharpe ratio under three different functional forms for the portfolio mean and variance, two of them depending on the Reserve Fund at the end-of-period target (about 100 years).Following the asset structure of the GPIF, we split the assets into cash and domestic and foreign bonds on one hand and domestic and foreign equity on the other.The first type of assets are assumed to be short memory while the second type are long memory.To obtain the optimal portfolio weights we rely on bootstrap.For the short memory returns, we use wild bootstrap (WB).Early work focuses on providing first-and second-order theoretical justification for the wild bootstrap in the classical linear regression model (see, e.g., Wu, 1986).Goncalves and Kilian (2004) show that WB is applicable for the linear regression model with conditional heteroscedastic such as stationary ARCH, GARCH, and stochastic volatility effects.For the long memory returns, we apply sieve bootstrap (SB).Bühlmann (1997) establishes consistency of the autoregressive sieve bootstrap.Assuming that the long memory process can be written as AR(∞) and MA(∞) processes, we estimate the long memory parameter by means of the Whittle's approximate likelihood (Beran, 1994).Given this estimator, the residuals are computed and resampled for the construction of the bootstrap samples, from which the optimal portfolio estimated weights are obtained.We study the usefulness of these procedures with an application to the GPIF assets.
Second, we consider the case when the returns are time dependent and follow a heavytailed.It is known that one of the stylized facts of financial returns are heavy tails.It is therefore reasonable to use the stable distribution, instead of the Gaussian, since it allows for skewness and fat tails.We couple this distribution with the conditional heteroskedasticity autoregressive nonlinear (CHARN) model that nests many well known time series models, such as ARMA and ARCH.We estimate the parameters and the optimal portfolio by means of empirical likelihood.
The article is organized as follows.Section 2 sets the Reserve Fund portfolio problem.Section 3 focuses on the first part, i.e. estimation in terms of the Sharp ratio and discusses the bootstrap procedure.Section 4 covers the CHARN model under stable innovations and the estimation by means of empirical likelihood.Section 5 concludes.

Reserve Funds Portfolio with end-of-period target
Let S i,t be the price of the i-th asset at time t (i = 1, . ..,k), and X i,t its log-return.Time runs from 0 to T .All over the article we consider that today is T 0 and T is the endof-period target.Hence the past and present observations run for t = 0, . . ., T 0 , and the future until the end-of-period target for t = T 0 + 1, . . ., T .The price S i,t can be written as where S i,0 is the initial price.Let F i,t denote the Reserve Fund on asset i at time t and defined by where c i,t denotes the maintenance cost at time t.By recursion F i,t can be written as where F i,0 = S i,0 .We gather the Reserve Funds in the vector F t = (F 1,t , . . ., F k,t ).Let F t (α) = α ′ F t be a portfolio form by the k Reserve Funds, which depend on the vector of weights α = (α 1 , . . ., α k ).The portfolio Reserve Fund can be expressed as a function of all past returns . (2.3) We are interested in maximizing F t (α) at the end-of-period target F T (α) )) . (2.4) It depends on the future returns, the maintenance cost and the portfolio weights.While the first two are assumed to be constant from T 0 to T (the constant return can be seen as the average return over the T − T 0 periods), we focus on the optimality of the weights that we denote by α opt .

Sharpe-ratio based optimal portfolios
In the first specification, the estimation of the optimal portfolio weights is based on the maximization of the Sharpe ratio: under different functional forms of the expectation µ(α) and the risk σ(α) of the portfolio.
We propose three functional forms, two of the them depending on the Reserve Fund.The first one is the traditional based on the returns: where E(X T ) and V (X T ) are the expectation and the covariance matrix of the returns at the end-of-period target.The second form for the portfolio expectation and risk is based on the vector of Reserve Funds: where E(F T ) and V (F T ) indicate the mean and covariance of the Reserve Funds at time T .Last, we consider the case where the portfolio risk depends on the lower partial moments of the Reserve Funs at the end-of-period target: where F is a given value.
Standard portfolio management rules are based on a mean-variance approach, for which risk is measured by the standard deviation of the future portfolio value.However, the variance often does not provide a correct assesment of risk under dependency and non-Gaussianity.To overcome this problem, various optimization models have been proposed such as mean-semivariance model, mean-absolute deviation model, mean-varianceskewness model ,mean-(C)VaR model and mean-lower partial moment model.The meanlower partial moment model is an appropriate model for reducing the influence of heavy tails.
The k returns are split into p and q-dimensional vectors {X S t ; t ∈ Z} and {X L t ; t ∈ Z}, where S and L stand for short and long memory respectively.The short memory returns correspond to cash and domestic and foreign bonds, which we generically denote by bonds.The long memory returns corresponds to domestic and foreign equity, which we denote as equity.
Cash and bonds follow the nonlinear model where µ S is a vector of constants, H : R mp → R p × R p is a positive definite matrix-valued measurable function, and ϵ S t = (ϵ S 1,t , . . ., ϵ S p,t ) are i.i.d.random vectors with mean 0 and covariance matrix Σ S .By contrast, equity returns follow a long memory nonlinear model where and We estimate the optimal portfolio weights by means of bootstrap.Let the superindexes (S, b) and (L, b) denote the bootstrapped samples for the bonds and equity respectively, and B the total number of bootstrapped samples.In the sequel we show the bootstrap procedure for both types of assets.

Bootstrap procedure for X (S,b) t
Step 1: Generate the i.i.d.sequences {ϵ Step 3: Compute {X , where ⊙ denotes the cellwise product.

Bootstrap procedure for X (L,b) t
Step 1: Estimate d from the observed returns by means of Whittle's approximate likelihood: d = arg min where Step 2: Compute {ε L t } for t = 1, . . ., T 0 : Step 3: Generate {ϵ Step 4: Generate {X We gather X And the bootstraped Reserve Fund portfolio is Finally, the estimated portfolio weights that give the optimal portfolio are αopt = arg max where µ (b) (α) and σ (b) (α) may taken any of the three forms introduced earlier, but evaluated in the bootstrapped returns or Reserve Funds.

An illustration
We consider monthly log-returns from January 31 1971 to October 31 2009 (466 observations) of the five types of assets considered earlier: Domestic Bond (DB), Domestic Equity (DE), Foreign Bond (FB), Foreign Equity (FE) and cash (cash).Cash and bonds are gathered in the short-memory panel ) and follow (3.2).Equity are gathered into the long-memory panel ) and follow (3.3). Figure 1 shows the five assets.Cash is virtually constant, and equity are significantly more volatile than bonds and with averages that are slightly higher than those of bonds.We estimate the optimal portfolio weights, αopt1 , αopt2 and αopt3 , corresponding to the three forms for the expectation and risk of the Sharpe ratio, and we compute the trajectory of the optimal Reserve Fund for t = T 0 + 1. . . ., T .Because of liquidity reasons, the portfolio weight for cash is kept constant to 5%.The target period is fixed to 100 years, and the maintenance cost is based on the 2004 Pension Reform.

Table 1 about here.
Table 1 shows the estimated optimal portfolio weights for the three different choices of the portfolio expectation and risk.The weight of domestic bonds is very high and clearly dominates over the other assets.Domestic bonds are low risk and medium return, which is in contrast with equity that show higher return but also higher risk, and with foreign bonds that show low return and risk.Therefore, in a sense, domestic bonds is a compromise between the characteristic of the four equities and bonds.
Figure 2 shows the trajectory of the future Reserve Fund for different values of the yearly return (assumed to be constant from T 0 + 1 to T ) ranging from 2.7% to 3.7%.Since the investment term is extremely long, 100 years, the Reserve Fund is quite sensitive to the choice of the yearly return.In the 2004 Pension Reform, authorities assumed a yearly return of the portfolio of 3.2%, which corresponds to the middle thick line of the Figure.

Optimal portfolio with time dependent returns and heavy tails
In this section we consider the second scenario where returns follow a dependent model with stable innovations.The theory of portfolio choice is mostly based on the assumption that investors maximize their expected utility.The most well-known utility is the Markowitz's mean-variance function that is optimal under Gaussianity.However, it is widely acknowledged that financial returns show fat tails and, frequently, skewness.Moreover, the variance may not always be the best risk measure.Since the purpose of GPIF is to avoid making a big loss at a certain point in future, risk measures that summarize the probability that the Reserve Fund is below the prescribed level at a certain future point, such as Value at Risk (VaR), are more appropriate (Shiohama et al. 2010).In addition, the traditional mean-variance approach considers that returns are i.i.d., which is not realistic as past information may help to explain today's distribution of returns.We need a specification that allows for heavy tails and skewness, and time dependencies.This calls for a general model with location and scale that are a function of past observations, and with innovations that are stable distributed.The location-scale model for the returns is the conditional heteroscedastic autoregressive nonlinear (CHARN), which is fairly general and it nests important models such as ARMA and ARCH.
Estimation of the parameters in a stable framework is not straightforward since the density does not have a closed form. 1 We rely on empirical likelihood, which is one of the nonparametric methods, as it has been already studied in this context (Ogata, 2009).Once the parameters are estimated, we simulate samples of the returns, which are used to compute the Reserve Fund at the end-of-period target, and estimate the optimal portfolio weights by means of minimizing the empirical VaR of the Reserve Fund at time T .
Suppose that the vector of returns X t follows the following CHARN model: where F µ : R kp → R k is a vector-valued measurable function with a parameter µ ∈ R p1 , H σ : R kq → R k × R k is a positive definite matrix-valued measurable function with a parameter σ ∈ R p2 .Each element of the vector of innovations ε t is standardized stable distributed: ∼ S(α i , β i , 1, 0).We set θ = (µ, σ, α, β) where α = (α 1 , . . ., α k ) and The stable distribution is often represented by its characteristic function: ) +iνδ where δ ∈ R is a location parameter, γ > 0 is a scale parameter, β ∈ [−1, 1] is a skewness parameter, and α ∈ (0, 2] is a characteristic exponent that captures the tail thickness of the distribution: the smaller α the heavier the tail.The distributions with α = 2 corresponds to the Gaussian.The existence of moments is given by α: moments of order higher than α do not exist, with the case of α = 2 being an exception, for which all moments exist.The lack of important moments may, in principle, render estimation by the method of moments difficult.However, instead of matching moments, it is fairly simple to match the theoretical and empirical characteristic function evaluated at a grid of frequencies (Ogata, 2009).Let be the residual of the CHARN model.If the parameters µ and σ are the true ones, the residuals ϵ i,t should be independently and identically distributed to S(α i , β i , 1, 0).So the aim is to find the estimated parameters such that the residuals are i.i.d. and stable distributed, meaning that their probability law is expressed by the above characteristic function.Or, in other words, estimate the parameters by matching the empirical and theoretical characteristic functions and minimizing their distance.Let J the number of frequencies at which we evaluate the characteristic function: ν 1 , . . ., ν J .That makes, in principle, a system of J matching equations.But since the characteristic function can be split into the real and imaginary parts, ϕ(ν) = E[cos(νε i,t )] + iE[sin(νε i,t )], we double the dimension of the system by matching these parts.Let Re ( ϕ(ν) ) and Im ( ϕ(ν) ) be the real and imaginary parts of the theoretical characteristic function, and cos(νϵ i,t ) and sin(νϵ i,t ) the empirical counterparts.The estimating functions are for each i = 1, . . ., k, and gather them into the vector The number of frequencies J, and the frequencies themselves, are chosen arbitrary.Feuerverger and McDunnough (1981) show that the asymptotic variance can be made arbitrarily close to the Cramér-Rao lower bound if the number of frequencies is sufficiently large and the grid is sufficiently fine and extended.Similarly Yu (section 2.1, 2004) argues that, from the viewpoint of the minimum asymptotic variance, many and fine frequencies is the appropriate.However, Carrasco and Florens (2002) show that too fine frequencies lead to a singular asymptotic variance matrix and we cannot calculate its inverse.Given the estimating functions (4.1), the natural estimator is constructed by GMM: where W is a weighting matrix defining the metric (its optimal choice is typically the inverse of the covariance matrix of ψ θ (ϵ t )), and the expectations are replaced by sample moments.GMM estimator can be generalized to the empirical likelihood estimator, which was originally proposed by Owen (1988) as a nonparametric methods of inference based on a data-driven likelihood ratio function (see also Owen, 2001, for a review and applications).It produces a better variance estimate in one step while, in general, the optimal GMM requires a preliminary step and a preliminary estimation of an optimal W matrix.We define the empirical likelihood ratio function for θ as where p = (p 1 , . . ., p T0 ), and the maximum empirical likelihood estimator is θ = arg max θ R(θ).Qin and Lawless (1994) show that this estimator is consistent, asymptotically Gaussian and with covariance matrix ) −1 where Once the parameters are estimated, we compute the optimal portfolio weights and the portfolio Reserve Fund at the end-of-period target.Because of a notational conflict, the weights are now denoted by a = (a 1 , . . ., a k ).And, for simplicity, we assume that there is no maintenance cost, so (2.4) simplifies to The procedure to estimate the optimal portfolio weights is as follows Step 1: For each asset i = 1, . . ., k, we simulate the innovation process εi,t based on the maximum empirical likelihood estimator ( αi , βi ).
Step 3: For a given portfolio weight a, we calculate the predicted values of fund at time T , F T (a), with (4.2).
Step 4: We repeat Step 1 -Step 3 M times and save F (a). Then we calculate the proportion that the predicted values fail below the prescribed level F , that is Step 5: Minimize g(a) with respect to a: a * = arg min a g(a).

An illustration
In this section, we apply the above procedure to real financial data.We consider the same monthly log-returns data in Section 3.1.Domestic Bond (DB), Domestic Equity (DE), Foreign Bond (FB) and Foreign Equity (FE) are assumed to follow the following ARCH(1) model respectively.Here b > 0 and ε t i.i.d.
∼ S(α, β, 1, 0).Cash is virtually constant so we assume the log-return of cash as 0, permanently.Set the present Reserve Fund F T0 = 1 and the target period is fixed to half years.

Table 2 about here.
Table 2 shows the estimated optimal portfolio weights for the different prescribed level F .The weights of domestic and foreign bonds tend to be high when F is small.Small F implies that we want to avoid the loss.On the contrary, the weights of equities become higher when F is larger.Large F implies that we do not want to miss the chance of big gain.This result seems to be natural because bonds are lower risk (less volatile) than equities.

Conclusions
In this article, we study the estimation of optimal portfolios for a Reserve Fund with an end-of-period target in two different setting.In the first setting, one assets are split into short memory (bonds) and long memory (equity), and the optimality of the portfolio is based on maximizing the Sharpe ratio.The simulation result shows that the portfolio weight of domestic bonds is quite high.The reason is that the investment term is extremely long (100 years).Because the investment risk for the Reserve Fund is exponentially amplified year by year, the portfolio selection problem for the Resurve Fund is quite sensitive to the year based portfolio risk.In the second setting, returns follow a conditional heteroskedasticity autoregressive nonlinear model and we study when the distribution of the innovation vector is heavy-tailed stable.Simulation studies show that we prefer the bonds when we want to avoid the big loss in the future.The result seems to be natural because the bonds are less volatile than the equities.

Figure 1
Figure 1 about here.

Figure 2
Figure 2 about here.