Accurate prediction of peak ground intensity measures is inevitably influenced by geotechnical variability. Variations in soil properties, subsurface conditions, and seismic inputs introduce complexities that challenge the reliability of predictions. This study introduces a Bayesian generalized linear model (GLM) to probabilistically predict peak ground acceleration (PGA) while accounting for uncertainties associated with geotechnical variability. Latin hypercube sampling (LHS) was employed to generate synthetic datasets of key geotechnical parameters, including plasticity index, shear wave velocity, soil thickness, input motion intensity, and unit weight of soil for hypothetical sites. Subsequently, a series of one-dimensional equivalent linear (1D-EQL) seismic site response analyses were performed, and PGA value at ground surface level were recorded for each analysis. The Bayesian GLM was then developed using these comprehensive datasets to probabilistically predict PGA. The performance and reliability of the developed model were evaluated on a separate test dataset. To benchmark its performance, a Bayesian neural network (BNN) was also developed and compared. In addition, a Shiny-based graphical user interface (GUI), named Bayes-PGA-predictor, was implemented to facilitate practical application. The findings demonstrate that the Bayesian GLM offers a robust and interpretable approach to predicting PGA while effectively quantifying uncertainty associated with geotechnical variability.