自助法直方圖
你打算在 5 月去 Cincinnati 度假,但你對 NO2 非常敏感。你擷取了 Cincinnati 過去幾年的 5 月空汙資料,並查看平均 NO2 濃度的自助法(bootstrap)估計結果。因為你只有單一估計值,要視覺化自助法結果最好的方式就是使用直方圖。
雖然你覺得單看自助法直方圖的直覺已經不錯,但和你一起去度假的夥伴更喜歡看到百分位區間。為了顧及他們,你決定將 95% 區間的範圍以底色標示出來。
本練習屬於課程
用 Python 改善你的資料視覺化
練習說明
- 針對
percentile()函式提供能得到 95% 區間所需的上下百分位數。 - 將圖表背景中屬於 95% 區間的範圍加上陰影。
- 以 100 個箱數繪製自助法平均值的直方圖。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
cinci_may_NO2 = pollution.query("city == 'Cincinnati' & month == 5").NO2
# Generate bootstrap samples
boot_means = bootstrap(cinci_may_NO2, 1000)
# Get lower and upper 95% interval bounds
lower, upper = np.percentile(boot_means, [____, ____])
# Plot shaded area for interval
plt.axvspan(____, ____, color = 'gray', alpha = 0.2)
# Draw histogram of bootstrap samples
sns.histplot(____, ____ = 100)
plt.show()