R/stratblockalloc.R

#' @export
stratblockalloc <-
function(n,blksize) {
    
    result = NULL
    numcompleteblocks = n %/% blksize
    if (numcompleteblocks > 0) {
      block = c( rep(1:numcompleteblocks,each=blksize),rep(numcompleteblocks+1,n%%blksize) )
      for (i in 1:numcompleteblocks) {
      result = c(result,sample(c(rep(0,blksize/2),rep(1,blksize/2)),blksize,replace=F))
    }}
    if (numcompleteblocks ==0) {
       block = rep(1,n)
    }
    if (max(block) > numcompleteblocks) {
       result = c(result,sample(c(rep(0,blksize/2),rep(1,blksize/2)),n%%blksize,replace=F))
    }
    result
}

Try the SeqAlloc package in your browser

Any scripts or data that you put into this service are public.

SeqAlloc documentation built on May 2, 2019, 3:14 p.m.