Measure changes in state across the phylogeny?

[maxgmarin]

Hello! First off, thank you for making this software.

I am interested in looking at changes in gene presence absence within a population. I have a phylogeny for my genome and I have a binary matrix that specifies per genome the presence (1) or absence (0) of a gene. I would like to calculate change in state for each gene at each node/branch of the phylogeny. Is there any functionality in the library for doing this?

I am interested in being able to identify the potential number of gain/loss events of a gene (and other features too).

(I could parse the ancestral states and the tree with ETE3 and calculate this my self, but I wanted to check if there is already implemented code.)

Thank you,
Max