##---- Should be DIRECTLY executable !! ----##-- ==> Define data, use random,##-- or do help(data=index) for the standard data sets.## The function is currently defined asfunction (n,nbin,r,s){x<-NAfor (iin1:n){pval<-rbeta(1,r,s)x[i]<-rbinom(1,nbin,pval)}x}