Hi
I'm wondering if you've ever run into this issue that your implementation deviates from R and scipy for extreme ORs with low p-values? Example is below. I did notice that scipy seems to use int64 and even states in a comment that int32 is not sufficient (though that is with their implementation and I honestly haven't looked hard enough yet to see how yours differs, scipy note about int32)
We need these low p-values because we are needing to do a multiplicity adjustment and we have a lot of p-values (hence why we need a fast implementation!) Let us know if you have any ideas.
Thanks,
Andrew
I think this is a limitation of the implementation here. If you need p-values that low then you must use scipy. I thought recent versions of scipy were quite fast.
Hi I'm wondering if you've ever run into this issue that your implementation deviates from R and scipy for extreme ORs with low p-values? Example is below. I did notice that scipy seems to use
int64
and even states in a comment thatint32
is not sufficient (though that is with their implementation and I honestly haven't looked hard enough yet to see how yours differs, scipy note about int32)We need these low p-values because we are needing to do a multiplicity adjustment and we have a lot of p-values (hence why we need a fast implementation!) Let us know if you have any ideas. Thanks, Andrew
scipy and brentp
R