library(Delaporte)
x <- c(0, seq_len(25))
plot(x, ddelap(x, lambda = 2, alpha=2, beta=1),type='b', col='blue', xlab='k', ylab=expression(P(X==k)))
par(new=TRUE)
plot(x, ddelap(x, lambda = 2, alpha=2, beta=2),type='b', col='red', axes=F, xlab='',ylab='')
par(new=TRUE)
plot(x, ddelap(x, lambda = 2, alpha=2, beta=3),type='b', col='green', axes=F, xlab='',ylab='')
par(new=TRUE)
plot(x, ddelap(x, lambda = 2, alpha=2, beta=4),type='b', col='black', axes=F, xlab='',ylab='')
legends=c(expression(lambda==2~~alpha==2~~beta==1),
expression(lambda==2~~alpha==2~~beta==2),
expression(lambda==2~~alpha==2~~beta==3),
expression(lambda==2~~alpha==2~~beta==4))
legend('topright', legends, col=c('blue','red','green','black'),
pch=1, lty=1, bty='n', cex=0.9)
##########################################################
x <- c(0, seq_len(50))
plot(x, ddelap(x, lambda = 1, alpha=2, beta=2),type='b', col='blue', xlab='k', ylab=expression(P(X==k)))
par(new=TRUE)
plot(x, ddelap(x, lambda = 5, alpha=2, beta=2),type='b', col='red', axes=F, xlab='',ylab='')
par(new=TRUE)
plot(x, ddelap(x, lambda = 10, alpha=2, beta=2),type='b', col='green', axes=F, xlab='',ylab='')
par(new=TRUE)
plot(x, ddelap(x, lambda = 25, alpha=2, beta=2),type='b', col='black', axes=F, xlab='',ylab='')
legends=c(expression(lambda==1~~alpha==2~~beta==2),
expression(lambda==5~~alpha==2~~beta==2),
expression(lambda==10~~alpha==2~~beta==2),
expression(lambda==25~~alpha==2~~beta==2))
legend('topright', legends, col=c('blue','red','green','black'),
pch=1, lty=1, bty='n', cex=0.9)
