Hi all, I've tried to plot a vector which has two peaks in the density. This link shows the figure.
http://docs.google.com/View?docid=dcvdrfrh_1dk9r2rc7 The red line is normal curve and green line is gamma curve. Notice that red line can correctly fit the histogram that has two peaks (i.e. red curve also has two peaks). But the gamma curve there only has one curve. Is there a way I can fit the gamma function such that it also yields two peaks? Below is my code that generate the graph: __BEGIN__ dgexp <- density(gexp.arr) gexp.ymax <- max(dgexp$y) # Histograms gehist <- hist(gexp.arr, col="blue", main=paste(genm,"-",desc, ", Mean=",ge_mean, ", SD=",ge_sd, ", Krt=",ge_kurt, ", AIC=",aic, ", BIC=",bic, sep=""), more = TRUE, xlab = "Exp Level", ylim = c(0, gexp.ymax), freq=FALSE ) # Normal Curve lines(density(gexp.arr), col = "red",lwd = 3) # Gamma Curve curve(dgamma(x,1.25,gexp.ymax),add=TRUE,col="green",lwd=3) __END__ - Gundala Viswanath Jakarta - Indonesia ______________________________________________ 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.
