Lecture 2
Lecture 2
Random Variable
A random variable x takes on a defined set of values with different probabilities.
For example, if you roll a die, the outcome is random (not fixed) and there are 6 possible outcomes, each of which occur with probability one-sixth. For example, if you poll people about their voting preferences, the percentage of the sample that responds Yes on Proposition 100 is a also a random variable (the percentage will be slightly differently every time you poll).
Roughly, probability is how frequently we expect different outcomes to occur if we repeat the experiment over and over (frequentist view)
Examples: blood pressure, weight, the speed of a car, the real numbers from 1 to 6.
Probability functions
A probability function maps the possible values of x against their respective probabilities of occurrence, p(x) p(x) is a number from 0 to 1.0. The area under a probability function is always 1.
1/6 x
P(x) ! 1
all x
1 2
5 6
Practice Problem:
The number of patients seen in the ER in any given hour is a random variable represented by x. The probability distribution for x is:
x P(x)
10 .4
11 .2
12 .2
13 .1
14 .1
Review Question 1
If you toss a die, whats the probability that you roll a 3 or less? a. b. c. d. e. 1/6 1/3 1/2 5/6 1.0
Review Question 1
If you toss a die, whats the probability that you roll a 3 or less? a. b. c. d. e. 1/6 1/3 1/2 5/6 1.0
Review Question 2
Two dice are rolled and the sum of the face values is six? What is the probability that at least one of the dice came up a 3? a. b. c. d. e. 1/5 2/3 1/2 5/6 1.0
Review Question 2
Two dice are rolled and the sum of the face values is six. What is the probability that at least one of the dice came up a 3? a. b. c. d. e. 1/5 2/3 1/2 5/6 1.0
How can you get a 6 on two dice? 1-5, 5-1, 2-4, 4-2, 3-3 One of these five has a 3. @1/5
Continuous case
The probability function that accompanies a continuous random variable is a continuous mathematical function that integrates to 1.
For example, recall the negative exponential function (in probability, this is called an exponential distribution ): f ( x) ! e x
g
e x ! e x
g 0
! 0 1!1
The probability that x is any exact particular value (such as 1.9976) is 0; we can only assign probabilities to possible ranges of x.
p(x)=e-x 1
x 1 2
P(1 e x e 2) ! e
1
x
! e
x
2 1
We can see its a probability distribution because it integrates to 1 (the area under the curve is 1): 1 1
1 ! x
0
!1 0 !1
P( uxu 0)=
All probability distributions are characterized by an expected value (mean) and a variance (standard deviation squared).
Expected value is just the average or mean () of random variable x. It s sometimes called a weighted average because more frequent values of X are weighted more highly in the average. It s also how we expect X to behave on-average over the long run ( frequentist view again).
E( X ) !
Continuous case:
x p(x )
i i all x
E( X ) !
xi p(xi )dx
all x
Symbol Interlude
E(X) =
x
X!
i !1
i !1
1 xi ( ) n
Expected Value
The Lottery (also known as a tax on people who are bad at math) A certain lottery works by picking 6 numbers from 1 to 49. It costs $1.00 to play the lottery, and if you win, you win $2 million after taxes. If you play the lottery once, what are your expected winnings or losses?
Lottery
Calculate the probability of winning in 1 try:
1 49 6 ! 1 1 ! ! 7.2 x 10 -8 49! 13,983,816 43!6!
Expected Value
The probability function
x$ -1 + 2 million p(x) .999999928 7.2 x 10--8
Expected Value E(X) = P(win)*$2,000,000 + P(lose)*-$1.00 = 2.0 x 106 * 7.2 x 10-8+ .999999928 (-1) = .144 - .999999928 = -$.86 Negative expected value is never good! You shouldnt play if you expect to lose money!
Expected Value
If you play the lottery every week for 10 years, what are your expected winnings or losses? 520 x (-.86) = -$447.20
Gambling (or how casinos can afford to give so many free drinks )
A roulette wheel has the numbers 1 through 36, as well as 0 and 00. If you bet $1 that an odd number comes up, you win or lose $1 according to whether or not that event occurs. If random variable X denotes your net gain, X=1 with probability 18/38 and X= -1 with probability 20/38. E(X) = 1(18/38) 1 (20/38) = -$.053
On average, the casino wins (and the player loses) 5 cents per game. The casino rakes in even more if the stakes are higher: E(X) = 10(18/38) 10 (20/38) = -$.53
If the cost is $10 per game, the casino wins an average of 53 cents per game. If 10,000 games are played in a night, that s a cool $5300.
Take the hit new show Deal or No Deal Everyone know the rules? Let s say you are down to two cases left. $1 and $400,000. The banker offers you $200,000. So, Deal or No Deal?
Deal or No Deal
x$ +$200,000
p(x) 1.0
x$ +$200,000
Q ! E ( X ) ! 200,000
p(x) 1.0
How to decide?
Variance! If you take the deal, the variance/standard deviation is 0. If you don t take the deal, what is average deviation from the mean? What s your gut guess?
Variance/standard deviation
W2=Var(x) =E(x-Q)2 The expected (or average) squared distance (or deviation) from the mean
W 2 ! Var ( x) ! E[( x Q ) 2 ] !
all x
( xi Q ) 2 p(xi )
Variance, continuous
Discrete case:
Var ( X ) !
(x
all x
Q ) p(xi )
2
Continuous case?:
Var ( X ) !
( xi Q ) p(xi )dx
2
all x
Symbol Interlude
Var(X)= W2 SD(X) = W
i !1
( xi x ) 2 n 1
1 ! ( xi x ) ( ) n 1 i !1
Division by n-1 reflects the fact that we have lost a degree of freedom (piece of information) because we had to estimate the sample mean before we could estimate the sample variance.
Variance
W !
2
(x
all x
Q ) p(xi )
2
W2 !
all x
( x i Q ) 2 p(xi ) !
2 2 2
Practice Problem
On the roulette wheel, X=1 with probability 18/38 and X= -1 with probability 20/38.
Answer
W !
2
(x Q )
i
p(xi )
! (1.053) 2 (18 / 38) (1 .053) 2 (20 / 38) ! (1.053) 2 (18 / 38) (.947) 2 (20 / 38) ! .997
W ! .997 ! .99
Standard deviation is $.99. Interpretation: On average, you re either 1 dollar above or 1 dollar below the mean, which is just under zero. Makes sense!
Review Question 3
The expected value and variance of a coin toss (H=1, T=0) are?
a. b. c. d.
Review Question 3
The expected value and variance of a coin toss are?
a. b. c. d.
e.g., 15 tosses of a coin; 20 patients; 1000 people surveyed e.g., head or tail in each toss of a coin; disease or no disease Generally called success and failure Probability of success is p, probability of failure is 1 p e.g., Probability of getting a tail is the same each time we toss the coin
A binary outcome
Binomial distribution
Take the example of 5 coin tosses. What s the probability that you flip exactly 3 heads in 5 coin tosses?
Binomial distribution
Solution: One way to get exactly 3 heads: HHHTT What s the probability of this exact arrangement? P(heads)xP(heads) xP(heads)xP(tails)xP(tails) =(1/2)3 x (1/2)2 Another way to get exactly 3 heads: THHHT Probability of this exact outcome = (1/2)1 x (1/2)3 x (1/2)1 = (1/2)3 x (1/2)2
Binomial distribution
In fact, (1/2)3 x (1/2)2 is the probability of each unique outcome that has exactly 3 heads and 2 tails. So, the overall probability of 3 heads and 2 tails is: (1/2)3 x (1/2)2 + (1/2)3 x (1/2)2 + (1/2)3 x (1/2)2 + .. for as many unique arrangements as there are but how many are there??
5 3
5C3
Outcome Probability THHHT (1/2)3 x (1/2)2 HHHTT (1/2)3 x (1/2)2 TTHHH (1/2)3 x (1/2)2 HTTHH (1/2)3 x (1/2)2 HHTTH (1/2)3 x (1/2)2 HTHHT (1/2)3 x (1/2)2 THTHH (1/2)3 x (1/2)2 HTHTH (1/2)3 x (1/2)2 HHTHT (1/2)3 x (1/2)2 THHTH (1/2)3 x (1/2)2 10 arrangements x (1/2)3 x (1/2)2
The probability of each unique outcome (note: they are all equal)
= 5!/3!2! = 10
5 3
x P(heads)3 x P(tails)2 =
10 x ()5=31.25%
number of heads
n X n X p (1 p ) X
X=# successes out of n trials p= probability of success
If I toss a coin 20 times, what s the probability of getting of getting 2 or fewer heads?
20! 20 0 20 (.5) 20 ! 9.5 x10 7 (.5) (.5) ! 20!0! 0 20! 20 1 19 (.5) 20 ! 20 x9.5 x10 7 ! 1.9 x10 5 (.5) (.5) ! 19! ! 1 1 20! 20 2 18 (.5) 20 ! 190 x9.5 x10 7 ! 1.8 x10 4 (.5) (.5) ! 18!2! 2 ! 1.8 x10 4
If X follows a binomial distribution with parameters n and p: X ~ Bin (n, p) Then: Note: the variance will E(X) = np always lie between 0*N-.25 *N Var (X) = np(1-p) p(1-p) reaches maximum at p=.5 np (1 p ) SD (X)= P(1-p)=.25
Practice Problem
1. You are performing a cohort study. If the probability of developing disease in the exposed group is .05 for the study duration, then if you (randomly) sample 500 exposed people, how many do you expect to develop the disease? Give a margin of error (+/- 1 standard deviation) for your estimate. 2. What s the probability that at most 10 exposed people develop the disease?
Answer
1. How many do you expect to develop the disease? Give a margin of error (+/- 1 standard deviation) for your estimate.
X ~ binomial (500, .05) E(X) = 500 (.05) = 25 Var(X) = 500 (.05) (.95) = 23.75 StdDev(X) = square root (23.75) = 4.87 @25 s 4.87
Answer
2. What s the probability that at most 10 exposed subjects develop the disease?
This is asking for a CUMULATIVE PROBABILITY: the probability of 0 getting the disease or 1 or 2 or 3 or 4 or up to 10. P(X10) = P(X=0) + P(X=1) + P(X=2) + P(X=3) + P(X=4)+.+ P(X=10)=
500 500 500 500 0 500 1 499 2 498 10 490 .01 (.05) (.95) (.05) (.95) (.05) (.95) ... (.05) (.95) 10 2 1 0
Practice Problem:
You are conducting a case-control study of smoking and lung cancer. If the probability of being a smoker among lung cancer cases is .6, whats the probability that in a group of 8 cases you have:
a. b. c.
Less than 2 smokers? More than 5? What are the expected value and variance of the number of smokers?
Answer
X 0 1 2 3 4 5 6 7 8 P(X) 8 1(.4) =.00065 1 7 8(.6) (.4) =.008 28(.6)2 (.4) 6 =.04 56(.6)3 (.4) 5 =.12 70(.6)4 (.4) 4 =.23 56(.6)5 (.4) 3 =.28 28(.6)6 (.4) 2 =.21 8(.6)7 (.4) 1=.090 1(.6)8 =.0168
0 1 2 3 4 5 6 7 8
Answer, continued
P(<2)=.00065 + .008 = .00865 P(>5)=.21+.09+.0168 = .3168
0 1 2 3 4 5 6 7 8
Review Question 4
In your case-control study of smoking and lung-cancer, 60% of cases are smokers versus only 10% of controls. What is the odds ratio between smoking and lung cancer?
a. b. c. d. e.
Review Question 4
In your case-control study of smoking and lung-cancer, 60% of cases are smokers versus only 10% of controls. What is the odds ratio between smoking and lung cancer?
a. b. c. d. e.
.6 .4 ! 3 x 9 ! 27 ! 13.5 .1 2 1 2 .9
Review Question 5
What s the probability of getting exactly 5 heads in 10 coin tosses?
a.
10 5 5 (.50) (.50) 0
10 b. (.50) 5 (.50) 5
c.
10 10 5 (.50) (.50) 5 10
d. 10 10 0 (.50) (.50)
Review Question 5
What s the probability of getting exactly 5 heads in 10 coin tosses?
a.
10 5 5 (.50) (.50) 0
10 b. (.50) 5 (.50) 5
c.
10 10 5 (.50) (.50) 5 10
d. 10 10 0 (.50) (.50)
Review Question 6
A coin toss can be thought of as an example of a binomial distribution with N=1 and p=.5. What are the expected value and variance of a coin toss?
a. b. c. d. e.
Review Question 6
A coin toss can be thought of as an example of a binomial distribution with N=1 and p=.5. What are the expected value and variance of a coin toss?
a. b. c. d. e.
Review Question 7
If I toss a coin 10 times, what is the expected value and variance of the number of heads?
a. b. c. d. e.
Review Question 7
If I toss a coin 10 times, what is the expected value and variance of the number of heads?
a. b. c. d. e.
Review Question 8
In a randomized trial with n=150, the goal is to randomize half to treatment and half to control. The number of people randomized to treatment is a random variable X. What is the probability distribution of X?
a. b. c. d. e.
Review Question 8
In a randomized trial with n=150, every subject has a 50% chance of being randomized to treatment. The number of people randomized to treatment is a random variable X. What is the probability distribution of X?
a. b. c. d. e.
Review Question 9
In the same RCT with n=150, if 69 end up in the treatment group and 81 in the control group, how far off is that from expected?
a. b. c. d.
Less than 1 standard deviation 1 standard deviation Between 1 and 2 standard deviations More than 2 standard deviations
Review Question 9
In the same RCT with n=150, if 69 end up in the treatment group and 81 in the control group, how far off is that from expected? Expected = 75
81 and 69 are both 6 away from the expected.
a. b. c. d.
Less than 1 standard deviation 1 standard deviation Between 1 and 2 standard deviations More than 2 standard deviations
Variance = 150(.25) = 37.5 Std Dev $ 6 Therefore, about 1 SD away from expected.
Proportions
The binomial distribution forms the basis of statistics for proportions. A proportion is just a binomial count divided by n.
For example, if we sample 200 cases and find 60 smokers, X=60 but the observed proportion=.30.
Statistics for proportions are similar to binomial counts, but differ by a factor of n.
Differs by a factor of n.
For proportion:
P-hat stands for sample proportion.
Qp ! p np (1 p) p(1 p ) Wp ! ! 2 n n p(1 p ) Wp ! n
2
Differs by a factor of n.
Statistics for proportions are based on a normal distribution, because the binomial can be approximated as normal if np>5
p(X)
Examples:
height weight age bone density IQ (mean=100; SD=15) SAT scores blood pressure ANYTHING YOU AVERAGE OVER A LARGE ENOUGH # (Central Limit Theorem)
f ( x) !
1 W 2T
1 xQ 2 ) ( 2 W e
This is a bell shaped curve with different centers and spreads depending on Q and W
Q W
1 W 2T 1 W 2T 1 W 2T
1 xQ 2 ( ) y e 2 W dx
Q 2W
Q 2W Q 3W
1 xQ 2 ( ) y e 2 W dx 1 xQ 2 ( ) y e 2 W dx
Q 3W
For example, whats the probability of getting a math SAT score below 575 if SAT scores are normally distributed with a mean of 500 and a std dev of 50??
575
@ P ( X e 575) !
200
(50)
1 2T
1 x 500 2 ( ) e 2 50 dx
Solve this?.Yikes!
X Q Z! W
For example, 575 in math SAT units translates to 1.5 standard deviations above the mean.
f (Z ) !
1 2T
500 0
575 1.5
X Z
(Q = 500, W = 50) (Q = 0, W = 1)
Example
So, Whats the probability of getting a math SAT score of 575 or less, Q=500 and W=50?
@ P ( X e 575) !
200
(50)
1 2T
1 x 500 2 ) ( 2 50 dx e
1.5
g
1 2T
1 Z2 e 2 dz
No need to do the integration! Just look up Z= 1.5 in standard normal chart no problem! = .9332
Z=1.50
a.
b.
Answer
a.
What is the chance of obtaining a birth weight of 141 oz or heavier when sampling birth records at random?
141 109 Z! ! 2.46 13
From the chart Z of 2.46 corresponds to a right tail (greater than) area of: P(Z2.46) = 1-(.9931)= .0069 or .69 %
Answer
b. What is the chance of obtaining a birth weight of 120 or lighter?
120 109 Z! ! .85 13
From the chart Z of .85 corresponds to a left tail area of: P(Z.85) = .8023= 80.23%
Not all continuous random variables are normally distributed!! It is important to evaluate how well the data are approximated by a normal distribution, because many statistics (ttest, ANOVA, linear regression) assume that the outcome variable is normally distributed
-1.6
17.8
-11.4
27.5
27
Review Problem 10
Which of the following about the normal distribution is NOT true?
a. Theoretically, the mean, median, and mode are the same. b. About 2/3 of the observations fall within 1 standard deviation from the mean. c. It is a discrete probability distribution. d. Its parameters are the mean, Q , and standard deviation, W.
Review Problem 10
Which of the following about the normal distribution is NOT true?
a. Theoretically, the mean, median, and mode are the same. b. About 2/3 of the observations fall within 1 standard deviation from the mean. c. It is a discrete probability distribution. d. Its parameters are the mean, Q , and standard deviation, W.
Review Problem 11
For some positive value of Z, the probability that a standard normal variable is between 0 and Z is 0.3770. The value of Z is:
a. b. c. d. 0.18. 0.81. 1.16 1.47.
Review Problem 11
For some positive value of Z, the probability that a standard normal variable is between 0 and Z is 0.3770. The value of Z is:
a. b. c. d. 0.18. 0.81. 1.16 1.47.
Review Problem 12
The probability that a standard normal variable Z is positive is ________.
a. b. c. d. 50% 100% 0% 95%
Review Problem 12
The probability that a standard normal variable Z is positive is ________.
a. b. c. d. 50% 100% 0% 95%
Review Problem 13
Suppose Z has a standard normal distribution with a mean of 0 and a standard deviation of 1. The probability that Z values are larger than __________ is 0.6985.
a. b. c. d. e. 1.0 0 -0.6 +0.6 -2.0
Review Problem 13
Suppose Z has a standard normal distribution with a mean of 0 and a standard deviation of 1. The probability that Z values are larger than __________ is 0.6985.
a. b. c. d. e. 1.0 0 -0.6 +0.6 -2.0
Homework