Hi.. I'm not sure why polygon returns an area above the standard normal curve.

z <- pretty(c(-3,3), 100)
ht <- dnorm(z)
data <- data.frame(z=z, ht=ht)
zc <- 1.645
plot(data, type="l")
lines(data)
t <- subset(data, z>zc)
polygon(t, col="red")

Thanks,
Lance


	[[alternative HTML version deleted]]

