每组python 1:1分层抽样

时间:2019-02-12 15:19:23

标签: python pandas multisampling subsampling

如何在python中进行1:1分层采样?

假设熊猫数据框df严重失衡。它包含一个二进制组和多列类别子组。

df = pd.DataFrame({'id':[1,2,3,4,5], 'group':[0,1,0,1,0], 'sub_category_1':[1,2,2,1,1], 'sub_category_2':[1,2,2,1,1], 'value':[1,2,3,1,2]})
display(df)
display(df[df.group == 1])
display(df[df.group == 0])
df.group.value_counts()

对于主要group==1的每个成员,我需要找到group==0与之的单个匹配项。

scikit-learn中的StratifiedShuffleSplit仅返回数据的随机部分,而不是1:1匹配。

1 个答案:

答案 0 :(得分:1)

如果我正确理解,可以使用np.random.permutation

import numpy as np
import pandas as pd

np.random.seed(42)

df = pd.DataFrame({'id': [1, 2, 3, 4, 5], 'group': [0, 1, 0, 1, 0], 'sub_category_1': [1, 2, 2, 1, 1],
                   'sub_category_2': [1, 2, 2, 1, 1], 'value': [1, 2, 3, 1, 2]})

# create new column with an identifier for a combination of categories
columns = ['sub_category_1', 'sub_category_2']
labels = df.loc[:, columns].apply(lambda x: ''.join(map(str, x.values)), axis=1)
values, keys = pd.factorize(labels)
df['label'] = labels.map(dict(zip(keys, values)))

# build distribution of sub-categories combinations
distribution = df[df.group == 1].label.value_counts().to_dict()

# select from group 0 only those rows that are in the same sub-categories combinations
mask = (df.group == 0) & (df.label.isin(distribution))

# do random sampling
selected = np.ravel([np.random.permutation(group.index)[:distribution[name]] for name, group in df.loc[mask].groupby(['label'])])

# display result
result = df.drop('label', axis=1).iloc[selected]
print(result)

输出

   group  id  sub_category_1  sub_category_2  value
4      0   5               1               1      2
2      0   3               2               2      3

请注意,此解决方案假定group 1的每个可能的sub_category组合的大小小于group 0中的相应子组的大小。一个更强大的版本涉及使用np.random.choice进行替换:

selected = np.ravel([np.random.choice(group.index, distribution[name], replace=True) for name, group in df.loc[mask].groupby(['label'])])

带有选择版本的版本与带有置换版本的假设不同,尽管每个子类别组合至少需要一个元素。

相关问题