lr_test function

Likelihood ratio test for nested models