1. Display the beta-binomial distribution cdf in png format for various values of alpha and beta and fixed n

require(VGAM)

create.beta.binomial.cdf <- function(N, alpha_v,beta_v,colour="black", pch=16) {

n <- max(length(N),length(alpha_v),length(colour),length(pch))
N <- rep(N,length=n); alpha_v <- rep(alpha_v,length=n)
colour <- rep(colour,length=n); pch <- rep(pch,length=n);labels=NULL 

add.one.series <- function(N, alpha_v,beta_v, colour, pch, maxN) {
  cdf <- pbetabin.ab(0:N, N, alpha_v,beta_v,log=FALSE)

  for (i in 1:N) lines (c(i-1,i), c(cdf[i], cdf[i]), type="s", col=colour)
  lines(c(N,maxN), c(1, 1), type="s", col=colour)
  points(0:N, cdf, col=colour, pch=pch)
}
par(lwd=1, cex=1, mar=c(5,4,2,2)+.1)
plot(x= c(0,max(N)),y=c(0,1),xlab=expression(italic(k)),
 ylab=expression(P(italic(X >= k))),type="n" ,main="")

 axis(1, at=0:(max(N)), lab=0:(max(N)),font.lab=5)
 
for (i in 1:n)
  add.one.series(N[i],alpha_v[i],beta_v[i],col=i,pch[i], max(N))
  

for(i in 1:n){ labels=c(labels,as.expression(bquote(paste(alpha==.(alpha_v[i]),", ", beta==.(beta_v[i])))))

}

legend("bottomright", inset=.05, title=paste("n=",max(N)),legend=labels

, pch=19,lty=1, col=1:n,bty="n") }

png("beta-binomial cdf.png") create.beta.binomial.cdf(10,c(.2,.7,2,600),c(.25,2,2,400)) dev.off()

我，本作品著作权人，特此采用以下许可协议发表本作品：

本文件采用知识共享署名-相同方式共享 3.0 未本地化版本许可协议授权。

您可以自由地：

• 共享 – 复制、发行并传播本作品
• 修改 – 改编作品

惟须遵守下列条件：

• 署名 – 您必须对作品进行署名，提供授权条款的链接，并说明是否对原始内容进行了更改。您可以用任何合理的方式来署名，但不得以任何方式表明许可人认可您或您的使用。
• 相同方式共享 – 如果您再混合、转换或者基于本作品进行创作，您必须以与原先许可协议相同或相兼容的许可协议分发您贡献的作品。

https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue

已授权您依据自由软件基金会发行的无固定段落及封面封底文字（Invariant Sections, Front-Cover Texts, and Back-Cover Texts）的GNU自由文件许可协议1.2版或任意后续版本的条款，复制、传播和／或修改本文件。该协议的副本请见“GNU Free Documentation License”。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

您可以选择您需要的许可协议。