Network models to improve robot advisory portfolios

Robot advisory services are rapidly expanding, responding to a growing interest people have in directly managing their savings. Robot-advisors may reduce costs and improve the quality of asset allocation services, making user’s involvement more transparent. Against this background, there exists the possibility that robot advisors underestimate market risks, especially during crisis times, when high order interconnections arise. This may lead to a mismatch between investors’ expected and actual risk. The aim of this paper is to overcome this issue, taking into account not only investors’ risk preference but also their attitude towards interconnectdness. To achieve this aim, we combine random matrix theory with correlation networks and extend the Markowitz’ optimisation problem to a third dimension. To demonstrate the practical advantage of our proposed approach we employ daily returns of a large set of Exchange Traded Funds, which are representative of the financial products employed by robot-advisors.


Introduction and literature review
Financial Technologies (FinTech) can be broadly defined as technologically enabled financial innovations that could result in new business models, applications, processes or products, with an associated material effect on financial markets, financial institutions, and on the provision of financial services (Carney 2017). In the last few years, FinTech innovations have increased exponentially, delivering new payments and lending methods, and penetrating the insurance and asset management sectors. The Financial Stability Board (FSB) in its two recent reports (FSB 2017a, b) has identified three common drivers for FinTech innovation; namely, the shifting of consumer preferences on the demand side, the evolution of data driven technologies and the changes in the financial regulation on the supply side.
Against this background, robot advisory services for automated investments are growing fast to address the need of directly managing savings. They are accessible via online platforms and, therefore, allow to act quickly and in the first person. According to Statista, the masses managed by automatic consultancy are estimated to be over 2552 billion in 2023. 1 The rapid growth of robot-advisors has determined the emergence of new financial risks. Robot-advisors that build personalised portfolios on the basis of automated algorithms have been suspected of underestimating investors' risk preferences. This may be the effect of the asset allocation models employed by robot advisors which, for the sake of transparency and user's engagement, are often over simplified, and do not take asset returns correlations properly into account.
To improve robot-advisory asset allocation, we propose to embed correlation networks into Markowitz' asset allocation, following the work of Clemente et al. (2019). The authors modify the objective function of the Markowitz minimum variance portfolio taking into account not only the volatility of individual assets but also their interconnectedness, expressed in terms of a clustering coefficient calculated on the correlation network. Their empirical results show that the resulting portfolios are more diversified and show better performances with respect to the classical Markowitz' portfolios. Similarly to Clemente et al. (2019), our proposal is based on the insertion of correlation network models into Markowitz' objective function. Our additional contribution is a parsimonious correlation network, based on random matrix theory filtering, and an interconnectdness measure that exploits the notion of network centrality. These methodological advancements are particularly useful in the context of robot advisory asset allocation, characterised by many assets, highly correlated with each other.
A related research work is Boginski et al. (2014) who, building on the work of Pattillo et al. (2013) and Boginski et al. (2006) exploit the concept of clique relaxations in weighted graphs to find profitable well-diversified portfolios. In their proposal, the weight of each asset corresponds to its return over the considered time period and each pair of assets is connected if the corresponding correlation exceeds a certain threshold value. This mechanism is able to ensure high returns of portfolios that are not guaranteed by the cliques themselves.
Another related paper is He and Zhou (2011), who modify Markowitz' problem by exploiting a different utility function, introducing a new measure of loss aversion for large payoffs, called the large-loss aversion degree (LLAD). The measure is applied to portfolio choice under the cumulative prospect theory in Tversky and Kahneman (1992).
A different stream of extensions of the Markowitz portfolio model considers multiobjective evolutionary algorithms (MOEAs), as in Metaxiotis and Liagkouras (2012). Within this framework, Cesarone et al. (2013) propose a heuristic solution based on a reformulation in terms of a Standard Quadratic Program to solve mean-variance portfolio issues related to the introduction of constraints based on cardinality (which limits the number of assets to be held in an efficient portfolios) and allocation shares (which determines the fraction of capital invested in each asset). Whereas Ehrgott et al. (2004) presents a method based on the application of four different heuristic solution techniques to test problems involving up to 1416 assets. Woodside-Oriakhi et al. (2011) consider the application of genetic algorithm, tabu search and simulated annealing meta heuristic approaches to find the cardinality constrained efficient frontier that arises in portfolio optimization. In addition, Doerner et al. (2004) introduces Pareto Ant Colony Optimization as an effective meta-heuristic solution, proposing a two-stage procedure that first identifies the solution space of all efficient portfolios and then locates the best solution within that space.
We contribute to the above literature by proposing a different objective function of Markowitz' portfolio allocation, that extends the classic formulation with an asset centrality term, function of the assets' similarity network.
Similarity networks among asset returns were introduced by Mantegna and Stanley (1999)), who expressed the distance between any pair of asset returns as a function of the pairwise correlations among the corresponding time series. They can reveal how assets are related in terms of the topology of a network (Newman 2018), and can allow to calculate centrality measures, which express their "importance" of each asset in the network (see e.g. Avdjiev et al. 2019). Among them, the eigenvector centrality (Bonacich 2007) assigns a relative score to all nodes in the network, based on the principle that connections to few high scoring nodes contribute more to the score of the node in question than equal connections to low scoring nodes.
To reduce the computational burden involved with the calculation of the eigenvector centrality on a fully connected matrix, Mantegna (1999) suggested to hierarchically cluster the assets, or groups of assets that are the "closest", leading to a parsimonious "Minimal Spanning Tree" (MST) representation. Tumminello et al. (2005) has extended Mantegna (1999) with a generalisation of the MST, the Planar Maximally Filtered Graph (PMFG), which retains the same hierarchical properties of the MST but adds more complex graph structures, such as loops and cliques.
Another important extension of Mantegna (1999) is Tola et al. (2008) who have shown how a MST calculated on a correlation matrix "filtered" from random noise through the Random Matrix Theory (RMT) approach can improve the performance of the optimal portfolios. Siimilar papers, based on the applications of RMT to asset management, are León et al. (2017), Raffinot (2017), Ren et al. (2017), Zhan et al. (2015), Bun et al. (2017) and Fraham and Jaekel (2005).
Our contribution to the above literature is twofold. From an applied viewpoint, we extend the application of the RMT approach in Tola et al. (2008) to Exchange Traded Fund returns (ETFs), investment funds that aim to replicate the index to which they refer (benchmark) through a totally passive management. From a methodological viewpoint, we propose a portfolio optimisation approach different from what proposed by Tola et al. (2008), as we include network centrality explicitly into the Markowitz objective function. Doing so, we do not rely only on the pairwise covariance between assets returns, but we also consider higher order information on assets' behavior.
In the paper we show that the inclusion of centrality measures into Markowitz model can improve portfolio returns. This because not only bivariate but also multivariate dependencies are taken into account. A further advantage is that investors' risk preferences receive a higher importance and, therefore, the matching between the expected and the actual risk profile improves.
The empirical findings obtained from the application of our proposed method confirm the validity of the proposed approach, which can thus become a new toolbox for robot-advisors. The results demonstrate that centrality measures can generate portfolio allocation strategies able to outperform the benchmark portfolios.
The structure of the paper is as follows: Sect. 2 presents our proposal, in an order that follows the data analysis flow: first we describe the Random Matrix Theory, to filter data from noise components; then the minimal spanning tree approach, to build a parsimonious correlation network among assets; and, finally, the proposed objective function that adds to Markowitz' the centrality measures calculated from the obtained correlation network. Section 3 presents the results of the application of the methodology to a database kindly provided by an anonymous robot advisor. Section 4 ends with some concluding remarks.
The intuition behind RMT is to extract the "systematic" part of a signal embedded in a correlation matrix, separating it from the "noise" component. To achieve this goal, RMT tests the eigenvalues of a correlation matrix: λ k < λ k+1 ; k = 1, . . . , N , against the null hypothesis that they are equal to the eigenvalues of a random Wishart matrix R = 1 T AA T of the same size, with A being a N × T matrix containing N time series of length T whose elements are independent and identically distributed standard Gaussian random variables.
Note that RMT does not require any assumption on the distribution of the asset returns. Rather, it compares the eigenvalues of the empirical covariance matrix with those that would be obtained if the returns were drawn from independent Gaussian distributions, leading to a Wishart matrix.
Let (x k =λ k ) < (x k+1 =λ k+1 ); k = 1, . . . , N be the observed sample eigenvalues. It can be shown (Marchenko and Pastur 1967) that, as N → ∞ and T → ∞, with a fixed ratio Q = T N ≥ 1, the density of each sample eigenvalue converges to: where x k ∈ (λ − , λ + ) and λ ± = 1 + 1 Q ± 2 1 Q . When x k > λ + the null hypothesis is rejected, for all the eigenvalues greater or equal than k. This implies that the relevant part of the signal contained in the correlation matrix of the returns can be obtained applying a singular value decomposition based only on the eigenvectors that correspond to eigenvalues that are greater than λ + . Doing so, RMT simplifies the correlation matrix into a filtered correlation matrix (Plerou et al. 2002;Eom et al. 2009).
More formally, let r i , for i = 1, . . . , N , be a time series of asset returns, computed, for any given time point t, as the difference between the logarithms of daily asset prices: (2) Given a set of N asset return series, a correlation coefficient between any two pairs can be defined as: where E(•) and σ (•) indicate, respectively, the mean and the standard deviation operators. Let C be the the correlation matrix.
According to the RMT the filtered correlation matrix is then given by: and V represents the matrix of the eigenvectors associated to the eigenvalues greater than λ + . We remark that the fact that the eigenvalues leaking out the "bulk" of the spectral density, the so called "spikes", contain the important information can also be obtained following Couillet (2015), who assume that the correlation matrix of the stock returns follows a spiked covariance matrix, in which all eigenvalues are equal to one except "r" spikes. In the random matrix theory approach the number of spikes is chosen comparing the sample eigenvalues with the theoretical ones obtained from a Wishart matrix.

Minimal spanning tree
When a large set of asset returns is considered, as in the case of Robot-Advisory, the (filtered) correlation matrix may be difficult to summarise. Its representation in terms of a correlation network (see e.g. Mantegna 1999) can help the task.
A correlation network can be obtained converting pairwise correlations in pairwise distances with the following function: where c i j are the elements of the filtered correlation matrix C and d i j is the Euclidean distance between return i and j. The set of all obtained pairwise distances can be organised into a distance matrix D = {d i j }. Then, a more parsimonious representation of the correlation matrix can be obtaned by means of the Minimal Spanning Tree method (MST, see e.g. (see Mantegna and Stanley 1999;Bonanno et al. 2003;Spelta and Araújo 2012)). The MST is obtained applying to the distance matrix D a single linkage clustering algorithm which associates each asset return to its closest neighbour, and avoiding loops. The term "minimal" refers to the fact that MST allows to reduce the number of links between asset returns from N (N −1) 2 to N − 1, the minimum number of links assuring connectivity of all nodes.
More formally, the MST algorithm proceeds in this way. Initially, it considers N clusters, corresponding to the N available asset returns (ETFs in our context). Then, at each subsequent step, two clusters l i and l j are merged into a single cluster if: with the distance between clusters being defined as: with p ∈ l i and q ∈ l j . The above steps are repeated until a single cluster emerges. 2 We remark that Raffinot (2017) extends the MST considering some clustering variants, such as complete linkage (CL), average linkage (AL) and Ward's Method (WM). He shows however that different algorithms differ in terms of grouping structures, but not in terms of performance.
We also remark that, to detect how financial relationships evolve over time we follow Spelta and Araújo (2012) by employing the residuality coefficient measure (R) that compares the relative strengths of the connections above and below a threshold value, ias follows: where L is the highest distance value that ensures the whole connectivity of the MST. We expect that during crisis phases, higher correlation patterns emerge, leading to a lower value of R, as the number of links increases, and vice-versa.

Eigenvector centrality
Having found a parsimonious representation of the return correlations, we now aim to summarise it, to understand which nodes (asset returns) act as hubs in the network. This is key for understanding how ETF returns behave in a multidimensional space, and to construct optimal portfolios that take the curse of dimensionality into account.
The research in network theory has dedicated much effort to develop measures aimed to detect the most important players in a network. The idea of "centrality" was initially proposed in the context of social systems, where a relationship between the location of a subject in the social network and its influence on the group processes was assumed.
Various measures of network centrality have been proposed, such as the count of the neighbors a node has: the degree centrality, which is a local centrality measure, or measures based on the spectral properties of the adjacency matrix (see Perra and Fortunato 2008), which are global measures. Examples of global centrality measures include the eigenvector centrality (Bonacich 2007), Katz's centrality (Katz 1953), PageRank (Brin and Page 1998), hub and authority centralities (Kleinberg 1999).
The eigenvector centrality measures the importance of a node by assigning relative scores to all nodes in the network, based on the principle that connections to few high scoring nodes contribute more to the score of the node in question than equal connections to low scoring nodes. More formally, for the i-th node, the centrality score is proportional to the sum of the scores of all nodes which are connected to it, as in the following: where x j is the score of node j,d i, j is the (i, j) element of the adjacency matrix of the network, λ is a constant and N is the number of nodes of the network. The previous equation can be rewritten in matrix notation: whereD is the adjacency matrix, λ is the eigenvalue of the matrixD, with an associated eigenvector x, an N-vector of scores (one for each node). Note that, in general, there will be many different eigenvalues λ for which a solution to the previous equation exists. However, the additional requirement that all the elements of the eigenvector be positive (a natural request in our context) implies (by the Perron-Frobenius theorem) that only the eigenvector corresponding to the largest eigenvalue provides the desired centrality measures. Therefore, once an estimate ofD is provided, network centrality scores can be obtained from the previous equation, as elements of the eigenvector associated to the largest eigenvalue.
We remark that, in the context of correlation networks that we are considering, the higher the centrality score associated to a node, the more the node is dissimilar from all other nodes in the network.

Portfolio construction
We now explain how centrality measures can be embedded in a portfolio optimisation framework, to improve portfolio performances.
Correlations between asset returns play a central role in investment theory and risk management, as key elements for optimisation problems as in Markowitz (1952) portfolio theory. It is natural that correlation networks, based on pairwise correlations, play an important role too.
Indeed, Onnela et al. (2003) have shown that the assets with the highest weights in Markowitz portfolios (Markowitz 1952) are always located in the outer nodes of a Minimal Spanning Tree, so that the optimal portfolios are mainly composed by assets in the periphery of the network, and not in its core. Pozzi et al. (2013) have shown that portfolios which include central assets are characterized by greater risks and lower returns with respect to portfolios which include peripheral assets. Giudici and Polinesi (2021), following Giudici and Abu-Hashish (2019) found a similar behaviour in crypto exchanges. Vỳrost et al. (2018) have suggested that network-based asset allocation strategies may improve risk/return trade-offs. Their work is based on the study of Peralta and Zareei (2016) which have found a negative relationships between asset return centralities and the optimal weights obtained under the Markowitz model.
Other authors have built on the above remarks by proposing novel portfolio optimisation strategies. For example, Plerou et al. (2002) and Conlon et al. (2007) have used the correlation matrix, filtered with the random matrix approach in the Markowitz model, and have shown that for the obtained portfolios the realized risk is closer to the expected one. Tola et al. (2008), combining MST with the RMT filtering, have shown performance improvement with respect to Markowitz portfolios. Finally, Tumminello et al. (2010) have demonstrated that the risk of the optimized portfolio obtained using a "filtered" correlation matrix is more stable than that associated with a "non filtered" matrix.
In line with the previous authors, we would like to to exploit centrality measures, based on the minimal spanning tree derived from the RMT filtered correlation matrix of asset returns to improve Markowitz portfolios. Differently from the previous authors, we extend Markowitz' approach using RMT and MST in the optimisation function itself, rather than applying Markowitz to the transformed (filtered and/or simplified) correlation matrix.
More formally, we propose to minimize the constrained objective function: where μ P indicates the mean return of the portfolio, obtained by averaging the mean return of each asset μ i ; the parameter γ represents a risk aversion coefficient, to be specified by investors, x i is the eigenvector centrality associated with the ETF return i and, finally, the The decision criteria specified by Eq. 9 can be better understood comparing it with the classical Markowitz approach. In the Markowitz model, efficient portfolios are found by solving the optimization problem: where λ is a parameter representing the investor's risk tolerance. If λ is large, 1 λ will be close to zero, meaning that the investor does not have much risk tolerance. Conversely, if λ is small, 1 λ will be large, placing more emphasis on returns.
We extend Markowitz model without modifying its underlying quadratic utility function. Specifically, we fix the second term and we minimise a modified version of the first term, which adds to the portfolio variance (first order risk) a new component that measures systemic risk (high order risk). To balance the two risks we introduce a parameter γ that defines their ratio, for a given level of return. When γ increases, investors prefer systemic risks over first order risks, and viceversa when γ decreases. The next Section contains a practical illustration of the implications of our proposal.
As a result, our optimal portfolio allocation formulation includes three dimensions: a first order risk, represented by the variance-covariance matrix, a high order (systemic) risk which depends on the network structure of ETF returns and, finally, the portfolio returns (fixed). We remark that the classical Markowitz efficient frontier is based only on two dimensions: the first order risk and the returns, and is therefore not comparable to our context. However, if we set γ equal to zero, our formulation is two dimensional and can, therefore, be placed on the efficient frontier and compared with Markowitz' solutions, as will be shown in the next Section.

Application
The data set we consider to illustrate the application of our proposal is composed by 92 ETFs returns' time series traded over the period January 2006-February 2018 (3173 daily observations). Table 1 shows the classification of the 92 ETFs in 11 asset classes, according to the classification provided by the Exchange where they are traded. From Table 1 note that the Emerging Market asset classes are the most frequent, followed by Corporate ETFs. Table 2 displays summary statistics for the considered asset classes and, specifically, the mean, variance and kurtosis of the returns' distribution, to describe their location and variability. From Table 2 note that the mean value of the returns is around 0 for each asset class, consistently with the efficient market hypothesis suggested by Malkiel and Fama (1970). Differently, the value assumed by the standard deviation depends on the considered asset class: Emerging Equity and Commodity classes are more volatile with respect to the Corporate classes. Moreover, the high values of the kurtosis confirm some known stylized facts: the distribution of most ETFs' returns tends to be non-Gaussian and heavy tailed.
To compare the behavior of the returns during financial crises and normal times the data set has been divided in two chronologically successive batches, from 2006 to 2012 (crisis)  and from 2013 to 2018 (post-crisis). Figure 1 provides temporal boxplots for ETFs' returns, grouped by their asset classes (as described in Table 1). Figure 1 shows that the volatility of the ETFs belonging to the Emerging Equity classes, as well as that of the Commodity asset class, is larger during the crisis period. This feature explains why their overall standard deviation, reported in 2, is the highest.

Transforming the correlation matrix: RMT and MST
In this subsection we show how the empirical correlation between ETF returns can be filtered and simplified, by means of the application of RMT and MST, described in the methodological section.
We first divide the data into consecutive overlapping time windows. The width of such windows has been set equal to T = 250 (12 trading months), with a window length of one month ( ∼ = 21 trading days) for a total of 140 overlapping windows. For each time window, we use 11 months ( ∼ = 229 trading days) of daily observations to build our model and the remained month to validate it. This means, in particular, that we calculate 140 correlation matrices between all 92 ETFs' returns, based on 11 months of data, to obtain the filtered correlation matrix applying the RMT approach, to derive the MST and the eigenvector centrality measure, and, finally, to derive the optimal portfolio, which is validated in an out-of-sample fashion using the twelfth month of each window. Figure 2 shows the ordered eigenvalue distribution of the empirical correlation matrix, for the last time window of the data set (March 2017-January 2018), compared with the theoretical Wishart correlation matrix that would be observed under random noise. Figure 2 shows that most of the eigenvalues' distribution lies between λ min and λ max , which are respectively equal to 0.16 and 2.71. This "bulk" may be considered as being generated by random fluctuations while the six deviating eigenvalues that are greater than λ max represent the effective characteristic dimension described by the correlation matrix.
Similar considerations can be made for other time windows, without loss of generality. As described in the methodological Section, if, for each time window, we reconstruct the correlation matrix using only the eigenvectors that correspond to the largest eigenvalues, we obtain a sequence of "filtered" correlation matrices which can be used to improve the Minimum Spanning Tree representation of the ETFs' returns. Figure 3 reports for both the filtered and the unfiltered correlation matrices and for each time window, the most central node, defined as the ETF return with the highest degree (the largest number of connected nodes) in the MST representation.
From Fig. 3 note that the application of the RMT filtering approach leads to different Minimal Spanning Tree configurations over time: the most central nodes are different and belong to different ETF classes. On the other hand, the Minimal Spanning Trees based on the unfiltered empirical correlation matrices do not seem to vary: the ETF labelled EIMI-IM, belonging to the Asia Emerging Market class, is for most of the time the most central node. To further evaluate how the MSTs dynamically change over time, we employ, as a summary measure, the Max link: the maximum distance value between two pairs of nodes used in the construction of the tree, and the residuality coefficient, which measures the ratio between links eliminated and maintained by the MST construction. Figure 4 shows the evolution of these two quantities over the considered period.
From Fig. 4 note that, during the 2008 financial crisis, the Max link sharply decreases, due to the decrease of most distances between ETF returns. This can be explained by the increased correlations between all returns, which synchronise during the crisis, consistently with the literature findings. While the Max link bounces back after the crisis, the residuality coefficient continues its decline until 2014. This may indicate the persistence of a set of strong connections in the market, that determine the relevance of a limited number of links.
To better understand the previous findings, Fig. 5 shows the MST topology during 2008, as representative of crisis times, and its topology during the last time window, taken as a reference period for a "business as usual" market phase. Figure 5 reflects how correlations increase during the crisis phase, leading to a high number of links in the network. Furthermore, in the crisis period, the MST reveals the importance of the Asian, American and World Emerging Market classes, which have the highest centralities. The importance of the American Emerging Market node declines post crisis, but the Asian class centrality remains high. This may explain the persistence of low values in the residuality coefficient, after the crisis phase.

Portfolio construction
We now present the application of our proposed portfolio strategy, in which the eigenvector centrality computed on the MST derived from the application of RMT to the empirical The optimal portfolio weights are obtained minimizing the constrained objective function in Eq. 9, for which the value of γ is set a priori, accordingly to the level of risk aversion of a   Cumulative returns for different portfolio strategies.The plot reports the cumulative Profit and Loss obtained using our proposal based on different values of γ , the "naive" strategy portfolio (orange line), the MSCI benchmark index (blue line) and, for robustness, our proposal based on an unfiltered correlation matrix, rather than RMT (Cov) and our proposal based on a Glasso regularisation, rather than MST (Glasso) hypothetical investor. A high value of γ indicates that, in the desired allocation, more central ETFs (such as the Emerging Markets ones) will have higher weights.
Once the optimal portfolio weights are derived, we can calculate portfolio returns, and the associated Profit and Loss, for the time windows described in the previous subsections. More precisely, we use rolling windows in each of which the last month acts as an out-of-sample month to predict. The remaining eleven months of observations are used as a build-up period, to apply RMT, MST, compute asset return centralities and obtain the consequent optimal portfolio weights. We then calculate the return which corresponds to the optimal portfolio over the next month, weighting each ETF with the obtained weights. Finally, we cumulate each monthly portfolio return, from December 2006 to February 2018, taking a re-balancing cost of 10 basis points into account.
According to the described computational procedure, Fig. 6 presents the cumulative returns obtained with different values of γ , using the model in Eq. 9. The figure also reports the portfolio Profit and Loss of a "naive" (equally weighted) strategy as well as the performance of a benchmark, the MSCI Index. We also compute the performances obtained by employing a non filtered correlation matrix and those obtained with a Glasso regularisation method of Friedman et al. (2008), rather than with the MST approach 3 The results in Fig. 6 are summarised in Table 3, which presents the annual Profit and Loss of each competing strategy. Figure 6 highlights that our proposed model performs better than the benchmark index, the "naive" portfolio strategy, the standard Markowitz portfolio (which correspond to γ = 0 and the portfolios obtained either without RMT or with the sparse Glasso regularisation. All of our strategies win in terms of end of sample cumulative returns, regardless the coefficient of the individual risk propension, ranging from γ = 0.05 to γ = 4. Note that the portfolio based on the non filtered covariance matrix produces the worst performance and, theefore, the RMT filter appears to be a fundamental condition for having adequate asset diversification in investment portfolios. Looking in more detail, during the crisis period (2007-2009) our strategy produces higher returns with respect to the competitor portfolios, as it correctly takes high order risks into account. However, it is not able to capture the growing rebound at the end of 2009. More generally, during non-crisis times, our strategy, despite producing positive returns, can not reach the performance of the other portfolios. This because, during normal times, it probably excessively overweights high order risks.
These results are indeed consistent with our proposed modification to the Markowitz' algorithm. We expect that, during crisis times, when high order interconnections among financial assets strongly increase, our strategy produces higher returns with respect to competitor portfolios, and that it will do especially so for higher values of γ . This intuition is confirmed: in the year 2008 profits increase, particularly when γ is high, and so does the annual Sharpe Ratio, as shown in Table 4.
To provide further insights on portfolio compositions, we report in Fig. 7 the dynamic of the portfolio weights for ETF classes, considering γ = 0.7. 4 From Fig. 7 it is clear that during crisis times the weight of the ETFs belonging to Emerging Equity classes is the highest. During non-crisis times, Emerging ETFs become less important, and their role is taken by other ETFs; in particular with Corporate ones.
To gain further insights about how portfolio performances change as market conditions mutate, the following tables report, as performance measures that take both risk and returns into consideration, the Sharpe Ratio (Sharpe 1994), the α of the Capital Asset Pricing Model (CAPM), the Value at Risk (VaR) and the Conditional VaR (CVaR). Table 4 specifically refers to the yearly Sharpe Ratio, defined as the ratio between the mean value of the excess returns and its standard deviation. Table 4 highlights how, during market crises (as in 2008), the Sharpe Ratio of our portfolio strategy is higher with respect to the "naive" one with respect to the Sharpe Ratio obtained with the Glasso regularisation method. The subsequent rebound of 2009 is not captured by our proposal and the lower Sharpe Ratio obtained under different vaues of γ reflect this feature. Notice that the worst values of the Sharpe Ratio are associated with the portfolio derived using a non-filtered covariance matrix: Table 3 Annual cumulative profits and losses     The table shows the annual Value at Risk for the portfolios related to the different strategies. All the values are expressed in absolute terms multiplied by a scale factor of 100   The table shows the Information Ratio of portfolios under different strategies. All measures are computed relatively to the benchmark strategy this explains, once more, the importance of the preliminary processing done by the RMT approach.
The value of the CAPM α measures the ability to choose potentially profitable assets, reflecting the expertise of asset managers in exploiting market signals and investing accordingly, thus generating positive extra-performances. Table 5 describes the α coefficient, which reflects portfolio extra/under performances with respect to the benchmark. Table 5 shows that our portfolios outperform the benchmark strategy, as they all report values greater than 0. They are also generally better than the "naive" and "Glasso" portfolios. Only during the growing rebound phase of 2009 our strategies seem to under perform.
Our portfolio allocation method can also be compared with other methods in terms of risk. Table 6 specifically refers to the Value at Risk. Table 6 highlights that our portfolio strategies, although becoming more risky during the crisis period (proportionally to risk aversion), has an overall risk that is lower than that of the benchmark portfolio and that of the "naive" one. Table 7 reports the values of the CVaR of the different portfolio strategies. This measure, introduced by Rockafellar et al. (2000), quantifies the potential extreme losses in the tail of the return distribution. The results are in line with those presented in Table 6: our strategies over perform in terms of expected losses the "naive" and benchmark portfolios, except for the last year considered in the analysis.
Finally, Table 8 reports the values of the information ratio which, differently from the Sharpe ratio, takes into account the tracking error. The results are in line with those obtained with the calculation of the Sharpe ratio.

Conclusions
In the paper we have shown how to improve robot-advisory portfolio allocation, typically based on several asset returns, highly correlated with each other, and especially during crisis times.
In particular, we have demonstrated how to build high performing portfolios by means of a mix of data processing strategies, that include: (i) filtering the correlation matrix among asset returns with the Random Matrix Theory approach; (ii) calculating the correlation network centrality of each asset return, after the application of the Minimal Spanning tree approach; (iii) selecting portfolio weights including in Markowitz' optimisation function the network centralities, thereby taking into account high order interconnection risks.
The application of our proposal to the observed returns of a set of Exchange Traded Funds (ETFs), which are representative of the assets traded by robot-advisors, shows that our proposal leads to higher returns, or to lower risks, when compared to standard portfolios, especially during crisis times.
We therefore believe that our proposal could be relevant, for robot-advisors aimed at improving their services, while maintaining the accessibility to a large audience of potential investors; but also for regulators and supervisors, aimed at measuring and preventing the under estimation of market risks coming from the adoption of robot-advisory financial consulting.
The research work could be extended in several directions, and we believe it would also be extremely important to apply what proposed to other datasets and robot-advisory settings and, in particular, to those concerning crypto assets (see e.g. Giudici and Pagnottoni 2020) and foreign exchanges (see e.g. . gramme" under the Grant Agreement No. 825215 (Topic: ICT-35-2018, Type of action: CSA). The authors thank the editor and the two anonymous referees for useful comments and suggestions. They also thank the company Moneyfarm, for having provided the data.
Funding Open access funding provided by Universitá Politecnica delle Marche within the CRUI-CARE Agreement.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.