This function compares two haplotypes and returns a summary of the differences.
diffHaplo(h, a = 1, b = 2)an object of class "haplotype".
two integers (or character strings) giving the indices (or labels) of the two haplotypes to be compared.
a data frame with three columns named pos (position of the
differences) and the labels of the two haplotypes compared.
This function prints the number of transitions and transversions between both sequences, and returns a data frame with three columns giving the positions of the differences and the nucleotides in each sequence at these positions.
# NOT RUN {
data(woodmouse)
h <- haplotype(woodmouse)
diffHaplo(h) # compares the 1st and 2nd haplotypes
diffHaplo(h, 1, 3)
diffHaplo(h, "I", "III") # same than above but using labels
# }
Run the code above in your browser using DataLab