R 列名和行名匹配时如何基于矩阵值获取百分比值
我有一个矩阵如下R 列名和行名匹配时如何基于矩阵值获取百分比值,r,R,我有一个矩阵如下 structure(c(4407, 0, 53, 46, 50, 0, 4407, 0, 0, 0, 53, 0, 742, 408, 411, 46, 0, 408, 688, 464, 50, 0, 411, 464, 672, 40, 0, 447, 359, 382, 72, 0, 439, 342, 356, 35, 0, 350, 321, 344, 45, 0, 365, 341, 364, 42, 0, 498, 372, 376, 102, 0,
structure(c(4407, 0, 53, 46, 50, 0, 4407, 0, 0, 0, 53, 0, 742,
408, 411, 46, 0, 408, 688, 464, 50, 0, 411, 464, 672, 40, 0,
447, 359, 382, 72, 0, 439, 342, 356, 35, 0, 350, 321, 344, 45,
0, 365, 341, 364, 42, 0, 498, 372, 376, 102, 0, 304, 333, 340,
52, 0, 396, 354, 343, 52, 0, 377, 371, 374, 80, 0, 458, 374,
385, 54, 0, 509, 367, 365, 52, 0, 416, 322, 352, 47, 0, 505,
368, 384, 85, 0, 338, 358, 374, 45, 0, 396, 328, 330, 45, 0,
351, 381, 358, 53, 0, 362, 293, 301, 52, 0, 396, 340, 346, 55,
0, 385, 343, 340, 52, 0, 493, 406, 406, 49, 0, 516, 384, 416,
29, 0, 344, 330, 335, 32, 0, 326, 299, 314, 50, 0, 520, 372,
376, 33, 0, 348, 336, 340, 33, 0, 380, 329, 358, 31, 0, 354,
337, 349, 33, 0, 320, 301, 336, 33, 0, 341, 316, 347, 52, 0,
340, 316, 329, 48, 0, 522, 388, 389, 56, 0, 300, 362, 334, 57,
0, 303, 327, 301, 44, 0, 487, 401, 400, 51, 0, 439, 362, 373,
38, 0, 361, 336, 353, 38, 0, 328, 315, 337, 52, 0, 561, 403,
402, 44, 0, 459, 386, 393, 37, 0, 381, 340, 337, 37, 0, 405,
344, 343, 27, 0, 335, 346, 347, 31, 0, 370, 370, 385, 31, 0,
369, 340, 348, 54, 0, 326, 305, 335, 39, 0, 368, 343, 357, 38,
0, 379, 347, 367, 42, 0, 527, 389, 393, 56, 0, 382, 504, 471,
56, 0, 382, 504, 471, 51, 0, 459, 379, 394, 41, 0, 507, 393,
381, 68, 0, 475, 396, 398, 30, 0, 386, 362, 371, 34, 0, 378,
363, 379, 24, 0, 341, 317, 350, 35, 0, 378, 348, 364, 34, 0,
334, 327, 345, 28, 0, 323, 311, 322, 47, 0, 486, 365, 378, 52,
0, 479, 395, 414, 40, 0, 391, 327, 350, 49, 0, 363, 344, 370,
22, 0, 342, 329, 338, 51, 0, 381, 463, 440, 46, 0, 389, 460,
476, 46, 0, 487, 371, 385, 35, 0, 331, 307, 330, 47, 0, 522,
388, 405, 39, 0, 359, 352, 365, 42, 0, 369, 342, 356, 27, 0,
364, 352, 364, 38, 0, 377, 327, 342, 29, 0, 369, 346, 380, 37,
0, 382, 330, 355, 29, 0, 291, 245, 275, 30, 0, 354, 332, 358,
31, 0, 338, 310, 326, 37, 0, 347, 325, 356, 46, 0, 505, 385,
377, 44, 0, 504, 411, 430, 34, 0, 315, 310, 333, 37, 0, 392,
367, 373, 60, 0, 339, 324, 352, 48, 0, 475, 385, 402, 48, 0,
376, 342, 314, 1, 0, 0, 0, 0), .Dim = c(5L, 91L), .Dimnames = list(
c("chr", "leftPos", "OC_AH_026C.res", "OC_AH_026C.1.res",
"OC_AH_026C.2.res"), c("chr", "leftPos", "OC_AH_026C.res",
"OC_AH_026C.1.res", "OC_AH_026C.2.res", "OC_AH_026T.res",
"OC_AH_058T.res", "OC_AH_084C.res", "OC_AH_084T.res", "OC_AH_086C.res",
"OC_AH_086C.1.res", "OC_AH_086C.2.res", "OC_AH_086C.3.res",
"OC_AH_086T.res", "OC_AH_088C.res", "OC_AH_088T.res", "OC_AH_094C.res",
"OC_AH_094C.1.res", "OC_AH_094C.2.res", "OC_AH_094C.3.res",
"OC_AH_094C.4.res", "OC_AH_094C.5.res", "OC_AH_094C.6.res",
"OC_AH_094C.7.res", "OC_AH_094T.res", "OC_AH_096C.res", "OC_AH_096T.res",
"OC_AH_100C.res", "OC_AH_100C.1.res", "OC_AH_100T.res", "OC_AH_127C.res",
"OC_AH_127T.res", "OC_AH_133C.res", "OC_AH_133T.res", "OC_ED_008C.res",
"OC_ED_008C.1.res", "OC_ED_008C.2.res", "OC_ED_008C.3.res",
"OC_ED_008T.res", "OC_ED_016C.res", "OC_ED_016T.res", "OC_ED_031C.res",
"OC_ED_031T.res", "OC_ED_036C.res", "OC_ED_036T.res", "OC_GS_001C.res",
"OC_GS_001T.res", "OC_QE_062C.res", "OC_QE_062T.res", "OC_RS_010C.res",
"OC_RS_010T.res", "OC_RS_027C.res", "OC_RS_027C.1.res", "OC_RS_027C.2.res",
"OC_RS_027T.res", "OC_SH_051C.res", "OC_SH_051T.res", "OC_ST_014C.res",
"OC_ST_014C.1.res", "OC_ST_014T.res", "OC_ST_016T.res", "OC_ST_020C.res",
"OC_ST_020T.res", "OC_ST_024C.res", "OC_ST_024T.res", "OC_ST_033C.res",
"OC_ST_033T.res", "OC_ST_034C.res", "OC_ST_034C.1.res", "OC_ST_034C.2.res",
"OC_ST_034T.res", "OC_ST_035C.res", "OC_ST_035T.res", "OC_ST_036C.res",
"OC_ST_036T.res", "OC_ST_037T.res", "OC_ST_040C.res", "OC_ST_040T.res",
"OC_WG_001T.res", "OC_WG_002C.res", "OC_WG_002T.res", "OC_WG_005C.res",
"OC_WG_005T.res", "OC_WG_006C.res", "OC_WG_006T.res", "OC_WG_009T.res",
"OC_WG_019C.res", "OC_WG_019T.res", "Means.res", "sd.res",
"count")))
NewMatrix <- t(sapply(row.names(res_matrix),
function(x) substr(x, 1, 6) == substr(colnames(res_matrix), 1, 6)))
NewMatrix您想要这样的东西吗：
library(dplyr)
library(tidyr)
long_data =
data %>%
as.data.frame %>%
mutate(rowname = rownames(.) ) %>%
gather(variable, value, OC_AH_026C.res:OC_WG_019T.res)
matches =
long_data %>%
filter(rowname == variable) %>%
select(rowname, control_value = value)
long_result =
long_data %>%
left_join(matches) %>%
mutate(ratio = value / control_value)
wide_result =
long_result %>%
select(-value, -control_value) %>%
spread(variable, ratio)
我在函数_list[[I]]（value）中遇到错误：找不到函数“add_rownames”您有最新版本的dplyr吗？我的版本是0.3.0.2请参阅编辑，这样可以避免使用add_rownames函数谢谢，但这是一个非常通用的答案。我不明白这是如何给出百分比或事实上匹配的。