R/branching.resample.R

In smcUtils: Utility functions for sequential Monte Carlo

branching.resample <-
function(weights,num.samples=length(weights), engine="R") {
  if (engine!="R") stop("Only engine='R' implemented for branching.resample")

  expected.num.samples = weights*num.samples
  deterministic.reps   = floor(expected.num.samples)
  random.reps          = rbinom(num.samples,1,expected.num.samples-deterministic.reps)
  
  return(rep2id(deterministic.reps+random.reps))
}