Abstract
We propose two numerical methods, namely the alternating block relaxation method and the alternating majorization method, for the problem of nearest correlation matrix with factor structure, which is highly nonconvex. In the block relaxation method, the subproblem is of the standard trust region problem, which is solved by Steighaug’s truncated conjugate gradient method or by the exact trust region method. In the majorization method, the subproblem has a closed-form solution. We then apply the majorization method to the case where nonnegative factors are required. The numerical results confirm that the proposed methods work quite well and are competitive against the best available methods.
Similar content being viewed by others
References
Albrecher, H., Ladoucette, S., Schoutens, W.: A generic one-factor Lévy model for pricing synthetic CDOs. In: Fu, M.C., Jarrow, R.A., Yen, J.-Y.J., Elliott, R.J. (eds.) Advances in Mathematical Finance. Applied and Numerical Harmonic Analysis, pp. 259–277. Birkhaüser, Boston (2007)
Birgin, E.G., Martínez, J.M., Raydan, M.: Nonmonotone spectral projected gradient methods on convex sets. SIAM J. Control Optim. 10, 1196–1211 (2000)
Birgin, E.G., Martínez, J.M., Raydan, M.: Algorithm 813: SPG-Software for convex-constrained optimization. ACM Trans. Math. Softw. 27, 340–349 (2001)
Borsdorf, R., Higham, N.J.: A preconditioned Newton algorithm for the nearest correlation matrix. IMA J. Numer. Anal. 94, 94–107 (2010)
Borsdorf, R., Higham, N.J., Raydan, M.: Computing a nearest correlation matrix with factor structure. SIAM J. Matrix Anal. Appl. 31, 2603–2622 (2010)
Boyd, S., Xiao, L.: Least-squares covariance matrix adjustment. SIAM J. Matrix Anal. Appl. 27, 532–546 (2005)
de Leeuw, J.: Application of convex analysis to multidimensional scaling. In: van Cutsem, B., et al. (eds.) Recent Advantages in Statistics. North Holland, Amsterdam (1977)
de Leeuw, J.: Block relaxation algorithms in statistics. In: Bock, H.H., et al. (eds.) Information Systems and Data Analysis, pp. 308–325. Springer, Berlin (1994)
Flury, B.: Common Principle Components and Related Multivariate Models. Wiley, New York (1988)
Gao, Y., Sun, D.F.: Calibrating least squares covariance matrix problems with equality and inequality constraints. SIAM J. Matrix Anal. Appl. 31, 1432–1457 (2009)
Golub, G.H., Van Loan, C.F.: Matrix Optimization, 3rd edn. John Hopkins University Press, Baltimore (1996)
Kiers, H.A.L.: Setting up alternating least squares and iterative majorization algorithms for solving various matrix optimization problems. Comput. Stat. Data Anal. 41, 150–170 (2002)
Heiser, W.J.: Convergent computation by iterative majorization: theory and applications in multidimensional data analysis. In: Krzanowski, W.J. (eds.) Recent Advances in Descriptive Multivariate Analysis, pp. 157–189. Oxford University Press, Oxford (1995)
Higham, N.J.: Computing the nearest correlation matrix—a problem from finance. IMA J. Numer. Anal. 22, 329–343 (2002)
Hull, J.C., White, A.: Forward rate volatilities, swap rate volatilities, and the implementation of the LIBOR market model. J. Fixed Income 10, 46–62 (2000)
Malick, J.: A dual approach to semidefinite least-squares problems. SIAM J. Matrix Anal. Appl. 26, 272–284 (2004)
Moré, J.J., Sorensen, D.C.: Computing a trust region step. SIAM J. Sci. Comput. 4, 553–572 (1983)
Nocedal, J., Wright, S.J.: Numerical Optimization. Springer, Berlin (1999)
Pietersz, R.: Pricing Models for Bermudan-Style Interest Rate Derivatives. ERIM PhD Series in Management. Erasmus University, Rotterdam (2005)
Pietersz, R., Grubišić, I.: Rank reduction of correlation matrices by majorization. Quant. Finance 4, 649–662 (2004)
Qi, H.D., Sun, D.F.: A quadratically convergent Newton method for computing the nearest correlation matrix. SIAM J. Matrix Anal. Appl. 28, 360–385 (2006)
Rojas, M., Santos, S.A., Sorensen, D.C.: Algorithm 873: LSTRS: MATLAB software for large-scale trust region subproblems and regularization. ACM Trans. Math. Softw. 34(2), 11 (2008)
Sonneveld, P., van Kan, J.J.I.M., Huang, X., Oosterlee, C.W.: Nonnegative matrix factorization of a correlation matrix. Linear Algebra Appl. 431, 334–349 (2009)
Steihaug, T.: The conjugate gradient method and trust regions in large scale optimization. SIAM J. Numer. Anal. 20, 626–637 (1983)
Toh, K.C.: An inexact path-following algorithm for convex quadratic SDP. Math. Program. 112, 221–254 (2008)
Toh, K.C., Tütüncü, R.H., Todd, M.J.: Inexact primal-dual path-following algorithms for a special class of convex quadratic SDP and related problems. Pac. J. Optim. 3, 135–164 (2007)
Yuan, Y.: On the truncated conjugate gradient method. Math. Program., Ser. A 87(3), 561–573 (2000)
Author information
Authors and Affiliations
Corresponding author
Additional information
Dedicated to Professor Liqun Qi on the occasion of his 65th birthday.
The work of N. Xiu was supported by the National Basic Research Program of China (2010CB732501).
Rights and permissions
About this article
Cite this article
Li, Q., Qi, H. & Xiu, N. Block relaxation and majorization methods for the nearest correlation matrix with factor structure. Comput Optim Appl 50, 327–349 (2011). https://doi.org/10.1007/s10589-010-9374-y
Received:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10589-010-9374-y