统计模型ztest:在proportions_ztest中使用了哪些方差？

Question

我在试着理解proportions_ztest是如何工作的

就像我跑步时所理解的那样

stat, pvalue = proportions_ztest(count = 40, nobs = 20000, value = 0.003, alternative = "smaller") 
pvalue

它应该给出与p值相同的结果。

stats.norm.cdf(40, 60, np.sqrt(20000* 0.003*(1-0.003)))

但事实并非如此。我的假设是，使用的是不同的差异，但我没有找到任何相关的文档。

Answer 1

statsmodels proportions_ztest默认使用基于估计比例的方差，而不是基于空比例的方差。

对此的一种解释是默认是Wald测试，而不是分数测试。

参见这里的答复和评论，https://stats.stackexchange.com/questions/447462/why-does-the-proportions-ztest-function-in-statsmodels-produce-different-values