I am comparing stats.ttest_ind() vs “manual” computation of the same test, and get different results. stats.ttest_ind() method: Out: Manual method: Out: We can see there’s a small difference. Why? Maybe because of how stats.ttest_ind() computes degrees of freedom? Any insight much appreciated. Answer The following works. It is your code from above, with only two rows changed. and it outputs