rcon: Constrain random samples

In ssokolen/metcourse: Tools for Metabolic Time-course Data

Constrain random samples

Description

Resamples supplied distribution to ensure generated values fall within given constraints.

Usage

rcon(n, rdist, con, max.iter = 100)

Arguments

n	Number of samples to draw.
rdist	User-defined function for generating random samples (must take n as only argument).
con	Vector of lower and upper constraints on generated samples.
max.iter	Maximum number of iterations before aborting.

Value

A vector of n values sampled from rdist such that the minimum sampled value is greater than or equal to con[1] and the the maximum is smaller than or equal to con[2].

Examples

# Using default distribution parameters
out <- rcon(10000, rnorm, c(-1, 1))
print(summary(out))
hist(out, 20, xlim=c(-2, 2), probability=TRUE, 
     main='', xlab='Random variable value')

# Modifying mean and standard deviation
out <- rcon(10000, function(n) {rnorm(n, 0.5, 0.5)}, c(-1, 1))
print(summary(out))
hist(out, 20, xlim = c(-2, 2), probability = TRUE, 
     main = '', xlab = 'Random variable value')

ssokolen/metcourse documentation built on Nov. 3, 2024, 11:22 p.m.