4

在 Pandas 中,我正在尝试手动编写卡方检验。我在下面的数据框中进行row 0比较。row 1

data
       2      3      5      10     30
0      3      0      6      5      0
1  33324  15833  58305  54402  38920

为此,我需要将每个单元格的预期单元格计数计算为cell(i,j) = rowSum(i)*colSum(j) / sumAll:在 R 中,我可以简单地通过获取outer()产品来做到这一点:

Exp_counts <- outer(rowSums(data), colSums(data), "*")/sum(data)    # Expected cell counts

我使用 numpy 的外积函数来模仿上述 R 代码的结果:

import numpy as np
pd.DataFrame(np.outer(data.sum(axis=1),data.sum(axis=0))/ (data.sum().sum()), index=data.index, columns=data.columns.values)
       2      3      5      10     30
0      2      1      4      3      2
1  33324  15831  58306  54403  38917

是否可以使用 Pandas 功能实现这一点?

4

1 回答 1

1

仅使用 Pandas 内置方法的完整解决方案:

def outer_product(row):
    numerator = df.sum(1).mul(row.sum(0))
    denominator = df.sum(0).sum(0)
    return (numerator.floordiv(denominator))

df.apply(outer_product)

图片

时序:对于 100 万行 DF。

在此处输入图像描述

于 2016-09-13T10:39:05.920 回答