One of the most important issues in nuclear astrophysics is the lithium problem, which is the disagreement between the theoretical and experimental values of the lithium abundance of the Big Bang nucleosynthesis. Solving this problem or establishing compatibility between the two mentioned values can be considered as the confirmation of one part of the Big Bang Theory. Bayesian statistics is one of the most recent methods used to improve the calculation accuracy of the abundance of primary nuclei after the Big Bang and their reaction rates. In this study, Bayesian statistics is used to fit the S-factor and also calculate the rate of astrophysical reactions d(p,γ)3He and 3He(d,p)4He using available experimental data. The results of this study are in good agreement with similar Bayesian researches, and can be considered as an improvement compared to the results of the other (non-Bayesian) methods.