A Bayesian method is proposed to predict ice loads according to loading conditions. Both the shape and the scale of load distribution vary according to the covariates. With increasing ice thickness ice loads increase in general. With increasing speed the maximum loads increase whereas the median loads decrease.