Python weightstats.ztest与手工计算

我正在比较我的z分数结果，从

weightstats.ztest

和手工计算得到的结果略有不同，z_package=-12.636和z_手册=-13.019

z_数据包由

z_数据包确定，p=weightstats.ztest（data1，data2）

样本量

data1

为574，而

data2

为1026

z_手册由以下代码确定：

y1 = np.array(data1).mean()
y2 = np.array(data2).mean()


s1 = np.array(data1)).std()
s2 = np.array(data2).std()

n1 = len(data1)
n2 = len(data2)

se = math.sqrt(s1**2/n1+s2**2/n2)

dy = y1-y2

z_handbook = (dy-0)/se

---

我曾试图阅读关于如何计算z分数的

weightstats.ztest

文章，但找不到任何好的解释。有人知道z分数的不同吗？

一项，但我不知道是否就这些。默认情况下，numpy std使用ddof=0。请尝试

.std（ddof=1）

第二期：

se

可以在假设差异相等的情况下合并，或者假设差异不相等，如您的示例所示。@Josef the

.std（ddof=1）

不幸地没有给出任何差异。是的，我的手工计算假设方差是非定量的，但我认为软件包假设它是相等的。这很烦人