MaskJointDensity (version 1.0)

rmulti: Simple way to generate noise

Description

Generates noise by selecting from a mixture distribution of normal densities with means, standard deviations and probabilities provided

Usage

rmulti(n, mean, sd, p)

Arguments

n

Size of the noise vector to be outputted

mean

Vector of means corresponding to each component density

sd

Vector of standard deviations corresponding to each component density

p

Probability of selecting from each component distribution

Value

Vector of random numbers

References

no references

Examples

# NOT RUN {
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (n, mean, sd, p) 
{
    x <- rnorm(n)
    u <- sample(1:length(mean), n, prob = p, replace = T)
    for (i in 1:length(mean)) x[u == i] = mean[i] + sd[i] * x[u == 
        i]
    return(x)
  }
# }