I am currently revamping small sample defaults of pyfixest (the goal is to match fixest SEs out of the box). Now some unit tests for Poisson comparing Python vs R fail for derived statistics (tstats, pvalues, CIs) because (I think) of restrictive tolerances.