python - MultiIndex Dataframe 将一个原始索引与其他索引进行比较
问题描述
我有一个以列表为值的数据框。
index=pd.MultiIndex.from_product([["file1", "file2", "file3"], ["a", "b"]])
index.names = ['file', 'name']
data = [
[[1,1],[0,0]],
[[],[]],
[[2,2,2],[7]],
[[],[]],
[[1],[4, 4]],
[[],[]],
]
df = pd.DataFrame(data, index=index, columns=['col1', 'col2'])
df
df
col1 col2
file name
file1 a [1, 1] [0, 0]
b [] []
file2 a [2, 2, 2] [7]
b [] []
file3 a [1] [4, 4]
b [] []
我想在每行和其他行的串联之间进行分组name
并运行 Kolmogorov-Smirnov 测试 ( )。scipy.stats.ks_2samp
名称的示例a
。{file1,a}
== [1,1]
。其他人的串联{file2,a}
+ {file3,a}
== [2,2,2]
+ [1]
== [2,2,2,1]
。他们之间的 KStest 是stats.ks_2samp([1,1], [2,2,2,1])
== 0.75
。我怎样才能得到下面的预期结果(手动完成)?
col1 col2
file name
file1 a 0.75 1.0
b NaN NaN
file2 a 1.0 1.0
b NaN NaN
file3 a 0.6 0.66
b NaN NaN
如果这个问题太临时,我很抱歉。
下面是我的尝试。我无法弄清楚如何优雅地将目标行从其他行中排除。
df.groupby(['name']).apply(
lambda per_name_df: per_name_df.apply(
lambda per_column: per_column.apply(
lambda cell: stats.ks_2samp(cell, np.concatenate(per_column.to_numpy())) if cell else cell)))
解决方案
...在单行和其他行的串联之间进行测试
由于您没有特别指定哪些行,我将给您一个示例来测试第一行和所有剩余行之间的情况:
from scipy.stats import ks_2samp
def ks(a, b):
b = [el for li in b for el in li]
if a and b:
return ks_2samp(a, b)[0]
df.groupby(df.index.get_level_values('name')).col1.apply(lambda x: ks(x[0],x[1:].to_list()))
结果:
name
a 0.75
b NaN
Name: col1, dtype: float64
已编辑问题的更新:
...在每一行和其他行的串联之间进行测试
def ks_all(a):
a = a.to_list()
return [ks(a[i],a[:i]+a[i+1:]) for i in range(0,len(a))]
df.groupby(df.index.get_level_values('name')).transform(ks_all)
结果:
col1 col2
file name
file1 a 0.75 1.000000
b NaN NaN
file2 a 1.00 1.000000
b NaN NaN
file3 a 0.60 0.666667
b NaN NaN
推荐阅读
- arrays - 为什么在 bash 中取消设置数组的第一个元素会擦除整个内容?
- javascript - 彩色方块 HTML/Javascript
- c++ - Events handling in QtScxml state machine
- ruby-on-rails-5 - 升级到 Rails 5 需要多长时间?
- javascript - 使用 C3.js 在仪表图中圆边
- python - 如何在 jupyter notebook 上找到并导入 excel 文件?
- angular - 如何从 Angular 默认错误处理程序中显示用户友好的错误消息?
- file - 我在将文件属性保存到 weblogic 时遇到问题:java.lang.IllegalArgumentException:URI 不是分层的
- php - 将平面数组按键分组为多维数组,如果我们不知道有多少层。PHP
- datetime - 按日期绘制