Skip to content
Prev 395 / 885 Next

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: