r - 密度関数データを見つけようとしています

Question

自分のデータに基づいて 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])

ただし、最後の行で dgumbel の代わりに qgumbel を使用すると、CDF (十分に近い) が得られます。CDF は便利ですが、PDF データ ポイントを取得したいと考えています。編集: データセットがかなり小さいため、rgumbel を実行しているため、500 ポイントを生成して曲線を滑らかにしたいと考えています。

ノート：

• PDF - 確率密度関数

• CDF - 累積分布関数