Hi, I'm trying to generate a mixture from 2 univariate normal distributions and would like to put a constrained on the means (Equal) of the two distributions. Here is my shot at the code: library(nor1mix) X <- norMix(mu=c(50, 60 ), sig2=c(5,10), w=c(0.5, 0.5)) mixData <- rnorMix(1000, X) plot(mixData, type="l") I did get the 1000 random numbers in equal proportion but how can I confirm if the means was constrained or not in the data generation process? Any help will be highly appreciated VD [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]]