resid.resamp: Residual Resampling

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

Usage

resid.resamp(p.wei, n = length(p.wei), nest.resamp = strat.resamp)

Arguments

p.wei
n
nest.resamp

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), nest.resamp = strat.resamp) 
{
    if ((sum(p.wei) > 0) & (n > 0)) {
        det.resamp <- floor(n * p.wei/sum(p.wei))
        p.idx <- c(rep(1:length(p.wei), det.resamp), nest.resamp(p.wei - 
            det.resamp * (sum(p.wei)/n), n = n - sum(det.resamp))$p.idx)
    }
    else {
        p.idx <- numeric(0)
    }
    list(p.idx = p.idx)
  }

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