Generating a genome-wide tree topology count

[MatteoSebastianelli]

Hello,

I have been generating tree sequences with tsinfer and would like to get a tree topology count across the entire genome. Would that be possible? I was thinking I should loop through my .trees files with the union function to merge them (it seems it cannot merge more than 2 files) and then use that output to count tree topologies. Does this sounds like a reasonable approach? I may be missing it, perhaps there is a functions that does this already?

Thanks in advance!

Matteo

[hyanwong]

Hi Matteo

If you mean that you want to summarize across multiple chromosomes, each inferred separately by tsinfer, then I think it may not be necessary to union the tree sequences: you can probably summarise the topology results across chromosomes by just adding the means (weighted by the chromosome length) and dividing by the total genome length. Note that you might want to remove centromeres and other "blank" regions from the length calculations.

Within a chromosome, however, there are fast methods for counting topologies: see tskit-dev/tskit#671 and the link from there to https://github.com/tskit-dev/tskit/blob/829fd1f42037ca404c44ea58eb6286df8d5d803d/python/tskit/combinatorics.py#L37. We should probably make this available as an API, if people are going to use it extensively.

[hyanwong]

P.s. I think this should be a tskit GitHub discussion (not a tsinfer GitHut issue), so I've moved it to tskit-dev/tskit#3082