Fit the model using fixed effects model and plm function.
For the model dplyr::selection, testing and validation see "./rhap/inst/extdata/model_testing.R"
fit_model(HIA_var)
Details on plm estimation can be found here: https://cran.r-project.org/web/packages/plm/plm.pdf#'
Regression model and coefficients for prediction