R - 从拟合数据到广义Pareto分布（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)

我希望得到类似于此解释的图表here