R 如何绘制广义贝塔分布图？

我试图通过以下方式使用R（软件）绘制beta-gumbel发行版：， genreal的想法是，在beta发行版的pdf中，我们使用gumbel的cdf，而不是插入x。但是我没有得到正确的情节

x <- seq(-3, 3, length=100)
Fx = pgumbel(x,loc=0,scale=1)
y = dbeta(Fx,shape1=0.5,shape2=0.5)
plot(x, y, type="l", lty=2, xlab="x value", ylab="Density",ylim=c(0,1))

---

x当你说你没有使用任何附加软件包时，我不相信你：
pgumbel（）
不在基本R中。
库（“sos”）；findFn（“pgumbel”）
在很多地方都能找到它，我使用了
evd
包

这里有几个小问题

library("evd")

最重要的是，您需要
length=100
而不是
by=100
（这将为您提供单个元素
x
向量）：

您需要更改
ylim
才能看到发生了什么。然而，我也认为您需要做一些事情来解释差异
dx
，以获得适当的密度函数（这更像是一个StackExchange而不是StackOverflow问题）

让我们知道您使用的是什么软件包。我没有安装任何外部软件包
pgumbel
来自何处？它只是一个嵌入在R中的函数，用于计算gumbel的cdf，请尝试在R中使用？pgumbel以获取帮助文件。gumbel和ordinal软件包中至少有
pgumbel
函数，你肯定得下载并安装这些。我试了好几次，但情节看起来不对劲。。。所以我只是输入了整个公式，结果成功了。
x <- seq(-3, 3, length=100)

Fx = pgumbel(x,loc=0,scale=1)
y = dbeta(Fx,shape1=0.5,shape2=0.5)

par(las=1,bty="l")  ## my personal preferences
plot(x, y, type="l", lty=2, xlab="x value", ylab="Density",ylim=c(0,40))