Как сделать массив диагональных матриц в R?

Question

Я знаю, как я могу сделать массив из N чисел р * р матриц:

m=array(x, c(p,p,N)) 

который x может быть вектор или скаляр. Я хочу создать массив диагональных матриц. Каждая матрица должна быть диагональной. Я попробовал несколько способов, но они не работают. Кто-нибудь знает об этом?

Answer 1

Если вы не возражаете, используя пакет, вы можете использовать abind:

library(abind) 
do.call(abind, c(rep(list(diag(3)), 2), along = 3)) 
#, , 1 
# 
#  [,1] [,2] [,3] 
#[1,] 1 0 0 
#[2,] 0 1 0 
#[3,] 0 0 1 
# 
#, , 2 
# 
#  [,1] [,2] [,3] 
#[1,] 1 0 0 
#[2,] 0 1 0 
#[3,] 0 0 1 

Answer 2

Это немного неуклюжим, но вы можете сделать это:

n = 3 
num_of_matrix = 2 

array(rep(c(diag(n)),num_of_matrix),c(n,n,num_of_matrix)) 
#, , 1 

#  [,1] [,2] [,3] 
#[1,] 1 0 0 
#[2,] 0 1 0 
#[3,] 0 0 1 

#, , 2 

#  [,1] [,2] [,3] 
#[1,] 1 0 0 
#[2,] 0 1 0 
#[3,] 0 0 1