Sample from the distribution given as a PMF object

Samples (with replacement) from the probability distribution specified by pmf.

Usage

PMF.sample(pmf, N_samples)

Arguments

pmf

the PMF object.

N_samples

number of samples.

Value

Samples drawn from the distribution specified by pmf.

See also

PMF.get_mean(), PMF.get_var(), PMF.get_quantile(), PMF.summary()

Examples

library(bayesRecon)

# Let's build the pmf of a Binomial distribution with parameters n and p
n <- 10
p <- 0.6 
pmf_binomial <- apply(matrix(seq(0,n)),MARGIN=1,FUN=function(x) dbinom(x,size=n,prob=p))

# Draw samples from the PMF object
set.seed(1)
samples <- PMF.sample(pmf=pmf_binomial,N_samples = 1e4)

# Plot the histogram computed with the samples and the true value of the PMF
hist(samples,breaks=seq(0,n),freq=FALSE)
points(seq(0,n)-0.5,pmf_binomial,pch=16)