版權聲明：本文為博主原創文章，遵循 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)) > MX Y X 1 0.337440 Y 2 0.1298763 1.2513904 -0.3824475 -0.451144

　

版權聲明：本文為CSDN博主「lili_wuwu」的原創文章，遵循 CC 4.0 BY-SA 版權協議，轉載請附上原文出處鏈接及本聲明。
原文鏈接：https://blog.csdn.net/lili_wuwu/article/details/51909229

轉載于:https://www.cnblogs.com/jiaxinwei/p/11503974.html

總結

以上是生活随笔為你收集整理的【转】R函数-diag()函数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。