Write an R function from scratch to conduct forward selection. Your function should allow the following user-specified options: selection based on adjusted R2 , AIC, as well as BIC. You can use built in functions in R like lm() or summary() but you cannot use any functions dedicated to model selection like step, etc. Need to handle continuous covariates. Have option for handling categorical covariates and hierarchically well-formulated interactions. Comment your code.