# R/CompoundDensPlot.r In GB2: Generalized Beta Distribution of the Second Kind: Properties, Likelihood, Estimation

```dplot.cgb2 <- function(x, shape1, scale, shape2, shape3, pl0, pl, w=rep(1,length(x)) ,decomp="r",
xmax = max(x)*(2/3), choicecol=1:3, kernel="epanechnikov", adjust=1, title=NULL, ylim=NULL){
para <- paste(" a=",round(shape1,2),", b=",round(scale),", p=",round(shape2,2),", q=",round(shape3,2))
sub0=paste("pl0 = (",round(pl0[1],3))
sub=paste("pl = (",round(pl[1],3))
pl1 <- length(pl0)-1
if (pl1 >= 2){
for (i in 2:pl1) {
sub0 <- paste(sub0,",", round(pl0[i],3))
sub <- paste(sub,",", round(pl[i],3))
}
}
sub0 <- paste(sub0,",",round(pl0[pl1+1],3),")")
sub <- paste(sub,",",round(pl[pl1+1],3),")")
fcgb2 <- function(x) dcgb2(x, shape1, scale, shape2, shape3, pl0, pl, decomp=decomp)
fgb2 <- function(x) dgb2(x, shape1, scale, shape2, shape3)
#    maxx <- max(x)*2/3   # change 28.04.2014: put as argument
if (is.null(ylim)) curve(fcgb2, col=choicecol[2], lwd=2, lty=1, from=0, to=xmax, ylab="Density")
else curve(fcgb2, col=choicecol[2], lwd=2, lty=1, from=0, to=xmax, ylab="Density", ylim=ylim)
curve(fgb2,col=choicecol[1], lwd=2, lty=2, add=TRUE)
if (is.null(title)) title <- "Comparison of densities"
title(title, sub = paste(para,"; ",sub0,";",sub),
cex.sub = 0.75, font.sub = 3)

## empirical counterparts
wk <- w/sum(w)
densk <- density(x, weights=wk, kernel= kernel, from=0, adjust=adjust)
lines(densk, col=choicecol[3], lwd=2, lty=3)
#     print("Please, place the cursor for the legend",quote = FALSE)
legend("topright",c("GB2 ","compound GB2 ","Kernel estimate"),
lwd=2,col=choicecol[c(1,2,3)], lty=c(2,1,3))
}
```

## Try the GB2 package in your browser

Any scripts or data that you put into this service are public.

GB2 documentation built on May 2, 2019, 5:53 a.m.