binom11<-rep(NA,5) binom12<-rep(NA,5) binom13<-rep(NA,5) binom14<-rep(NA,5) binom15<-rep(NA,5) binom16<-rep(NA,5) binom17<-rep(NA,5) binom18<-rep(NA,5) binom19<-rep(NA,5) binom11[1]<-pbinom(0, 5, 0.1) binom12[1]<-pbinom(0, 5, 0.2) binom13[1]<-pbinom(0, 5, 0.3) binom14[1]<-pbinom(0, 5, 0.4) binom15[1]<-pbinom(0, 5, 0.5) binom16[1]<-pbinom(0, 5, 0.6) binom17[1]<-pbinom(0, 5, 0.7) binom18[1]<-pbinom(0, 5, 0.8) binom19[1]<-pbinom(0, 5, 0.9) for(i in 2:5) { binom11[i]<-pbinom(i-1, 5, 0.1) } for(i in 2:5) { binom12[i]<-pbinom(i-1, 5, 0.2) } for(i in 2:5) { binom13[i]<-pbinom(i-1, 5, 0.3) } for(i in 2:5) { binom14[i]<-pbinom(i-1, 5, 0.4) } for(i in 2:5) { binom15[i]<-pbinom(i-1, 5, 0.5) } for(i in 2:5) { binom16[i]<-pbinom(i-1, 5, 0.6) } for(i in 2:5) { binom17[i]<-pbinom(i-1, 5, 0.7) } for(i in 2:5) { binom18[i]<-pbinom(i-1, 5, 0.8) } for(i in 2:5) { binom19[i]<-pbinom(i-1, 5, 0.9) } round((binom<-cbind(binom11,binom12, binom13,binom14,binom15, binom16,binom17,binom18, binom19)),digits=3)