I'm having trouble understanding the concept of comparing models of evolution with LRT. I have following models: M0 (neutrality), M1A (negative selection) and M2a (positive selection) with following results from Maximum Likelihood:
M0: w=0.81; lnL=-5702.55 M1a: wo=0.50, po=0.65; w1=1.0, p1=0.35; lnL=-5650.21 M2a: wo=0.46, po=0.64; w1=1.0, p1=0.25; w2=2.38, p2=0.11; lnL=-5632.10.
I want to do the Likelihood Ratio Test to compare the models and make assumptions about the selective pressures on the analyzed sequences. So, first I compare M0 and M1a and calculate LRT=2(lnL1-lnL0)=104.68.
Now I want to take a look at a Chi table, but what am I using there? Am I supposed to used the po value for the M1a model? And if that value exceeds the value of the p-value in the table then the simpler model is chosen?