hide long namesshow long names
hide short namesshow short names
Integer type:  int32  int64  nag_int  show int32  show int32  show int64  show int64  show nag_int  show nag_int

PDF version (NAG web site, 64-bit version, 64-bit version)
Chapter Contents
Chapter Introduction
NAG Toolbox

NAG Toolbox: nag_correg_corrmat_nearest_kfactor (g02ae)

Purpose

nag_correg_corrmat_nearest_kfactor (g02ae) computes the factor loading matrix associated with the nearest correlation matrix with kk-factor structure, in the Frobenius norm, to a given square, input matrix.

Syntax

[g, x, iter, feval, nrmpgd, ifail] = g02ae(g, k, 'n', n, 'errtol', errtol, 'maxit', maxit)
[g, x, iter, feval, nrmpgd, ifail] = nag_correg_corrmat_nearest_kfactor(g, k, 'n', n, 'errtol', errtol, 'maxit', maxit)

Description

A correlation matrix CC with kk-factor structure may be characterised as a real square matrix that is symmetric, has a unit diagonal, is positive semidefinite and can be written as C = XXT + diag(IXXT)C=XXT+diag(I-XXT), where II is the identity matrix and XX has nn rows and kk columns. XX is often referred to as the factor loading matrix.
nag_correg_corrmat_nearest_kfactor (g02ae) applies a spectral projected gradient method to the modified problem min GXXT + diag(XXTI)FminG-XXT+diag(XXT-I)F such that xiT21xiT21, for i = 1,2,,ni=1,2,,n, where xixi is the iith row of the factor loading matrix, XX, which gives us the solution.

References

Birgin E G, Martínez J M and Raydan M (2001) Algorithm 813: SPG–software for convex-constrained optimization ACM Trans. Math. Software 27 340–349
Borsdorf R, Higham N J and Raydan M (2010) Computing a nearest correlation matrix with factor structure. SIAM J. Matrix Anal. Appl. 31(5) 2603–2622

Parameters

Compulsory Input Parameters

1:     g(ldg,n) – double array
ldg, the first dimension of the array, must satisfy the constraint ldgnldgn.
GG, the initial matrix.
2:     k – int64int32nag_int scalar
kk, the number of factors and columns of XX.
Constraint: 0 < kn0<kn.

Optional Input Parameters

1:     n – int64int32nag_int scalar
Default: The first dimension of the array g and the second dimension of the array g. (An error is raised if these dimensions are not equal.)
nn, the order of the matrix GG.
Constraint: n > 0n>0.
2:     errtol – double scalar
The termination tolerance for the projected gradient norm. See references for further details. If errtol0.0errtol0.0 then 0.010.01 is used. This is often a suitable default value.
Default: 0.00.0 
3:     maxit – int64int32nag_int scalar
Specifies the maximum number of iterations in the spectral projected gradient method.
If maxit0maxit0, 4000040000 is used.
Default: 00 

Input Parameters Omitted from the MATLAB Interface

ldg ldx

Output Parameters

1:     g(ldg,n) – double array
ldgnldgn.
A symmetric matrix (1/2)(G + GT)12(G+GT) with the diagonal elements set to unity.
2:     x(ldx,k) – double array
ldxnldxn.
Contains the matrix XX.
3:     iter – int64int32nag_int scalar
The number of steps taken in the spectral projected gradient method.
4:     feval – int64int32nag_int scalar
The number of evaluations min GXXT + diag(XXTI)FminG-XXT+diag(XXT-I)F.
5:     nrmpgd – double scalar
The norm of the projected gradient at the final iteration.
6:     ifail – int64int32nag_int scalar
ifail = 0ifail=0 unless the function detects an error (see [Error Indicators and Warnings]).

Error Indicators and Warnings

Errors or warnings detected by the function:
  ifail = 1ifail=1
Constraint: 0 < kn0<kn.
Constraint: ldgnldgn.
Constraint: ldxnldxn.
Constraint: n > 0n>0.
  ifail = 2ifail=2
Spectral gradient method fails to converge in __ iterations.
  ifail = 999ifail=-999
Dynamic memory allocation failed.

Accuracy

The returned accuracy is controlled by errtol and limited by machine precision.

Further Comments

Arrays are internally allocated by nag_correg_corrmat_nearest_kfactor (g02ae). The total size of these arrays is n × n + 4 × n × k + (nb + 3) × n + n + 50n×n+4×n×k+(nb+3)×n+n+50 double elements and 6 × n6×n integer elements. Here nbnb is the block size required for optimal performance by nag_lapack_dsytrd (f08fe) and nag_lapack_dormtr (f08fg) which are called internally. All allocated memory is freed before return of nag_correg_corrmat_nearest_kfactor (g02ae).
See nag_mv_factor (g03ca) for constructing the factor loading matrix from a known correlation matrix.

Example

function nag_correg_corrmat_nearest_kfactor_example
g = [2, -1,  0,  0;
    -1,  2, -1,  0;
     0, -1,  2, -1;
     0,  0, -1,  2];
k = int64(2);
% Calculate nearest correlation matrix
[g, x, iter, feval, nrmpgd, ifail] = nag_correg_corrmat_nearest_kfactor(g, k);

if (ifail == 0)
  fprintf('\n Factor Loading Matrix x:\n');
  disp(x);
  fprintf('\n Number of Newton steps taken:   %d\n', iter);
  fprintf(' Number of function evaluations: %d\n', feval);

  % Generate Nearest k factor correlation matrix

  fprintf('\n Nearest Correlation Matrix:\n');
  disp(x*transpose(x) + diag(diag(eye(4)-x*transpose(x))));
end
 

 Factor Loading Matrix x:
    0.7665   -0.6271
   -0.4250    0.9052
   -0.4250   -0.9052
    0.7665    0.6271


 Number of Newton steps taken:   5
 Number of function evaluations: 6

 Nearest Correlation Matrix:
    1.0000   -0.8935    0.2419    0.1943
   -0.8935    1.0000   -0.6388    0.2419
    0.2419   -0.6388    1.0000   -0.8935
    0.1943    0.2419   -0.8935    1.0000


function g02ae_example
g = [2, -1,  0,  0;
    -1,  2, -1,  0;
     0, -1,  2, -1;
     0,  0, -1,  2];
k = int64(2);
% Calculate nearest correlation matrix
[g, x, iter, feval, nrmpgd, ifail] = g02ae(g, k);

if (ifail == 0)
  fprintf('\n Factor Loading Matrix x:\n');
  disp(x);
  fprintf('\n Number of Newton steps taken:   %d\n', iter);
  fprintf(' Number of function evaluations: %d\n', feval);

  % Generate Nearest k factor correlation matrix

  fprintf('\n Nearest Correlation Matrix:\n');
  disp(x*transpose(x) + diag(diag(eye(4)-x*transpose(x))));
end
 

 Factor Loading Matrix x:
    0.7665   -0.6271
   -0.4250    0.9052
   -0.4250   -0.9052
    0.7665    0.6271


 Number of Newton steps taken:   5
 Number of function evaluations: 6

 Nearest Correlation Matrix:
    1.0000   -0.8935    0.2419    0.1943
   -0.8935    1.0000   -0.6388    0.2419
    0.2419   -0.6388    1.0000   -0.8935
    0.1943    0.2419   -0.8935    1.0000



PDF version (NAG web site, 64-bit version, 64-bit version)
Chapter Contents
Chapter Introduction
NAG Toolbox

© The Numerical Algorithms Group Ltd, Oxford, UK. 2009–2013