shape1 = 1; shape2 = 2; y = 0:30
proby = dbetageom(y, shape1, shape2, log=FALSE)
plot(y, proby, type="h", col="blue", ylab="P[Y=y]", main=paste(
"Y ~ Beta-geometric(shape1=",shape1,", shape2=",shape2,")", sep=""))
sum(proby)
Run the code above in your browser using DataLab