假设你有序列AAA和ATA,怎么用R比较它们的差异,即第二个字符,并返回差异的位点与字符?...我用谷歌搜索这个问题时发现stackoverflow上有类似的提问,但不完全一致,基本就是问找出差异的字符,并没有我想要的这么全。...do.call(setdiff, strsplit(c(a, b), split = ""))
# 或者
Reduce(setdiff, strsplit(c(a, b), split = ""))
a,b是两个字符串...AAA"), split = ""))
[1] "T"
> Reduce(setdiff, strsplit(c("ATA", "AAA"), split = ""))
[1] "T"
神奇的是,如果你将两个序列呼唤...na.omit(diff.info)
if(only.position){
diff.info$position
}else diff.info
}
这个函数 可以同时记录位置和原始序列