We are developing a new method based on the Bayesian inference technique for the ground motion computations. The ultimate goal of this project is to have a better understanding of earthquake distribution. We are interested in accurately computing quantities of interests, such as spectral acceleration and peak ground acceleration. These quantities will then be used to produce engineering graphs which are employed to design seismic-resistant buildings and structures.