如何在直方图的一部分下找到区域？ muhats=replicate（200000，平均值（rexp（18，比率=3））*exp（（-1/2）*平均值（rexp（18，比率=3））） myhist

您可以使用此函数，该函数使用梯形规则对密度进行数值积分，并返回超过特定值的第一个值：

muhats = replicate(200000,mean(rexp(18, rate = 3))*exp(((-1/2)*mean(rexp(18, rate = 3)))) )
myhist <- hist(muhats)
plot(myhist)

因此，x轴上该点左侧的密度面积为0.09

但请注意，这给出的答案与刚才所做的大致相同：

density_area_exceeds(muhat, 0.09)
#> [1] 0.1961847

---

你的意思是

myhist你肯定想要直方图下的区域吗？或者你是想从你的数据中估计一个分位数？估计分位数的方法比看直方图下的面积更好。@Dason试图估计一个分位数是
g
和
muhats
同样的事情？
分位数（g，0.09）
是否更合适？@r2evans我认为它更像
ecdf
的倒数。如果你做了ecdf（g）（分位数（g，x））

你得到了

x

density_area_exceeds(muhat, 0.09)
#> [1] 0.1961847

quantile(g, 0.09)
#>        9% 
#> 0.1971049