如何提取符合R中高斯分布的值？

我有一个数据框X，有两列a和b，a是类字符，b是类数字。 我使用fitdist（fitdistrplus包）函数在b上拟合高斯分布

data.fit <- fitdist(x$b,"norm", "mle")

data.fit通过调用
unclass（data.fit）
可以看到构成
data.fit
对象的所有部分，包括：

$estimate
     mean        sd 
0.1125554 1.2724377 

这意味着您可以通过以下方式访问估计的平均值和标准偏差：

data.fit$estimate['sd']
data.fit$estimate['mean']

要计算拟合分布的上5个百分位数，可以使用
qnorm（）
函数（q表示分位数，顺便说一句）：

threshold <- 
    qnorm(p = 0.95,
          mean=data.fit$estimate['mean'],
          sd=data.fit$estimate['sd'])

x[x$b > threshold,# an indicator of the rows to return
  'a']# the column to return