extarct_f2_subset does not copy block_lengths file about admixtools HOT 2 OPEN

Early on, when I wrote the extract_f2_subset() function, I thought it would be a good idea to compute f2-statistics for a large number of populations, and then create subsets of that. Because the missing data problem tends to get worse for more populations, I abandoned that idea, and I haven't updated the extract_f2_subset() in a while. It might work if you copy the rds files with the block length information into the new directory.

Another thing you could try is to pass the pops argument to f2_from_precomp():
f2_blocks = f2_from_precomp('./f2sub', pops = c('2','3','5','8','9'))
I have never tested this with numbers as populations names. If you run into problems with that, try using names that start with letters.

Please let me know if this works, or if there is a reason why extract_f2_subset() would be better in your case than passing the pops argument to f2_from_precomp()!

from admixtools.

Hello
Copying the rds files works OK. However your suggestion for:
f2_blocks = f2_from_precomp('./f2', pops = c('2','3','5','8','9'))
seems more reasonable. I just did not pay attention to this option. Working with numbers as characters works fine for me.
Thank you

from admixtools.