Question

我在一个名为x的字符串向量中有一组核苷酸序列。

我想检查x中是否存在一些（比如10个）图案。我想生成一个数据框或表，其中行是X中的序列，列是模式/主题位于向量sdseqs中。

sdframe <- data.frame
sdseqs = c("AGGAG.+ATG", 
"AGAAG.+ATG","AAAGG.+ATG","GGAGG.+ATG","GAAGA.+ATG",
"GGAGA.+ATG","AAGGT.+ATG","AGGAA.+ATG","AAGGA.+ATG","GTGGA.+ATG")
for (i in 1:10) {
sdframe <- cbind(sdframe,(grepl(sdseqs[i], x)))
}

此代码工作正常，但数据框的第一列将为空，带有问号。其他列填充了true和false - 这就是我想要的。

我试图在开头的循环外定义一个空数据框。我是R的新手，我来自Perl。这就是我在Perl中经常做的事情：你定义了在外部循环中使用的变量。我怎么能在R中这样做？

另外，一个可行的选择是从我的数据框中删除第一列，但这对我来说似乎并不那么简单。

感谢任何帮助。

输出我现在使用我的代码：

  sdframe                                                            
[1,] ?       TRUE  FALSE TRUE  TRUE  FALSE TRUE  TRUE  TRUE  TRUE  FALSE
[2,] ?       FALSE TRUE  FALSE FALSE FALSE FALSE FALSE FALSE TRUE  TRUE 
[3,] ?       FALSE FALSE TRUE  FALSE TRUE  FALSE TRUE  TRUE  TRUE  TRUE 
[4,] ?       TRUE  FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[5,] ?       FALSE TRUE  FALSE FALSE TRUE  FALSE FALSE FALSE FALSE FALSE
[6,] ?       FALSE FALSE FALSE TRUE  FALSE FALSE FALSE TRUE  FALSE TRUE 
[7,] ?       FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE  FALSE FALSE
[8,] ?       FALSE FALSE TRUE  FALSE FALSE TRUE  FALSE FALSE TRUE  FALSE
[9,] ?       FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[10,] ?       FALSE FALSE FALSE FALSE TRUE  FALSE FALSE FALSE FALSE FALSE
[11,] ?       FALSE FALSE TRUE  FALSE FALSE FALSE FALSE FALSE FALSE FALSE

我想要相同但没有第一列？注意我的x有11个序列，我检查的主题是列（10列，11列第一个？）

Answer 1

常见的R解决方案将使用apply系列中的函数在向量上应用函数。

sdseqs = c(
  "AGGAG.+ATG",
  "AGAAG.+ATG",
  "AAAGG.+ATG",
  "GGAGG.+ATG",
  "GAAGA.+ATG",
  "GGAGA.+ATG",
  "AAGGT.+ATG",
  "AGGAA.+ATG",
  "AAGGA.+ATG",
  "GTGGA.+ATG"
)

sdframe <- sapply(sdseqs, function(one.motif) {
  grepl(one.motif, x = x)
})

sdframe

     AGGAG.+ATG AGAAG.+ATG AAAGG.+ATG GGAGG.+ATG GAAGA.+ATG GGAGA.+ATG AAGGT.+ATG AGGAA.+ATG AAGGA.+ATG GTGGA.+ATG
[1,]      FALSE       TRUE      FALSE      FALSE       TRUE       TRUE       TRUE      FALSE       TRUE      FALSE
[2,]      FALSE       TRUE      FALSE      FALSE       TRUE       TRUE       TRUE      FALSE       TRUE      FALSE
[3,]      FALSE       TRUE      FALSE      FALSE       TRUE       TRUE       TRUE      FALSE       TRUE      FALSE

sdframe.t <- t(sdframe)

sdframe.t

            [,1]  [,2]  [,3]
AGGAG.+ATG FALSE FALSE FALSE
AGAAG.+ATG  TRUE  TRUE  TRUE
AAAGG.+ATG FALSE FALSE FALSE
GGAGG.+ATG FALSE FALSE FALSE
GAAGA.+ATG  TRUE  TRUE  TRUE
GGAGA.+ATG  TRUE  TRUE  TRUE
AAGGT.+ATG  TRUE  TRUE  TRUE
AGGAA.+ATG FALSE FALSE FALSE
AAGGA.+ATG  TRUE  TRUE  TRUE
GTGGA.+ATG FALSE FALSE FALSE

Answer 2

事实上，在第一行中您不会创建data.frame。所以你的输出是list。

而不是cbind，您需要rbind来添加行：

sdframe <- data.frame()
sdseqs = c("AGGAG.+ATG", 
       "AGAAG.+ATG","AAAGG.+ATG","GGAGG.+ATG","GAAGA.+ATG",
       "GGAGA.+ATG","AAGGT.+ATG","AGGAA.+ATG","AAGGA.+ATG","GTGGA.+ATG")
for (i in 1:10) {
sdframe <- rbind(sdframe,(grepl(sdseqs[i], x)))
}

在数据帧输出的序列向量中搜索图案矢量

2 个答案: