hello all,
using ape3.0-8, I seem to get pairwise genetic distances that are not
consistent with the IUPAC code.
> dist.dna( as.DNAbin( rbind( c("a","a"),c("s","a") ) ), model="raw" )
1
2 0
S means C or G, so I would expect a distance of 0.5 ?
Am I doing anything wrong ?
Thank you,
Oliver
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
