R/importanceResamplingWithoutReplacementS2.R

In rohan-shah/particleMethodsBernoulli: Monte Carlo methods for a simple Bernoulli example

importanceResamplingWithoutReplacementS2 <- function(nBernoullis, lowerBound, probability, probabilities, n, seed, k)
{
	if(missing(probabilities))
	{
		if(missing(probability) || missing(nBernoullis))
		{
			stop("If input probabilities is missing, then inputs nBernoullis and probability must be given")
		}
		probabilities <- rep(probability, nBernoullis)
	}
	else
	{
		if(!missing(probability) || !missing(nBernoullis))
		{
			stop("If input probabilities is used, then inputs nBernoullis and probability must be missing")
		}
	}
	start <- Sys.time()
	result <- .Call("importanceResamplingWithoutReplacementS2", lowerBound, probabilities, n, seed, k, PACKAGE="particleMethodsBernoulli")
	end <- Sys.time()
	return(new("importanceResamplingResult", estimate = result$estimate, call = match.call(), start = start, end = end))
}