我有一些fasta序列,由数十万个字母组成“a”“g”“c”“t”
library(seqinr)
mydata<-read.fasta(file="sequence.fasta")
mydata1<-mydata[[1]]
我想把它切成相等长度的矢量,比方说大小为100.手动它会看起来像
vec1<-mydata1[c(1:100)]
vec2<-mydata1[c(101:200)]
等。我无法访问任何其他编码程序,因此必须在r中完成。我正在考虑某种类型的for循环,但我不知道如何实现它。这项任务在r中是否可行?
编辑:这个问题与突出显示的问题不同,因为我没有数字向量,但DNA序列由字母组成。
答案 0 :(得分:1)
如果我们需要用长度为100的向量进行分割
lst <- split(my.data, as.integer(gl(length(my.data), 100, length(my.data))))