strat.resamp: Stratified Resampling

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

Usage

strat.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)) 
{
    vec.length <- length(p.wei)
    cum.wei <- cumsum(p.wei)
    if (cum.wei[vec.length] > 0) {
        cum.wei <- cum.wei/cum.wei[vec.length]
        samp.vec <- seq(0, (n - 1)/n, 1/n) + runif(n, 0, 1/n)
        p.idx <- findInterval(samp.vec, 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.