This function returns an object having the following elements:
"loglik.Modr" maximized log-likelihood of modr; "loglik.Modf" maximized
log-likelihood of modf;
"dif.loglik" difference in log-likelihood between both models, and
"dif.df" difference in degrees of freedong of both models, and
"p-value" is the p-value for the LRT.