Output fitted theta parameter files

[LuZhangH]

Hi,

I have enjoyed applying RFMix in my population genetics project! I have one question about the output files: whether it is possible to output the fitted parameter files, so that I can customize my inference procedures after plugging in the theta parameters to the linear chain CRF, instead of directly using the marginal probabilities by the forward-backward algorithm (fb.tsv files) or the most likely subpopulation assignment by the Viterbi algorithm (msp.tsv files). A quick example will be, having the fitted theta parameters, I can sample the local ancestries given the haplotypes jointly from the linear chain CRF.

Thank you a lot for the help!

Best,
Lu