swh:1:snp:ffdd0a7d2c8ea15ad41d45b3b178f668bd942287
Tip revision: 41940cb670c6adeaf624dd1015755490189dd568 authored by Derek Young on 07 February 2020, 11:20:02 UTC
version 1.2.0
version 1.2.0
Tip revision: 41940cb
rmvnormmix.R
# Note: normmixrm.sim is here for backwards compatibility
rmvnormmix <- normmixrm.sim <- function(n,lambda=1,mu=0,sigma=1) {
m <- length(lambda) # nb of components
mu <- matrix(mu, nrow=m)
sigma <- matrix(sigma, nrow=m)
if ((r <- NCOL(mu)) != NCOL(sigma)) {
stop("mu and sigma must have the same number of columns", call.=FALSE)
}
z <- sample(m,n,replace=TRUE,prob=lambda) # component
matrix(rnorm(n*r,mean=as.vector(mu[z,]),sd=as.vector(sigma[z,])),n,r)
}