如何摆脱“”对象名称并循环工作？

Question

我不知道如何精确定义这个问题；但我想不通

    genx <- list(scaffold_1 = c("AAATTTTTATAT"),scaffold_2 = c("AAATTTTTATAT"),
    scaffold_3 = c("AAATTTTTATAT"),scaffold_4 = c("AAATTTTTATAT"),
    scaffold_5 = c("AAATTTTTATATA"),scaffold_6 = c("AAATTTTTATAT"),
scaffold_7 = c("AAATTTTTATAT"),scaffold_8 = c("AAATTTTTATATA"))


TATA = "TATA"   
      myobs <- paste("genx$scaffold_", 1:8, sep = "")

我想将以下函数应用于 myobs 元素（是对象）的每个元素：

source("http://www.bioconductor.org/biocLite.R")

biocLite("Biostrings")
require((Biostrings)

countPattern (TATA, genx$scaffold_1, max.mismatch = 1)
[1] 3  

当我使用以下内容时：

countPattern (TATA, myobs[1], max.mismatch = 1)

不工作，我认为解释为：

countPattern (TATA, "genx$scaffold_1", max.mismatch = 1)
 [1] 0 

这和上面的不一样。如何摆脱“”并创建一个循环来执行这项工作，感谢您的建议：

Answer 1

您可以直接使用列表并执行sapply，它作用于列表的每个元素。这是一些示例代码

sapply(genx, countPattern, pattern = TATA, max.mismatch = 1)

Answer 2

你可以使用get函数

a  <- 12
get("a") # returns 12

Answer 3

这个呢：

sapply(genx, function (x) { countPattern(TATA, x, max.mismatch = 1) })

Answer 4

有问题：

get("genx$scaffold_1")

因为 R 认为您正在寻找具有该全名的对象，而不是 'genx' 的 'scaffold_1' 组件。

应该起作用的是：

eval(parse(text="genx$scaffold_1"))

但请看：

fortune(106)

更多麻烦点可以在“The R Inferno”中找到：http://www.burns-stat.com/pages/Tutor/R_inferno.pdf

Answer 5

我建议你将genx$scaffold定义为向量，然后你就可以轻松使用apply函数了：

genx = data.frame(scaffold = c(1, 2, 3, 4)) # genx$scaffold is a vector

现在您可以轻松地在每个 genx$scaffold 的元素上运行 countPattern() 函数

apply(as.matrix(genx$scaffold), 1, 
    function (x) { countPattern(TATA, x, max.mismatch = 1) })

您可以轻松扩展解决方案，以便 genx$scaffold 也可以是矩阵等。