coin3 <- data.frame( N_Heads = c(0, 1, 2, 3), C_Heads = c(1, 3, 3, 1) ) attach(coin3) coin3$P_Heads <- C_Heads/sum(C_Heads)

normal_dist <- function(x, b){ mu = mean(x) sigma = sd(x) x1 <- seq(mu - 6*sigma, mu + 6*sigma, length = 100) # Normal curve fun <- dnorm(x1, mean = mu, sd = sigma) # Histogram options(repr.plot.width=7, repr.plot.height=4) hist(x, prob = TRUE, col = "white", breaks = b, xlim = c(mu - 6*sigma, mu + 6*sigma), ylim = c(0, max(fun)), main = "Histogram overlayed with normal curve") lines(x1, fun, col = 2, lwd = 2) }

data(chickwts) attach(chickwts) normal_dist(chickwts$weight, 10)

shapiro.test(chickwts$weight)

chickwts$zscore = (weight-mean(weight))/sd(weight) tail(chickwts)