1

我正在尝试根据我的数据获取 PDF 数据点。它曾经工作到昨晚，但我现在得到的数字非常少。我在这里做错了什么？

dgumbel <- function(x,a,b){ # PDF
  exp((a - x)/b - exp((a - x)/b))/b
}

pgumbel <- function(q,a,b){ # CDF
  exp(-exp(-((q - a)/b)))
}
qgumbel <- function(p, a, b){ # quantile function
  a-b*log(-log(p))
}
gfit <- fitdist(year[,1], "gumbel", start=list(a=5, b=5), method="mle")
#year[1] is the dataset I will be using
para = rgumbel(500,gfit$estimate[1],gfit$estimate[2])
d<-dgumbel(para,gfit$estimate[1],gfit$estimate[2])

但是，如果我在最后一行使用 qgumbel 而不是 dgumbel，我会得到一个 CDF（足够接近）。CDF 很有帮助，但我希望获得 PDF 数据点。编辑：我正在做 rgumbel 因为数据集是一个相当小的数据集，因此我希望通过产生 500 个点来平滑曲线。

笔记：

  • PDF - 概率密度函数

  • CDF - 累积分布函数

4

0 回答 0