Individual weights of a feral strain of cat are normally distributed with a mean of 2400 g and variance of 642 g2.
Calculate the following probabilities:
meanCats <- 2400
varCats <- 642
sdCats <- sqrt(varCats)
pnorm(2385, meanCats, sdCats)
## [1] 0.2769244
pnorm(2385, meanCats, sdCats)^2
## [1] 0.0766871
mean10 <- meanCats
sd10 <- sdCats/sqrt(10)
pnorm(2385, mean10, sd10)
## [1] 0.03059796
meanDiff <- 0
varDiff <- varCats * 2
pnorm(23.85, meanDiff, sqrt(varDiff), lower.tail = FALSE) * 2
## [1] 0.5056741
# Simulation
simSize <- 1e6
cat1 <- rnorm(simSize, meanCats, sdCats)
cat2 <- rnorm(simSize, meanCats, sdCats)
diffs <- cat2 - cat1
length(diffs[diffs > 23.85 | diffs < -23.85])/simSize
## [1] 0.505616