版權聲明：本文為博主原創文章，遵循 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
)
