Statistics of Partitioning method applied to the computation of coupled electronic-nuclear dynamics in LiH induced by strong pulses including photoionization