Add support to more complex regressions in beta_gen

[wleoncio]

beta_gen() currently calculates regression parameters based on simple linear regression. It could be expanded to support a formula argument so it could take more complex functions and interactions between covariates.