получить 9меры из сиквенса пептида по заданной позиции. В моём коде мне вылетает warning, и в колонку записывается только первое число. Как мне сохранить весь лист из 9 элементов?
l <- 9
data.missence.prot.pos <- mydf %>%
mutate(pos = gsub("\\D", "", c(mutation.aa))) %>%
mutate(peps = substring(peptide,
(as.numeric(pos)-l+1):as.numeric(pos),
as.numeric(pos):(as.numeric(pos)+l-1)))
Тут хорошо бы пример исходных данных. А то трудно код просто так читать.
Обсуждают сегодня