lrt function

Computes a likelihood ratio test between a reduced model and a full model