我设法写了一个for loop
来比较以下向量中的字母:
bases <- c("G","C","A","T")
test <- sample(bases, replace=T, 20)
test
将返回
[1] "T" "G" "T" "G" "C" "A" "A" "G" "A" "C" "A" "T" "T" "T" "T" "C" "A" "G" "G" "C"
使用该功能,Comp()
我可以检查一个字母是否与下一个字母匹配
Comp <- function(data)
{
output <- vector()
for(i in 1:(length(data)-1))
{
if(data[i]==data[i+1])
{
output[i] <-1
}
else
{
output[i] <-0
}
}
return(output)
}
导致;
> Comp(test)
[1] 0 0 0 0 0 1 0 0 0 0 0 1 1 1 0 0 0 1 0
这是有效的,但是数量很大时速度很慢。为此我尝试了sapply()
Comp <- function(x,i) if(x[i]==x[i+1]) 1 else 0
unlist(lapply(test, Comp, test))
不幸的是它不起作用......(Error in i + 1 : non-numeric argument to binary operator
)我无法弄清楚如何访问向量中的前一个字母来比较它。此外length(data)-1
,“不比较”最后一个字母的 , 可能会成为一个问题。
谢谢大家的帮助!
干杯幸运