Question

我正在努力使用翻译功能：

我有一个序列矩阵，我无法弄清楚为什么翻译功能不起作用。

这是我的剧本：

 head(myseq)
     [,1]                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             
[1,] "CAATTAAGATGCAGTTACTTCGCTGTTTTTCAATATTTTCTGTTATTGCTAGCGTTTTAGCAGCTGGAGAAGTCCAGTTGGTTGAATCCGGTGGCGGGTTGGTTCAACCTGGCCGTTCCCTGCGCCTTTCTTGTGCCGCTTCAGGCTTTACCTTCGACGATTATGCCATGCATTGGGTACGCCAAGCTCCGGGCAAAGGTCTGGAATGGGTGAGTGCCATTACCTGGAATAGCGGCCACATTGACTATGCGGATAGCGTCGAAGGGCGCTTTACCATTAGCCGTGATAACGCGAAGAACAGCCTCTATCTGGACATGAACAGCTTACGTGCGGAAGATACCGCAGTGTACTATTGCGCTAAGGTCAGCTATCTGAGTACTGCCAGTTCCCTGGATTATTGGGGTCAAGGAACCCTGGTGACCGTGTCAAGCGGCGGAGGCGGTTCTGGTGGTGGAGGCAGTBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB"
[2,] "CAATTAAGATGCAGTTACTTCGCTGTTTTTCAATATTTTCTGTTATTGCTAGCGTTTTAGCAGCTGGAGAAGTCCAGTTGGTTGAATCCGGTGGCGGGTTGGTTCAACCTGGCCGTTCCCTGCGCCTTTCTTGTGCCGCTTCAGGCTTTACCTTCGACGATTATGCCATGCATTGGGTACGCCAAGCTCCGGGCAAAGGTCTGGAATGGGTGTCTGCGGCTAGTTGGAGTGGTGGCCACATCGATTATGCCGATAGCGTCGAAGGGCGCTTTACCATTAGCCGTGATAACGCGAAGAACAGCCTCTATCTGGACATGAACAGCTTACGTGCGGAAGATACCGCAGTGTACTATTGCGCTAAGGTCAGCTATCTGAGTACTGCCAGTTCCCTGGATTATTGGGGTCAAGGAACCCTGGTGACCGTGTCAAGCGGCGGAGGCGGTTCTGGTGGTGGAGGCAGTBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB"

translate(myseq[1,],frame=2,sens=F)

Error in s2n(seq, levels = s2c("tcag")) :sequence is not a vector of chars

我检查了我的参数：

class(myseq[1,])

[1] "character"

is.vector(myseq[1,])

[1] TRUE

所以我认为我错过了一些东西，但我无法弄清楚是什么。

Answer 1

<强>解决方案：

如果您使用的是seqinr::translate，那么您可以尝试用s2c包裹DNA字符串。

myString <- "AAACCCGGGTTT"
library(seqinr)

# Without s2c
translate(myString)
Error in s2n(seq, levels = s2c("tcag")) : 
  sequence is not a vector of chars

# With s2c
translate(s2c(myString))
[1] "K" "P" "G" "F"

<强>解释

translate()所需字符的向量是"A" "A" "A" "C" "C" "C" "G" "G" "G" "T" "T" "T"而不是"AAACCCGGGTTT" 您可以使用strsplit()生成这样的矢量：

strsplit(myString, "")
[[1]]
 [1] "A" "A" "A" "C" "C" "C" "G" "G" "G" "T" "T" "T"

all(strsplit(myString, "")[[1]] == s2c(myString))
[1] TRUE

翻译功能：“错误：序列不是字符的矢量”

1 个答案: