如何提取矩阵的上三角矩阵:
例如:已知矩阵,如何求除对角线值以外均值。
问题描述
提示:求矩阵最大值、最小值、均值;可以值提取上/下三角形矩阵;借助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包,如何提取矩阵的上/下三角矩阵内容啦,快去试一试吧。