npar.boot.bin<-function(n.one,nsize,nboots=1000,conf.coef=.95) {
    data.vec <- c(rep(0,nsize-n.one),rep(1,n.one))    # observed data
        phat <- n.one/nsize
      pboots <- matrix(0,nrow=nboots,ncol=1)
    for (ii in 1:nboots) {
       data.boot.vec <- sample(data.vec,replace=T)
        pboots[ii,1] <- sum(data.boot.vec)/nsize
   }
    alpha <- (1-conf.coef)
    list(phat=phat,confid.coef=conf.coef,
        CI=quantile(pboots[,1],c(alpha/2,1-alpha/2)))
 }

 npar.boot.bin(n.one=10,nsize=20,nboots=1000,conf.coef=.95)