Here is my function for a two-sample pooled binomial z-test based on NIST Engineering Statistics Handbook: (all the p’s are actually p^hat’s)