SPEI infinity values

[phyvera]

Has anyone found a solution to this problem of infinite SPEI values? I calculate the global SPEI using CESM-LE data and get an infinite value (-Inf or -Inf). I already consider na.rm =TRUE even though I have no missing value in my data.
Should I replace -Inf with -3 and Inf with 3, does this make sense?

I have already calculated SPEI ( PET with THornthwaite, Hargreave and Penman) taking different distributions such as Logistic, Gamma, PearsonIII, but I still have these values Infinite and NaN.

Any help and advice is welcome

[sbegueria]

There’s no solution to this, as it is not really a bug. Simply, it’s not possible to fit this distributions to certain data. Look at the data that produce this values, probably it’s way out of range, or contains mostly zeros, or similar.