如何繪制廣義beta分布圖？

Question

我正在嘗試通過以下方式使用R（軟件）繪制beta-gumbel分布。一般的想法是，在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))

Answer 1

當您說您未使用任何附加軟件包時，我不相信您： pgumbel()不在基本R中pgumbel() library("sos"); findFn("pgumbel") library("sos"); findFn("pgumbel")在很多地方找到它，我使用了evd包。

這里有幾個小問題。

library("evd")

最主要的是您希望length=100而不是by=100 （這給您一個單元素x向量）：

x <- seq(-3, 3, length=100)

實際的計算是可以的：

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

您需要更改ylim才能看到正在發生的事情。 但是，我也認為您需要采取一些措施來解決差分dx以獲得適當的密度函數（盡管這比StackOverflow問題更像是StackExchange，而不是StackOverflow問題）。

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