Stochastic portfolio theory explained

Stochastic portfolio theory (SPT) is a mathematical theory for analyzing stock market structure and portfolio behavior introduced by E. Robert Fernholz in 2002. It is descriptive as opposed to normative, and is consistent with the observed behavior of actual markets. Normative assumptions, which serve as a basis for earlier theories like modern portfolio theory (MPT) and the capital asset pricing model (CAPM), are absent from SPT.

SPT uses continuous-time random processes (in particular, continuous semi-martingales) to represent the prices of individual securities. Processes with discontinuities, such as jumps, have also been incorporated* into the theory (*unverifiable claim due to missing citation!).

Stocks, portfolios and markets

SPT considers stocks and stock markets, but its methods can be applied to other classes of assets as well. A stock is represented by its price process, usually in the logarithmic representation. In the case the market is a collection of stock-price processes

X _{i,} for i =1 , ... , n, each defined by a continuous semimartingaled log X _{i(t)} = \gamma _{i(t)} dt + \sum ^{d} _{\nu =1} \xi _{i} (t) dW _{\nu} (t) where W : = (W _{1,} ... , W _{d)} is an n -dimensional Brownian motion (Wiener) process with d \geq n , and the processes \gamma _{i} and \xi _{i} are progressively measurable with respect to the Brownian filtration\{ l{F} _{t\}} = \{ l{F} ^{W} _{t\}} . In this representation \gamma _{i(t)} is called the (compound) growth rate of X _{i,} and the covariance between log X _{i} and log X _{j} is \sigma _{ij} (t) = \sum ^{d} _{\nu =1} \xi _{i} (t) \xi _{j} (t). It is frequently assumed that, for all i, the process \xi ^{2 (t)} _{i, 1} + … + \xi ^{2 (t)} _{i d} is positive, locally square-integrable, and does not grow too rapidly as t → inf ty.

The logarithmic representation is equivalent to the classical arithmetic representation which uses the rate of return

\alpha _{i(t),} however the growth rate can be a meaningful indicator of long-term performance of a financial asset, whereas the rate of return has an upward bias. The relation between the rate of return and the growth rate is\alpha _{i} (t) = \gamma _{i(t)} + (\sigma _{ii} (t))/(2) The usual convention in SPT is to assume that each stock has a single share outstanding, so X _{i(t)} represents the total capitalization of the i -th stock at time t, and X(t) = X _{1(t)} + … + X _{n(t)} is the total capitalization of the market. Dividends can be included in this representation, but are omitted here for simplicity.

An investment strategy

\pi = (\pi _{1} , … , \pi _{n)} is a vector of bounded, progressively measurableprocesses; the quantity \pi _{i(t)} represents the proportion of total wealth invested in the i -th stock attime t , and \pi _{0(t)} : = 1 - \sum ^{n} _{i =1} \pi _{i(t)} is the proportion hoarded (invested in a money market with zero interest rate). Negative weights correspond to short positions. The cash strategy \kappa \equiv 0 (\kappa _{0} \equiv 1 ) keeps all wealth in the money market. A strategy \pi is called portfolio, if it is fully invested in the stock market, that is \pi _{1(t)} + … + \pi _{n} (t) = 1 holds, at all times.

The value process

Z _{\pi} of a strategy \pi is always positive and satisfiesd log Z _{\pi} (t) = \sum ^{n} _{i =1} \pi _{i(t)} d log X _{i(t)} + \gamma ^{*(t)} _{\pi} dt

where the process

\gamma ^{*} _{\pi} is called the excess growth rate process and is given by\gamma ^{*(t)} _{\pi} : = (1)/(2) \sum ^{n} _{i =1} \pi _{i(t)} \sigma _{ii} (t) - (1)/(2) \sum ^{n} _{i,j =1} \pi _{i(t)} \pi _{j(t)} \sigma _{ij} (t)

This expression is non-negative for a portfolio with non-negative weights

\pi _{i(t)} and has been usedin quadratic optimization of stock portfolios, a special case of which is optimization with respect to the logarithmic utility function.

The market weight processes,

\mu _{i(t)} : = (X _{i(t)})/(X _{1(t)} + … + X _{n(t)}) where i =1 , ... , n define the market portfolio \mu . With the initial condition Z _{\mu(0)} = X( 0 ), the associated value process will satisfy Z _{\mu} (t) = X(t) for all t.

A number of conditions can be imposed on a market, sometimes to model actual markets and sometimes to emphasize certain types of hypothetical market behavior. Some commonly invoked conditions are:

(\sigma _{ij} (t)) _{1} are bounded away from zero. It has bounded variance if the eigenvalues are bounded.

1. A market is coherent if

\operatorname{lim} _{t → infty} t ^{-1} log (\mu _{i(t))} = 0 for all i = 1 , ... , n.

1. A market is diverse on

[ 0 , T] if there exists \varepsilon > 0 such that \mu _{max} (t) \leq 1 - \varepsilon for t \in [ 0 , T].

1. A market is weakly diverse on

[ 0 , T] if there exists \varepsilon > 0 such that(1)/(T) \int ^{T} _{0} \mu _{max} (t) dt \leq 1 - \varepsilon

Diversity and weak diversity are rather weak conditions, and markets are generally far more diverse than would be tested by these extremes. A measure of market diversity is market entropy, defined by

S(\mu(t)) = - \sum ^{n} _{i =1} \mu _{i(t)} log (\mu _{i(t)).}

Stochastic stability

We consider the vector process

(\mu _{(1)} (t), ... , \mu _{(n)} (t)), with 0 \leq t < inf ty of ranked market weightsmax _{1\leq} \mu _{i(t)} = : \mu _{(1)} (t) \geq \mu _{(2)} (t) \geq … \mu _{(n)} (t) : = min _{1\leq} \mu _{i(t)} where ties are resolved “lexicographically”, always in favor of the lowest index. The log-gapsG ^{(k,k+1)} (t) : = log (\mu _{(k)} (t)/\mu _{(k+1)} (t)), where 0 \leq t < inf ty and k =1 , ... , n -1 are continuous, non-negative semimartingales; we denote by Λ ^{(k,k+1)} (t) = L ^{G ^{(k,}} (t ; 0 ) their local times at the origin. These quantities measure the amount of turnover between ranks k and k +1 during the time-interval [ 0 , t] .

A market is called stochastically stable, if

(\mu _{(1)} (t), … , \mu _{(n)} (t)) converges in distribution as t → inf ty to a random vector (M _{(1)} , … , M _{(n)} ) with values in the Weyl chamber \{(x _{1,} ... , x _{n} ) \mid x _{1} > x _{2} > ... > x _{n} and \sum ^{n} _{i =1} x _{i} = 1 \} of the unit simplex, and if the strong law of large numbers\lim _{t} (Λ ^{(k,k+1)} (t))/(t) = λ ^{(k,k+1)} > 0 holds for suitable real constants λ ^{(1,2)} , ... , λ ^{(n-1,} .

Arbitrage and the numeraire property

Given any two investment strategies

\pi,\rho and a real number T > 0 , we say that \pi is arbitrage relative to \rho over the time-horizon [ 0 ,T] , if P (Z _{\pi(T)} \geq Z _{\rho} (T)) \geq 1 and P (Z _{\pi(T)} > Z _{\rho} (T)) > 0 both hold; this relative arbitrage is called “strong” if P (Z _{\pi} (T) > Z _{\rho} (T)) = 1 . When \rho is \kappa \equiv 0 , we recover the usual definition of arbitrage relative to cash.We say that a given strategy \nu has the numeraire property, if for any strategy \pi the ratio Z _{\pi/Z} _{\nu} is a P −supermartingale. In such a case, the process 1 /Z _{\nu} is called a “deflator” for the market.

No arbitrage is possible, over any given time horizon, relative to a strategy

\nu that has the numeraire property (either with respect to the underlying probability measure P , or with respect to any other probability measure which is equivalent to P ). A strategy \nu with the numeraire property maximizes the asymptotic growth rate from investment, in the sense that\limsup _{T → infty} (1)/(log \left( T) (Z _{\pi(T)})/(Z _{\nu(T)}) \right) \leq 0 holds for any strategy \pi ; it also maximizes the expected log-utility from investment, in the sense that for any strategy \pi and real number T > 0 we haveE [ log (Z _{\pi(T)]} \leq E [ log (Z _{\nu(T))].} If the vector \alpha(t) = (\alpha _{1(t),} … , \alpha _{n(t))'} of instantaneous rates of return, and the matrix \sigma(t) = (\sigma(t)) _{1\leq} of instantaneous covariances, are known, then the strategy\nu(t) = \arg max _{p \in R ^{n}} (p' \alpha(t) - \tfrac{ 1 }{ 2 }p' \alpha(t) p) forall 0 \leq t < inf ty has the numeraire property whenever the indicated maximum is attained.

The study of the numeraire portfolio links SPT to the so-called Benchmark approach to Mathematical Finance, which takes such a numeraire portfolio as given and provides a way to price contingent claims, without any further assumptions.

A probability measure

Q is called equivalent martingale measure (EMM) on a given time-horizon [ 0 ,T] , if it has the same null sets as P on l{F} _{T} , and if the processes X _{1(t),} ... , X _{n(t)} with 0 \leq t \leq T are all Q −martingales. Assuming that such an EMM exists, arbitrage is not possible on [ 0 ,T] relative to either cash \kappa or to the market portfolio \mu (or more generally, relative to anystrategy \rho whose wealth process Z _{\rho} is a martingale under some EMM). Conversely, if \pi, \rho are portfolios and one of them is arbitrage relative to the other on [ 0 ,T] then no EMM can exist on this horizon.

Functionally-generated portfolios

Suppose we are given a smooth function

G : U → ( 0 , inf ty) on some neighborhood U of the unit simplex in R ^{n} . We call\pi ^{G} _{i} (t) : = \mu _{i(t)} \left( D _{i} log ( G (\mu(t))) + 1 - \sum ^{n} _{j =1} \mu _{j(t)} D _{j} log ( G (\mu(t))) \right) for 1 \leq i \leq n the portfolio generated by the function G . It can be shown that all the weights of this portfolio are non-negative, if its generating function G is concave. Under mild conditions, the relative performance of this functionally-generated portfolio \pi _{G} with respect to the market portfolio \mu , is given by the F-G decompositionlog \left( (Z _{\pi ^{G}} (T))/(Z _{\mu} (T)) \right) =log \left( (G (\mu(T)))/(G (\mu( 0 ))) \right) + \int ^{T} _{0} g(t) dt which involves no stochastic integrals. Here the expression g(t) : = (-1)/(2 G (\mu(t))) \sum ^{n} _{i =1} \sum ^{2} _{ij} G (\mu(t)) \mu _{i(t)} \mu _{j(t)} \tau ^{\mu(t)} _{ij} is called the drift process of the portfolio (and it is a non-negative quantity if the generating function G is concave); and the quantities \tau ^{\mu(t)} _{ij} : = \sum ^{n} _{\nu =1} (\xi _{i\nu} (t) - \xi ^{\mu(t))(\xi} _{j\nu} (t) - \xi ^{\mu(t)),} _{\nu} \xi _{i\nu} (t) : = \sum ^{n} _{i =1} \mu _{i(t)} \xi _{i\nu} (t)

with

1 \leq i, j \leq n are called the relative covariances between log (X _{i)} and log (X _{j)} with respect to the market.

Examples

1. The constant function

G : = w > 0 generates the market portfolio \mu ,

1. The geometric mean function

H (x) : = (x _{1} … x ^{(1)/(n)} _{n)} generates the equal-weighted portfolio \varphi _{i(n)} = (1)/(n) for all 1 \leq i \leq n ,

1. The modified entropy function

S ^{c(x)} = c - \sum ^{n} _{i =1} x _{i} ⋅ log (x _{i)} for any c > 0 generates the modified entropy-weighted portfolio,

1. The function

D ^{(p)} (x) : = (\sum ^{n} _{i =1} x ^{p)} _{i} ^{(1)/(p)} with 0 < p < 1 generates the diversity-weighted portfolio \delta ^{(p)} _{i} (t) = ((\mu ^{p} _{i(t))})/(\sum ^{p} _{i(t))}) with drift process ( 1- p)\gamma ^{*(t)} _{\delta ^{(p)}} .

Arbitrage relative to the market

The excess growth rate of the market portfolio admitsthe representation

2 \gamma ^{*(t)} _{\mu} = \sum ^{n} _{i =1} \mu _{i(t)} \tau ^{\mu(t)} _{ii} as a capitalization-weighted average relative stockvariance. This quantity is nonnegative; if it happens to be bounded away from zero, namely\gamma ^{*(t)} _{\mu} = (1)/(2) \sum ^{n} _{i =1} \mu _{i(t)} \tau ^{\mu(t)} _{ii} \geq h > 0 , for all 0 \leq t < inf ty for some real constant h , then it can be shown using the F-G decomposition that, for every T > S (\mu( 0 ))/h, there exists a constant c > 0 for which the modified entropic portfolio \Theta ^{(c)} is strict arbitrage relative to the market \mu over [ 0 , T] ; see Fernholz and Karatzas (2005) for details. It is anopen question, whether such arbitrage exists over arbitrary time horizons (for two special cases, inwhich the answer to this question turns out to be affirmative, please see the paragraph below andthe next section).

If the eigenvalues of the covariance matrix

(\sigma _{ij} (t)) _{1\leq} are bounded away from both zero and infinity, the condition \gamma ^{*\geq} _{\mu} h > 0 can be shown to be equivalent to diversity, namely \mu _{max} \leq 1 - \varepsilon for a suitable \varepsilon \in ( 0 , 1 ). Then the diversity-weighted portfolio \delta ^{(p)} leads to strict arbitragerelative to the market portfolio over sufficiently long time horizons; whereas, suitable modificationsof this diversity-weighted portfolio realize such strict arbitrage over arbitrary time horizons.

An example: volatility-stabilized markets

We consider the example of a system of stochastic differential equations

d log (X _{i(t))} = (\alpha)/(2 \mu _{i(t)}) dt + (\sigma)/(\mu _{i(t)}) dW _{i(t)} with 1 \leq i \leq n given real constants \alpha \geq 0 and an n -dimensional Brownian motion (W _{1} , ... , W _{n).} It follows from the work of Bass and Perkins (2002) that this system has a weak solution, which is unique in distribution. Fernholz and Karatzas (2005) show how to construct this solution in terms of scaled and time-changed squared Bessel processes, and prove that the resulting system is coherent.

The total market capitalization

X behaves here as geometric Brownian motion with drift, and has the same constant growth rate as the largest stock; whereas the excess growth rate of the marketportfolio is a positive constant. On the other hand, the relative market weights \mu _{i} with 1 \leq i \leq n have the dynamics of multi-allele Wright-Fisher processes. This model is an example of a non-diverse market with unbounded variances, in which strong arbitrage opportunities with respect to the market portfolio \mu exist over arbitrary time horizons, as was shown by Banner and Fernholz (2008). Moreover, Pal (2012) derived the joint density of market weights at fixed times and at certain stopping times.

Rank-based portfolios

We fix an integer

m \in \{ 2 , ... , n - 1 \} and construct two capitalization-weighted portfolios: one consisting of the top m stocks, denoted \zeta , and one consisting of the bottom n - m stocks, denoted η . More specifically,\zeta _{i(t)} = (\sum ^{m} _{k =1} \mu _{(k)} (t) 1 _{\{\mu _{i(t)} = \mu _{(k)} (t)\})/() }\qquad \text\eta_i(t) = \fracfor 1 \leq i \leq n. Fernholz (1999), (2002) showed that the relative performance of the large-stock portfolio with respect to the market is given aslog \left( (Z _{\zeta(t)})/(Z _{\mu(t)}) \right) =log \left( (\mu _{(1)} (T) + … + \mu _{(m)} (T))/(\mu _{(1)} ( 0 ) + … + \mu _{(m)} ( 0 )) \right) - (1)/(2) \int ^{T} _{0} (\mu _{(m)} (t))/(\mu _{(1)} (t) + … + \mu _{(m)} (t)) d Λ ^{(m,} (t). Indeed, if there is no turnover at the mth rank during the interval [ 0 , T] , the fortunes of \zeta relativeto the market are determined solely on the basis of how the total capitalization of this sub-universeof the m largest stocks fares, at time T versus time 0; whenever there is turnover at the m -th rank,though, \zeta has to sell at a loss a stock that gets “relegated” to the lower league, and buy a stockthat has risen in value and been promoted. This accounts for the “leakage” that is evident in thelast term, an integral with respect to the cumulative turnover process Λ ^{(m,m+1)} of the relative weight in the large-cap portfolio \zeta of the stock that occupies the mth rank.

The reverse situation prevails with the portfolio

η of small stocks, which gets to sell at a profit stocks that are being promoted to the “upper capitalization” league, and buy relatively cheaply stocks that are being relegated:log \left( (Z _{η(t)})/(Z _{\mu(t)}) \right) =log \left( (\mu _{(m+1)} (T) + … + \mu _{(n)} (T))/(\mu _{(m+1)} ( 0 ) + … + \mu _{(n)} ( 0 )) \right) + (1)/(2) \int ^{T} _{0} (\mu _{(m+1)} (t))/(\mu _{(m+1)} (t) + … + \mu _{(n)} (t)) . It is clear from these two expressions that, in a coherent and stochastically stable market, the small-stock cap-weighted portfolio \zeta will tend to outperform its large-stock counterpart η , at least overlarge time horizons and; in particular, we have under those conditions\lim _{T → infty} (1)/(log \left( T) (Z _{η(t)})/(Z _{\mu(t)}) \right) = λ ^{(m,} E \left( (M _{(1)})/(M _{(1)} + … + M _{(m)}) + (M _{(m+1)})/(M _{(m+1)} + … + M _{(n)}) \right) > 0 . This quantifies the so-called size effect. In Fernholz (1999, 2002), constructions such as these are generalized to include functionally generated portfolios based on ranked market weights.

First- and second-order models

First- and second-order models are hybrid Atlas models that reproduce some of the structure of real stock markets. First-order models have only rank-based parameters, and second-order models have both rank-based and name-based parameters.

Suppose that

X _{1,\ldots,X} _{n} is a coherent market, and that the limits\sigma ^{2} _{k} = \lim _{t\toinfty} t ^{-1} \langle log \mu _{(k)} \rangle (t)

and

g _{k} = \lim _{T\toinfty} (1)/(T) \int ^{n} _{i =1} 1 _{\{r _{t(i)=k\}} }\,d\log\mu_i(t)

exist for

k =1 ,\ldots,n , where r _{t(i)} is the rank of X _{i(t)} . Then the Atlas model {\widehat X} _{1,\ldots,{\widehat} X} _{n} defined byd log {\widehat X} _{i(t)} = \sum ^{n} _{k =1} g _{k1} _{\{{\hat} _{t(i)=k\}}dt +} \sum ^{n \sigma} _{k} 1 _{\{{\hat} _{t(i)=k\}}dW} _{i(t),}

where

{\hat r} _{t(i)} is the rank of \widehat X _{i(t)} and (W _{1,\ldots,W} _{n)} is an n -dimensional Brownian motion process, is the first-order model for the original market, X _{1,\ldots,X} _{n} .

Under reasonable conditions, the capital distribution curve for a first-order model will be close to that of the original market. However, a first-order model is ergodic in the sense that each stock asymptotically spends

( 1 /n) -th of its time at each rank, a property that is not present in actual markets. In order to vary the proportion of time that a stock spends at each rank, it is necessary to use some form of hybrid Atlas model with parameters that depend on both rank and name. An effort in this direction was made by Fernholz, Ichiba, and Karatzas (2013), who introduced a second-order model for the market with rank- and name-based growth parameters, and variance parameters that depended on rank alone.

References

• Fernholz, E.R. (2002). Stochastic Portfolio Theory. New York: Springer-Verlag.