- library(graphics)
- xdnorm <- function(x) {return(dnorm(x, mean = 3)) }
- plot(dnorm, from = -3, to = 6, ylab = "Freq", ylim = c(0, 0.45))
- plot(xdnorm, from = -3, to =6, add = TRUE)
- vec1 <- seq(-3, 1.5, by = 0.05)
- vec2 <- seq(1.5, 6, by = 0.05)
- polygon(c(vec1, 1.5), c(xdnorm(vec1), 0), col = "grey50")
- polygon(c(1.5, vec2), c(0, dnorm(vec2)), col = "grey70")
- segments(0, 0, 0, dnorm(0), lty = 3, lwd = 2)
- segments(3, 0, 3, xdnorm(3), lty = 3, lwd = 2)
- text(x = 0, y = 0.42, labels = "G1")
- text(x = 3, y = 0.42, labels = "G2")
- text(x = 0, y = -0.01, labels = "Mu1")
- text(x = 3, y = -0.01, labels = "Mu2")
复制代码
|