Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Build a block diagonal matrix from a known matrix

Tags:

r

matrix

diagonal

I want to build a block diagonal matrix (A) from a known matrix (B) by putting B in diagonal positions of A.

Let's say my B:

> matrix(c(1,3,4,5),nrow=2)
      [,1] [,2]
[1,]    1    4
[2,]    3    5 

I am looking for a function like this: function(B,3) (3 is just a random number) which returns matrix A like this:

[1,] 1 4 . . . .
[2,] 3 5 . . . .
[3,] . . 1 4 . .
[4,] . . 3 5 . .
[5,] . . . . 1 4
[6,] . . . . 3 5

Really appreciate any help

like image 844
Square9627 Avatar asked Dec 23 '17 18:12

Square9627


1 Answers

The desired matrix is the kronecker product of an identity matrix and B. No packages used.

B <- matrix(c(1,3,4,5),nrow=2)
diag(3) %x% B

giving:

     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    1    4    0    0    0    0
[2,]    3    5    0    0    0    0
[3,]    0    0    1    4    0    0
[4,]    0    0    3    5    0    0
[5,]    0    0    0    0    1    4
[6,]    0    0    0    0    3    5
like image 171
G. Grothendieck Avatar answered Sep 22 '22 05:09

G. Grothendieck