Hi all,
With the caveat that RevGadgets is its own animal, I wonder if anyone has some tips here.
I have an ASE tree inferred from a DEC+J model and am using proccessAncStates() and plotAncStatesPie() to make a labeled tree. My issue is getting appropriate states to be plotted.
I have 6 areas and therefore 64 states, of which only 20 were recovered in end_states_1, _2, or _3. These states are not entirely consecutive, and some are state numbers greater than 23. RevGadgets wants to plot states 1:23, regardless of "occupancy". This leaves me with several unused labels/colors, while omitting several states that should be shown.
I feel this must be a relatively doable reconfiguration, but I'm getting lost in all the ggtree statements.
Any thoughts greatly appreciated!
-Ziv