Question

我的Amino_acid突变数据有一个问题。例如，

p.K303R   
p.?
p.R1450*

我希望这样outopt

AA_mutation wt_residue  position    mt_residue
p.K303R       K          303           R
p.?         
p.R1450*      R           1450         *

我想删除“p。”，“？”从数据中又分成三个新变量。我已经在excel上管理它但不在R平台上管理它。有人能帮我请R. 亲切的问候 Fakhrul

Answer 1

使用dplyr和extract（有some reading here），我们可以使用：

library(dplyr)
df <- data.frame(AA_mutation = c("p.K303R", "p.?", "p.R1450*"))
df <- df %>%
  extract(AA_mutation, 
          into = c("wt_residue", "position", "mt_residue"), 
          regex = "p\\.([A-Z])?(\\d+)?([A-Z*])?",
          remove = FALSE)
df

屈服

  AA_mutation wt_residue position mt_residue
1     p.K303R          K      303          R
2         p.?       <NA>     <NA>       <NA>
3    p.R1450*          R     1450          *

Answer 2

以下是base R

的选项

lst <- regmatches(df1[[1]], gregexpr("([a-z]+)|([A-Z*]+)|[0-9]+", df1[[1]], perl = TRUE))
res <- do.call(rbind.data.frame, lapply(lst, `length<-`, max(lengths(lst))))
names(res) <- c("AA_mutation", "wt_residue",  "position",    "mt_residue")
res
#  AA_mutation wt_residue position mt_residue
#1           p          K      303          R
#2           p       <NA>     <NA>       <NA>
#3           p          R     1450          *

在R中减少并分成几列

2 个答案: