Lecture 04: Statistical Models: Dr. Nguyen Tai Hung
Lecture 04: Statistical Models: Dr. Nguyen Tai Hung
Lecture 04: Statistical Models: Dr. Nguyen Tai Hung
09 October 2010
Goals for Today
2
2
Topics
3
3
Discrete and Continuous
Random Variables
4
Discrete versus Continuous Random Variables
Discrete Random Variable Continuous Random Variable
3. f ( x) 0, if x is not in RX
p X x x
x
x
p (x i ) p X x f t dt 0
i
b
p a X b f x dx
a
5
Expectation
E ( x) xi p( xi )
all i
E ( x) xf ( x)dx
2
n n n
V X x i 2 p xi xi
p x i x i p x i
2
x
i 1 i 1 i 1
V X x x f x dx x 2 f x dx
2
x2
6
Statistical Models: Properties
7
Discrete Probability Distributions
Bernoulli Trials
Binomial Distribution
Geometric Distribution
Poisson Distribution
Poisson Process
8
Modeling of Random Events with
Two-States
Bernoulli Trials
Binomial Distribution
9
Bernoulli Trials
Context: Random events with two possible values
Two events: Yes/No, True/False, Success/Failure
Two possible values: 1 for success, 0 for failure.
Example: Tossing a coin, Packet Transmission Status,
An experiment comprises n trial.
Probability Mass Function (PMF): Probability in one trial
p, x j 1, j 1, 2,..., n
PMF: p one trial p (x j ) 1 p q , x j 0 , j 1, 2 ,..., n
0, otherwise
n k n k n n!
P X k p 1 p with
k k k ! n k !
where
k = 0, 1, 2, ......., n k: number of success
n = 1, 2, 3, ....... n: number of trials
0<p<1 p = success probability
11
Binomial Distribution
Probability of k=2 successes in n=3 trials?
E1: 1 1 0 P(E1)= p(1 and 1 and 0) = p p (1-p) = p2 (1-p)
or
E2: 1 0 1 P(E2)= p(1 and 0 and 1) = p (1-p) p = p2 (1-p)
or
E3: 0 1 1 P(E3)= p(1 and 1 and 0) = p p (1-p) = p2 (1-p)
Administrative issues
• Groups Formation
13
Modeling of Discrete Random
Time
Geometric Distribution
14
Geometric Distribution
Context: the number of Bernoulli trials until achieving the
FIRST SUCCESS.
It is used to represent random time until a first transition occurs.
PMF
q k 1 p , k 0,1, 2,..., n
PMF: p (X k )
0, otherwise
CDF: F X p X k 1 1 p
k
1
Expected Value : E X
p k
q 1 p
Variance :V X 2
2
p p2
15
Modeling of Random Number of
Arrivals/Events
Poisson Distribution
Poisson Process
16
Poisson Distribution
Context: number of events occurring in a fixed period of time
Events occur with a known average rate and are independent
Possion distribution is characterized by the average rate
The average number of arrival in the fixed time period.
Examples
The number of cars passing a fixed point in a 5 minute
interval. Average rate: = 3 cars/5 minutes
The number of calls received by a switchboard during a
given period of time. Average rate: =3 call/minutes
The number of message coming to a router per second
The number of travelers arriving to the airport for flight
registration
17
Poisson Distribution
The Poisson distribution with the average rate parameter
Variance: V X
CDF
18
Example: Poisson Distribution
The probability of having exactly 15 cars entering the parking in one hour:
2015
p 15 P X 15 exp 20 0.051649
15!
or
p 15 F 15 F 14 0.156513 0.104864 0.051649
The probability of having more than 3 cars entering the parking in one hour:
p X 3 1 p X 3 1 F 3
1 p 0 p 1 p 2 p 3
USE EXCEL/MATLAB
0.9999967 FOR COMPUTATIONS
19
Example: Poisson Distribution
20k k
20i
pX k exp 20
k! Fk pX k
i 0
i!
exp 20
20
Five Minutes Break
You are free to discuss with your classmates about the previous slides, or to
refresh a bit, or to ask questions.
Administrative issues
• Groups Formation
21
Modeling of Random Number of
Arrivals/Events
Poisson Distribution
Poisson Process
22
Poisson Process
Process N(t)
random variable N that depends on time t
Poisson Process: a Process that has a Poisson Distribution
N(t) is called a counting poisson process
There are two types:
Homogenous Poisson Process: constant average rate
Non-Homogenous Poissoon Process: variable average rate
23
What is “Stochastic Process”?
State Space = {SUNNY,
RAINNY}
X day i " S " or " R " : RANDOM VARIABLE that varies with the DAY
X day 2 " S " X day 4 " S " X day 6 " S "
X day 1 " S " X day 3 " R " X day 5 " R " X day 7 "S "
Day
Day 1 Day 2 Day 3 Day 4 Day 5 Day 6 Day 7
THU FRI SAT SUN MON TUE WED
25
(Homogenous) Poisson Process
N t N t
describes the number of events in time interval t ,t
The mean and the variance are equal
E N V N
26
(Homogenous) Poisson Process
Properties of Poisson process
Arrivals occur one at a time (not simultaneous)
N , 0has stationary increments, which means
N t N s N t s
The number of arrivals in time s to t is also Poisson-
distributed with mean
has independent t sincrements
N , 0
27
Inter-Arrival Times of a Poisson Process
Inter-arrival time: time between two consecutive arrivals
The inter-arrival times of a Poisson process are random.
What is its distribution?
Consider the inter-arrival times of a Poisson process (A1, A2, …), where Ai is
the elapsed time between arrival i and arrival i+1
The first arrival occurs after time t MEANS that there are no
arrivals in the interval [0,t], As a consequence:
p A1 t p N t 0 exp t
p A1 t 1 p A1 t 1 exp t
28
The Inter-arrival times of a Poisson process are 28
exponentially distributed and independent with mean 1/
Splitting and Pooling
Splitting
p N1(t) ~ Poi[p]
N(t) ~ Poi()
Pooling
N1(t) ~ Poi[] 1
N(t) ~ Poi(12)
N2(t) ~ Poi[] 2
29
Modeling of Random Number of
Arrivals/Events
Poisson Distribution
Non Homogenous Poisson Process
30
Non Homogenous (Non-stationary) Poisson Process
(NSPP)
The number of cars that cross the intersection of King Fahd Road and Al-Ourouba
Road is distributed according to a non homogenous Poisson process with a mean
(t) defined as follows:
80 cars/mn if 8 am t 9am
60 cars/mn if 9am t 11pm
t
50 car/mn if 11am t 15 pm
70 car/mn if 15 pm t 17 pm
Let us consider the time 8 am as t=0.
Q1. Compute the average arrival number of cars at 11H30?
Q2. Determine the equation that gives the probability of having only 10000 car
arrivals between 12 pm and 16 pm.
Q3. What is the distribution and the average (in seconds) of the inter-arrival time
of two cars between 8 am and 9 am?
32
Example: Non-stationary Poisson Process
(NSPP)
You are free to discuss with your classmates about the previous slides, or to
refresh a bit, or to ask questions.
Administrative issues
• Groups Formation
34
Continuous Probability Distributions
Uniform Distribution
Exponential Distribution
Normal (Gaussian) Distribution
35
Uniform Distribution
36
Continuous Uniform Distribution
The continuous uniform distribution is a family of probability
distributions such that for each member of the family, all intervals of
the same length on the distribution's support are equally probable
A random variable X is uniformly distributed on the interval [a,b],
U(a,b), if its PDF and CDF are:
0, x a
1 x a
, a x b
PDF: f (x ) b a CDF: F ( x ) , a x b
0, b a
otherwise
1, x b
Expected value: E X
a b a b 2
Variance: V X
2 12
37
Uniform Distribution
Properties
p x 1 X x 2is proportional to the
length of the interval
X 2 X1
F X 2 F X 1
b a
CDF
Special case: a standard uniform
distribution U(0,1).
Very useful for random number
generators in simulators
38
Exponential Distribution
39
Exponential Distribution
0, x 0 0, x 0
CDF: F (x ) x t
F (x ) x
0 e dt 1 e , x 0
x
1 exp , x 0
1 1
Expected value: E X Variance: V X 2
2 40
Exponential Distribution
µ=20 µ=20
1 x 0, x 0
exp , x 0
f (x ) 20 20 F (x ) x
0, 1 exp , x 0
otherwise 20
41
Exponential Distribution
The memoryless property: In probability theory, memoryless is a property of
certain probability distributions: the exponential distributions and the
geometric distributions, wherein any derived probability from a set of random
samples is distinct and has no information (i.e. "memory") of earlier samples.
Formally, the memoryless property is:
For all s and t greater or equal to 0:
p X s t | X s p X t
This means that the future event do not depend on the past event, but only on
the present event
The fact that Pr(X > 40 | X > 30) = Pr(X > 10) does not mean that the events X > 40 and X >
30 are independent; i.e. it does not mean that
Pr(X > 40 | X > 30) = Pr(X > 40).
42
Exponential Distribution
p X s t | X s p X t
The fact that Pr(X > 40 | X > 30) = Pr(X > 10) does not mean that the events
X > 40 and X > 30 are independent; i.e. it does not mean that
Pr(X > 40 | X > 30) = Pr(X > 40).
43
Example: Exponential Distribution
The probability that the repair time lasts for another hour given it has been
operating for 2.5 hours:
Using the memoryless property of the exponential distribution, we have:
1
p X 2.5 1 | X 2.5 p X 1 1 p X 1 exp 0.717
3
44
Normal (Gaussian) Distribution
45
Normal Distribution
46
Normal or Gaussian Distribution
48
Normal Distribution
Evaluating the distribution:
Independent of and using the standard normal distribution:
Z ~ N 0,1 z 1 t 2 / 2
, where ( z ) e dt
2
X
Transformation of variables: let Z
x
F ( x ) P X x P Z
( x ) / 1 z2 / 2
e dz
2
( x ) /
( z )dz ( x )
49
Normal Distribution
Example: The time required to load an oceangoing vessel, X, is
distributed as N(12,4)
The probability that the vessel is loaded in less than 10 hours:
10 12
F (10) (1) 0.1587
2
Using the symmetry property, (1) is the complement of (-1)
50
Empirical Distribution
51
Empirical Distributions
52
Empirical Distributions
In statistics, an empirical distribution function is a cumulative
probability distribution function that concentrates probability 1/n
at each of the n numbers in a sample.
Let X1, X2, …, Xn be iid random variables in with the CDF equal to
F(x).
The empirical distribution function Fn(x) based on sample X1, X2, …,
Xn is a step function defined by
n
number of element in the sample x 1
Fn x
n
n I X
i 1
i x
1 if X i x
where I(A) is the indicator of event A. I Xi x
0 otherwise
For a fixed value x, I(Xi≤x) is a Bernoulli random variable with
parameter p=F(x), hence nFn(x) is a binomial random variable with
mean nF(x) and variance nF(x)(1-F(x)).
53
End of Chapter
4
54