## Time-stamp: <2013-09-16 14:11 Koistinen Petri P M>

##
## 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:
y <- pnorm(x)
hist(y, prob = TRUE)
ugrid <- seq(0,1)
lines(ugrid, dunif(ugrid), 'l', col = 'red')


## Local Variables:
## coding: utf-8
## End: