This paper considers the Bayesian model averaging of inverse Gaussian regression models for regression analysis in situations that the response observations are positive and right-skewed. The computational challenges related to computing the essential quantities for executing of this methodology and their dominating ways are discussed. Providing closed form expressions for the interested posterior quantities by considering suitable prior distributions is an attractive aspect of the proposed methodology. The proposed approach has been evaluated via a simulation study and its applicability is expressed by using a real example related to the seismic studies.