获取每行的字符串的特定部分

时间:2017-09-11 20:29:14

标签: r regex bioinformatics

我有一个数据,我想采取它的具体部分

    DoseResponse_Curves/drCurve_AAATT.pdf
    DoseResponse_Curves/drCurve_AGMK1.pdf
    DoseResponse_Curves/drCurve_AGU.pdf
    DoseResponse_Curves/drCurve_ALH1L2.pdf
    DoseResponse_Curves/drCurve_ALKB1.pdf
    DoseResponse_Curves/drCurve_AS2.pdf
    DoseResponse_Curves/drCurve_ANK1.pdf
    DoseResponse_Curves/drCurve_ANKRD54.pdf

我只想接受第二个_之前和之后的任何事情。这意味着输出看起来像这样

AAATT 
AGMK1
AGU
ALH1L2
ALKB1
AS2
ANK1
ANKRD54

注意:因为我们正在处理基因名称,所以它们可能包含c(".", "-")等字符。

2 个答案:

答案 0 :(得分:2)

您可以使用sub和正则表达式执行此操作。

Files = c(
    'DoseResponse_Curves/drCurve_AAATT.pdf',
    'DoseResponse_Curves/drCurve_AGMK1.pdf',
    'DoseResponse_Curves/drCurve_AGU.pdf',
    'DoseResponse_Curves/drCurve_ALH1L2.pdf',
    'DoseResponse_Curves/drCurve_ALKB1.pdf',
    'DoseResponse_Curves/drCurve_AS2.pdf',
    'DoseResponse_Curves/drCurve_ANK1.pdf',
    'DoseResponse_Curves/drCurve_ANKRD54.pdf')

sub(".*?_.*?_(.*?)\\..*", "\\1", Files)
[1] "AAATT"   "AGMK1"   "AGU"     "ALH1L2"  "ALKB1"   "AS2"     "ANK1"   
[8] "ANKRD54"

答案 1 :(得分:0)

有很多方法:

# Example data with gene names with dots and dashes
Files = c('DoseResponse_Curves/drCurve_ALKB1.pdf',
          'DoseResponse_Curves/drCurve_BAC05914.1.pdf',
          'DoseResponse_Curves/drCurve_ALDH1L1-AS1.pdf',
          'DoseResponse_Curves/drCurve_AL953854.2-002.pdf')

# as parts of path are all same, we can sub with "":
gsub("DoseResponse_Curves/drCurve_|.pdf", "", Files)
# [1] "ALKB1"          "BAC05914.1"     "ALDH1L1-AS1"    "AL953854.2-002"

# Or, as we are working with path and filenames maybe:
gsub("drCurve_", "", tools::file_path_sans_ext(basename(Files)))
# [1] "ALKB1"          "BAC05914.1"     "ALDH1L1-AS1"    "AL953854.2-002"


# @G5W answer doesn't handle extra dots in gene names
sub(".*?_.*?_(.*?)\\..*", "\\1", Files)
# [1] "ALKB1"       "BAC05914"    "ALDH1L1-AS1" "AL953854"