How to calculate LRT in GREML?