r - R - 从拟合数据到广义帕累托分布 (GPD) 的双对数图

Question

我想将数据集拟合到 GPD，然后绘制一个双对数图表，显示所有使用数据点，趋势线显示渐近 GPD。我有代码可以使用 VGAM 包来拟合我的数据，但不知道如何获得所需的绘图。代码是：

library("VGAM")
gdata <- data.frame(V1 = rexp(n = 3000, rate = 2))
fit <- vglm(V1 ~ 1, gpd(), data = gdata, trace = TRUE)
Coef(fit)

我想得到一个类似于这里解释的图表