Re: [R-sig-phylo] PEGAS: assignment to haplotype when missing information

2020-02-26 Thread Emmanuel Paradis
(Cc. to r-sig-genetics) I'm going to modify the algorithm in haplotype.DNAbin() as follows: 1. Find the sequences that are exactly identical, so that, eg, the 3 sequences: A- AR AA would be treated as different at this step. 2. Substitute the leading and trailing "-" for N (thus keeping the al

Re: [R-sig-phylo] PEGAS: assignment to haplotype when missing information

2020-02-26 Thread Emmanuel Paradis
Hi Hirra, The assignment is not random, it follows the order of the sequences in the data: - Seqs. A and B are compared and found to be identical so they are both assigned to haplotype I. - Seq. C is compared to haplotype I (effectively seq. A) and found to be different so it is assigned to hap

[R-sig-phylo] PEGAS: assignment to haplotype when missing information

2020-02-25 Thread Hirra Farooq
Hello, I am using the pegas R package to assign sequences into haplotypes. I recently tried out a test examples with 4 sequences. 2 of the sequences (A and B) are identical, 1 sequence (Seq C) differs from these at only one position (pos 648). The 4th sequence (Seq D) is identical to all but sh