matrixcalc (version 1.0-3)

K.matrix: K Matrix

Description

This function returns a square matrix of order p = r * c that, for an r by c matrix A, transforms vec(A) to vec(A') where prime denotes transpose.

Usage

K.matrix(r, c = r)

Arguments

r

a positive integer row dimension

c

a positive integer column dimension

Value

An order \(\left( {r\;c} \right)\) matrix.

Details

The \(r \times c\) matrices \({\bf{H}}{}_{i,j}\) constructed by the function H.matrices are combined using direct product to generate the commutation product with the formula \({{\bf{K}}_{r,c}} = \sum\limits_{i = 1}^r {\sum\limits_{j = 1}^c {\left( {{{\bf{H}}_{i,j}} \otimes {{{\bf{H'}}}_{i,j}}} \right)} }\)

References

Magnus, J. R. and H. Neudecker (1979). The commutation matrix: some properties and applications, The Annals of Statistics, 7(2), 381-394.

Magnus, J. R. and H. Neudecker (1999) Matrix Differential Calculus with Applications in Statistics and Econometrics, Second Edition, John Wiley.

See Also

H.matrices

Examples

Run this code
# NOT RUN {
K <- K.matrix( 3, 4 )
A <- matrix( seq( 1, 12, 1 ), nrow=3, byrow=TRUE )
vecA <- vec( A )
vecAt <- vec( t( A ) )
y <- K %*% vecA
print( y )
print( vecAt )
# }

Run the code above in your browser using DataLab