I usually go for as high a precision as possible first, and then reduce it a bit to take account of cross machine differences. But that's often for more complex code (in statsmodels) where differences to other packages can be anywhere from 4 decimals to floating point precision.

So far we haven't tried to go precision hunting in corner cases in scipy.stats. Although precision problems show up in some issues.