数组的bootstrapped后的置信区间问题？

下面我们来分析这段代码

a=np.random.randint(0,5,(2,5))
b=np.random.randn(5)
print(a)
print(b)
means=b[a].mean(axis=0) 
print(means)
confint=np.percentile(means,[25,50]) 
print(confint)

这里有一个统计学上的bootstrap方法，本来是要做b这个array的均值统计量估计，怎么弄，第一行代码里的那个“2”就是做两次，最好呢，是做1000次，好这里简化做两次，a是一个array，行是两行，列数跟b的维度一样。b[a]这个就是取数，取的b的数，维度跟a一致，axis=1是考虑行的，axis=0是考虑列的，mean()带空参数就是整个flatted plane都算，这样means出来的就多个平均，取四分位，中位数的值来表示最后参数估计值。

只是不明白跟置信区间什么关系呢？百分位就对应置信区间。我上面的例子就是25的置信区间？还是感觉结果有些奇怪？

参考

这100道练习，带你玩转Numpy。Url：https://www.kesci.com/home/project/59f29f67c5f3f5119527a2cc

原文：https://www.cnblogs.com/zwydx/p/12327402.html