fit_ancova function

ANCOVA Analysis