# NOT RUN {
x <- 1:7
GA <- c(1,3,1.5,8); DE <- c(1,3,4,2.5)
par(mfrow=c(2,2))
for (j in 1:4) {
plot(x,dgBetaD(x,gam=GA[j],del=DE[j]),type="h",
panel.first=points(x,dgBetaD(x,gam=GA[j],del=DE[j]),pch=19),
main=paste("gamma=",GA[j]," delta=",DE[j],sep=""),ylim=c(0,.6),
ylab="dgBetaD(x)")
}
# }
Run the code above in your browser using DataLab