How can I obtain per site log-likelihoods for a topology against an alignment without optimising branch lengths? (which software can I use to do this)
Per site log-likelihoods can be used for various tests of whether a topology fits the data in a multiple sequence alignment well using e.g. the Shimodaira-Hasagawa test, the AU test etc.
RAxML, Garli and FastTree and other software packages can calculate per site log-likelihoods given a tree, alignment and substitution model but all will automatically optimise branch lengths first.
I would like to test a tree in which a pair of sequences (A and B) are monophyletic, where in the maximum likelihood topology they are not. My problem is that when I infer a tree with a constraint requiring A and B to be monophyletic, the edge length is set to the minimum within the allowed precision. But having an edge length so close to zero does not make biological sense if the tree is describing a monophyletic clade for A and B: the minimum length of the edge must be equivalent to one substitution or character change, else A and B are in a polytomy containing clade with other sequences which is not the biological scenario I'd like to test (in that case there would be ambiguity regarding the relationship of A and B to other groups emerging from that polytomy).
So . . . I'd like to edit my constrained tree and increase the edge length to the A and B clade to the minimum one substitution length, then perform the AU test against the per site log likelihoods for that topology. But the above software optimises the edge lengths, resetting the edge length to A and B to effectively 0!
thanks