qoiqpwqr 发表于 2014-5-4 12:11 
我不懂DNA,所以没看懂。但是不管怎么说,用类似
x[x %in% c("AG")]
install.packages("stringr")
library(stringr)
sub("--","NA",original[1:2181])
fun2<-function(x){
a<-substr(x[2],start=1,stop=1)
for(i in 2:2181){
x
<-str_count(x,a)
}
return(x) #自编函数
}
my.2<-apply(gene,1,fun2)
这个是我自己写的函数,但是不能处理含缺失值的数据?能帮我看下吗?谢谢了