#>source("c:/Work07/ShortCourse07/Programs/NormalMixture.txt",print.eval=TRUE)#
#This gives the distribution of the mixture of two normals#
e<-.25; nsim<-5000;m<-2;s<-1;
u<-(runif(nsim)<e);z<-rnorm(nsim);z1<-rnorm(nsim,mean=m,sd=s);
#This plots histogram and density#
hist(u*z+(1-u)*z1,xlab="x",xlim=c(-5,5),freq=F,col="green",breaks=50,)
mix<-function(x)e*dnorm(x)+(1-e)*dnorm(x,mean=m,sd=s)
xplot<-c(-50:50)/10
par(new=T)
plot(xplot,mix(xplot), xlim=c(-5,5),type="l",yaxt="n",ylab="")