r - 将列向量转换为扩展对角矩阵

Question

考虑以下列向量:

vec <- rbind(c(0.5),c(0.6))

我想将其转换为以下 4x4 对角矩阵:

0.5   0   0    0 
  0 0.6   0    0 
  0   0 0.5    0
  0   0   0  0.6

我知道我可以通过以下代码做到这一点:

dia <- diag(c(vec,vec))

但是如果我想把它转换成一个 1000x1000 的对角矩阵怎么办。那么上面的代码就是这么高效。也许我可以用 rep ，但我不完全确定该怎么做。我怎样才能更有效地做到这一点？

Answer 1

这是使用回收的另一种方式:

diag(c(vec), length(vec)*2)