I'm trying to draw the density function of a mixed normal distribution in the form of: .6*N(.4,.1)+ .4*N(.8,.1) At first I generate a random sample with size 200 by the below code: means = c(.4,.8) sds = sqrt(c(.1,.1)) ind = sample(1:2, n, replace=TRUE, prob=c(.6,.4)) x=rnorm(n,mean=means[ind],sd=sds[ind])
Then I use the below code for drawing the graph: plot(density(x)) The plot doesn't seem to be belonging to the desired distribution, because there is just one mode in it (I've seen the real graph of this mixed normal in a paper, it has two clear distinct modes). Even the hist() doesn't draw a plot similar to the real graph. I think the generation code isn't correct. Is it? (I've asked the generation code here!) -- Free e-mail accounts at http://zworg.com ______________________________________________ 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.
