#'@rdname CR_multinomial
#'@title Coupled Resampling: multinomial
#'@description This function performs independent multinomial resampling on each system
#'@return Two vectors of ancestors, column-binded in a matrix.
#'@export
CR_multinomial <- function(xparticles1, xparticles2, normweights1, normweights2, ...){
nparticles <- ncol(xparticles1)
ancestors1 <- multinomial_resampling_n(normweights1, nparticles)
ancestors2 <- multinomial_resampling_n(normweights2, nparticles)
# ancestors1 <- sample(x = 1:nparticles, size = nparticles, replace = TRUE, prob = normweights1)
# ancestors2 <- sample(x = 1:nparticles, size = nparticles, replace = TRUE, prob = normweights2)
return(cbind(ancestors1, ancestors2))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.