如何将Data Frame转换为带空白/NA角的非对称矩阵
将长格式DataFrame转换为非对称矩阵
你可以用R的tidyr或reshape2包,把现有长格式的群体对比数据转成行对应pop_1、列对应pop_2的非对称矩阵,缺失单元格自动填充NA(需要空白可后续替换),具体操作如下:
方法一:使用tidyr包的pivot_wider
# 先安装并加载tidyr包(未安装过的话执行第一行) install.packages("tidyr") library(tidyr) # 长格式转宽格式 asym_matrix <- pivot_wider(table, names_from = pop_2, values_from = value) # 将pop_1设为行名,移除原pop_1列 rownames(asym_matrix) <- asym_matrix$pop_1 asym_matrix <- asym_matrix[, -1] # 查看最终矩阵 print(asym_matrix)
方法二:使用reshape2包的dcast
# 安装并加载reshape2包 install.packages("reshape2") library(reshape2) # 转换为宽格式 asym_matrix_dcast <- dcast(table, pop_1 ~ pop_2, value.var = "value") # 设置行名并移除pop_1列 rownames(asym_matrix_dcast) <- asym_matrix_dcast$pop_1 asym_matrix_dcast <- asym_matrix_dcast[, -1]
转换后的矩阵会保留所有pop_1与pop_2的组合,无对应值的位置用NA填充。如果需要把NA替换成空白,可执行以下代码(注意这会将矩阵转为字符型):
asym_matrix_char <- replace(asym_matrix, is.na(asym_matrix), "")
内容的提问来源于stack exchange,提问作者trilisser
相关产品推荐
相关产品推荐

