This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to calculate the p-value (?) using LRT from codeml output ?

HI

I constructed LRT

ΔLRT = 2×(lnL1 - lnL0) 
         = (2×(-1114.093613 − (-1116.287406)) = 4.387586.

The degree of freedom is 1 (np1 - np0 = 20 - 19)  = 1

How do I calculate the p-value ?

Little confused,

Suggestions please.

alignment dna gene rna
• 1,837 views • link
written by sunnykevin97 sunnykevin97

Hey check out this tutorial..

https://evosite3d.blogspot.com/2011/09/identifying-positive-selection-in.html
• 0 views • link
written by pinn pinn

Seems good tutorial, I'm not so clear ... Clear explanation is needed in understanding/calculating the p-value (?) using the LRT and df score.

• 0 views • link
written by sunnykevin97 sunnykevin97

Calculating p-value using dchisq() in R.

dchisq(4.387586,1) 
p-value 0.02123465
• 0 views • link
written by sunnykevin97 sunnykevin97

0 answers

No answers yet.

Log in to answer this question.

More posts like this