# Time-stamp: <2012-09-19 15:28 Petri Koistinen>
# Havainnollistetaan kertymäfunktiomuunnosta

# jakauma on standardinormaalijakauma N(0,1)

n <- 10^5
x <- rnorm(n)

hist(x, prob = TRUE, ylim = c(0, 0.5))
rug(x[1:100])
xgrid <- seq(-3, 3, len = 501)
# dnorm() laskee N(0,1)-jakauman tiheysfunktion arvot:
lines(xgrid, dnorm(xgrid), 'l', col = 'red')

# pnorm() laskee N(0,1)-jakauman kertymäfunktion arvot:
hist(pnorm(x), prob = TRUE)
ugrid <- seq(0,1)
lines(ugrid, dunif(ugrid), 'l', col = 'red')