Two partitioning schemes were tested, one unlinking the model parameters between nuclear and plastid regions, the other unlinking all markers (i.e three partitions). The best-fit model was chosen by Bayes factor test based on the harmonic mean of the respective long likelihoods .