通過刪除部分字符串來重命名樹提示標簽

Edit: clarification

我想在R中重命名我的nwk樹的tip.labels。當前標簽的示例是:

AB177299.1 Uncultured bacterium gene for 16S rRNA, clone: ODP1251B1.3

但我只是想要開頭的數字ID,所以AB177299.1。我的第一個傾向是在tree$tip.label上使用gsub(),但我的問題是每個數字ID和下面的文本對于每個提示都是不同的。然后,如果我能夠使用order()準確地替換每個ID,我創建了一個所有數字ID的字符向量,這些ID可以用作提示名稱。

tip.ids <- rownames(taxonomy)
tree$tip.label[order(tree$tip.label) %in%
                         tip.ids] <- tip.ids

這并沒有改變我樹上的任何東西,我不知所措。

為了再現性。taxonomy不包括在內,因為您只需要tip.ids

> dput(tip.ids)
c("AB109878.1", "AB109879.1", "AB109880.1", "AB109881.1", "AB109882.1", 
"AB109883.1", "AB109884.1", "AB109885.1")

> dput(tree)
structure(list(edge = structure(c(9L, 10L, 10L, 9L, 11L, 11L, 
12L, 13L, 13L, 12L, 14L, 14L, 15L, 15L, 10L, 1L, 2L, 11L, 3L, 
12L, 13L, 4L, 5L, 14L, 6L, 15L, 7L, 8L), dim = c(14L, 2L)), edge.length = c(0.0341921975, 
5e-09, 0.12821348, 0.000367458500000008, 0.027617765, 0.037677039, 
0.028633124, 0.014468092, 5e-09, 0.009763081, 0.078168769, 0.021640684, 
0.341568464, 0.092957415), Nnode = 7L, node.label = c("root", 
"0.917", "", "0.929", "0.921", "0.302", "0.692"), tip.label = c("'AB109881.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-4'", 
"'AB109880.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-3'", 
"'AB109883.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-6'", 
"'AB109879.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-2'", 
"'AB109884.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-7'", 
"'AB109878.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-1'", 
"'AB109882.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-5'", 
"'AB109885.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-8'"
)), class = "phylo", order = "cladewise")
? 最佳回答:

也許其中一個選項對你有用。

library(tidyverse, quietly = TRUE)                                                                                                                                                               

text <- "'AB109880.1 Uncultured archaeon gene for 16S rRNA, partial sequence, clone:pMLA-3'"

# Option 1
paste0("AB", parse_number(text))
# Option 2
str_split(text, " ")[[1]][1]
主站蜘蛛池模板: 无码国产伦一区二区三区视频| 国产精品毛片一区二区| 日韩在线观看一区二区三区| 怡红院美国分院一区二区| 一区二区三区四区在线播放 | 人妻免费一区二区三区最新| 国产福利电影一区二区三区久久老子无码午夜伦不 | 无码乱码av天堂一区二区| 综合激情区视频一区视频二区| 亚洲国产老鸭窝一区二区三区| 国产精品资源一区二区| AV鲁丝一区鲁丝二区鲁丝三区| 三级韩国一区久久二区综合| 久久久国产精品亚洲一区| 中文字幕在线观看一区二区三区| 精品国产乱子伦一区二区三区| 久久精品国产一区二区三| 亚洲国产精品成人一区| 日韩十八禁一区二区久久| 亚洲码欧美码一区二区三区| 亚洲一区二区三区免费观看| 乱精品一区字幕二区| 台湾无码一区二区| 极品人妻少妇一区二区三区| 无码国产精品一区二区免费3p | 日韩精品一区二区三区老鸭窝| 亚洲bt加勒比一区二区| 久久久久成人精品一区二区| 久久精品国产一区| 3d动漫精品成人一区二区三| 亚洲日本久久一区二区va| 亚洲AV无码一区二区三区牲色| 国产精品电影一区二区三区| 国产成人av一区二区三区不卡| 日本精品一区二区三本中文| av无码一区二区三区| 人妻夜夜爽天天爽爽一区| 无码少妇一区二区浪潮免费| jizz免费一区二区三区| 精品一区二区三区波多野结衣 | 国产亚洲日韩一区二区三区|