I am calculating 2*(modelLR - nestedModelLR) and then in R doing pchisq(value, df, lower.tail=FALSE). I decided on lower.tail=FALSE because my hypothesis is that the model with positive selection will be a better fit (so a one-direction prediction).
My questions are:
(1) Is this the correct chi-squared test?
(2) Do I have the tails correct in my chi-squared test?
(3) Do I use the same statistical method for branch, site, and branch-site models?
Thank you in advance, and apologies for the basic question--I haven't found a good explanation for this step.
Best,
Kim