【转】R函数-diag()函数

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lili_wuwu/article/details/51909229

diag()提取或替换一个矩阵的对角线,或构造一个对角矩阵。
语法：
diag(x = 1, nrow, ncol)
diag(x) <- value
解析：
x：一个矩阵,向量或一维数组,或不填写。
nrow, ncol：可选 行列。
value ：对角线的值，可以是一个值或一个向量。

例子

> x <- matrix(1:16,nrow=4,ncol=4)
> x
     [,1] [,2] [,3] [,4]
[1,]    1    5    9   13
[2,]    2    6   10   14
[3,]    3    7   11   15
[4,]    4    8   12   16
> diag(x)
[1]  1  6 11 16

#展示x的对角线数据

> diag(diag(x))
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    6    0    0
[3,]    0    0   11    0
[4,]    0    0    0   16

#由x的对角线数据生成对角矩阵

> diag(4)
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    1    0    0
[3,]    0    0    1    0
[4,]    0    0    0    1

#对角线为4个元素每个元素默认为1的对角矩阵

> diag(10,3,4)
     [,1] [,2] [,3] [,4]
[1,]   10    0    0    0
[2,]    0   10    0    0
[3,]    0    0   10    0

#对角线数据位10,三行四列的矩阵

> a<-diag(4)+1
> a
     [,1] [,2] [,3] [,4]
[1,]    2    1    1    1
[2,]    1    2    1    1
[3,]    1    1    2    1
[4,]    1    1    1    2

　　

看几个较复杂的例子
例1

> all(diag(1:3) == {m <- matrix(0,3,3); diag(m) <- 1:3; m})
[1] TRUE

　

例2

> diag(var(M <- cbind(X = 1:5, Y = stats::rnorm(5))))
       X        Y
2.500000 0.472455

　　

例3

> rownames(M) <- c(colnames(M),rep("",3))
> M
  X         Y
X 1  0.337440
Y 2  0.129876
  3  1.251390
  4 -0.382447
  5 -0.451144

　

版权声明：本文为CSDN博主「lili_wuwu」的原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/lili_wuwu/article/details/51909229