Likelihood-based haplotype frequency modeling using variable-order Markov chains