## Paul Johnson
## 2013-02-18

## response to Fan's question about assingment.

myMu <- rep(0, 1000)


getNormalMean <- function(mu, sigma, N){
    e <- rnorm(N, m = mu, sd = sigma)
    mean(e)
}

##first arg is "mu" from myMu.
myMeans1 <- sapply(myMu, getNormalMean, sigma=23, N = 15)

hist(myMeans1)

plot of chunk unnamed-chunk-1

getPoissonMean <- function(lambda, N){
    e <- rpois(N, l = lambda)
    mean(e)
}



myLambda <- rep(0.7, 1000)

myMeans2 <- sapply(myLambda, getPoissonMean, N = 10)



myLambda <- rep(0.7, 1000)

myMeans2 <- sapply(myLambda, getPoissonMean, N = 1500)

hist(myMeans2, prob = TRUE)
lines(density(myMeans2))

plot of chunk unnamed-chunk-1