1
$\begingroup$

I am looking at 59 mitochondrial CO1 DNA sequences of length 563bp, but it is a haplotype file and it has 9 haplotypes within.

When I use .arp files in arlequin to run mismatch distribution it says that can’t converge beyond 2000 steps and also no graphs are produced. Can someone suggest what can I do to overcome this problem?

$\endgroup$
0

1 Answer 1

0
$\begingroup$

What you appear to have is a Bayesian MCMC for 2000 "steps" which isn't a lot and normally the term 'steps' is used for parsimony, so I'm a bit confused.

However, if a Bayesian analysis does not reach convergence, i.e. duplicate runs do not achieve the same results then the data space has not been solved so you get no result.

9 haplotypes does not sound much to me, it is possible there are too many sequences, it is also possible there is simply not enough variation (9 haplotypes ... how many alleles?).

Suggestions

  • Arlequin should provide a default convergence value which you can lower to force a result
  • Decrease the number of sequences but maintain the number of haplotypes
  • Try their example data set

Do you get any associated file such as posterior probability, any output file with a minus sign in front of it?

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.