In this paper we propose Bayesian methods to fit econometric regression models, including those where the variability is assumed to follow a regression structure. We formulate the main functions of the statistical R-package BSPADATA, developed according to the proposed methods to obtain posteriori parameter inferences. After that, we include results of simulated studies to illustrate the use of this package and the performance of the proposed methods. Finally, we provide studies to illustrate the applications of the models and compare our results with that obtained by maximum likelihood.