This function will fit the full and reduced models and calculate LRT raw p-value and adjusted p-value based on BH Method
Useful links