system.resamp: Systematic Resampling

In mpoll/scale: The Scaleable Langevin Exact Algorithm -- Reference Implementation

Usage

system.resamp(p.wei, n = length(p.wei))

Arguments

p.wei
n

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (p.wei, n = length(p.wei)) 
{
    if ((sum(p.wei) > 0) & (n > 0)) {
        cum.wei <- cumsum(p.wei)/sum(p.wei)
        samp.vec <- seq(runif(1, 0, 1/n), 1, 1/n)
        p.idx <- numeric(0)
        for (i in 1:n) {
            p.idx <- c(p.idx, length(cum.wei[samp.vec[i] >= cum.wei]) + 
                1)
        }
    }
    else {
        p.idx <- numeric(0)
    }
    list(p.idx = p.idx)
  }

mpoll/scale documentation built on Dec. 9, 2019, 7:15 a.m.