Example(s) of Biomodal distributions
Would this work? m1 <- 3.5 m2 <- 6.5 s1 <- 0.8 s2 <- 0.8 dist1 <- dnorm(x, mean=m1, sd=s1) dist2 <- dnorm(x, mean=m2, sd=s2) biModalDist <- dist1 + dist2 q <- qnorm(.9,0,1) z1 <- q*s1 z2 <- q*s2 zLow <- m1 - z1 zUp <- m2 + z xShade <- seq(zLow, zUp, 0.1) d1Shade <- dnorm(xShade, mean=m1, sd=s1) d2Shade <- dnorm(xShade, mean=m2, sd=s1) yShade <- d1Shade + d2Shade windows() plot(x, biModalDist, type="l") polygon(x=c(zUp, zLow, xShade), y=c(0,0,yShade), col="gray")
On 1/15/2012 2:41 PM, AbouEl-Makarim Aboueissa wrote:
Dear R users:
I am currently teaching a course in Statistics. Can someone give an R code(s) to create a biomodal curve(s) with shaded area of 90% and with 5% in each tail
With many thanks
abou
==========================
AbouEl-Makarim Aboueissa, Ph.D.
Associate Professor of Statistics
Graduate Program Coordinator
Department of Mathematics& Statistics
University of Southern Maine
96 Falmouth Street
P.O. Box 9300
Portland, ME 04104-9300
USA
Tel: (207) 228-8389
Fax: (207) 780-5607
Email: aaboueissa at usm.maine.edu
aboueiss at yahoo.com
Office: 301C Payson Smith
[[alternative HTML version deleted]]
_______________________________________________ R-sig-teaching at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-teaching ----- No virus found in this message. Checked by AVG - www.avg.com Version: 2012.0.1901 / Virus Database: 2109/4744 - Release Date: 01/15/12