Harvard University Statistics Lab Homework Lab Homework #save your plots as an image file and upload separately (export > save as image) #save them as pn
Harvard University Statistics Lab Homework Lab Homework
#save your plots as an image file and upload separately (export > save as image)
#save them as png with filename your_name_Lab3_plot
#Exercise 1
#The rain probability for the month of July in the city of Portland is 72%
#(a) What is the probability that it rains exactly 11 days
#(b) What is the probability that it rains 11 or fewer days?
#(c) What is the probability that it rains 27 or more days?
#(d) What is the probability that it rains more than 20 but less than 25 days?
#Exercise 2
#Suppose the winnings of gamblers at Las Vegas are normally distributed with
#mean $670 and standard deviation $38.
#(a) What is the probability of winning exactly $500?
#(b) What is the probability of winning $500 or less?
#(c) What is the probability of winning more than $800?
#(d) What is the probability of winning between $500 to $875?
#(e) Generate 10000 random numbers using rnorm function. Show how this
#distribution would look like. Add in the plot the mean, median and 20%
#trimmed mean (provide these values).
#Exercise 3
#In the previous problem, how much someone has to win to be in the top 5%?
——————————————————————-
LAB 3 EXAMPLE
#Lab 3-Contents
#1. Binomial Probability Distribution
#2. Normal Probability Distribution
#————————————–
# 1. Binomial Probability Distribution
#————————————–
# The binomial distribution is one of the most important distributions in statisitcs.
# Like it’s name suggests, bi-nomial, indicates that we are dealing with
# numbers that have only two possibilites (eg. Yes/No, Dead/Alive etc.)
#There are three main types of problems dealing with binomial probabilites:
#A) Finding the probability of EXACTLY x successes (Yes responses, deaths, boys)
#B) Finding the probability of < x successes OR > x successes
#C) Finding the probability of < x1 successes AND > x2 successes
# Let’s try one by semi-hand (we’ll have R do some computations)
# and then learn how to use R to solve these kinds of problems.
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXAMPLE 1: Let’s say the probability for divorce is 45%. Out of 20 married couples,
#what is the probability of EXACTLY 5 of them getting divorced?
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
N=20 #n
K=5 # x
p=0.45
q=1-p
Nf=factorial(N) #n!
Kf=factorial(K) # x!
Nf/(Kf*factorial(N-K)) * p^K * (1-p)^(N-K)
#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^#
#Binomial Probability (EXACT): dbinom(k, n, p)
#Where k: number of successes, n: number of trials,
#p: probability of success
#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^#
#Let’s use the dbinom() command from R to answer this question
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXERCISE 1-1: Use the dbinom() command to solve the problem in Example 1
#
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
dbinom(5, 20, 0.45)
# What if we wanted to know the probability of 5 OR fewer divorces?
# One strategey is to compute the probability for 5, 4, 3, 2, 1, 0 divorces and sum the probabilities
p5=dbinom(5, 20, 0.45)
p4=dbinom(4, 20, 0.45)
p3=dbinom(3, 20, 0.45)
p2=dbinom(2, 20, 0.45)
p1=dbinom(1, 20, 0.45)
p0=dbinom(0, 20, 0.45)
p0+p1+p2+p3+p4+p5
# This is a bit tedious, so let’s just use a new command called pbinom()
#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^#
#Binomial Probability (x <= k): pbinom(k, n, p) #Where k: <= number of successes, n: number of trials, #p: probability of success #^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^# #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# #EXERCISE 1-2: Use the pbinom() command to solve the problem of computing the probability for # 5 or fewer divorces. # #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# pbinom(5, 20, 0.45) #??????????????????????????????????????????????????????????????# #Thought Question 1: How would the command change if the question # was worded as "Find the probabilty for less than 5 divorces"? pbinom(4, 20, 0.45) #??????????????????????????????????????????????????????????????# # The pbinom() command is useful, but only gives us probabilities for <= k # What if we wanted to know the probability of having >= k responses or > k?
# Because the probabilities have to add to exactly 1, we can use some simple math to figure this out
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXAMPLE 2: The probability for divorce is 45%. Out of 20 married couples,
#what is the probability of >= 6 of them getting divorced?
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
# We know from above that the probability of <= 5 couples divorcing is 0.05533419 # If we think about it, then the probability of >=6 is just 1-0.05533419 = 0.9446658
# We could even use the pbinom() command with this by putting the 1 minus in front
1-pbinom(5, 20, 0.45)
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXERCISE 1-3: Use the pbinom() command to solve the problem of computing the probability for
# > 7 divorces.
#
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
1-pbinom(7, 20, 0.45)
# Thus far we’ve covered how to answer 2 out of the 3 possible binomial probability probelms (A and B).
# Let’s think about how we could answer a question asking about a range of possible sucesses.
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXAMPLE 3: The probability for divorce is 45%. Out of 20 married couples,
#what is the probability of >= 7 and <= 10 divorces? # 7=< p<=10 ? #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# # This problem has 2 parts. #1) Compute the probability for <=10 divorces pbinom(10, 20, 0.45) #2) Compute the probability for <= 6 divorces (YES, 6 not 7) pbinom(6, 20, 0.45) # The difference between these two probabilites is the probability between 7 and 10 pbinom(10, 20, 0.45) - pbinom(6, 20, 0.45) #Note: Probability will always be positive, so make sure you do the subtraction in the correct order #??????????????????????????????????????????????????????????????# #Thought Question 2: How would computing the probability # in example 3 change if the direction of the signs were reversed? # eg. <=7 & >=10
#p<=7 and p>=10
#??????????????????????????????????????????????????????????????#
pbinom(7, 20, 0.45)
1-pbinom(9, 20, 0.45)
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXERCISE 1-4: Use the pbinom() command to solve the problem of computing the probability for
# > 4 AND <= 15 divorces. 4
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
1-pnorm(1.96, 0, 1)
#We can generalize this to normal distributions that do not have a mean of 0 and SD of 1
#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#
#EXAMPLE 5: What is the probability of having a bmi < 25.2 if the mean bmi is 35 and SD = 5 # #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# #We can use pnorm() in this situation as well #1) pnorm(25.2, 35, 5) # 2) Z=(25.2-35)/5 pnorm(Z,0,1) #Notice that the answer we get: 0.0249979 is the same as the answer we got in Exercise 2-1. #??????????????????????????????????????????????????????????????# #Thought Question 3: Why do you think the answers are the same in # example 5 and exercise 2-1? #??????????????????????????????????????????????????????????????# (25.2-35)/5 #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# #EXERCISE 2-2: The average age for incoming college freshman is 18.2 years with SD=0.5. # What is the probability that a freshman will be between 18 and 19 years old? #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# pnorm(19, 18.2, .5) - pnorm(18, 18.2, .5) # We can also do these problems in reverse, where we are given the probability and have to find the values # at which that probability exists. When we do this, we are calculating the quantile of the normal distribution. #^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^# #Normal Dist Probability (Quantile): qnorm(p, mean, sd) #^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^# # Let's think about the problem we did in example 5 # We found the probability of a bmi < 25.2 when the mean of the distribution was 35, SD=5 # The probability was 0.0249979 # Let's see if qnorm() returns the correct result: qnorm(0.0249979, 35, 5) #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# #EXERCISE 2-3: The average age for incoming college freshman is 18.2 years with SD=0.5. # A researcher tells you that the 65% of the freshman were less than... # He can't seem to remember what age it was at that 65% of the freshman were less than, help him out. #*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*#*# qnorm(0.65, 18.2, .5) #????????????????????????????????????????????????????????????????????????????????????# #Thought Question 4: What if instead of knowing the age for the bottom 65% # we wanted to know how old a person needed to be in the top 10% of age in their class #?????????????????????????????????????????????????????????????????????????????????????????# qnorm(0.9, 18.2, .5)
