Question

我正在使用Python中的pandas库。假设我有以下方式从正态分布中抽取四个随机样本：

np.random.seed(12345)

df = pd.DataFrame([np.random.normal(32000,20000,3650), 
                   np.random.normal(43000,10000,3650), 
                   np.random.normal(43500,14000,3650), 
                   np.random.normal(48000,7000,3650)], 
                  index=[1992,1993,1994,1995])

我希望每个样本获得95％的置信区间，因此我计算：

mean_value=df.mean(axis=1)
std_value=df.std(axis=1,ddof=0)
lower_bound=mean_value-1.96*std_value
upper_bound=mean_value+1.96*std_value
diff = upper_bound-lower_bound

对于每个置信区间，我想将其切成11个等距间隔。我有类似以下的想法：

low=lower_bound.values[1]
high=upper_bound.values[1]
diff=0.09*diff.values[1]
np.arange(low,high,diff)

这不太有效，因为切割间隔不会在置信区间的上端结束。我怎样才能获得等距间隔？

Answer 1

我不确定你想要什么，但是使用NumPy的linspace函数获得相等间隔的间隔非常容易。以下是第一次分发的11个区间。

np.linspace(lower_bound.values[0], upper_bound.values[0], 12)
array([ -7.18705879e+03,  -3.82825067e+01,   7.11049377e+03,
         1.42592701e+04,   2.14080463e+04,   2.85568226e+04,
         3.57055989e+04,   4.28543752e+04,   5.00031514e+04,
         5.71519277e+04,   6.43007040e+04,   7.14494803e+04])

多个置信区间

1 个答案: