일반화 된 베타 배포판을 만드는 방법은 무엇입니까?

Question

R (소프트웨어)를 사용하여 베타 - 검버섯 분포를 다음과 같이 플롯하려고합니다. 베타 배포판의 PDF에서 장황한 생각은 x를 연결하는 대신 gumbel 대신 cdf를 사용한다는 것입니다. 그러나 나는 올바른 음모를 얻을 수 없었다.

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()는 library("sos"); findFn("pgumbel")는 다양한 장소에서 발견 기본 R.에 있지, 나는 evd 패키지를 사용했다.

여기에는 몇 가지 작은 문제가 있습니다.

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

실제 계산은 OK입니다 :

library("evd") 

중요한 것은 당신이 (당신에게 단일 요소 x 벡터를 제공하는) length=100보다는 by=100 할 것입니다

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

무슨 일이 일어나고 있는지 보려면 ylim을 변경해야합니다. 그러나, 나는 또한 적절한 밀도 함수를 얻으려면 차이점 dx을 설명하기 위해 무언가를해야한다고 생각합니다. StackOverflow보다 더 많은 StackExchange가 있습니다.

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