Next Article in Journal
Evaluating Volatility Using an ANFIS Model for Financial Time Series Prediction
Previous Article in Journal
Mapping Financial Connections: Market Integration in Emerging Economies through Graph Theory
Previous Article in Special Issue
Quick Introduction into the General Framework of Portfolio Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modifying Sequential Monte Carlo Optimisation for Index Tracking to Allow for Transaction Costs

by
Leila Hamilton-Russell
1,
Thomas Malan O’Callaghan
1,
Dmitrii Savin
2 and
Erik Schlögl
1,3,4,*
1
The African Institute of Financial Markets and Risk Management (AIFMRM), University of Cape Town, Cape Town 7701, South Africa
2
Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
3
School of Mathematical and Physical Sciences, University of Technology Sydney, Ultimo, NSW 2007, Australia
4
Faculty of Science, Department of Statistics, University of Johannesburg, Johannesburg 2006, South Africa
*
Author to whom correspondence should be addressed.
Risks 2024, 12(10), 155; https://doi.org/10.3390/risks12100155
Submission received: 14 March 2024 / Revised: 17 September 2024 / Accepted: 18 September 2024 / Published: 30 September 2024
(This article belongs to the Special Issue Portfolio Theory, Financial Risk Analysis and Applications)

Abstract

:
Managing a portfolio whose value closely tracks an index by trading only in a subset of the index constituents involves an NP-hard optimisation problem. In the prior literature, it has been suggested that this problem be solved using sequential Monte Carlo (SMC, also known as particle filter) methods. However, this literature does not take transaction costs into account, although transaction costs are the primary motivation for attempting to replicate the index by trading in a subset, rather than the full set of index constituents. This paper modifies the SMC approach to index tracking to allow for proportional transaction costs and implements this extended method on empirical data for a variety stock indices. In addition to providing a more practically useful tracking strategy by allowing for transaction costs, we find that including a penalty for transaction costs in the optimisation objective can actually lead to better tracking performance.

1. Introduction

Index tracking, also known as index replication, is a passive investment strategy that has gained popularity due to the higher costs associated with active investing and scepticism regarding the ability of active funds to outperform the market (see e.g., Rompotis (2013), Barber and Odean (2000)). Financial indices track the performance of a collection of financial assets or an overall market. Indices are not tradable instruments themselves; however, replication can be achieved by creating a portfolio of assets that closely match the performance of an index. To perfectly replicate the performance of an index, one would need to invest in all its constituent stocks in the appropriate proportions. This may require continuously trading in hundreds of assets, which is both costly and often infeasible. More specifically, index replication’s effectiveness can be influenced by the constituent assets’ liquidity. Furthermore, continuous trading can lead to high transaction costs due to the costs incurred with each trade and the adjustments made during the rebalancing of the portfolio’s composition. The goal of sparse index tracking is to create a portfolio that closely replicates the performance of a benchmark index while avoiding the need to invest in all its constituent assets. The desired solution is a portfolio of a relatively small subset of assets that replicates an index at a lower cost.
Satpathy and Shah (2022) proposed a Sequential Monte Carlo (SMC) algorithm to address the sparse index tracking problem. The key idea is to interpret the index-tracking problem as a maximisation problem over the set of all possible (fixed-size) subsets of stocks composing some index and then using SMC as a global optimisation technique. Despite the novelty and the success of this approach when compared to other methods, such as Tibshirani (1996), Xu et al. (2015), Benidis et al. (2018), it is not clear whether this technique is an optimal index-tracking strategy, as it falls short in addressing the equally significant objective of reducing transaction costs. Indeed, when rebalancing the tracking portfolio, the algorithm does not explicitly minimise transaction costs. It focuses solely on minimising the tracking error between the portfolio and the index returns. This approach may lead to marginal improvements in the tracking error, but the cost associated with the portfolio rebalancing could outweigh any benefit gained.
The purpose of this work is to improve the Sequential Monte Carlo algorithm proposed by Satpathy and Shah (2022) by accounting for transaction costs in portfolio rebalancing. To this end, we redefine the dynamics of the index-tracking portfolio in such a way as to account for (proportional) transaction costs. Relying on this, we first reformulate the problem by adding a penalisation term to the objective function of the index-tracking problem and then study this modified problem by adopting similar techniques.
The remainder of the report is structured as follows. The general formulation of the index tracking problem, as proposed by Satpathy and Shah (2022), is presented in Section 2. In Section 3, a description of how to adjust the index tracking problem for transaction costs is given. Section 4 details the methods used for model implementation and testing, followed by the results in Section 5. Finally, the report concludes in Section 6 with a brief summary of the outcomes as well as a short description of possible further studies.

2. Preliminaries

2.1. The Index Tracking Problem

In this section, we recall the general formulation of the index–tracking problem. Let us consider an index I comprising n stocks. Denote by r I ( t ) R T the column vector of daily returns of the index over a period of T days before day t. For each constituent stock i, let r i ( t ) R T be the corresponding returns and set R ( t ) : = [ r 1 ( t ) , , r n ( t ) ] R T × n . The index tracking problem at time t consists of finding β * ( t ) R n , where β * ( t ) = ( β 1 * ( t ) , , β n * ( t ) ) , such that
β * ( t ) : = arg min β R n | | r I ( t ) R ( t ) β | | 2 2
s . t . 0 β 1
β T 1 = 1
| | β | | 0 = p .
Here, time t is the rebalancing point. The composition of the tracking portfolio is based on market data up until time t, with a predefined lookback period. Thus, the optimisation problem (1) uses only the subset of the data specific to that lookback period. Constraint (3) represents a full budget constraint, while (2) imposes a constraint on β that limits investors to taking only long positions. The cardinality constraint (4) restricts the number of stocks used to track the index to a fixed number p N . An increase in p might lead to a reduction in tracking errors, but possibly incurring larger transaction costs. Thus, p allows us to control the trade-off between tracking error and trading costs.
The cardinality constraint for β is non-convex and gives rise to an NP-hard problem. A wealth of methods exist in the literature to solve this. Popular methods include applying L 1 / 2 regularisation, using a Lagrangian version with an approximated L 0 -norm, replacing the L 0 -norm with an L 1 -norm (the LASSO technique), employing Mixed-Integer Quadratic optimisation and various genetic algorithm methods. We refer to Satpathy and Shah (2022) and the references therein for a more general discussion on this problem as well as an overview of these methods.

Optimisation Problem as Sampling Problem

In Satpathy and Shah (2022), a Sequential Monte Carlo (SMC) sampler is used to solve the optimisation problem (1). There, the key idea is to interpret the L 0 optimisation problem as a maximisation problem over the set of all possible (fixed-size) permutations of stocks composing some index and then using SMC sampling as a global optimisation technique. Hereafter, we briefly recall the main idea in Satpathy and Shah (2022); see also Duan (2019).
The index tracking problem with a L 0 -norm penalty is a regression problem where the objective function in Equation (1) is minimised. Similarly, one can consider the following maximisation problem:
β * ( t ) : = arg max β R n exp ( | | r I ( t ) R ( t ) β | | 2 2 ) s . t . 0 β 1 β T 1 = 1 | | β | | 0 = p .
Introducing the L 0 -norm penalty translates to finding the optimal permutation P * ( t ) of p stocks (a fixed number of regressors) that solves
P * ( t ) : = arg max P exp ( | | r I ( t ) R P ( t ) β ^ P | | 2 2 ) ,
where R P is the submatrix of R corresponding to columns in R for p regressors and β ^ P is given via1
β ^ P ( t ) : = arg min β p R p | | r I ( t ) R P ( t ) β P | | 2 2 s . t . 0 β P 1 β P T 1 = 1 .
Moreover, since for every fixed permutation P ( t ) , exp ( | | r I ( t ) R P ( t ) β ^ P ( t ) | | 2 2 ) > 0 and bounded above by 1, if a normalisation constant η > 0 is introduced, we can interpret this as a discrete probability distribution function over the permutation space. The target distribution to be maximised is then expressed as
T ( P ( t ) ) : = exp ( | | r I ( t ) R P ( t ) ( t ) β ^ P ( t ) ( t ) | | 2 2 ) η .
The optimal composition of p stocks is then the permutation that maximises this T ( P ( t ) ) .

2.2. Sequential Monte Carlo Sampling

A naive approach to solve problem (6) would comprise drawing samples from T ( P ) and choosing the optimal P that maximises it. However, it is not easy to sample from T ( P ) . For such a case, T ( P ) can be empirically approximated via SMC methods. SMC methods sample particles2 from an initial proposal distribution, after which a sequence of (sequential) importance sampling, resampling, and support boosting steps are employed. A collection of weighted particles is then obtained, which then represents the approximation of the target distribution. In what follows, we describe in more detail the steps within the SMC algorithm as employed in Satpathy and Shah (2022). See also Del Moral et al. (2006) for a more general discussion of the SMC method.

2.2.1. Sequential Importance Sampling

The first step in an SMC algorithm is importance sampling. It can be described as follows: N particles, denoted by ( P i ) i = 1 N , are sampled from an easy-to-sample distribution I ( P ) , also called proposal distribution. Then, for each i = 1 , , N , normalised importance weights w i are computed:
w i : = T ( P i ) I ( P i ) j = 1 N T ( P j ) I ( P j ) 1 .
The pair ( P i , w i ) i = 1 N completely characterises the target distribution approximation, and all quantities of interest from the distribution can be computed.
The quality of the sample drawn depends on the choice of I ( P ) , in particular on its closeness to the target distribution. Moreover, a necessary condition for I ( P ) to be a good initial choice is that its support includes the support of T ( P ) .
A tool to gauge the quality of the sample (and thus of the proposal distribution) is given by the computation of its Effective Sample Size defined by
ESS : = 1 i = 1 N w i 2 .
Since the weights w i are normalised, the ESS value will range from 1 to N, where we remember that N is the number of samples drawn from the proposal distribution. An ESS value of 1 indicates a poor choice of I ( P ) , contrary to a value equal to N, which stands for a good choice of I ( P ) . Indeed, this happens if, for instance, each importance weight is equal to 1 N .
However, when approximating a target distribution that is high-dimensional and complex, the simple importance sampling method tends to yield poor approximations (Satpathy and Shah (2022)). To account for this, Sequential Importance Sampling (SIS) methods are used. They aim to gradually approximate the target distribution with a sequence of approximating distributions.
SIS methods are commonly employed in filtering problems, where a sequence of new information arrives sequentially. However, in the context of this report, a sequence of new information used to approximate T ( P ) does not occur in the same manner. To account for this, density tempering is introduced.
Density tempering is a process proposed by Del Moral et al. (2006) which allows for sampling from target distributions with higher dimensions.
Its main ideas can be summarised as follows: a collection of synthetic, intermediate target distributions ( T γ i ( P ) ) i = 0 M are introduced, where for all i = 1 , , M , γ i [ 0 , 1 ] with 0 = γ 0 < γ 1 < γ 2 < γ M = 1 . The sequence of γ -modulated distribution is computed as follows:
γ i , T γ i ( P ) = I ( P ) T ( P ) I ( P ) γ i .
Notice that T γ 0 = I ( P ) and T γ M = T ( P ) , corresponding to the proposal and (true) target distribution. Similarly, the importance weights between each intermediate step are calculated according to the following formula:
w γ j ( P ) : = w γ i ( P ) T ( P ) I ( P ) γ j γ i , γ i < γ j .
The initial proposal distribution I ( P ) suggested by Satpathy and Shah (2022) is described as follows. Consider the regression coefficient of determination R j 2 for the jth stock returns (vis-à-vis the index to be replicated). A high R j 2 would imply that stock j is more likely to appear in the final P that maximises T ( P ) . A sampling strategy without replacement is considered, and the probability of the first jth stock to be chosen is described as
q j = R j 2 j = 1 n R j 2
Similarly, the probability of choosing the next mth stock in a set of n 1 stocks is q m 1 q j . Following this logic, I ( P ) for a given permutation P of p stocks is computed as follows:
I ( P ) = q 1 q 2 q p ( 1 q 1 ) ( 1 ( q 1 + q 2 ) ) ( 1 i = 1 p 1 q i )
Notice that I ( P ) described above is sequence-dependent. The order in which stocks are chosen to compute I ( P ) matters. However, this is not the case for the regression solution that we consider. Indeed, the value of T ( P ) is the same for each combination of p stocks despite the order in which they appear within the permutation. We are only concerned with the choice of stocks used to track the index and their proportion of the current portfolio wealth (their respective weights).

2.2.2. Resampling

As T γ ( P ) evolves, and with each re-weighting process, more variability is introduced, and fewer samples retain significant weights. This causes weight degeneracy and consequently a low ESS value. To overcome this phenomenon, one gets rid of samples with low-importance weights and repeatedly samples those with high-importance weights. This is in practice carried out by applying some resampling methods (see e.g., Chen (2003), Speekenbrink (2016), Gellert and Schlögl (2021)) whenever the ESS value falls below a certain threshold (commonly set to N/2). In Satpathy and Shah (2022), a multinomial resampling technique is implemented. Here, we instead consider the systematic resampling technique,3 which can be described as follows. Fix the number of particles to be resampled to r N . For k { 1 , , r } , and let
u k : = ( k 1 ) + u ˜ r with u ˜ U [ 0 , 1 ) .
Consider the m-th resampling in the algorithm. Replace particle k with particle i satisfying the following condition:
u k j = 1 i 1 w j ( m ) , j = 1 i w j ( m ) ,
where w j ( m ) denotes the weight corresponding to the particle j on the m-th resampling.

2.2.3. Support Boosting and Metropolis–Hastings Algorithm

Since, in the resampling step, samples with high weights are repeatedly resampled, this method leads to a loss in particle diversity, causing the support of the sequence of distributions represented by the particles and their weights to shrink. To account for this problem, an additional step is added with the intention of boosting the support. To this end, several moves of the Metropolis–Hastings (MH) algorithm are performed in Satpathy and Shah (2022). We do not employ the support-boosting step, since it results in a prohibitive and stochastic run time; however, we briefly describe the method as follows: after each resampling step, the MH algorithm is used for sampling from the current intermediate target distribution T γ ( P ) by using a new distribution for proposing a new sample P and then accepting or rejecting the new sample with some probability.
The MH probability of acceptance, that is, the probability of replacing the current particle P with the new sample P , is computed as follows:
α γ ( P P ) = min 1 , T γ ( P ) T γ ( P ) h w ( P | P ) h w ( P | P ) ,
where h ω ( P | P ) denotes the proposal density. Notice that if the ratio T γ i ( P ) T γ i ( P ) is high, there is a greater chance of moving from P to P .
The acceptance or rejection step in the MH algorithm is applied until the accumulative acceptance rate reaches 500% to ensure the support of the distribution is sufficiently boosted.
The choice of the proposal distribution4 h is performed as follows. Initially, a count-based probability Q γ i ( P ) is considered, meaning that the probability of choosing a stock to be sampled in a permutation is proportional to the number of times the stock appears in the current sample of permutations. Notice that this choice reflects the relative importance of particles after the SMC algorithm has reached the stage indicated by the current γ . Then, the proposal distribution h is defined as a weighted sum of I ( P ) , the initial proposal distribution described in Section 2.2.1, and Q γ i ( P ) :
h ω ( P | P ) = ω h Q ( P | P ) + ( 1 ω ) h I ( P | P ) ,
where ω [ 0 , 1 ] . Moreover, the distribution employed in Satpathy and Shah (2022) h is defined in such a way that only a given subset A of the permutation is replaced. More precisely, to describe the computation of h Q ( P | P ) and h I ( P | P ) as in Satpathy and Shah (2022), consider a subset A of the permutation P that we wish to replace in order to transform P to P . Since only a subset of P is replaced, for every P , it holds that P A = P A . Thus, we sample A from the set of stocks that exclude stocks in P A . In particular, if the set of all stocks is S , we sample A from S P A . We then define h Q ( P | P ) as the probability of sampling A from S P A given the count-based probability Q γ i ( P ) described above. Similarly, h I ( P | P ) can be computed based on I ( P ) applied to the set of S P A stocks.
Finally, Satpathy and Shah (2022) set the number of stocks for constructing an index tracking portfolio as a constant, p, representing 20% of the total number of stocks in the index. We instead use Principal Component Analysis (PCA) to set the number of stocks chosen for the index tracking portfolio. The variance explained by the principal components is set to 95%. The variable p is then assigned a value equal to the number of principal components computed under this constraint. The heuristic (which of course represents an approximation in our context) behind this approach is that indices are linear combinations of stock prices, and if stock price dynamics were fully spanned by p linear factors, p stocks would be sufficient to replicate the index (assuming constant index weights).

Summary of SMC Algorithm

  • Fix the number of stocks used to replicate the index to p (using PCA).
  • Draw n permutations of size p from T γ = 0 ( P ) = I ( P ) , for I ( P ) computed as in Equation (9). Set the importance weights for each permutation to 1 N .
  • Choose δ and increment γ such that γ = γ + δ .
  • At γ + δ , compute the importance weights for each permutation using recursive formulas in Section 2.2.1.
  • Compute ESS
  • While γ 1 ,
    (a)
    If ESS N 2 proceed to Step 3.
    (b)
    If ESS < N 2 . Resample the particles (Satpathy and Shah (2022)) suggest that if support boosting is applied, initialise MH with an accumulative acceptance rate of 500%).
    (c)
    Set the weights of new samples computed in MH to 1 N and proceed to Step 3.
  • At γ = 1 , samples from T ( P ) are obtained, represented as pairs of permutations and their respective weights.
  • Resample once more to obtain new permutations with equal weights.
  • Compute T ( P ) for the new sample of permutations to obtain the optimal P * that maximises the T ( P ) given by Equation (8).

3. Transaction Costs

As briefly mentioned in Section 1, in Satpathy and Shah (2022), transaction costs have not been considered when rebalancing the index tracking portfolio. This can potentially cause stocks to enter and exit the tracking portfolio with minimal benefit but at an increased cost. In order to consider this aspect, in this section, we reformulate the index-tracking problem by adding an additional term to the objective function, which results in a more gradual rebalancing of the index-tracking portfolio.

3.1. Dynamics of the Wealth Process Given Proportional Transaction Costs

We start by describing the evolution of the value of the index-tracking portfolio over rebalancing times. Consider the rebalancing time point t. Let us also denote the money invested in the j-th stock at time t as ψ j ( t ) . If transaction costs are not considered, the value of the index tracking (self-financing) portfolio X, for an index comprising n stocks, is given by
X ( t ) = X ( t ) = j = 1 n ψ j ( t ) ,
where t stands for the moment just before the rebalancing.
If transaction costs are instead allowed, we model the evolution of X ( t ) as follows:
X ( t ) = X ( t ) · C
where C ( 0 , 1 ) is the proportion of our wealth spent on the transaction costs at t. Since the portfolio wealth is lower due to transaction costs, we take into account that there is now less money to invest in each stock. Here, the portfolio is constructed in such a way that after paying the transaction costs, the portfolio weights are strictly equal to the target portfolio weights: ( β 1 ( t ) , , β n ( t ) ) . Thus, the following equation represents the value of the portfolio after rebalancing:
X ( t ) = X ( t ) · C j = 1 n β j ( t ) = j = 1 n ψ j ( t ) = C j = 1 n ψ j ( t ) .
Recall that only the p of all β ’s are not zero. Now recall that we pay a transaction cost for each unit of money we buy or sell. Let ε be the percentage we pay per unit of money we buy or sell, i.e., the rate of transaction costs. We assume here that the transaction costs for buying and selling are the same. We will incur a transaction cost of ε enforced on the difference in the amount of money invested in each stock before and after the rebalancing. In some sense, it means that instead of buying a stock for USD 1, we buy it for USD (1 + ε ).
It was demonstrated in Ruf and Xie (2019) that knowing the target portfolio weights, ( β 1 ( t ) , , β n ( t ) ) , and the amount of money invested in each stock prior to rebalancing, ( ψ 1 ( t ) , , ψ n ( t ) ) , is sufficient to compute C. Refer to Ruf and Xie (2019) for a detailed approach on the computation of C. Here, we introduce notations and ideas from Ruf and Xie (2019) necessary to present the exact formula for C.
Recall that the portfolio is self-financing; that is, the amount of currency used to buy extra stocks should be exactly the amount of currency obtained from selling redundant stocks, and hence, the following should be satisfied:
( 1 + ε ) j = 1 n ( ψ j ( t ) ψ j ( t ) ) + = ( 1 ε ) j = 1 n ( ψ j ( t ) ψ j ( t ) ) + ,
where dividends are not taken into account. We introduce D ^ and c j as follows:
D ^ = ( 1 ε ) j = 1 n ψ j ( t ) 1 β j ( t ) = 0 X ( t ) ,
c j = β j ( t ) β j ( t ) 1 β j ( t ) > 0 .
If we now divide both parts of (12) by X ( t ) we will obtain the following equation:
( 1 + ε ) j = 1 n ( C c j ) + β j ( t ) = ( 1 ε ) j = 1 n ( c j C ) + β j ( t ) + D ^ ,
where the LHS of (15) is a continuous function of C and strictly increasing from 0 to 1, as C changes from min i { 1 , , n } c i to . Conversely, the RHS of (15) is a continuous function of C strictly decreasing from to D ^ 0 , as C changes from to max i { 1 , , n } c i , and equals to D ^ as C changes from max i { 1 , , n } c i to . Therefore, both sides of Equation (15), treated as functions of C, must intersect at a unique point, indicating the existence of a unique solution for Equation (15).
Given the above, an expression for C now follows:
C ( β ) = ( 1 + ε ) j = 1 n c j β j ( t ) 1 c j c ^ + ( 1 ε ) j = 1 n β j ( t ) 1 c j > c ^ + D ^ ( 1 + ε ) j = 1 n β j ( t ) 1 c j c ^ + ( 1 ε ) j = 1 n β j ( t ) 1 c j > c ^ .
Here, c ^ = c k , where k = arg max k ( 1 , , n ) { D ^ k ; D ^ k D ^ } with D ^ k defined as follows:
D ^ k = ( 1 + ε ) j = 1 n ( c k c j ) + β j ( t ) ( 1 ε ) j = 1 n ( c j c k ) + β j ( t ) .
Using these formulas, we will be able to update the wealth process after each rebalancing step. The formula was given in (11).

3.2. Index Tracking Problem with Proportional Transaction Costs

Our goal is to minimise the transaction costs incurred at each rebalancing step, i.e., to minimise the money lost due to transaction costs, ( 1 C ) X ( t ) . In order to do that, we will introduce the updated version of the regression problem from Satpathy and Shah (2022). The objective function is the following:
β * ( t ) = arg min β ( t ) | | r I ( t ) R ( t ) β ( t ) | | 2 2 + λ | | β ( t ) β ( t ) | | 2 2 ,
where λ > 0 is a flexible parameter needed to model the sensitivity of the algorithm to transaction costs. The intuition behind this formula is that we do not want to change the weights too abruptly or quickly as it implies higher transaction costs.
Given all these considerations, we introduce the (simplified) index tracking problem (allowing for transaction costs) as follows:
β * ( t ) : = arg min β ( t ) | | r I ( t ) R ( t ) β ( t ) | | 2 2 + λ | | β ( t ) β ( t ) | | 2 2 s . t . 0 β ( t ) 1 β ( t ) T 1 = 1 , | | β ( t ) | | 0 = p .
Similarly to the approach in Satpathy and Shah (2022), we introduce distribution on the space of permutations as a means to solve this problem
T ( P ( t ) ) = exp | | r I ( t ) R P ( t ) ( t ) β P ( t ) ( t ) | | 2 2 λ | | β P ( t ) ( t ) β P ( t ) ( t ) | 2 2 η .
Then, applying similar reasoning as in the case of the optimisation problem without transaction costs, in order to find the optimal permutation at time t, T ( P ( t ) ) is approximated via SMC methods using density tempering.

4. Methodology and Data

The index-tracking algorithm utilised in this study builds upon the framework introduced in Section 2, with further details provided in Section 3. Here, we delve into the methodological framework used for implementing index-tracking strategies, emphasising their reliance solely on historical data. To clarify, at each rebalancing point (t), the composition of the tracking portfolio is determined based solely on market data available up to time t, within a predefined lookback period of 30 trading days. This approach ensures that no future data beyond time t influence our backtesting methodology. Furthermore, we have defined the temporal scope to encompass 20 portfolio rebalancing instances for each strategy under investigation. The rebalancing frequency is set to every 60 observations, corresponding to a span of 60 trading days.

4.1. Model Cases

The model cases examined in this analysis primarily differ in their hyperparameter settings. To provide a concise overview of the diverse models formulated, Table 1 summarises the hyperparameters employed for each case.
In order to assess the performance of the updated model that incorporates transaction costs compared to the initial model, we maintain the hyperparameters of both models at the values specified in Table 1. When λ = 0 , this corresponds to the initial model without accounting for transaction costs.
Following an initial examination of the efficacy of the R 2 proposal distribution in identifying accurate stock selections, it became evident that its capacity to accurately pinpoint assets relevant to the experimental framework was inadequate. Consequently, an alternative proposal distribution rooted in multiple linear regression was employed. Specifically, this approach entails evaluating all assets within the designated universe at a given time point and computing the coefficients of multiple linear regression for these stock prices relative to index prices. Subsequently, the absolute values of these regression coefficients are obtained, normalised to ensure their cumulative sum equals one, thereby constituting the basis of our prior/proposal distribution. An exploratory analysis unveiled the superiority of this proposal distribution over the R 2 variant originally posited by Satpathy and Shah (2022).
The enhanced efficacy of the revised proposal distribution, coupled with considerations regarding numerical efficiency, culminated in the decision to forgo the Metropolis–Hastings support boosting step. The proposal distribution gives a good sampling of the state space and is close enough to the posterior distribution to make additional steps such as support boosting redundant or even counterproductive. Due to this, we have chosen to not use the support boosting step in our version of the particle filter algorithm.

4.2. Constraints and Optimisation

In our numerical experiments, we implemented the index-tracking algorithm with the inclusion of no-short-selling constraints and full-budget constraints. It is worth noting the implications of introducing these constraints. When incorporating no-short-selling and full-budget constraints, an explicit solution of the β parameters cannot be obtained, necessitating the use of numerical optimisation techniques. This introduces a critical consideration, as the numerical optimisation algorithm needs to be run multiple times during the execution of the index-tracking algorithm, particularly when the Effective Sample Size (ESS) falls below the threshold, requiring Metropolis–Hastings steps. Consequently, this can significantly increase the overall runtime and potentially limit the algorithm’s performance and efficiency.
Fortunately, the problem at hand benefits from the presence of a closed-form expression for the Jacobian, which provides valuable gradient information. The closed-form expression for the Jacobian, obtained by simply differentiating Equation (19), is given by
( R T R ) β ( t ) R T r I + 2 λ ( β ( t ) β ( t ) ) .
Numerical optimisers that utilise this gradient information can leverage it to enhance the convergence speed and improve overall efficiency. In our study, we employed the Sequential Least Squares Quadratic Programming (SLSQP) algorithm while incorporating the provided Jacobian. This approach allowed us to reduce the computation time by a substantial factor of approximately 7 (compared to algorithms not using gradient information).

4.3. Data

4.3.1. Simulated Data

In the process of algorithm development, we employed simulated index and stock data where the true parameters were known. To generate the stock price paths for N stocks up to some future time T > 0 , we employed the Standard Geometric Brownian Motion (GBM) model. Each stock price path follows the process
d S i ( t ) = μ i S i ( t ) d t + σ i S i ( t ) d W i ( t )
where μ i and σ i denote the drift and volatility of the i-th stock, respectively, and { W i ( t ) : t 0 } denotes a standard Brownian Motion process. The Brownian Motion processes are correlated such that E ( d W i ( t ) · d W j ( t ) ) = ρ i , j · d t and ρ i , i = 1 . This leads to an explicit form to determine the stock price paths for each stock i = 1 , , n using
S i ( t k ) = S i ( 0 ) exp j = 1 k ( μ i 1 2 σ i 2 ) Δ t j + Δ t j σ i X i
with time points 0 = t 0 , t 1 , , t n = T and time increments Δ t j = t j t j 1 . X i corresponds to the i-th element in a multivariate normal vector X N n ( 0 , ρ ) , where ρ = C orr ( X , X ) denotes the correlation matrix. The multivariate normal random vectors can be calculated as X = LZ where Z N n ( 0 , I ) and L denotes the Cholesky decomposition of the correlation matrix ρ .
To account for the dynamic nature of index constituents, we established a fixed maximum number of stocks (n) in the asset universe and a maximum number of stocks that could compose the index ( n I n ). From the asset universe S = { 1 , , n } , a random selection of stocks was made to form the index, along with random time points indicating when these stocks entered or exited the index. This gives the set of stocks that the index comprises at each time point as I t S . This simulation setup mirrors the dynamic nature of real-world indices, where stocks are added or removed based on specific criteria and timing.
In order to construct the simulated index time series, we take the weighted average of stock returns in the index at each time point, i.e.,
r I ( t ) = i I t β i ( t ) R i ( t )
where at each time point, we have that
i I t β i ( t ) = 1
The simulated data spanned a predetermined temporal window, with a daily frequency. The parameters for drift, volatility, and initial stock prices for the constituent stocks in the index were uniformly sampled within appropriate ranges. The rationale behind uniformly sampling GBM parameters was primarily to introduce random variation among the simulated stocks. By incorporating this variation, we aimed to capture the inherent diversity and unpredictability observed in real-world stock markets. This approach ensured that the simulated data reflected a realistic scenario, allowing us to evaluate the algorithm’s performance under more representative conditions.
To evaluate the accuracy and effectiveness of the index-tracking algorithm, we calculated index returns as the average of the individual stock returns that constituted it (i.e., β i ( t ) = 1 | I t | ). By aggregating these returns, we obtained the index value, which served as a benchmark for evaluating the algorithm’s ability to accurately select the constituent stocks and their weights in the tracking portfolio.
Throughout the simulation experiments, the index-tracking algorithm consistently demonstrated the capability to effectively identify the stocks forming the index and their constituent weights. This successful performance suggests that the algorithm holds promise in real-world scenarios, where the true underlying stock dynamics are not known. The key test, however, is the algorithm’s performance on such real-world data, which will be considered next.

4.3.2. Empirical Data

Daily data for various indices were obtained from Bloomberg for an empirical study. This resulted in a diverse set of empirical datasets, encompassing the number of constituent stocks and economy sectors. Additionally, data points were excluded from the analysis when the timelines for the indices did not align with those of the individual stocks. The results presented in this paper are based on the indices from Table 2 and the period spanning 18 October 2018 to 5 July 2023. The presented number of stocks in each index is approximate as stocks enter and exit the indices from time to time.

4.4. Metrics

In this section, we outline the metrics of interest. These metrics fall into two categories: “Performance metrics,” which directly impact the decision-making process for selecting investment strategies and parameters, and “Research metrics,” which we explore in order to understand the nature of the algorithms more deeply.

4.4.1. Tracking Error

The primary performance measure we utilise is the tracking error, as defined in Satpathy and Shah (2022). For actual index returns r I ( t ) and tracking portfolio returns β ( t ) R ( t ) at times t = 1 , , T , the tracking error can be calculated as
TE = 1 T 1 t = 1 T ( r I ( t ) β ( t ) R ( t ) ) 2
The tracking error quantifies the level of deviation between the returns of the index tracking portfolio and the target index. It can be calculated as the standard deviation of the return residuals, which are the differences between the returns of the index and the returns of the index tracking portfolio.

4.4.2. Total Transaction Costs

The second performance metric we looked at was the total transaction costs. In the context of transaction cost-conscious investing, monitoring the evolution of the wealth process over time becomes crucial. If transaction costs are defined as in Section 3, then the total transaction costs can be calculated as
TC = i = 1 n ( 1 C ( β ( t ) ) ) X ( t )
Recall that t is the moment just before the rebalancing at t. Tracking the total transaction costs allows us to assess the overall cost-effectiveness of the tracking strategy. These costs are calculated at the end of each rebalancing period, representing the expenses incurred during the adjustment of the portfolio’s composition.
The interpretation of the results regarding transaction costs should take into account the specific investment objectives and constraints of the investor. Lower total transaction costs are generally desired, as they indicate a more efficient allocation of resources. However, it is important to balance these cost considerations with the tracking error. Intuitively, there is often a trade-off between minimising the tracking error and minimising transaction costs. Achieving a lower tracking error may require more frequent portfolio rebalancing, leading to higher transaction costs. Conversely, reducing transaction costs may involve a less frequent rebalancing strategy, potentially resulting in a higher tracking error.

4.4.3. Wealth Error

The assessment of tracking error conventionally involves the examination of the standard deviation of error terms derived from the difference between the simple returns of an index and those of a tracking portfolio. While this method provides valuable insights into tracking performance, we sought to devise a composite metric that not only evaluates tracking accuracy but also integrates considerations of cost efficiency inherent to a tracking strategy. To this end, we have opted to employ a metric that combines both aspects by computing the average absolute differences between the cumulative returns of the portfolio wealth and the benchmark returns. Formally, our proposed metric, denoted as the Wealth Error (WE), is defined as:
WE = 1 T t = 1 T u = 1 t ( 1 + r I ( u ) ) X ( t ) X ( 1 )
where T represents the total number of observations, r I ( u ) denotes the return of the index at time u, and X ( t ) represents the portfolio wealth at time t. Furthermore, u = 1 t ( 1 + r I ( u ) ) represents the cumulative return of the index, whereas X ( t ) X ( 1 ) represents the cumulative return of the portfolio. This formulation encapsulates both the magnitude and direction of disparities between the portfolio’s cumulative returns and those of the benchmark index, thereby providing a comprehensive measure of tracking performance while considering the cost efficiency of the tracking strategy.
By evaluating both the tracking error and the total transaction costs, we can comprehensively assess the performance of the index tracking strategy. It allows us to understand the trade-offs between achieving accurate index replication and managing transaction costs, enabling investors to make informed decisions based on their specific investment goals and preferences.
The following metrics relates to the “Research metrics” category.

4.4.4. Other Cost Metrics

To better understand the interplay and distribution of costs in a tracking strategy, we additionally investigate the minimum, mean and maximum transaction costs associated with a rebalancing over the investment horizon.

4.4.5. Retention

This metric indicates the proportion of stocks retained in the tracking portfolio after each rebalancing. For instance, if initially there are 20 stocks in the portfolio and after rebalancing 12 of them are replaced with others from the index, the retention is 40%.
We aggregate the values of Min, Mean, and Max Retention over all the rebalancing time points in the experiment.

4.4.6. Maximum Weight

Lastly, in order to assess the effect that our objective function has on portfolio concentration, we will look at the maximum weight metric, which is defined as the overall maximum weight achieved by a single asset in the tracking portfolio, over the entire investment horizon.

5. Results and Discussion

In this section, we present an analysis of the experimental results. Our primary objective is to investigate the various patterns exhibited by the metrics concerning changes in the number of particles, denoted as N, and the cost aversion parameter, represented by λ . The stochastic nature of the optimisation algorithm means that there is some random variation in the outcomes, which we also need to quantify. Thus, we consider each of these metrics for different choices of N and λ , where for each ( N , λ ) pair, we run 40 (independent) replications of the algorithm. The averages for each metric and ( N , λ ) choice are reported in Table 3, Table 4, Table 5, Table 6 and Table 7, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).

5.1. Tracking Error

Ex ante, one would intuitively expect there to be a trade–off between transaction costs and tracking error, as including transaction costs in the optimisation objective reduces the amount of rebalancing in the tracking portfolio, therefore leading one to expect tracking performance to be worse. However, the tracking errors reported in Table 3, Table 4 and Table 5 do not bear this out. The performance of the tracking strategy, in terms of tracking error, remains largely unaffected by the size of the cost aversion parameter λ . Plotting the tracking portfolio performance against the target indices (see Appendix A.1 for graphs of the 40 tracking portfolio runs against the target index in each scenario) illustrates that the introduction of cost aversion has minimal impact on tracking performance and in some cases even improves this performance. Consequently, it is possible to reduce total transaction costs without significantly affecting tracking performance.
In the tracking plots in Appendix A.1, one does observe a few instances where increased cost aversion leads to poor tracking performance in some portfolio wealth trajectories. This is a consequence of the random element in the optimisation algorithm, it can be avoided by conducting multiple runs of the algorithm and eliminating “outlier” results. However, an interesting observation is that for some indices, when cost aversion increases and tracking performance deteriorates, the tracking portfolio actually tends to outperform the index: This is evident in Figure A2, Figure A6, Figure A7, Figure A8, Figure A10, Figure A11 and Figure A14.
While it might seem counterintuitive that a larger transaction cost weighting λ does not imply a larger tracking error, this may be because the future return behaviour of real-world data may differ from past behaviour, and in particular, correlations can change. Thus, the optimal tracking portfolio based on past returns is not necessarily optimal going forward. As including transaction costs tends to stabilise the tracking portfolio weights, this may result in better tracking performance, compared to ignoring transaction costs and rebalancing fully to reflect any change in return correlations in the lookback estimation window, which may be partly due to random noise.
Similarly, one might expect that increasing the number of particles would improve tracking error, given that a larger particle set should enable the particle filter algorithm to achieve more accurate index tracking. However, the results presented in Table 3, Table 4 and Table 5 show that the increase in particle numbers has minimal impact on tracking error. In some instances, employing only ten particles results in better tracking performance compared to a higher number of particles.
From a computational efficiency perspective, using fewer particles significantly reduces the execution time of the particle filter algorithm. This highlights the potential for optimising particle numbers to balance between computational demands and tracking performance.

5.2. Total Cost

The investigation into total cost (Figure 1) shows that increasing the weighting λ of transaction costs in the optimisation has the desired effect of reducing the cumulative transaction costs. Unlike when considering tracking performance, the number of particles used in the algorithm does have a substantial impact here: The variability of cumulative transaction costs across the 40 runs of the algorithm is reduced as the number of particles increases, and cumulative transaction costs are substantially reduced when this is part of the optimisation objective (i.e., for λ = 1000 and λ = 10,000).

5.3. Retention

Table 6 and Table 7 offer insights into the behaviour of portfolio retention under varying levels of λ , representing transaction cost sensitivity. When λ is set to 0, indicating no consideration for transaction costs, there are notable instances where the entire portfolio experiences complete turnover during rebalancing.
However, as λ increases, reaching values as high as 10,000, the frequency of rebalancing events resulting in no stock replacements becomes more pronounced. This trend suggests a shift towards a conservative approach, where transaction costs heavily influence decision-making processes, leading to fewer adjustments in the portfolio composition.
Further examination of the mean retention reveals additional nuances. On average across all indices, around 25% of stocks persist in the portfolio post-rebalancing when λ = 0 . This relatively low retention rate underscores the high level of churn as a consequence of ignoring transaction costs.
In contrast, with increasing λ sets, the average retention rate increases, signifying a more stable portfolio composition. This implies that higher λ values effectively curb the turnover of stocks, promoting greater stability.
For an exemplary run of the algorithm on DAX data, this is visualised in Figure 2 (additional examples are given in Appendix A.3). This chart vividly showcases the dynamic evolution of the tracking portfolio’s composition, illustrating how larger values of λ lead to a more stable portfolio, whereas smaller values of λ permit frequent changes in the portfolio composition at each rebalancing step.
Overall, the observed patterns suggest that the adjustment of λ plays a crucial role in determining the level of portfolio turnover, and thus including transaction costs in the index tracking algorithm acts as a deterrent to excessive stock hopping and fosters a more stable investment strategy.
This is also reflected by the maximum weight that the algorithm assigns to any single asset over the life of the tracking strategy, as reported in Table 7: Increasing the aversion to transaction costs reduces this maximum weight, while ignoring transaction costs allows the algorithm to “hop” frequently between large allocations to individual assets in an inherently unstable manner, exposing the tracking portfolio to potentially large idiosyncratic risk.
Unlike tracking performance (but like transaction costs), the number of particles used in the algorithm does impact the asset retention results: The algorithm achieves the objective of minimising transaction costs better when a larger number of particles is used, and this is reflected by higher rates of asset retention. This can also be seen in the distribution of mean asset retention rates across the 40 runs in each ( N , λ ) scenario for each index, as represented by the histograms in Appendix A.4: Increasing the number of particles tends to shift the distribution of asset retention to the right. This shift suggests that a higher number of particles results in portfolios where dominant assets are retained for longer periods over the tracking strategy’s time horizon. This observation aligns with the expectation that more particles enhance the exploration of possible subsets of assets used to track the index, thereby allowing the algorithm to find better solutions in terms of the objective of minimising transaction costs.

6. Conclusions

Overall, the incorporation of the cost aversion parameter into the algorithm serves as a mechanism for managing transaction costs by enhancing asset retention across rebalancing periods. This mitigates excessive rebalancing due to estimation noise. Increased asset retention leads to portfolio strategies characterised by not only greater stability but also more appropriate diversification. Notably, including transaction costs in the tracking objective has little impact on tracking performance, suggesting that the benefit of cost control and portfolio stability can be achieved without any substantial performance trade–off.
However, it is important to recognise that the nature and intrinsic properties of the index being analysed, such as its size and the maturity of the market that it represents, can also influence the outcomes of the strategy. For instance, an analysis of the S&P 500 index reveals that the strategy consistently performs better compared to other indices, especially those representing emerging markets. This indicates that the algorithm may be more effective in more liquid markets.
These conclusions emphasise the complex interplay between tracking error minimisation and the incorporation of transaction costs into the objective function. The study reveals that there are additional factors at play beyond the straightforward trade-off between minimising tracking error and integrating transaction costs. This highlights the need for a deeper understanding of the underlying mechanisms that govern the relationship between these variables.

Opportunities for Further Study

Future research to be explored could include the following:
  • Considering the market impact, as it can also lead to additional losses when working with large volumes. Hence, investors should come up with smart order routing algorithms to optimise trade execution and reduce market impact. These algorithms can help split large orders into smaller ones and execute them in a way that minimises the price impact.
  • Considering liquidity screening to avoid investing in assets that may have high transaction costs or are illiquid, as these can significantly impact the performance of your tracking strategy.
  • Considering short sales so one can also track “short” indices, as well as considering the short position, i.e., negative β ’s, which will relax the constraints in the proposed method.
  • Considering other types of transaction costs: tiered commissions, fixed-plus-percentage commissions, flat-fee commissions.

Author Contributions

Conceptualisation, E.S.; Formal analysis, T.M.O., L.H.-R. and D.S.; Investigation, T.M.O., L.H.-R. and D.S.; Methodology, T.M.O., L.H.-R., D.S. and E.S.; Software, T.M.O., L.H.-R. and D.S.; Supervision, E.S.; Visualisation, T.M.O., L.H.-R. and D.S.; Writing—original draft, T.M.O., L.H.-R. and D.S.; Writing—review and editing, T.M.O., L.H.-R., D.S. and E.S. Authors are listed in alphabetical order in the article citation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Market data (history of closing prices for indices and component stocks) that support the findings of this study are available from Bloomberg. Commercial restrictions apply to the availability of this data, which were used under license for this study.

Acknowledgments

The initial research for this paper was conducted as part of the Financial Mathematics Team Challenge (FMTC) at AIFMRM, University of Cape Town. The authors thank Francesca Primavera for her contributions to this work as Team Leader at the FMTC and for subsequent helpful discussions about the formulation of the index tracking problem under transaction costs in Section 2.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

This part provides additional plots obtained during the experiments.

Appendix A.1. Tracking Performance

In this subsection, you will find plots that illustrate the returns of the constructed portfolio compared to the real index returns across various values of λ -s and numbers of particles. The red line indicates the index returns and the blue indicates the tracking portfolio’s returns, with one blue line for each of the 40 independent runs of the algorithm.
Figure A1. DAX Tracking portfolio return (blue) versus actual index returns (red).
Figure A1. DAX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a1
Figure A2. HSI Tracking portfolio return (blue) versus actual index returns (red).
Figure A2. HSI Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a2
Figure A3. JALSH Tracking portfolio return (blue) versus actual index returns (red).
Figure A3. JALSH Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a3
Figure A4. SENSEX Tracking portfolio return (blue) versus actual index returns (red).
Figure A4. SENSEX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a4
Figure A5. SMI Tracking portfolio return (blue) versus actual index returns (red).
Figure A5. SMI Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a5
Figure A6. SPTSX Tracking portfolio return (blue) versus actual index returns (red).
Figure A6. SPTSX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a6
Figure A7. SPX Tracking portfolio return (blue) versus actual index returns (red).
Figure A7. SPX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a7
Figure A8. UKX Tracking portfolio return (blue) versus actual index returns (red).
Figure A8. UKX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a8
Figure A9. CAC Tracking portfolio return (blue) versus actual index returns (red).
Figure A9. CAC Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a9
Figure A10. IBEX Tracking portfolio return (blue) versus actual index returns (red).
Figure A10. IBEX Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a10
Figure A11. IBOV Tracking portfolio return (blue) versus actual index returns (red).
Figure A11. IBOV Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a11
Figure A12. IPSA Tracking portfolio return (blue) versus actual index returns (red).
Figure A12. IPSA Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a12
Figure A13. NKY Tracking portfolio return (blue) versus actual index returns (red).
Figure A13. NKY Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a13
Figure A14. SHSZ300 Tracking portfolio return (blue) versus actual index returns (red).
Figure A14. SHSZ300 Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a14
Figure A15. TWSE Tracking portfolio return (blue) versus actual index returns (red).
Figure A15. TWSE Tracking portfolio return (blue) versus actual index returns (red).
Risks 12 00155 g0a15

Appendix A.2. Transaction Costs

The plots below illustrate the cumulative transaction costs resulting from running the algorithm various numbers of particles.
Figure A16. Cumulative transaction costs paid. Here, green is for λ = 0 , orange is λ = 1000 and red is λ = 10,000.
Figure A16. Cumulative transaction costs paid. Here, green is for λ = 0 , orange is λ = 1000 and red is λ = 10,000.
Risks 12 00155 g0a16

Appendix A.3. Composition Plots

These plots illustrate how the composition of the tracking portfolio was changing during the experiments for various values of λ .
Figure A17. Changes in the composition of the tracking portfolio for the HSI index.
Figure A17. Changes in the composition of the tracking portfolio for the HSI index.
Risks 12 00155 g0a17
Figure A18. Changes in the composition of the tracking portfolio for the JALSH index.
Figure A18. Changes in the composition of the tracking portfolio for the JALSH index.
Risks 12 00155 g0a18
Figure A19. Changes in the composition of the tracking portfolio for the SENSEX index.
Figure A19. Changes in the composition of the tracking portfolio for the SENSEX index.
Risks 12 00155 g0a19
Figure A20. Changes in the composition of the tracking portfolio for the SMI index.
Figure A20. Changes in the composition of the tracking portfolio for the SMI index.
Risks 12 00155 g0a20
Figure A21. Changes in the composition of the tracking portfolio for the SPTSX index.
Figure A21. Changes in the composition of the tracking portfolio for the SPTSX index.
Risks 12 00155 g0a21
Figure A22. Changes in the composition of the tracking portfolio for the SPX index.
Figure A22. Changes in the composition of the tracking portfolio for the SPX index.
Risks 12 00155 g0a22
Figure A23. Changes in the composition of the tracking portfolio for the UKX index.
Figure A23. Changes in the composition of the tracking portfolio for the UKX index.
Risks 12 00155 g0a23

Appendix A.4. Mean Assets Retention

The plots below show the distribution (across 40 independent runs of the algorithm) of the average (over the life of the strategy) proportion of the portfolio that was maintained between rebalancing periods for various values of the number of particles N and the cost aversion coefficient λ . The red line indicates the mean of the distribution.
Figure A24. CAC Mean Assets retention.
Figure A24. CAC Mean Assets retention.
Risks 12 00155 g0a24
Figure A25. HSI Mean Assets retention.
Figure A25. HSI Mean Assets retention.
Risks 12 00155 g0a25
Figure A26. JALSH Mean Assets retention.
Figure A26. JALSH Mean Assets retention.
Risks 12 00155 g0a26
Figure A27. SENSEX Mean Assets retention.
Figure A27. SENSEX Mean Assets retention.
Risks 12 00155 g0a27
Figure A28. SMI Mean Assets retention.
Figure A28. SMI Mean Assets retention.
Risks 12 00155 g0a28
Figure A29. SPTSX Mean Assets retention.
Figure A29. SPTSX Mean Assets retention.
Risks 12 00155 g0a29
Figure A30. SPX Mean Assets retention.
Figure A30. SPX Mean Assets retention.
Risks 12 00155 g0a30
Figure A31. UKX Mean Assets retention.
Figure A31. UKX Mean Assets retention.
Risks 12 00155 g0a31
Figure A32. DAX Mean Assets retention.
Figure A32. DAX Mean Assets retention.
Risks 12 00155 g0a32
Figure A33. IBEX Mean Assets retention.
Figure A33. IBEX Mean Assets retention.
Risks 12 00155 g0a33
Figure A34. IBOV Mean Assets retention.
Figure A34. IBOV Mean Assets retention.
Risks 12 00155 g0a34
Figure A35. IPSA Mean Assets retention.
Figure A35. IPSA Mean Assets retention.
Risks 12 00155 g0a35
Figure A36. NKY Mean Assets retention.
Figure A36. NKY Mean Assets retention.
Risks 12 00155 g0a36
Figure A37. SHSZ300 Mean Assets retention.
Figure A37. SHSZ300 Mean Assets retention.
Risks 12 00155 g0a37
Figure A38. TWSE Mean Assets retention.
Figure A38. TWSE Mean Assets retention.
Risks 12 00155 g0a38

Notes

1
In Satpathy and Shah (2022) the sparse index problem is solved without considering any constraints on the β ^ P . In such a case the optimal solution admits an explicit form.
2
Going forward particles, samples and permutations will be used interchangeably, i.e., each particle represents a particular subset of p of the n stocks comprising the index.
3
Results in the literature suggest that multinomial resampling is not a good choice in this context. For example, Hol et al. (2006) note that “theoretical analysis shows that the resampling quality can be improved by using a different algorithm than multinomial resampling.” In particular, residual and stratified resampling have lower variances, and Hol et al. (2006) advocate systematic resampling over stratified resampling as it has the lowest discrepancy.
4
The proposal in MH algorithm has no relation to the proposal in the Importance Sampling step.

References

  1. Barber, Brad, and Terrance Odean. 2000. Trading is hazardous to your wealth: The common stock investment performance of individual investors. Journal of Finance 55: 773–806. [Google Scholar] [CrossRef]
  2. Benidis, Konstantinos, Yiyong Feng, and Daniel P. Palomar. 2018. Optimization methods for financial index tracking: From theory to practice. Foundations and Trends® in Optimization 3: 171–279. [Google Scholar] [CrossRef]
  3. Chen, Zhe. 2003. Bayesian filtering: From Kalman filters to particle filters, and beyond. Statistics 182: 1–69. [Google Scholar] [CrossRef]
  4. Del Moral, Pierre, Arnaud Doucet, and Ajay Jasra. 2006. Sequential Monte Carlo Samplers. Journal of the Royal Statistical Society Series B: Statistical Methodology 68: 411–36. [Google Scholar] [CrossRef]
  5. Duan, Jin-Chuan. 2019. Variable selection with big data based on zero norm and via sequential Monte Carlo. SSRN. [Google Scholar] [CrossRef]
  6. Gellert, Karol, and Erik Schlögl. 2021. Parameter Learning and Change Detection Using a Particle Filter with Accelerated Adaptation. Risks 9: 228. [Google Scholar] [CrossRef]
  7. Hol, Jeroen D., Thomas B. Schon, and Fredrik Gustafsson. 2006. On resampling algorithms for particle filters. Paper presented at the 2006 IEEE Nonlinear Statistical Signal Processing Workshop, Cambridge, UK, September 13–15; pp. 79–82. [Google Scholar]
  8. Rompotis, Gerasimos G. 2013. Actively vs. passively managed exchange traded funds. Aestimatio 6: 116–35. [Google Scholar]
  9. Ruf, Johannes, and Kangjianan Xie. 2019. The impact of proportional transaction costs on systematically generated portfolios. SIAM Journal on Financial Mathematics 11: 881–96. [Google Scholar] [CrossRef]
  10. Satpathy, Tanmay, and Rushabh Shah. 2022. Sparse index tracking using sequential Monte Carlo. Quantitative Finance 22: 1579–92. [Google Scholar] [CrossRef]
  11. Speekenbrink, Maarten. 2016. A tutorial on particle filters. Journal of Mathematical Psychology 73: 140–52. [Google Scholar] [CrossRef]
  12. Tibshirani, Robert. 1996. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society Series B: Statistical Methodology 58: 267–88. [Google Scholar] [CrossRef]
  13. Xu, Kelvin, Jimmy Ba, Ryan Kiros, Kyunghyun Cho, Aaron Courville, Ruslan Salakhudinov, Rich Zemel, and Yoshua Bengio. 2015. Show, attend and tell: Neural image caption generation with visual attention. In International Conference on Machine Learning. Birmingham: PMLR, pp. 2048–57. [Google Scholar]
Figure 1. Total cost for each index and all cases of λ . Here, the green line represents λ = 0 , the orange line represents λ = 1000 and the red line represents λ = 10,000. Each line corresponds to one independent run of the algorithm.
Figure 1. Total cost for each index and all cases of λ . Here, the green line represents λ = 0 , the orange line represents λ = 1000 and the red line represents λ = 10,000. Each line corresponds to one independent run of the algorithm.
Risks 12 00155 g001
Figure 2. Changes in the composition of the tracking portfolio for the DAX index.
Figure 2. Changes in the composition of the tracking portfolio for the DAX index.
Risks 12 00155 g002
Table 1. Hyperparameters chosen for implementation and comparison of model cases on real-world data.
Table 1. Hyperparameters chosen for implementation and comparison of model cases on real-world data.
Hyper-ParametersValues
Step Size ( δ )0.2
Lookback Window Size (T)30
ESS Threshold N 2
Number of Particles (n) { 10 , 100 , 1000 }
Variance Explained95%
Rebalance Frequency (Observations)60
Transaction Cost per Unit Nominal ( ϵ )0.001
Transaction Cost Penalty ( λ ) { 0 , 1000 , 10,000 }
Table 2. Equity indices used to test the performance of the algorithm, ordered by the number of stocks.
Table 2. Equity indices used to test the performance of the algorithm, ordered by the number of stocks.
IndexLocationNumber of StocksSectors
SMISwitzerland20Financials, industrials, pharmaceuticals, and technology
MervalArgentina20Financials, energy, consumer staples, and industrials.
IBEXSpain35Financials, telecommunications, utilities, and energy.
CACFrance40Financials, industrials, consumer goods, and healthcare.
HSIHong Kong50Financials, real estate, technology, and consumer goods.
DAXGermany30Automotive, financial services, industrials, and technology.
IPSAChile30Utilities, financials, industrials, and consumer goods.
SENSEXIndia30Information technology, financial services, energy, and consumer goods.
IBOVBrazil73Financials, energy, consumer discretionary, and materials.
UKXUnited Kingdom100Financials, energy, consumer goods, and pharmaceuticals.
NKYJapan225Industrials, technology, consumer goods, and financials.
JALSHSouth Africa164Financials, mining, consumer goods, and telecommunications.
SHSZ300China300Financials, industrials, technology, and healthcare.
SPTSXCanada230Financials, energy, materials, and industrials.
SPXUnited States500Financials, technology, healthcare, and consumer discretionary.
TWSETaiwan900+Technology, manufacturing, financials, and telecommunications.
Table 3. Tracking error, Wealth Error, and transaction costs observed for 10 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Table 3. Tracking error, Wealth Error, and transaction costs observed for 10 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Mean Abs Wealth ErrorTotal CostTracking Error
N = 10 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.111 ± 0.0180.128 ± 0.0250.145 ± 0.02437,461.2 ± 1300.021,823.1 ± 1300.014,280.5 ± 1600.03.52 × 10 5 ± 2.2 × 10 6 4.93 × 10 5 ± 4.6 × 10 6 6.78 × 10 5 ± 8.3 × 10 6
DAX0.15 ± 0.020.149 ± 0.0290.144 ± 0.02237,367.5 ± 1200.023,546.9 ± 1200.019,775.7 ± 1600.03.89 × 10 5 ± 1.6 × 10 6 4.31 × 10 5 ± 3.1 × 10 6 5.33 × 10 5 ± 5.1 × 10 6
HSI0.162 ± 0.0250.123 ± 0.0310.113 ± 0.02237,738.1 ± 1100.026,915.6 ± 1100.022,143.1 ± 1200.08.11 × 10 5 ± 5.4 × 10 6 4.51 × 10 5 ± 2.5 × 10 6 4.98 × 10 5 ± 4.5 × 10 6
IBEX0.087 ± 0.0170.128 ± 0.0290.15 ± 0.04432,016.2 ± 1100.015,890.0 ± 1100.010,220.6 ± 1100.05.74 × 10 5 ± 3.6 × 10 6 7.19 × 10 5 ± 6.8 × 10 6 8.30 × 10 5 ± 1.1 × 10 5
IBOV0.172 ± 0.0330.183 ± 0.0310.186 ± 0.03543,602.5 ± 1900.033,271.0 ± 1400.029,226.9 ± 1700.05.53 × 10 5 ± 4.1 × 10 6 4.60 × 10 5 ± 4.5 × 10 6 5.34 × 10 5 ± 5.9 × 10 6
IPSA0.174 ± 0.0190.137 ± 0.0180.183 ± 0.04725,533.6 ± 750.017,580.6 ± 1300.015,129.8 ± 1400.09.84 × 10 5 ± 4.7 × 10 6 8.05 × 10 5 ± 5.3 × 10 6 1.26 × 10 4 ± 3.1 × 10 5
JALSH0.221 ± 0.0280.179 ± 0.0250.175 ± 0.02734,301.3 ± 1000.026,258.0 ± 960.025,863.5 ± 1200.04.82 × 10 5 ± 1.8 × 10 6 4.56 × 10 5 ± 1.9 × 10 6 4.76 × 10 5 ± 2.7 × 10 6
MERVAL0.77 ± 0.130.84 ± 0.190.84 ± 0.2179,791.8 ± 8500.046,439.5 ± 6000.028,360.7 ± 4300.02.03 × 10 4 ± 1.4 × 10 5 1.66 × 10 4 ± 1.2 × 10 5 1.79 × 10 4 ± 1.6 × 10 5
NKY0.125 ± 0.0150.102 ± 0.0130.085 ± 0.01333,800.9 ± 520.028,365.3 ± 530.028,252.4 ± 530.01.8 × 10 5 ± 1 × 10 6 1.46 × 10 5 ± 6.5 × 10 7 1.42 × 10 5 ± 6.4 × 10 7
SENSEX0.189 ± 0.0530.187 ± 0.0380.166 ± 0.03648,238.7 ± 2300.020,878.6 ± 1800.013,487.9 ± 1700.04.31 × 10 5 ± 2.4 × 10 6 5.01 × 10 5 ± 2.9 × 10 6 5.83 × 10 5 ± 4.1 × 10 6
SHSZ3000.256 ± 0.0450.3 ± 0.0350.399 ± 0.04457,871.9 ± 2200.049,309.1 ± 1200.052,331.2 ± 1300.08.15 × 10 5 ± 5.4 × 10 6 3.9 × 10 5 ± 2 × 10 6 3.13 × 10 5 ± 1.5 × 10 6
SMI0.134 ± 0.0240.08 ± 0.0150.074 ± 0.01337,641.7 ± 1200.08376.9 ± 990.04133.0 ± 740.02.41 × 10 5 ± 1.5 × 10 6 3.12 × 10 5 ± 2.5 × 10 6 3.23 × 10 5 ± 2.6 × 10 6
SPTSX0.226 ± 0.0390.17 ± 0.020.16 ± 0.02348,127.1 ± 1600.039,464.1 ± 790.037,829.3 ± 820.03.52 × 10 5 ± 2.4 × 10 6 2.31 × 10 5 ± 1.3 × 10 6 1.92 × 10 5 ± 1.4 × 10 6
SPX0.09 ± 0.0140.075 ± 0.0120.082 ± 0.01245,785.0 ± 1000.040,077.8 ± 670.038,728.4 ± 810.02.82 × 10 5 ± 3.5 × 10 6 1.85 × 10 5 ± 9.6 × 10 7 1.42 × 10 5 ± 8.3 × 10 7
TWSE0.18 ± 0.020.152 ± 0.0220.135 ± 0.01643,858.6 ± 900.040,046.5 ± 1300.038,459.0 ± 840.05.72 × 10 5 ± 1.8 × 10 6 5.88 × 10 5 ± 1.9 × 10 6 6.00 × 10 5 ± 1.4 × 10 6
UKX0.159 ± 0.0310.119 ± 0.0230.11 ± 0.0238,613.5 ± 1100.026,332.3 ± 1300.024,699.3 ± 1200.03.19 × 10 5 ± 1.9 × 10 6 2.75 × 10 5 ± 2.2 × 10 6 2.62 × 10 5 ± 2.1 × 10 6
Table 4. Tracking error, Wealth Error and transaction costs observed for 100 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Table 4. Tracking error, Wealth Error and transaction costs observed for 100 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Mean Abs Wealth ErrorTotal CostTracking Error
N = 100 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.119 ± 0.0260.213 ± 0.0340.166 ± 0.03138,713.0 ± 1400.011,502.7 ± 1500.08424.6 ± 1300.04.08 × 10 5 ± 3.8 × 10 6 6.08 × 10 5 ± 5.8 × 10 6 5.82 × 10 5 ± 7.5 × 10 6
DAX0.145 ± 0.0240.143 ± 0.0210.149 ± 0.02338,583.7 ± 1500.011,689.6 ± 750.09031.4 ± 710.04.64 × 10 5 ± 1.8 × 10 6 6.90 × 10 5 ± 7.5 × 10 6 7.11 × 10 5 ± 6.5 × 10 6
HSI0.178 ± 0.0290.121 ± 0.0240.151 ± 0.03638,144.5 ± 1200.016,605.8 ± 1000.013,927.5 ± 1200.07.95 × 10 5 ± 3.2 × 10 6 7.10 × 10 5 ± 4.6 × 10 6 8.18 × 10 5 ± 1.1 × 10 5
IBEX0.103 ± 0.0190.187 ± 0.0350.148 ± 0.03233,051.8 ± 1200.09396.5 ± 910.05140.7 ± 530.06.51 × 10 5 ± 4.1 × 10 6 6.68 × 10 5 ± 6.8 × 10 6 5.98 × 10 5 ± 8.3 × 10 6
IBOV0.158 ± 0.0330.243 ± 0.0540.256 ± 0.05243,072.5 ± 1700.023,722.2 ± 1600.021,695.1 ± 1600.06.64 × 10 5 ± 5.8 × 10 6 6.04 × 10 5 ± 5.8 × 10 6 7.73 × 10 5 ± 8.8 × 10 6
IPSA0.159 ± 0.0210.155 ± 0.0180.178 ± 0.02426,124.1 ± 1000.08361.2 ± 780.05599.9 ± 600.09.38 × 10 5 ± 2.8 × 10 6 1.11 × 10 4 ± 1.4 × 10 5 1.20 × 10 4 ± 1.5 × 10 5
JALSH0.198 ± 0.0210.173 ± 0.0290.205 ± 0.05635,285.6 ± 730.023,295.5 ± 1400.021,555.3 ± 2000.04.95 × 10 5 ± 1.3 × 10 6 5.64 × 10 5 ± 4.7 × 10 6 6.47 × 10 5 ± 6.8 × 10 6
MERVAL0.711 ± 0.0791.04 ± 0.210.84 ± 0.1578,626.3 ± 5700.033,692.0 ± 4700.017,779.8 ± 2500.03.007 × 10 4 ± 8.9 × 10 6 1.62 × 10 4 ± 1.1 × 10 5 1.61 × 10 4 ± 1.3 × 10 5
NKY0.137 ± 0.0130.095 ± 0.0130.089 ± 0.01333,699.9 ± 510.025,742.7 ± 570.024,273.2 ± 660.01.98 × 10 5 ± 7.8 × 10 7 1.651 × 10 5 ± 9 × 10 7 1.65 × 10 5 ± 1.5 × 10 6
SENSEX0.145 ± 0.0240.16 ± 0.0290.154 ± 0.02645,068.3 ± 1600.011,596.8 ± 1400.07147.5 ± 1000.04.58 × 10 5 ± 1.8 × 10 6 5.76 × 10 5 ± 4.3 × 10 6 6.20 × 10 5 ± 5.9 × 10 6
SHSZ3000.257 ± 0.0550.32 ± 0.0420.417 ± 0.05658,082.6 ± 2700.044,531.1 ± 1200.045,820.7 ± 1700.07.43 × 10 5 ± 3.6 × 10 6 4.19 × 10 5 ± 3.8 × 10 6 3.76 × 10 5 ± 4.8 × 10 6
SMI0.152 ± 0.0260.076 ± 0.0110.0712 ± 0.008438,880.2 ± 1200.03686.9 ± 370.01305.6 ± 190.02.8 × 10 5 ± 1 × 10 6 2.61 × 10 5 ± 1.7 × 10 6 2.25 × 10 5 ± 1.3 × 10 6
SPTSX0.238 ± 0.0390.199 ± 0.0380.18 ± 0.0448,773.5 ± 1600.034,953.4 ± 1200.032,363.8 ± 1300.04.32 × 10 5 ± 7.5 × 10 6 2.58 × 10 5 ± 3.1 × 10 6 2.76 × 10 5 ± 3.6 × 10 6
SPX0.094 ± 0.0140.088 ± 0.0170.092 ± 0.01845,003.7 ± 1000.035,540.5 ± 800.033,498.3 ± 730.03.12 × 10 5 ± 4.1 × 10 6 2.05 × 10 5 ± 2.1 × 10 6 1.80 × 10 5 ± 2.4 × 10 6
TWSE0.154 ± 0.0210.151 ± 0.0210.157 ± 0.01845,331.9 ± 970.036,982.7 ± 1200.035,035.3 ± 1200.05.76 × 10 5 ± 8.8 × 10 7 6.11 × 10 5 ± 2.7 × 10 6 6.32 × 10 5 ± 2.3 × 10 6
UKX0.153 ± 0.0290.134 ± 0.0260.107 ± 0.01838,755.3 ± 1200.019,185.3 ± 1000.016,099.2 ± 710.03.29 × 10 5 ± 1.9 × 10 6 3.20 × 10 5 ± 3.1 × 10 6 3.34 × 10 5 ± 2.9 × 10 6
Table 5. Tracking error, Wealth Error and transaction costs observed for 1000 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Table 5. Tracking error, Wealth Error and transaction costs observed for 1000 particles (N) and all cost aversions ( λ ) for each index observed. Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Mean Abs Wealth ErrorTotal CostTracking Error
N = 1000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.154 ± 0.0210.175 ± 0.0240.17 ± 0.0236,766.1 ± 1300.06915.6 ± 910.05066.1 ± 760.06.06 × 10 5 ± 5.8 × 10 6 5.48 × 10 5 ± 7.2 × 10 6 4.78 × 10 5 ± 5.6 × 10 6
DAX0.123 ± 0.0140.149 ± 0.0220.16 ± 0.02438,991.0 ± 1100.07031.1 ± 560.05010.0 ± 360.05.43 × 10 5 ± 2.1 × 10 6 5.95 × 10 5 ± 5.7 × 10 6 6.74 × 10 5 ± 5.1 × 10 6
HSI0.194 ± 0.0260.158 ± 0.0340.154 ± 0.03838,762.8 ± 1100.012,179.1 ± 860.09883.1 ± 830.07.83 × 10 5 ± 2.2 × 10 6 9.19 × 10 5 ± 7.7 × 10 6 1.04 × 10 4 ± 1.1 × 10 5
IBEX0.098 ± 0.0190.175 ± 0.0340.152 ± 0.02733,193.2 ± 1100.05741.9 ± 530.03089.9 ± 320.06.99 × 10 5 ± 3.9 × 10 6 5.66 × 10 5 ± 4.8 × 10 6 4.90 × 10 5 ± 7.2 × 10 6
IBOV0.165 ± 0.0370.235 ± 0.0510.29 ± 0.0643,649.0 ± 2000.017,291.7 ± 1100.015,714.1 ± 1500.07.46 × 10 5 ± 6.2 × 10 6 7.21 × 10 5 ± 7.4 × 10 6 7.99 × 10 5 ± 6.8 × 10 6
IPSA0.16 ± 0.020.167 ± 0.0150.173 ± 0.01925,937.2 ± 880.04367.5 ± 310.03349.6 ± 330.09.43 × 10 5 ± 2.3 × 10 6 1.02 × 10 4 ± 1.8 × 10 5 1.08 × 10 4 ± 1.4 × 10 5
JALSH0.202 ± 0.0290.207 ± 0.0310.3 ± 0.0735,465.8 ± 990.022,845.9 ± 1800.022,747.9 ± 2200.05.32 × 10 5 ± 1.5 × 10 6 7.66 × 10 5 ± 9.3 × 10 6 9.19 × 10 5 ± 7.6 × 10 6
MERVAL0.62 ± 0.0941.04 ± 0.231.06 ± 0.2378,527.7 ± 4800.023,606.4 ± 3500.017,217.0 ± 2400.03.024 × 10 4 ± 3.7 × 10 6 1.39 × 10 4 ± 7.9 × 10 6 1.46 × 10 4 ± 1.3 × 10 5
NKY0.147 ± 0.0130.088 ± 0.0130.083 ± 0.01333,843.6 ± 490.024,069.7 ± 680.022,391.7 ± 580.02.03 × 10 5 ± 6.4 × 10 7 1.98 × 10 5 ± 1.3 × 10 6 1.85 × 10 5 ± 1.3 × 10 6
SENSEX0.135 ± 0.0260.167 ± 0.0270.158 ± 0.03244,861.1 ± 1500.06759.4 ± 670.04589.2 ± 580.05.03 × 10 5 ± 1.6 × 10 6 6.12 × 10 5 ± 5.7 × 10 6 6.37 × 10 5 ± 6.1 × 10 6
SHSZ3000.365 ± 0.0640.37 ± 0.060.407 ± 0.06462,652.2 ± 2700.041,541.1 ± 1600.040,719.1 ± 1900.07.82 × 10 5 ± 4.4 × 10 6 4.59 × 10 5 ± 5.2 × 10 6 4.08 × 10 5 ± 6.3 × 10 6
SMI0.18 ± 0.030.065 ± 0.0110.058 ± 0.00638,389.3 ± 1200.02128.6 ± 340.0968.2 ± 190.02.97 × 10 5 ± 8.7 × 10 7 2.16 × 10 5 ± 2.2 × 10 6 2.17 × 10 5 ± 2.7 × 10 6
SPTSX0.249 ± 0.0410.245 ± 0.0550.18 ± 0.0449,544.5 ± 1600.032,659.8 ± 1500.028,436.5 ± 1200.03.9 × 10 5 ± 3 × 10 6 3.94 × 10 5 ± 7.6 × 10 6 2.85 × 10 5 ± 4.3 × 10 6
SPX0.08 ± 0.0120.1 ± 0.020.1 ± 0.0246,020.7 ± 900.033,298.2 ± 780.031,047.2 ± 750.03.45 × 10 5 ± 4.5 × 10 6 2.27 × 10 5 ± 1.7 × 10 6 1.97 × 10 4 ± 2 × 10 6
TWSE0.144 ± 0.0170.198 ± 0.0350.17 ± 0.0445,756.3 ± 810.033,299.1 ± 1800.031,386.4 ± 1700.05.84 × 10 5 ± 1.6 × 10 6 6.3 × 10 5 ± 4 × 10 6 6.60 × 10 5 ± 3.6 × 10 6
UKX0.177 ± 0.0260.113 ± 0.0270.105 ± 0.02639,896.9 ± 950.014,552.5 ± 930.011,522.6 ± 800.03.56 × 10 5 ± 2.2 × 10 6 3.11 × 10 5 ± 2.9 × 10 6 3.41 × 10 5 ± 2.9 × 10 6
Table 6. Minimum and maximum portfolio retention under varying levels of the transaction cost penalty λ and number of particles (N). Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Table 6. Minimum and maximum portfolio retention under varying levels of the transaction cost penalty λ and number of particles (N). Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Min RetentionMax Retention
N = 10 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.019 ± 0.0140.142 ± 0.0120.188 ± 0.0210.501 ± 0.0260.576 ± 0.0190.635 ± 0.022
DAX0.0031 ± 0.00610.136 ± 0.0110.16 ± 0.0150.465 ± 0.0230.576 ± 0.0270.571 ± 0.022
HSI0.049 ± 0.0120.102 ± 0.0120.138 ± 0.0130.388 ± 0.0160.455 ± 0.0150.459 ± 0.019
IBEX0.013 ± 0.0120.149 ± 0.0170.2 ± 0.020.546 ± 0.0280.649 ± 0.0250.684 ± 0.021
IBOV0.05 ± 0.010.109 ± 0.0130.131 ± 0.0140.389 ± 0.0160.445 ± 0.0170.463 ± 0.012
IPSA0 ± 00.136 ± 0.0220.1661 ± 0.00920.548 ± 0.0310.611 ± 0.0250.654 ± 0.015
JALSH0.102 ± 0.0110.14 ± 0.010.151 ± 0.0110.372 ± 0.0130.402 ± 0.0120.398 ± 0.011
MERVAL0 ± 00.152 ± 0.0280.24 ± 0.020.671 ± 0.0440.749 ± 0.0310.796 ± 0.025
NKY0.1329 ± 0.00680.1611 ± 0.00840.1667 ± 0.00960.3314 ± 0.00810.354 ± 0.0090.355 ± 0.0072
SENSEX0.011 ± 0.0120.18 ± 0.020.247 ± 0.0210.56 ± 0.0260.702 ± 0.0260.706 ± 0.022
SHSZ3000.1115 ± 0.00850.1299 ± 0.00740.141 ± 0.0070.302 ± 0.00920.3266 ± 0.00770.3362 ± 0.0088
SMI0.042 ± 0.0230.255 ± 0.0260.339 ± 0.0220.721 ± 0.0360.866 ± 0.0350.912 ± 0.031
SPTSX0.1154 ± 0.00760.1429 ± 0.00860.1551 ± 0.00740.324 ± 0.0110.343 ± 0.00740.351 ± 0.008
SPX0.1271 ± 0.00770.143 ± 0.0070.1492 ± 0.00680.3066 ± 0.00850.3226 ± 0.00780.3361 ± 0.0072
TWSE0.1214 ± 0.00710.1511 ± 0.00840.155 ± 0.00840.3138 ± 0.00930.328 ± 0.00750.3354 ± 0.0087
UKX0.07 ± 0.0120.131 ± 0.0120.152 ± 0.0110.392 ± 0.0140.419 ± 0.0140.426 ± 0.016
Min RetentionMax Retention
N = 100 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10000
CAC0.013 ± 0.0120.231 ± 0.0160.3 ± 0.020.501 ± 0.0260.658 ± 0.0230.727 ± 0.026
DAX0.0071 ± 0.00980.211 ± 0.0210.264 ± 0.0170.479 ± 0.0220.635 ± 0.0240.682 ± 0.021
HSI0.051 ± 0.0110.156 ± 0.0130.186 ± 0.0110.398 ± 0.0190.511 ± 0.0210.508 ± 0.013
IBEX0.01 ± 0.010.25 ± 0.020.332 ± 0.0220.543 ± 0.0310.717 ± 0.0270.795 ± 0.028
IBOV0.0538 ± 0.00970.169 ± 0.0130.18 ± 0.010.372 ± 0.0140.494 ± 0.0170.503 ± 0.018
IPSA0 ± 00.212 ± 0.0210.29 ± 0.020.545 ± 0.0310.726 ± 0.0250.773 ± 0.026
JALSH0.09 ± 0.010.1754 ± 0.00960.181 ± 0.0110.355 ± 0.0130.427 ± 0.0140.435 ± 0.012
MERVAL0 ± 00.247 ± 0.0270.399 ± 0.0380.603 ± 0.0380.863 ± 0.0360.964 ± 0.025
NKY0.1184 ± 0.00860.1767 ± 0.00710.1744 ± 0.00890.3191 ± 0.00930.357 ± 0.0110.3661 ± 0.0084
SENSEX0.011 ± 0.0120.285 ± 0.0190.359 ± 0.0240.55 ± 0.0290.74 ± 0.030.8 ± 0.02
SHSZ3000.115 ± 0.0070.1508 ± 0.00690.1586 ± 0.00770.308 ± 0.0090.3414 ± 0.00960.3442 ± 0.0087
SMI0.008 ± 0.0110.403 ± 0.0190.491 ± 0.0130.74 ± 0.040.933 ± 0.0310.981 ± 0.018
SPTSX0.109 ± 0.0080.1505 ± 0.00950.1699 ± 0.00670.317 ± 0.00740.363 ± 0.0110.3599 ± 0.0091
SPX0.1136 ± 0.00980.161 ± 0.00840.1704 ± 0.00680.2989 ± 0.00810.3304 ± 0.00660.3442 ± 0.0087
TWSE0.1086 ± 0.00780.1651 ± 0.00750.1657 ± 0.00940.3066 ± 0.00820.3462 ± 0.00950.35 ± 0.01
UKX0.084 ± 0.0110.171 ± 0.0130.183 ± 0.0110.383 ± 0.0190.46 ± 0.0120.474 ± 0.015
Min RetentionMax Retention
N = 1000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.013 ± 0.0120.329 ± 0.0180.38 ± 0.020.457 ± 0.0230.693 ± 0.0190.789 ± 0.025
DAX0.0031 ± 0.00610.294 ± 0.0150.361 ± 0.0220.514 ± 0.0280.7 ± 0.0280.747 ± 0.019
HSI0.049 ± 0.0120.208 ± 0.0120.233 ± 0.0140.405 ± 0.0160.519 ± 0.0140.535 ± 0.016
IBEX0.01 ± 0.0110.371 ± 0.0220.442 ± 0.0260.539 ± 0.0230.787 ± 0.0280.86 ± 0.018
IBOV0.05 ± 0.0120.199 ± 0.0120.231 ± 0.0130.394 ± 0.0180.513 ± 0.0130.517 ± 0.013
IPSA0 ± 00.334 ± 0.0180.402 ± 0.0220.626 ± 0.0460.802 ± 0.0180.827 ± 0.025
JALSH0.09 ± 0.010.182 ± 0.0110.2 ± 0.010.346 ± 0.0130.429 ± 0.0120.459 ± 0.012
MERVAL0 ± 00.434 ± 0.0280.508 ± 0.0270.599 ± 0.0360.989 ± 0.0150.989 ± 0.015
NKY0.1185 ± 0.00860.1861 ± 0.00860.1822 ± 0.00720.32 ± 0.010.3678 ± 0.00780.3778 ± 0.0095
SENSEX0.033 ± 0.0190.364 ± 0.0280.451 ± 0.0260.583 ± 0.0320.8 ± 0.0260.864 ± 0.023
SHSZ3000.1135 ± 0.00790.1629 ± 0.00690.171 ± 0.0080.2982 ± 0.00760.3465 ± 0.00720.36 ± 0.01
SMI0.033 ± 0.0210.495 ± 0.00680.51 ± 0.00940.785 ± 0.0310.9958 ± 0.00820 ± 0
SPTSX0.1114 ± 0.00920.1578 ± 0.00850.1744 ± 0.00870.31 ± 0.010.3615 ± 0.00810.378 ± 0.011
SPX0.1213 ± 0.00790.1699 ± 0.00810.1733 ± 0.00790.3 ± 0.00880.3506 ± 0.00770.3536 ± 0.0081
TWSE0.1174 ± 0.00690.169 ± 0.0080.1798 ± 0.00650.31 ± 0.010.3515 ± 0.00730.3586 ± 0.0089
UKX0.078 ± 0.0130.1968 ± 0.00910.2294 ± 0.00970.372 ± 0.0130.488 ± 0.0180.497 ± 0.014
Table 7. Max weight and mean portfolio retention under varying levels of the transaction cost penalty λ and number of particles (N). Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Table 7. Max weight and mean portfolio retention under varying levels of the transaction cost penalty λ and number of particles (N). Values are averages based on 40 independent runs of the algorithm, along with a plus/minus value of 1.96 standard deviations of these averages (i.e., corresponding to the upper and lower bounds of a 95% confidence interval).
Max WeightMean Retention
N = 10 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.567 ± 0.0270.622 ± 0.0230.708 ± 0.0270.229 ± 0.0110.3462 ± 0.00770.4 ± 0.011
DAX0.709 ± 0.0520.613 ± 0.0350.699 ± 0.0320.197 ± 0.0110.33 ± 0.010.3547 ± 0.0094
HSI0.767 ± 0.0380.36 ± 0.0210.475 ± 0.0430.2038 ± 0.00720.2687 ± 0.00660.2887 ± 0.0063
IBEX0.704 ± 0.0390.706 ± 0.0340.778 ± 0.0290.25 ± 0.010.39 ± 0.010.4419 ± 0.0091
IBOV0.48 ± 0.0330.399 ± 0.0280.465 ± 0.0410.2072 ± 0.00650.271 ± 0.00630.2931 ± 0.0054
IPSA0.968 ± 0.0240.64 ± 0.0310.71 ± 0.030.225 ± 0.0120.358 ± 0.0110.4048 ± 0.0082
JALSH0.288 ± 0.0410.335 ± 0.0230.362 ± 0.0270.2249 ± 0.00610.2671 ± 0.00560.2765 ± 0.0051
MERVAL0.971 ± 0.0130.809 ± 0.0270.82 ± 0.020.267 ± 0.0150.433 ± 0.0110.512 ± 0.014
NKY0.275 ± 0.0340.179 ± 0.0150.206 ± 0.0140.2273 ± 0.00390.255 ± 0.0040.2646 ± 0.0045
SENSEX0.604 ± 0.0310.666 ± 0.0330.73 ± 0.030.26 ± 0.010.42 ± 0.010.47 ± 0.011
SHSZ3000.813 ± 0.0140.1371 ± 0.00660.139 ± 0.0120.2056 ± 0.00410.2313 ± 0.00330.235 ± 0.0033
SMI0.763 ± 0.0380.804 ± 0.0180.726 ± 0.0360.361 ± 0.0120.537 ± 0.0130.6113 ± 0.0096
SPTSX0.31 ± 0.0230.15 ± 0.010.176 ± 0.0180.2124 ± 0.00410.2414 ± 0.00340.2493 ± 0.0034
SPX0.322 ± 0.0270.138 ± 0.0080.149 ± 0.0120.2116 ± 0.00330.231 ± 0.00320.2412 ± 0.0031
TWSE0.276 ± 0.0250.1287 ± 0.00750.1336 ± 0.00850.2163 ± 0.00350.2372 ± 0.00320.243 ± 0.004
UKX0.438 ± 0.0210.391 ± 0.0370.41 ± 0.040.222 ± 0.0050.274 ± 0.0050.286 ± 0.006
Max WeightMean Retention
N = 100 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.586 ± 0.0360.634 ± 0.0220.61 ± 0.0270.2273 ± 0.00970.4294 ± 0.00870.5 ± 0.01
DAX0.847 ± 0.0530.717 ± 0.0230.735 ± 0.0260.214 ± 0.0120.41 ± 0.010.47 ± 0.01
HSI0.809 ± 0.0310.49 ± 0.020.55 ± 0.020.208 ± 0.00620.3172 ± 0.00670.3433 ± 0.0066
IBEX0.785 ± 0.0350.659 ± 0.0310.624 ± 0.0210.231 ± 0.0120.4794 ± 0.00990.557 ± 0.009
IBOV0.497 ± 0.0360.477 ± 0.0230.49 ± 0.0220.2106 ± 0.00630.3213 ± 0.00590.3365 ± 0.0057
IPSA0.9991 ± 0.00160.748 ± 0.0290.747 ± 0.0280.234 ± 0.0130.453 ± 0.0110.5179 ± 0.0092
JALSH0.293 ± 0.0350.345 ± 0.0170.396 ± 0.0190.218 ± 0.0050.297 ± 0.00560.3074 ± 0.0045
MERVAL0.96 ± 0.0080.694 ± 0.0270.685 ± 0.0260.243 ± 0.0140.582 ± 0.0140.677 ± 0.015
NKY0.34 ± 0.0350.218 ± 0.0150.256 ± 0.0130.219 ± 0.0040.2649 ± 0.00330.2747 ± 0.0041
SENSEX0.622 ± 0.0250.618 ± 0.0250.58 ± 0.020.26 ± 0.010.494 ± 0.0120.578 ± 0.012
SHSZ3000.819 ± 0.0170.198 ± 0.0150.214 ± 0.0180.2072 ± 0.00360.241 ± 0.0030.251 ± 0.0037
SMI0.804 ± 0.0250.599 ± 0.0270.501 ± 0.0250.357 ± 0.0120.6483 ± 0.00920.738 ± 0.011
SPTSX0.348 ± 0.0360.223 ± 0.0170.262 ± 0.0140.2099 ± 0.00440.2539 ± 0.00340.2608 ± 0.0039
SPX0.385 ± 0.0440.209 ± 0.0170.236 ± 0.0150.2047 ± 0.00340.2458 ± 0.00330.2541 ± 0.0048
TWSE0.248 ± 0.0250.159 ± 0.0110.197 ± 0.0160.2054 ± 0.00310.249 ± 0.00370.2523 ± 0.0039
UKX0.476 ± 0.0370.43 ± 0.020.452 ± 0.0140.2132 ± 0.00550.3158 ± 0.00550.3249 ± 0.0057
Max WeightMean Retention
N = 1000 λ = 0 λ = 1000 λ = 10,000 λ = 0 λ = 1000 λ = 10,000
CAC0.711 ± 0.0370.524 ± 0.0260.502 ± 0.0240.2208 ± 0.00970.5016 ± 0.00980.5695 ± 0.0095
DAX0.93 ± 0.040.568 ± 0.0270.605 ± 0.0270.2195 ± 0.00970.49 ± 0.0080.5488 ± 0.0092
HSI0.793 ± 0.0260.462 ± 0.0170.47 ± 0.020.209 ± 0.0060.3572 ± 0.00510.38 ± 0.006
IBEX0.832 ± 0.0270.551 ± 0.0260.557 ± 0.0330.24 ± 0.0110.56 ± 0.010.6378 ± 0.0093
IBOV0.523 ± 0.0240.463 ± 0.0180.449 ± 0.0180.2071 ± 0.00670.3478 ± 0.00570.3732 ± 0.0067
IPSA0.9972 ± 0.00560.606 ± 0.0270.63 ± 0.0270.26 ± 0.0110.568 ± 0.00890.61 ± 0.01
JALSH0.271 ± 0.0320.35 ± 0.020.38 ± 0.0140.2129 ± 0.00590.303 ± 0.0050.3204 ± 0.0053
MERVAL0.9597 ± 0.00340.613 ± 0.0230.583 ± 0.0210.259 ± 0.0120.714 ± 0.0130.773 ± 0.015
NKY0.324 ± 0.0330.232 ± 0.0120.2529 ± 0.00960.2161 ± 0.00470.273 ± 0.0030.2789 ± 0.0034
SENSEX0.635 ± 0.0260.533 ± 0.0180.484 ± 0.0210.29 ± 0.010.587 ± 0.0110.659 ± 0.013
SHSZ3000.812 ± 0.0170.219 ± 0.0130.237 ± 0.0150.2048 ± 0.00410.2527 ± 0.00330.2603 ± 0.0031
SMI0.75 ± 0.0290.435 ± 0.0230.437 ± 0.0190.37 ± 0.00990.75 ± 0.0110.837 ± 0.012
SPTSX0.38 ± 0.0470.262 ± 0.0210.256 ± 0.0110.2086 ± 0.00430.261 ± 0.00440.2728 ± 0.0044
SPX0.44 ± 0.0570.236 ± 0.0190.25 ± 0.0150.205 ± 0.00470.2583 ± 0.00370.2647 ± 0.0032
TWSE0.277 ± 0.0250.195 ± 0.0120.226 ± 0.0110.2089 ± 0.00430.2633 ± 0.00440.2682 ± 0.0041
UKX0.51 ± 0.0310.375 ± 0.0170.385 ± 0.0130.2154 ± 0.00560.3268 ± 0.00610.3546 ± 0.0051
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hamilton-Russell, L.; O’Callaghan, T.M.; Savin, D.; Schlögl, E. Modifying Sequential Monte Carlo Optimisation for Index Tracking to Allow for Transaction Costs. Risks 2024, 12, 155. https://doi.org/10.3390/risks12100155

AMA Style

Hamilton-Russell L, O’Callaghan TM, Savin D, Schlögl E. Modifying Sequential Monte Carlo Optimisation for Index Tracking to Allow for Transaction Costs. Risks. 2024; 12(10):155. https://doi.org/10.3390/risks12100155

Chicago/Turabian Style

Hamilton-Russell, Leila, Thomas Malan O’Callaghan, Dmitrii Savin, and Erik Schlögl. 2024. "Modifying Sequential Monte Carlo Optimisation for Index Tracking to Allow for Transaction Costs" Risks 12, no. 10: 155. https://doi.org/10.3390/risks12100155

APA Style

Hamilton-Russell, L., O’Callaghan, T. M., Savin, D., & Schlögl, E. (2024). Modifying Sequential Monte Carlo Optimisation for Index Tracking to Allow for Transaction Costs. Risks, 12(10), 155. https://doi.org/10.3390/risks12100155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop