Ship operation and ice loading in floe ice fields have received considerable interest during recent years. There have been several numerical simulators developed by different institutes which can simulate ship navigation through floe ice fields and estimate ship performance and local ice loads. However, public data obtained from full-scale measurement covering comprehensively ship performance and ice loads under various ice thicknesses, concentrations and floe sizes are rare. The 2018/19 Antarctic voyage of the Polar Supply and Research Vessel (PSRV) S.A. Agulhas II gathered considerable data of the ship in floe ice fields under various thicknesses, concentrations, and floe sizes. The aim of this paper is to carry out statistical analysis to seek suitable probability distributions which adequately fit the measured ice load and therefore suitable to be used as parent distributions for long-term estimation. For this aim, three categories of probability distributions, namely standard distributions, truncated distributions and mixture distributions are tested. It is found that truncated distributions can fit the load data better than standard distributions bounded at the threshold. In addition, mixture distributions are shown to have promising features, which fit the data well and are able to separate distribution components. Subsequentially, the well-performed distributions are used as parent distributions to make long-term load estimations. The estimation results demonstrate that long-term estimations are sensitive to the selection of parent distribution, which addresses the importance of finding correct distribution to model short-term ice loads. The data of ten selected cases will be published for the use of other researchers.