Statistics for A path integral Monte Carlo method for obtaining thermodynamic properties of nonadiabatic systems