R/coupledresampling-systematic-given.R

In pierrejacob/CoupledPF: Coupling of Particle Filters

#'@rdname CR_systematic_given
#'@title Coupled Resampling: systematic, given vector of ancestor
#'@description This function performs systematic resampling on the second system
#'@return A vector of ancestors
#'@export
CR_systematic_given <- function(xparticles1, xparticles2, normweights1, normweights2, 
                                              parameters = list(),
                                              ancestors_ref, uniforms){
  # u <- runif(1)
  nparticles <- ncol(xparticles1)
  ancestors2 <- systematic_resampling_n(normweights2, nparticles, uniforms[1])
  return(ancestors2)
}