Use same functions for fitting model for main and interaction effects