我有如下矩陣
prob.matx
:
prob.matx <- structure(c(0.182, 0.212, 0.364, 0.242, 0.242, 0.152, 0.394,
0.212, 0.364, 0.091, 0.273, 0.273, 0.333, 0.242, 0.364, 0.061,
0.273, 0.333, 0.212, 0.182, 0.303, 0.333, 0.182, 0.182),
dim = c(4L, 6L),
dimnames = list(c("A", "C", "G", "T"),
c("V1", "V2", "V3", "V4", "V5", "V6")))
看起来像这样:
V1 V2 V3 V4 V5 V6
A 0.182 0.242 0.364 0.333 0.273 0.303
C 0.212 0.152 0.091 0.242 0.333 0.333
G 0.364 0.394 0.273 0.364 0.212 0.182
T 0.242 0.212 0.273 0.061 0.182 0.182
还有一个向量DNA
:
DNA <- c("A", "C", "C", "C", "C", "A")
我想根据 DNA
向量中的元素和位置提取矩阵值,并返回一个向量。也就是说:
c(prob.matx[DNA[1], 1], prob.matx[DNA[2], 2], prob.matx[DNA[3], 3],
prob.matx[DNA[4], 4], prob.matx[DNA[5], 5], prob.matx[DNA[6], 6])
[1] 0.182 0.152 0.091 0.242 0.333 0.303
这似乎很简单,但我很难找到一个一步函数来完成这个任务,不使用apply
或for loop。
DNA
向量来说可能会占用过多内存。 - ThomasIsCoding