【R语言基础】如何提取矩阵的上三角矩阵

如何提取矩阵的上三角矩阵：

例如：已知矩阵，如何求除对角线值以外均值。

---

问题描述

提示：求矩阵最大值、最小值、均值；可以值提取上/下三角形矩阵；借助R包

实现方法

借助R包gdata,函数upper.tri 与函数 lower.tri

library(openxlsx)
library(gdata)

# 读取原始上下三角形矩阵
G=read.xlsx(xlsxFile ="G.xlsx",sheet=1,rowNames=T,colName=T)

# 提取上三角矩阵，即下三角矩阵值替换为空 / 0  / NA
G[!upper.tri(G, diag = TRUE)] <- " " 
G[!upper.tri(G, diag = TRUE)] <- 0
G[!upper.tri(G, diag = TRUE)] <- NA
# diag参数如果为TRUE，表示包含矩阵对角线值

# 同理生成下三角矩阵
G[!lower.tri(G, diag = TRUE)] <- " "

#也可以生成按行排列的结果 
G_byrow=upperTriangle(G,byrow = T)

# 求最大值、最小值
> signif(mean(G_byrow,na.rm = T),4) #忽略缺失值
[1] -0.009783
> signif(max(G_byrow,na.rm = T),4)
[1] 0.9994
> signif(min(G_byrow,na.rm = T),4)
[1] -0.2613

---

以上就是借助R包，如何提取矩阵的上/下三角矩阵内容啦，快去试一试吧。