pstocRs<-function(X,Y,B) { n<-length(X) u<-unique(sort(X)) g<-length(u) Y<-as.matrix(Y) k<-dim(Y)[2] XX<-matrix(0,n,g-1) for (i in u[-g]) { XX[1:n,i]<-(X<=i)*1 } mB<-apply(matrix(1:n,n,B),2,sample) mB[,1]<-1:n Z<-as.matrix(Y) Z.star<-Z T<-array(NA,c(B,k,(g-1))) for (b in 1:B){ Z.star<-as.matrix(Z[mB[,b],]) for (d in 1:(g-1)){ T[b,,d]<-apply(Z.star,2,function(col){sum(col[XX[,d]==0])}) } } return(T=T) }