# R Notes from October 8th class


pnorm(15)
pnorm(15, mean=0, sd=10)
pnorm(1.5)

heights <- c(180.34, 170.18, 175.26, 177.8, 172.72, 160.02, 172.72, 182.88,
			 177.8, 177.8, 167.64, 180.34, 180.34, 172.72, 165.1, 154.94,
			 180.34, 172.72, 165.1, 167.64, 182.88, 175.26, 182.88, 177.8,
			 175.26, 185.42, 175.26, 167.64, 187.96, 175.26, 180.34, 175.26,
			 198.12, 177.8, 185.42, 175.26, 180.34, 187.96, 182.88, 187.96,
			 177.8, 182.88, 187.96, 170.18, 182.88, 182.88, 175.26, 170.18,
			 182.88, 180.34, 180.34, 170.18, 180.34, 187.96, 193.04, 175.26,
			 193.04, 182.88, 177.8, 167.64, 170.18, 160.02, 172.72, 193.04,
			 187.96, 190.5, 172.72, 175.26, 193.04, 180.34, 162.56, 187.96,
			 182.88, 180.34, 177.8, 172.72, 185.42, 180.34, 180.34, 182.88,
			 185.42, 180.34, 195.58, 185.42, 170.18, 170.18, 172.72, 180.34,
			 190.5, 172.72, 182.88, 170.18, 177.8, 175.26, 162.56, 162.56,
			 175.26, 167.64, 170.18, 177.8)/2.54

ggplot(data.frame(heights = heights), aes(heights)) + geom_histogram(bins = 10)
ggplot(data.frame(heights = heights), aes(heights)) + geom_density()

DATA606::qqnormsim(heights)

.65^9 * .35

n <- 4
p <- .35
dbinom(0, n, p)
dbinom(1, n, p)
dbinom(2, n, p)
dbinom(3, n, p)
dbinom(4, n, p)