BayesLinRegDTR.model.fit function

Main function for fitting a Bayesian likelihood-based linear regression model