SAGA06
拟合GPD,以下是程序,运行总不出,谢谢
curve(dgpd(x, shape=0), lwd=3, xlim=c(-.1,2), ylim=c(0,2),ylab="", main="density of the Generalized Pereto Distribution”)
for (s in seq(-2,2,by=.1)) {
curve(dgpd(x, shape=s), add=T)
}
curve(dgpd(x, shape=-1), lwd=3, add=T, col=’red’)
curve(dgpd(x, shape=1), lwd=3, add=T, col=’blue’)