I am calculating SPI values for different timescales using the SPEI package in R with the PearsonIII distribution instead of the Gamma distribution like it was originally st