Python Mann-Whitney置信区间

import numpy as np
from scipy.stats import mannwhitneyu

# 两组样本数据
group1 = np.array([1, 2, 3, 4, 5])
group2 = np.array([6, 7, 8, 9, 10])

# 自助法重复抽样次数
n_bootstrap = 1000

# 存储Mann-Whitney U检验统计量的数组
u_values = np.zeros(n_bootstrap)

# 进行自助法重复抽样并计算Mann-Whitney U检验统计量
for i in range(n_bootstrap):
    resampled_group1 = np.random.choice(group1, size=len(group1), replace=True)
    resampled_group2 = np.random.choice(group2, size=len(group2), replace=True)
    u_values[i] = mannwhitneyu(resampled_group1, resampled_group2, alternative='two-sided').statistic

# 计算置信区间
confidence_interval = np.percentile(u_values, [2.5, 97.5])

print("Mann-Whitney U检验置信区间:", confidence_interval)