2

我想使用 R 和 sub 提取符号周围的字符。我尝试了很多正则表达式，但没有得到我想要的。

我的载体：

c("G>GA", "T>A", "G>A", "G>A", "A>T", "CT>C", "T>C", "T>C", "A>T", "T>C", "T>A", "A>G", "CCGCCGCGGCCGCCGTCTTCCACCAACAACATGGCGGA>C", "C>T", "T>A", "T>C", "T>G", "G>C", "T>G", "T>A", "G>A")

我只需要在 . 之前和之后的一个字符>。

我最好的尝试是：

sub("(.*?)>", ">", aa, perl = TRUE)
4

3 回答 3

9

您需要在正则表达式中使用捕获组：

vec <- c("G>GA", "T>A", "G>A", "G>A", "A>T", "CT>C", "T>C", "T>C", "A>T", "T>C", "T>A", "A>G", "CCGCCGCGGCCGCCGTCTTCCACCAACAACATGGCGGA>C", "C>T", "T>A", "T>C", "T>G", "G>C", "T>G", "T>A", "G>A")
> sub(".*(.)>(.).*","\\1\\2",vec)
 [1] "GG" "TA" "GA" "GA" "AT" "TC" "TC" "TC" "AT" "TC" "TA" "AG" "AC" "CT" "TA"
[16] "TC" "TG" "GC" "TG" "TA" "GA"

换句话说，正则表达式匹配任何零次或多次，.*然后捕获下一个字符，(.)然后匹配大于号>，然后捕获下一个字符(.)，然后在最后匹配任何零次或多次.*。用捕获的两个字符替换所有这些\\1\\2。

于 2014-01-13T15:59:37.153 回答
5

提供一个可重现的例子

> x = c("A>G", "AT>GC")

找到您感兴趣的符号的索引（使用fixed=TRUE，因为您实际上并不是在寻找正则表达式）。

> i = regexpr(">", x, fixed=TRUE)

然后提取前面和/或后面的字符

> substr(x, i-1, i-1)
[1] "A" "T"
> substr(x, i+1, i+1)
[1] "G" "G"

或获取序列

> substr(x, i-1, i+1)
[1] "A>G" "T>G"

也许您的可重现示例包括边缘情况

> x = c("A>G", "AT>GC", "", ">G", "A>", ">", NA)

然后需要更多处理？

于 2014-01-13T16:07:24.320 回答
0

看起来您正在尝试获取参考和替代等位基因？只寻找一个字符表明您只对 SNP 感兴趣？您可以使用 strsplit 生成 ref 和 alt 等位基因的数据框。

test <- c("G>GA", "T>A", "G>A", "G>A", "A>T", "CT>C", "T>C", "T>C", "A>T", "T>C", "T>A", "A>G", "CCGCCGCGGCCGCCGTCTTCCACCAACAACATGGCGGA>C", "C>T", "T>A", "T>C", "T>G", "G>C", "T>G", "T>A", "G>A")
Alleles <- data.frame(t(data.frame(sapply(test, function(x)   strsplit(x,split=">")))),row.names=NULL,stringsAsFactors=F)
colnames(Alleles) <- c("Ref","Alt")
Alleles$bases <- apply(Alleles,1,function(x) sum(length(unlist(strsplit(x[1],split=""))),length(unlist(strsplit(x[2],split="")))))
SNPs <- Alleles[Alleles$bases == 2,]

仅在替换 (>) 的任一侧取一个碱基就会给你错误的遗传信息。变体“CCGCCGCGGCCGCCGTCTTCCACCAACAACATGGCGGA>C”将简化为“A>C”——它看起来像一个简单的 SNP，但与删除最后 38 个碱基“CGCCGCGGCCGCCGTCTTCCACCAACAACATGGCGGA>-”相同。

这就是你所追求的吗？

于 2014-01-14T02:35:02.213 回答