如何将R中矩阵的对角线元素转换为缺失值？

首先我们需要了解的是，对角元素仅在具有方矩阵的情况下才有用，否则设置对角元素将毫无意义，这几乎是所有数学家都知道的，但是有些新生可能会感到困惑，因为我们可以在对角元素中创建对角元素。非正方形矩阵，不应称为对角线。在R中，我们可以使用diag函数将矩阵的对角元素设置为缺失值/ NA。

例1

> M1<-matrix(1:16,nrow=4)
> M1

输出结果

  [,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(M1)<-NA
> M1

输出结果

    [,1] [,2] [,3] [,4]
[1,] NA    5    9   13
[2,]  2   NA   10   14
[3,]  3    7   NA   15
[4,]  4    8   12   NA

例2

> M2<-matrix(rpois(100,10),nrow=10)
> M2

输出结果

    [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 13   18    9    5    15   7   14    7    9    14
[2,] 11   10    12   6    11   13  11   12    6    11
[3,] 12   13     6   6    11   18   9   12    9     6
[4,] 13 14 11 17 17 6 5 10 8 11
[5,] 10 9 15 11 10 14 6 5 5 15
[6,] 8 15 6 10 8 8 10 7 11 7
[7,] 10 12 8 9 12 15 19 9 10 15
[8,] 15 9 4 13 4 13 10 9 11 11
[9,] 7 6 11 12 3 8 12 8 8 11
[10,] 12 13 13 9 11 11 6 6 7 10

示例

> diag(M2)<-NA
> M2

输出结果

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] NA 18 9 5 15 7 14 7 9 14
[2,] 11 NA 12 6 11 13 11 12 6 11
[3,] 12 13 NA 6 11 18 9 12 9 6
[4,] 13 14 11 NA 17 6 5 10 8 11
[5,] 10 9 15 11 NA 14 6 5 5 15
[6,] 8 15 6 10 8 NA 10 7 11 7
[7,] 10 12 8 9 12 15 NA 9 10 15
[8,] 15 9 4 13 4 13 10 NA 11 11
[9,] 7 6 11 12 3 8 12 8 NA 11
[10,] 12 13 13 9 11 11 6 6 7 NA

例子3

> M3<-matrix(rpois(25,3),nrow=5)
> M3

输出结果

[,1] [,2] [,3] [,4] [,5]
[1,] 2 3 4 2 4
[2,] 2 2 3 5 3
[3,] 3 7 6 1 5
[4,] 1 1 4 1 4
[5,] 1 2 0 1 5

示例

> diag(M3)<-NA
> M3

输出结果

[,1] [,2] [,3] [,4] [,5]
[1,] NA 3 4 2 4
[2,] 2 NA 3 5 3
[3,] 3 7 NA 1 5
[4,] 1 1 4 NA 4
[5,] 1 2 0 1 NA

例子4

> M4<-matrix(rnorm(36,5,2.5),nrow=6)
> M4

输出结果

[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 7.449650 2.977026 3.631003 1.588073 6.4641515 4.9833353
[2,] 4.326318 4.327728 6.790520 2.367960 0.9577471 5.4084005
[3,] 3.143312 6.766010 5.364526 2.544719 7.9256179 -0.8116725
[4,] 5.983452 5.340530 2.530320 6.830628 4.0030207 7.1645111
[5,] 4.073458 1.795408 -2.391053 6.320859 4.5389839 1.0296674
[6,] 3.739790 6.036844 6.171213 6.901320 5.5595449 4.9644731

示例

> diag(M4)<-NA
> M4

输出结果

[,1] [,2] [,3] [,4] [,5] [,6]
[1,] NA 2.977026 3.631003 1.588073 6.4641515 4.9833353
[2,] 4.326318 NA 6.790520 2.367960 0.9577471 5.4084005
[3,] 3.143312 6.766010 NA 2.544719 7.9256179 -0.8116725
[4,] 5.983452 5.340530 2.530320 NA 4.0030207 7.1645111
[5,] 4.073458 1.795408 -2.391053 6.320859 NA 1.0296674
[6,] 3.739790 6.036844 6.171213 6.901320 5.5595449 NA

范例5

> M5<-matrix(runif(25,2,5),nrow=5)
> M5

输出结果

[,1] [,2] [,3] [,4] [,5]
[1,] 2.753600 3.017767 2.923270 2.532730 3.425423
[2,] 2.742611 2.398763 4.906797 4.273442 3.416288
[3,] 2.744358 4.227648 3.690739 3.508798 3.423980
[4,] 3.124375 2.363392 3.336010 4.770364 3.294046
[5,] 3.468065 3.424600 2.667044 2.623021 4.055215

示例

> diag(M5)<-NA
> M5

输出结果

   [,1] [,2] [,3] [,4] [,5]
[1,] NA 3.017767 2.923270 2.532730 3.425423
[2,] 2.742611 NA 4.906797 4.273442 3.416288
[3,] 2.744358 4.227648 NA 3.508798 3.423980
[4,] 3.124375 2.363392 3.336010 NA 3.294046
[5,] 3.468065 3.424600 2.667044 2.623021 NA