Continuous Random Variables: Density, PMF F (X) Density, PDF F (Y) Y/6, 2 y 4

Download as pdf or txt
Download as pdf or txt
You are on page 1of 28

Chapter 3

Continuous Random Variables

3.1 Introduction
Rather than summing probabilities related to discrete random variables, here for
continuous random variables, the density curve is integrated to determine probability.

Exercise 3.1 (Introduction)


Patient’s number of visits, X, and duration of visit, Y .

density, pmf f(x) density, pdf f(y) = y/6, 2 < y <


_4
probability less than 1.5 =
1 sum of probability 1
at specific values
P(X < 1.5) = P(X = 0) + P(X = 1) probability less than 3
0.75 2/3 = area under curve,
= 0.25 + 0.50 = 0.75
P(Y < 3) = 5/12
0.50 1/2 probability at 3,
P(X = 2) = 0.25 P(Y = 3) = 0
0.25 1/3

0 1 2 0 1 2 3 4 x

probability (distribution): cdf F(x)


1 1
probability =
0.75 value of function,
0.75 F(3) = P(Y < 3) = 5/12
0.50 probability less than 1.5 =
value of function
0.25 F(1.5 ) = P(X < 1.5) = 0.75
0.25

0 1 2 0 1 2 3 4 x

Figure 3.1: Comparing discrete and continuous distributions

73
74 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

1. Number of visits, X is a (i) discrete (ii) continuous random variable,


and duration of visit, Y is a (i) discrete (ii) continuous random variable.

2. Discrete

(a) P (X = 2) = (i) 0 (ii) 0.25 (iii) 0.50 (iv) 0.75


(b) P (X ≤ 1.5) = P (X ≤ 1) = F (1) = 0.25 + 0.50 = 0.75
requires (i) summation (ii) integration and is a value of a
(i) probability mass function (ii) cumulative distribution function
which is a (i) stepwise (ii) smooth increasing function
P R
(c) E(X) = (i) xf (x) (ii) xf (x) dx
(d) V ar(X) = (i) E(X 2 ) − µ2 (ii) E(Y 2 ) − µ2
 
(e) M (t) = (i) E etX (ii) E etY
(f) Examples of discrete densities (distributions) include (choose one or more)
(i) uniform
(ii) geometric
(iii) hypergeometric
(iv) binomial (Bernoulli)
(v) Poisson

3. Continuous

(a) P (Y = 3) = (i) 0 (ii) 0.25 (iii) 0.50 (iv) 0.75


R3 2
ix=3
32 22
(b) P (Y ≤ 3) = F (3) = 2 x6 dx = x12 = 12 − 12 5
= 12
x=2
requires (i) summation (ii) integration and is a value of a
(i) probability density function (ii) cumulative distribution func-
tion
which is a (i) stepwise (ii) smooth increasing function
P R
(c) E(Y ) = (i) yf (y) (ii) yf (y) dy
(d) V ar(Y ) = (i) E(X 2 ) − µ2 (ii) E(Y 2 ) − µ2
 
(e) M (t) = (i) E etX (ii) E etY
(f) Examples of continuous densities (distributions) include (choose one or
more)
(i) uniform
(ii) exponential
(iii) normal (Gaussian)
(iv) Gamma
(v) chi-square
(vi) student-t
(vii) F
Section 2. Definitions (LECTURE NOTES 5) 75

3.2 Definitions
Random variable X is continuous if probability density function (pdf) f is continuous
at all but a finite number of points and possesses the following properties:
• f (x) ≥ 0, for all x,
R∞
• −∞ f (x) dx = 1,
Rb
• P (a < X ≤ b) = a f (x) dx
The (cumulative) distribution function (cdf) for random variable X is
Z x
F (x) = P (X ≤ x) = f (t) dt,
−∞

and has properties


• limx→−∞ F (x) = 0,
• limx→∞ F (x) = 1,
• if x1 < x2 , then F (x1 ) ≤ F (x2 ); that is, F is nondecreasing,
Rb
• P (a ≤ X ≤ b) = P (X ≤ b) − P (X ≤ a) = F (b) − F (a) = a f (x) dx,
Rx
• F 0 (x) = dx
d
−∞
f (t) dt = f (x).

f(x) is positive
density, f(x)
probability probability
_ a)
cdf F(a) = P(X < total area = 1
_ b) = F(b) - F(a)
P(a < X <

a x a b

Figure 3.2: Continuous distribution

The expected value or mean of random variable X is given by


Z ∞
µ = E(X) = xf (x) dx,
−∞

the variance is

σ 2 = V ar(X) = E[(X − µ)2 ] = E(X 2 ) − [E(X)]2 = E(X 2 ) − µ2


76 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)


with associated standard deviation, σ = σ2.
The moment-generating function is
Z ∞
M (t) = E etX = etX f (x) dx
 
−∞

for values of t for which this integral exists.


Expected value, assuming it exists, of a function u of X is
Z ∞
E[u(X)] = u(x)f (x) dx
−∞

The (100p)th percentile is a value of X denoted πp where


Z πp
p= f (x) dx = F (πp )
−∞

and where πp is also called the quantile of order p. The 25th, 50th, 75th percentiles
are also called first, second, third quartiles, denoted p1 = π0.25 , p2 = π0.50 , p3 = π0.75
where also 50th percentile is called the median and denoted m = p2 . The mode is the
value x where f is maximum.

Exercise 4.2 (Definitions)

1. Waiting time.
Let the time waiting in line, in minutes, be described by the random variable
X which has the following pdf,
 1
6
x, 2 < x ≤ 4,
f (x) =
0, otherwise.

density, pdf f(x) probability, cdf F(x)


1 1
probability =
probability less than 3 value of function,
2/3 = area under curve, F(3) = P(X < 3) = 5/12 0.75
P(X < 3) = 5/12
1/2 probability at 3,
P(X = 3) = 0
1/3 0.25

0 1 2 3 4 x 0 1 2 3 4 x

Figure 3.3: f(x) and F(x)


Section 2. Definitions (LECTURE NOTES 5) 77

(a) Verify function f (x) satisfies the second property of pdfs,


Z ∞ Z 4 x=4
1 x2 42 22 12
f (x) dx = x dx = = − = =
−∞ 2 6 12 x=2 12 12 12
(i) 0 (ii) 0.15 (iii) 0.5 (iv) 1
(b) P (2 < X ≤ 3) =
3 x=3
x2 32 22
Z
1
x dx = = − =
2 6 12 x=2 12 12
5 9
(i) 0 (ii) 12
(iii) 12
(iv) 1
32 32
(c) P (X = 3) = P (3− < X ≤ 3) = 12
− 12
= 0 6= f (3) = 1
6
· 3 = 0.5
(i) True (ii) False
1
So the pdf f (x) = 6
x determined at some value of x does not determine probability.

(d) P (2 < X ≤ 3) = P (2 < X < 3) = P (2 ≤ X ≤ 3) = P (2 ≤ X < 3)


(i) True (ii) False
because P (X = 3) = 0 and P (X = 2) = 0

(e) P (0 < X ≤ 3) =
3 x=3
x2 32 22
Z
1
x dx = = − =
2 6 12 x=2 12 12
5 9
(i) 0 (ii) 12
(iii) 12 (iv) 1
Why integrate from 2 to 3 and not 0 to 3?

(f) P (X ≤ 3) =
3 x=3
x2 32 22
Z
1
x dx = = − =
2 6 12 x=2 12 12
5 9
(i) 0 (ii) 12
(iii) 12
(iv) 1
(g) Determine cdf (not pdf) F (3).
3 x=3
x2 32 22
Z
1
F (3) = P (X ≤ 3) = x dx = = − =
2 6 12 x=2 12 12
5 9
(i) 0 (ii) 12
(iii) 12
(iv) 1
(h) Determine F (3) − F (2).
3 x=3
x2 32 22
Z
1
F (3)−F (2) = P (X ≤ 3)−P (X ≤ 2) = P (2 ≤ X ≤ 3) = x dx = = − =
2 6 12 x=2 12 12
5 9
(i) 0 (ii) 12
(iii) 12 (iv) 1
because everything left of (below) 3 subtract everything left of 2 equals what is between 2 and 3
78 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

(i) The general distribution function (cdf) is


 Rx
 R−∞ 0 dt = 0,it=x
 x ≤ 2,
x t t2 x2 4
F (x) = 2 6
dt = 12 = 12
− 12
, 2 < x ≤ 4,

 t=2
1, x > 4.
2 2
In other words, F (x) = x12 − 12
4
= x 12−4 on (2, 4]. Both pdf density and cdf
distribution are given in the figure above.
(i) True (ii) False
22 4 5 9
(j) F (2) = 12
− 12
= (i) 0 (ii) 12
(iii) 12
(iv) 1.
32 4 5 9
(k) F (3) = 12
− 12
= (i) 0 (ii) 12
(iii) 12
(iv) 1.
   
3.52 4 2.52 4
(l) P (2.5 < X < 3.5) = F (3.5) − F (2.5) = 12
− 12
− 12
− 12
=
(i) 0 (ii) 0.25 (iii) 0.50 (iv) 1.
 
3.52 4
(m) P (X > 3.5) = 1 − P (X ≤ 3.5) = 1 − F (3.5) = 1 − 12
− 12
=
(i) 0.3125 (ii) 0.4425 (iii) 0.7650 (iv) 1.
(n) Expected value. The average wait time is
∞ 4 4 4
x2 x3 43 23
Z Z   Z
1
µ = E(X) = xf (x) dx = x x dx = dx = = − =
−∞ 2 6 2 6 18 2 18 18
23
(i) 9
(ii) 28
9
(iii) 31
9
(iv) 35
9
.
library(MASS) # INSTALL once, RUN once (per session) library MASS
mean <- 4^3/18 - 2^3/18; mean; fractions(mean) # turn decimal to fraction

[1] 3.111111
[1] 28/9

(o) Expected value of function u = x2 .


∞ 4 4 4
x3 x4 44 24
Z Z   Z
2
 2 2 1
E X = x f (x) dx = x x dx = dx = = − =
−∞ 2 6 2 6 24 2 24 24

(i) 9 (ii) 10 (iii) 11 (iv) 12.


(p) Variance, method 1. Variance in wait time is
 2
2 2
 2 28
σ = V ar(X) = E X − µ = 10 − =
9
23 26
(i) 81
(ii) 81 (iii) 31
81
35
(iv) 81 .
library(MASS) # INSTALL once, RUN once (per session) library MASS
sigma2 <- 10 - (28/9)^2 ; fractions(sigma2) # turn decimal to fraction
Section 2. Definitions (LECTURE NOTES 5) 79

[1] 26/81

(q) Variance, method 2. Variance in wait time is


Z ∞
2 2
σ = V ar(X) = E[(X − µ) ] = (x − µ)2 f (x) dx
−∞
Z 2  
4
28 1
= x− x dx
2 9 6
Z 4 3
56x2 784x

x
= − + dx
2 6 54 486
4
x4 56x3 784x2
= − + =
24 162 972 2
23 26
(i) 81
(ii) 81 (iii) 31
81
35
(iv) 81 .
library(MASS) # INSTALL once, RUN once (per session) library MASS
sigma2 <- (4^4/24 - 56*4^3/162 + 784*4^2/972) - (2^4/24 - 56*2^3/162 + 784*2^2/972)
fractions(sigma2) # turn decimal to fraction

[1] 26/81

(r) Standard deviation. Standard deviation in time spent on the call is,


r
26
σ = σ2 = ≈
81

(i) 0.57 (ii) 0.61 (iii) 0.67 (iv) 0.73.


(s) Moment generating function.
Z ∞
tX
M (t) = E[e ] = etx f (x) dx
−∞
Z 4  
tx 1
= e x dx
2 6
1 4 tx
Z
= xe dx
6 2
 4
1 tx x 1
= e − 2
6 t t
  x=2  
1 4t 4 1 1 2t 2 1
= e − − e − , t 6= 0.
6 t t2 6 t t2

(i) True (ii) False


v du where u = x, dv = etx dx.
R R
use integration by parts for t 6= 0 case: u dv = uv −
80 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

x2 −4
(t) Median. Since the distribution function F (x) = 12
, then median m =
π0.50 occurs when

m2 − 4 1
F (m) = P (X ≤ m) = =
12 2
so r
12
m = π0.50 = +4≈
2
(i) 1.16 (ii) 2.16 (iii) 3.16
negative answer m ≈ −3.16 is not in 2 < x ≤ 4, so m 6= −3.16

2. Density with unknown a. Find a such that


 2
ax , 1 < x ≤ 5,
f (x) =
0, otherwise.

density, pdf f(x) probability, cdf F(x) = P(X < x)


1 1
0.75 0.75
0.50 0.50

0.25 0.25

0 1 2 3 4 5 x 0 1 2 3 4 5 x

Figure 3.4: f(x) and F(x)

(a) Find constant a. Since the cdf F (5) = 1, and


x t=x
t3 ax3 a a(x3 − 1)
Z
2
F (x) = at dt = a = − =
1 3 t=1 3 3 3

then
a(53 − 1) 124a
F (5) = = = 1,
3 3
3 2 1
so a = (i) 124
(ii) 124
(iii) 124
(b) Determine pdf.
f (x) = ax2 =
3 2 1
(i) 124
x2 (ii) 124
x2 (iii) 124
x2
Section 2. Definitions (LECTURE NOTES 5) 81

(c) Determine cdf.


a(x3 − 1) 3 x3 − 1
F (x) = = × =
3 124 3
3 2 1
(i) 124
(x3 − 1) (ii) 124
(x3 − 1) (iii) 124
(x3 − 1)
(d) Determine P (X ≥ 2).
1 3
P (X ≥ 2) = 1 − P (X < 2) = 1 − F (2) = 1 − (2 − 1) =
124
61 117 61
(i) 124
(ii) 124
(iii) 117
(e) Determine P (X ≥ 4).
1 3
P (X ≥ 4) = 1 − P (X < 4) = 1 − F (4) = 1 − (4 − 1) =
124
61 117 61
(i) 124
(ii) 124
(iii) 117
(f) Determine P (X ≥ 4|X ≥ 2).
P (X ≥ 4 ∩ X ≥ 2) P (X ≥ 4) 61/124
P (X ≥ 4|X ≥ 2) = = = =
P (X ≥ 2) P (X ≥ 2) 117/124
116 117 61
(i) 124
(ii) 124
(iii) 117
(g) Expected value.
Z ∞ 5 5 5
3x3 3x4 3 · 54 3 · 14
Z   Z
3 2
µ = E(X) = xf (x) dx = x x dx = dx = = − =
−∞ 1 124 1 124 496 1 496 496
114 115 116 117
(i) 31
(ii) 31
(iii) 31
(iv) 31
.
2
(h) Expected value of function u = x .
Z ∞ Z 5  Z 5 4 5
3x5 3 · 55 3 · 15

2 2 2 3 2 3x
E(X ) = x f (x) dx = x x dx = dx = = − ≈
−∞ 1 124 1 124 620 1 620 620
(i) 15.12 (ii) 16.12 (iii) 17.12 (iv) 18.12.
(i) Determine 25th percentile (first quartile), p1 = π0.25 .

1
Since F (x) = 124
(x3 − 1), then,
1 3 1
F (π0.25 ) = (π0.25 − 1) = 0.25 = ,
124 4
and so r
3 124
π0.25 = +1≈
4
(i) 2.38 (ii) 3.17 (iii) 4.17
The π0.25 is that value where 25% of probability is at or below (to the left of) this value.
82 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

(j) Determine π0.01 .


1 3
F (π0.01 ) = (π0.01 − 1) = 0.01,
124
and so √
3
π0.01 = 0.01 · 124 + 1 ≈
(i) 1.31 (ii) 3.17 (iii) 4.17
The π0.01 is that value where 1% of probability is at or below (to the left of) this value.

3. Piecewise pdf. Let random variable X have pdf



 x, 0 < x ≤ 1,
f (x) = 2 − x, 1 < x ≤ 2,
0, elsewhere.

density, pdf f(x) probability, cdf F(x) = P(X < x)


1 1

0.75 0.75

0.50 0.50

0.25 0.25

0 1 2 x 1 x
0 2 3

Figure 3.5: f(x) and F(x)

(a) Expected value.


Z ∞
µ = E(X) = xf (x) dx
−∞
Z 1 Z 2
= x (x) dx + x (2 − x) dx
0 1
Z 1 Z 2
2
2x − x2 dx

= x dx +
0 1
 3 1  2 2
x 2x x3
= + −
3 0 2 3 1
 3 3
23 13
    
1 0 2 2
= − + 2 − − 1 − =
3 3 3 3
(i) 1 (ii) 2 (iii) 3 (iv) 4.
Section 2. Definitions (LECTURE NOTES 5) 83

(b) E (X 2 ).
Z 1 Z 2
2 2
x2 (2 − x) dx

E X = x (x) dx +
Z0 1 Z 2 1
x3 dx + 2x2 − x3 dx

=
0 1
 4 1  3 2
x 2x x4
= + −
4 0 3 4 1
 4 4
2(2)3 24 2(1)3 14
    
1 0
= − + − − − =
4 4 3 4 3 4
4 5 6
(i) 6
(ii) 6
(iii) 6
(iv) 76 .

(c) Variance.
7
σ 2 = Var(X) = E X 2 − µ2 = − 12 =

6
1 1 1 1
(i) 3
(ii) 4
(iii) 5 (iv) 6 .

(d) Standard deviation.



r
1
σ= σ2 = ≈
6
(i) 0.27 (ii) 0.31 (iii) 0.41 (iv) 0.53.

(e) Median. Since the distribution function is




 0 ix x≤0
 x t dt = t2 2
= x2 ,
 R
0 < x ≤ 1,

0 2
F (x) = t=0 ix
Rx t2 1 x2 12

 1
(2 − t) dt + F (1) = 2t − 2 + 2
= 2x − 2
− (2(1) − 2
) + 12 , 1 < x ≤ 2,
t=1


1, x > 2.

1 1
F (1) = 2
, so add 2
to F (x) for 1 < x ≤ 2
then median m occurs when
 m2
 2 = 12 , 0 < x ≤ 1,
2
F (x) = 2m − m2 − 1 = 12 , 1 < x ≤ 2,

1, x > 2.

so for both 0 < x ≤ 1 and 1 < x ≤ 2, m = (i) 1 (ii) 1.5 (iii) 2


84 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

3.3 The Uniform and Exponential Distributions


Two special probability density functions are discussed: uniform and exponential.
The continuous uniform (rectangular) distribution of random variable X has density
 1
b−a
, a ≤ x ≤ b,
f (x) =
0, elsewhere,

where “X is U (a, b)” means “X is uniform over [a,b]”, distribution function,



 0 x < a,
x−a
F (x) = a ≤ x ≤ b,
 b−a
1 x > b,

and its expected value (mean), variance and standard deviation are,

a+b 2 (b − a)2 p
µ = E(X) = , σ = V ar(X) = , σ= V ar(X),
2 12
and moment-generating function is
  etb − eta
M (t) = E etX = , t 6= 0.
t(b − a)
The continuous exponential random variable X has density
 −λx
λe , x ≥ 0,
f (x) =
0, elsewhere,
distribution function, 
0 x < 0,
F (x) =
1 − e−λx x ≥ 0,
and its expected value (mean), variance and standard deviation are,
1 1 1
µ = E(X) = , σ 2 = V (Y ) = , σ= ,
λ λ2 λ
and moment-generating function is
λ
M (t) = E etX =
 
, t < λ.
λ−t
1
Notice if θ = λ
in the exponential distribution,
1 x x θ
f (x) = e− θ , F (x) = 1 − e− θ , µ = σ = θ, M (t) = .
θ 1 − tθ

Exercise 3.3 (The Uniform and Exponential Distributions)


Section 3. The Uniform and Exponential Distributions (LECTURE NOTES 5) 85

1. Uniform: potato weights. An automated process fills one bag after another
with Idaho potatoes. Although each filled bag should weigh 50 pounds, in fact,
because of the differing shapes and weights of each potato, each bag weight, X,
is anywhere from 49 pounds to 51 pounds, with uniform density:

0.5, 49 ≤ x ≤ 51,
f (x) =
0, elsewhere.

(a) Since a = 49 and b = 51, the distribution is



 0 x < 49,
x−49
F (x) = 49 ≤ x ≤ 51,
 51−49
1 x > 52,

and so graphs of density and distribution are given in the figure.


(i) True (ii) False

density, pdf f(x) probability, cdf F(x) = P(X < x)


1 1

0.75 0.75

0.50 0.50

0.25 0.25

x 0 49 50 x
0 49 50 51 51

Figure 3.6: Distribution function: continuous uniform

(b) Determine P (49.5 < X < 51) by integrating pdf.


Z 51
1 x i51 51 49.5 1.5
P (49.5 < X < 51) = dx = = − =
49.5 2 2 x=49.5 2 2 2

(i) 0.25 (ii) 0.50 (iii) 0.75 (iv) 1.


(c) Determine P (49.5 < X < 51) using cdf.
51 − 49 49.5 − 49
P (49.5 < X < 51) = F (51) − F (49.5) = − =
51 − 49 51 − 49
(i) 0.25 (ii) 0.50 (iii) 0.75 (iv) 1.
1 - punif(49.5,49,51) # uniform P(49.5 < X < 51) = 1 - P(X < 49.5), 49 < x < 51

[1] 0.75
86 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

(d) The chance the bags weight more than 49.5


49.5 − 49
P (X > 49.5) = 1 − P (X ≤ 49.5) = 1 − F (49.5) = 1 − =
51 − 49
(i) 0.25 (ii) 0.50 (iii) 0.75 (iv) 1.
1 - punif(49.5,49,51) # uniform P(X > 49.5) = 1 - P(X < 49.5), 49 < x < 51

[1] 0.75

(e) What is the mean weight of a bag of potatoes?


a+b 49 + 51
µ = E(X) = = =
2 2
(i) 49 (ii) 50 (iii) 51 (iv) 52.
(f) What is the standard deviation in the weight of a bag of potatoes?
s s
2
(b − a) (51 − 49)2
σ= = ≈
12 12

(i) 0.44 (ii) 0.51 (iii) 0.55 (iv) 0.58.


(g) Determine probability within 1 standard deviation of mean.

P (µ − σ < X < µ + σ) ≈ P (50 − 0.58 < X < 50 + 0.58)


= P (49.42 < X < 50.58)
Z 50.58
1 x i50.58
= dx = =
49.42 2 2 x=49.42

(i) 0.44 (ii) 0.51 (iii) 0.55 (iv) 0.58.


(h) Function of X. If it costs $0.0556 (5.56 cents) per pound of potatoes,
then the cost of X pounds of potatoes is Y = 0.0556X. Determine the
probability a bag of potatoes chosen at random costs at least $2.78.
Z 51
1 x i51
P (Y ≥ 2.78) = P (0.0556X ≥ 2.78) = P (X ≥ 50) = dx = =
50 2 2 x=50

(i) 0.25 (ii) 0.50 (iii) 0.75 (iv) 1.

2. Exponential: atomic decay. Assume atomic time-to-decay obeys exponential


pdf with (inverse) mean rate of decay λ = 3,
 −3x
3e , x > 0,
f (x) =
0, otherwise.
Section 3. The Uniform and Exponential Distributions (LECTURE NOTES 5) 87

density, pdf f(x) probability, cdf F(x) = P(X < x)


4 4
3 3
2 2
1 1
x x
0 1 2 3 4 5 0 1 2 3 4 5

Figure 3.7: f(x) and F(x)

(a) Find F (x). By substituting u = −3x, then du = −3 dx and so


Z Z Z
−3x −3x
F (x) = 3e dx = − e (−3) dx = − eu du = −eu =

(i) −e−3x (ii) e−3x (iii) −3e−3x


then t=x
F (x) = −e−3t t=0 = −e−3x − (−e−3(0) ) =
 

(i) e−3x − 1 (ii) e−3x + 1 (iii) 1 − e−3x


(b) Chance atomic decay takes at least 2 microseconds, P (X ≥ 2), is

P (X ≥ 2) = 1 − P (X < 2) = 1 − F (2) = 1 − 1 − e−3(2) ≈




(i) 0.002 (ii) 0.003 (iii) 0.004


1 - pexp(2,3) # exponential P(X > 2), lambda = 3

[1] 0.002478752

(c) Chance atomic decay is between 1.13 and 1.62 microseconds is

P (1.13 < X < 1.62) = F (1.62) − F (1.13)


= 1 − e−(3)(1.62) − 1 − e−(3)(1.13)
 

= e−(3)(1.13) − e−(3)(1.62) ≈
(i) 0.014 (ii) 0.026 (iii) 0.034 (iv) 0.054.
pexp(1.62,3) - pexp(1.13,3) # exponential P(1.13 < X < 1.62), lambda = 3

[1] 0.02595819

(d) Calculate probability atomic rate between mean and 1 standard deviation
88 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

below mean. Since λ = 3, so µ = σ = 31 , so


 
1 1 1
P (µ − σ < X < µ) = P − <X<
3 3 3
 
1
= P 0<X<
3
 
1
= F − F (0)
3
 1

= 1 − e−(3) 3 − 1 − e−(3)(0) ≈


(i) 0.33 (ii) 0.43 (iii) 0.53 (iv) 0.63.


(e) Check if f (x) is a pdf.
Z ∞ Z b
−3x
x=b
3e−3x dx = lim −e−3x x=0 = lim −e−3b − (−ex(0) ) = 1
  
3e dx = lim
0 b→∞ 0 b→∞ b→∞

or, equivalently,
lim F (b) = lim 1 − e−3b = 1
b→∞ b→∞

(i) True (ii) False


0
(f) Show F (x) = f (x).
d d
F 0 (x) = 1 − e−3x =

F (x) =
dx dx
(i) 3e−3x (ii) e−3x + 1 (iii) 1 − e−3x
λ
(g) Determine µ using mgf. Since M (t) = λ−t ,
 
0 d λ λ
M (t) = = ,
dt λ − t (λ − t)2
1
so µ = M 0 (0) = (i) λ (ii) 1 (iii) 1 − λ (iv) λ
.
0 0 (λ−t)(0)−λ(−1)
u λ vu −uv λ 1
use quotient rule, since f = v
= λ−t
, f0 = v2
= (λ−t)2
= (λ−t)2
= λ
if t = 0

(h) Memoryless property of exponential. Chance atomic decay lasts at least 10


microseconds is
P (X > 10) = 1 − F (10) = 1 − (1 − e−3(10) ) =
(i) e−10 (ii) e−20 (iii) e−30 (iv) e−40 .

and chance atomic decay lasts at least 15 microseconds, given it


has already lasted at least 5 microseconds is
P (X > 15, Y > 5) P (X > 15) 1 − (1 − e−3(15) )
P (X > 15|X > 5) = = = =
P (X > 5) P (X > 5) 1 − (1 − e−3(5) )
Section 3. The Uniform and Exponential Distributions (LECTURE NOTES 5) 89

(i) e−10 (ii) e−20 (iii) e−30 (iv) e−40 .

so, in other words,

P (X > 15) P (X > 10 + 5)


P (X > 15|X > 5) = = = P (X > 10)
P (X > 5) P (X > 5)
or
P (X > 10 + 5) = P (X > 10) · P (X > 5)
This is an example of the “memoryless” property of the exponential, it implies time intervals are
independent of one another. Chance of decay after 15 microsecond, given decay after 5 microseconds,
same as chance of decay after 10 seconds; it is as though first 5 microseconds “forgotten”.

3. Verifying if experimental decay is exponential. In an atomic decay study, let


random variable X be time-to-decay where x ≥ 0. The relative frequency
distribution table and histogram for a sample of 20 atomic decays (measured
in microseconds) from this study are given below. Does this data follow an
exponential distribution?

0.60, 0.07, 0.66, 0.17, 0.06, 0.14, 0.15, 0.19, 0.07, 0.36
0.85, 0.44, 0.71, 1.02, 0.07, 0.21, 0.16, 0.16, 0.01, 0.88

0.40
0.35
0.30 exponential density curve, pdf, f(x)
relative frequency

0.25
0.20
relative frequency
0.15 histogram approximates
0.10 density curve
0.05
0
0 0.2 0.4 0.6 0.8 1.0 1.2
time-to-decay (microseconds), x

Figure 3.8: Histogram of exponential decay times


90 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

bin frequency relative exponential


frequency probability
11
(0,0.2] 11 20
= 0.55 0.44
(0.2,0.4] 2 0.10 0.25
(0.4,0.6] 2 0.10 0.14
(0.6,0.8] 2 0.10 0.08
(0.8,1] 2 0.10 0.04
(1,1.2] 1 0.05 0.02

(a) Sample mean time-to-decay.


0.60 + 0.07 + · · · + 0.88
x̄ = =
20
(i) 0.330 (ii) 0.349 (iii) 0.532 (iv) 0.631.
x <- c(0.60,0.07,0.66,0.17,0.06,0.14,0.15,0.19,0.07,0.36,0.85,0.44,0.71,1.02,0.07,0.21,0.16,0.16,0.01,0.88)
mean(x)

[1] 0.349

(b) Approximate λ for exponential. Since approximate mean time-to-decay is

µ ≈ x̄ = 0.349,

then, if time-to-decay is exponential, approximate value of λ parameter,


since µ = λ1 , is
1 1
λ= ≈ ≈
µ x̄
(i) 2.9 (ii) 3.0 (iii) 3.1 (iv) 3.2

and so possible model of data might be exponential density


(i) 2.9e−3x (ii) 3e−2.9x (iii) 2.9e−2.9x (iv) 3e−3x .
(c) Approximate probabilities.
Z 0.2
P (0 ≤ X ≤ 0.2) = 2.9e−2.9x dx ≈
0

(i) 0.44 (ii) 0.25 (iii) 0.14 (iv) 0.08.


remaining probabilities in table calculated in a similar way
pexp(0.2,2.9) - pexp(0,2.9) # exponential P(0 < X < 0.2), lambda approx 2.9

[1] 0.4401016

(d) Since sample relative frequencies do not closely match exponential with
λ = 2.9 in table, this distribution (i) is (ii) is not a good fit to the data.
Section 4. The Normal Distribution (LECTURE NOTES 5) 91

3.4 The Normal Distribution


The continuous normal distribution of random variable X, defined on the interval
(−∞, ∞), has pdf with parameters µ and σ, that is, “X is N (µ, σ 2 )”,
1 2
f (x) = √ e−(1/2)[(x−µ)/σ] ,
σ 2π
and cdf Z x
1 2
F (x) = P (X ≤ x) = √ e−(1/2)[(t−µ)/σ] dt,
−∞ σ 2π

and its expected value (mean), variance and standard deviation are,
p
E(X) = µ, V ar(X) = σ 2 , σ = V ar(X),
and mgf is
σ 2 t2
 
M (t) = exp µt + .
2
A normal random variable, X, may be transformed to a standard normal, Z,
1 2
f (z) = √ e−z /2 ,

where “Z is N (0, 1)” and
Z z
1 2
Φ(z) = P (Z ≤ z) = √ e−t /2 dt,
−∞ 2π
t2 /2
where µ = 0, σ = 1 and M (t) = e using the following equation,
X −µ
Z= .
σ
The distribution of this density does not have a closed–form expression and so must be solved using numerical
integration methods. We will use both R and the tables to obtain approximate numerical answers.

Exercise 3.4 (The Normal Distribution)


1. Nonstandard normal: IQ scores.
It has been found that IQ scores, Y , can be distributed by a normal distribution.
Densities of IQ scores for 16 year olds, X1 , and 20 year olds, X2 , are given by
1 2
f (x1 ) = √ e−(1/2)[(y−100)/16] ,
16 2π
1 2
f (x2 ) = √ e−(1/2)[(y−120)/20] .
20 2π
A graph of these two densities is given in the figure.
92 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

f(x)

16 year old IQs

σ = 16
f(x)
20 year old IQs

σ = 20

x
µ = 100 µ = 120

Figure 3.9: Normal distributions: IQ scores

(a) Mean IQ score for 20 year olds is


µ = (i) 100 (ii) 120 (iii) 124 (iv) 136.
(b) Average (or mean) IQ scores for 16 year olds is
µ = (i) 100 (ii) 120 (iii) 124 (iv) 136.
(c) Standard deviation in IQ scores for 20 year olds
σ = (i) 16 (ii) 20 (iii) 24 (iv) 36.
(d) Standard deviation in IQ scores for 16 year olds is
σ = (i) 16 (ii) 20 (iii) 24 (iv) 36.
(e) Normal density for 20 year old IQ scores is
(i) broader than normal density for 16 year old IQ scores.
(ii) as wide as normal density for 16 year old IQ scores.
(iii) narrower than normal density for 16 year old IQ scores.
(f) Normal density for the 20 year old IQ scores is
(i) shorter than normal density for 16 year old IQ scores.
(ii) as tall as normal density for 16 year old IQ scores.
(iii) taller than normal density for 16 year old IQ scores.
(g) Total area (probability) under normal density for 20’s IQ scores is
(i) smaller than area under normal density for 16’s IQ scores.
(ii) the same as area under normal density for 16’s IQ scores.
(iii) larger than area under normal density for 16’s IQ scores.
(h) Number of different normal densities:
(i) one (ii) two (iii) three (iv) infinity.
2. Percentages: IQ scores.
Densities of IQ scores for 16 year olds, X1 , and 20 year olds, X2 , are given by
1 2
f (x1 ) = √ e−(1/2)[(y−100)/16] ,
16 2π
1 2
f (x2 ) = √ e−(1/2)[(y−120)/20] .
20 2π
Section 4. The Normal Distribution (LECTURE NOTES 5) 93

f(x) N(100,162 ) P(96 < X < 120) = ? f(x)


P(X > 84) = ? 2
N(100,16 )

84 x 100 120 x
100 96
(a) (b)

P(X > 84) = ? f(x) 2


P(96 < X < 120) = ? f(x)
2
N(120,20 ) N(120,20 )

x x
84 120 96 120

(c) (d)

Figure 3.10: Normal probabilities: IQ scores

(a) For the sixteen year old normal distribution, where µ = 100 and σ = 16,
Z ∞
1 2
P (X1 > 84) = √ e−(1/2)[(y−100)/16] dx1 ≈
84 16 2π

(i) 0.4931 (ii) 0.9641 (iii) 0.8413 (iv) 0.3849.


1 - pnorm(84,100,16) # normal P(X > 84), mean = 100, SD = 16

[1] 0.8413447

(b) For sixteen year old normal distribution, where µ = 100 and σ = 16,
P (96 < X1 < 120) ≈ (i) 0.4931 (ii) 0.9641 (iii) 0.8413 (iv) 0.3849.
pnorm(120,100,16) - pnorm(96,100,16) # normal P(96 < X < 120), mean = 100, SD = 16

[1] 0.4930566

(c) For twenty year old, where µ = 120 and σ = 20,


P (X2 > 84) ≈ (i) 0.4931 (ii) 0.9641 (iii) 0.8413 (iv) 0.3849.
1 - pnorm(84,120,20) # normal P(X > 84), mean = 120, SD = 20

[1] 0.9640697

(d) For twenty year old, where µ = 120 and σ = 20,


P (96 < X2 < 120) ≈ (i) 0.4931 (ii) 0.9641 (iii) 0.8413 (iv) 0.3849.
pnorm(120,120,20) - pnorm(96,120,20) # normal P(96 < X < 120), mean = 120, SD = 20

[1] 0.3849303
94 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

(e) Empirical Rule (68-95-99.7 rule). If X is N (120, 202 ), find probability


within 1, 2 and 3 SDs of mean.

P (µ−σ < X < µ+σ) = P (120 − 20 < X < 120 + 20) = P (100 < X < 140) ≈

(i) 0.683 (ii) 0.954 (iii) 0.997 (iv) 1.

P (µ − 2σ < X < µ + 2σ) = P (80 < X < 160) ≈

(i) 0.683 (ii) 0.954 (iii) 0.997 (iv) 1.

P (µ − 3σ < X < µ + 3σ) = P (60 < X < 180) ≈


(i) 0.683 (ii) 0.954 (iii) 0.997 (iv) 1.
Empirical rule is true for all X which are N (µ, σ 2 ).
pnorm(140,120,20) - pnorm(100,120,20); pnorm(160,120,20) - pnorm(80,120,20); pnorm(180,120,20) - pnorm(60,120,20)

[1] 0.6826895 [1] 0.9544997 [1] 0.9973002

3. Standard normal.
Normal densities of IQ scores for 16 year olds, X1 , and 20 year olds, X2 ,
1 2
f (x1 ) = √ e−(1/2)[(y−100)/16] ,
16 2π
1 2
f (x2 ) = √ e−(1/2)[(y−120)/20] .
20 2π
Both densities may be transformed to a standard normal with µ = 0 and σ = 1
using the following equation,
Y −µ
Z= .
σ

(a) Since µ = 100 and σ = 16, a 16 year old who has an IQ of 132 is
z = 132−100
16
= (i) 0 (ii) 1 (iii) 2 (iv) 3
standard deviations above the mean IQ, µ = 100.
(b) A 16 year old who has an IQ of 84 is
z = 84−100
16
= (i) −2 (ii) −1 (iii) 0 (iv) 1
standard deviations below the mean IQ, µ = 100.
(c) Since µ = 120 and σ = 20, a 20 year old who has an IQ of 180 is
z = 180−120
20
= (i) 0 (ii) 1 (iii) 2 (iv) 3
standard deviations above the mean IQ, µ = 120.
(d) A 20 year old who has an IQ of 100 is
z = 100−120
20
= (i) −3 (ii) −2 (iii) −1 (iv) 0
standard deviations below the mean IQ, µ = 120.
Section 4. The Normal Distribution (LECTURE NOTES 5) 95

16 year olds
mean 100, SD 16

52 68 84 100 116 132 148 X, nonstandard


-3 -2 -1 0 1 2 3 Z, standard

110

20 year olds
mean 120, SD 20

60 80 100 120 140 160 180 X, nonstandard


-3 -2 -1 0 1 2 3 Z, standard

Figure 3.11: Standard normal and (nonstandard) normal

(e) Although both the 20 year old and 16 year old scored the same, 110, on
an IQ test, the 16 year old is clearly brighter relative to his/her age group
than is the 20 year old relative his/her age group because
110 − 100 110 − 120
z1 = = 0.625 > z2 = = −0.5.
16 20
(i) True (ii) False
(f) The probability a 20 year old has an IQ greater than 90 is
 
90 − 120
P (X2 > 90) = P Z2 > = P (Z2 > −1.5) ≈
20
(i) 0.93 (ii) 0.95 (iii) 0.97 (iv) 0.99.
1 - pnorm(90,120,20) # normal P(X > 90), mean = 120, SD = 20
1 - pnorm(-1.5) # standard normal P(Z > -1.5), defaults to mean = 0, SD = 1

[1] 0.9331928

Or, using Table C.1 from the text,


P (Z2 > −1.5) = 1 − Φ(−1.5) ≈ 1 − 0.0668 ≈
(i) 0.93 (ii) 0.95 (iii) 0.97 (iv) 0.99.
Table C.1 can only used on standard normal questions, where Z is N (0, 1).

(g) The probability a 20 year old has an IQ between 125 and 135 is
 
125 − 120 135 − 120
P (125 < X2 < 135) = P < Z2 < = P (0.25 < Z2 < 0.75) =
20 20
(i) 0.13 (ii) 0.17 (iii) 0.27 (iv) 0.31.
96 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

pnorm(135,120,20) - pnorm(125,120,20) # normal P(125 < X < 135), mean = 120, SD = 20


pnorm(0.75) - pnorm(0.25) # standard normal P(0.25 < Z < 0.75)

[1] 0.1746663

Or, using Table C.1 from the text,

P (0.25 < Z2 < 0.75) = Φ(0.75) − Φ(0.25) ≈ 0.7734 − 0.5987 ≈

(i) 0.13 (ii) 0.17 (iii) 0.27 (iv) 0.31.


(h) If a normal random variable X with mean µ and standard deviation σ
can be transformed to a standard one Z with mean µ = 0 and standard
deviation σ = 1 using
X −µ
Z= ,
σ
then Z can be transformed to X using

X = µ + σZ.

(i) True (ii) False


(i) A 16 year old who has an IQ which is three (3) standards above the mean
IQ has an IQ of x1 = 100 + 3(16) =
(i) 116 (ii) 125 (iii) 132 (iv) 148.
(j) A 20 year old who has an IQ which is two (2) standards below the mean
IQ has an IQ of x2 = 120 − 2(20) =
(i) 60 (ii) 80 (iii) 100 (iv) 110.
(k) A 20 year old who has an IQ which is 1.5 standards below the mean IQ
has an IQ of x2 = 120 − 1.5(20) =
(i) 60 (ii) 80 (iii) 90 (iv) 95.
(l) The probability a 20 year old has an IQ greater than one (1) standard
deviation above the mean is

P (Z2 > 1) = P (X2 > 120 + 1(20)) = P (X2 > 140) =

(i) 0.11 (ii) 0.13 (iii) 0.16 (iv) 0.18.


1 - pnorm(1) # standard normal P(Z > 1), defaults to mean = 0, SD = 1
1 - pnorm(140,120,20) # normal P(X > 140), mean = 120, SD = 20

[1] 0.1586553

(m) Percentile. Determine 1st percentile, π0.01 , when X is N (120, 202 ).


Z π0.01
1 2
F (π0.01 ) = √ e−(1/2)[(y−120)/20] = 0.01,
−∞ 20 2π
and so π0.01 ≈ (i) 72.47 (ii) 73.47 (iii) 75.47
Section 5. Functions of Continuous Random Variables (LECTURE NOTES 5) 97

qnorm(0.01,120,20) # normal 1st percentile, mean = 120, SD = 20

[1] 73.47304

Or, using Table C.1 “backwards”, F (π0.01 ) = 0.01 when

Z≈

(i) −2.33 (ii) −2.34 (iii) −2.35

but since X is N (120, 202 ), 1st percentile

π0.01 = X = µ + σZ ≈ 120 + 20(−2.33) =

(i) 73.4 (ii) 74.1 (iii) 75.2


The π0.01 is that value where 1% of probability is at or below (to the left of) this value.

3.5 Functions of Continuous Random Variables


Assuming cdf F (x) is strictly increasing (as opposed to just monotonically increasing)
on a < x < b, one method to determine the pdf of a function, Y = U (X), of
random variable X, requires first determining the cdf of X, F (x), then using F (X)
to determine the cdf of Y , F (Y ), and finally differentiating the result,

dF
f (y) = F 0 (x) = .
du
The second fundamental theorem of calculus is sometimes used in this method,
Z u2 (y)
d du2 du1
fX (x) dx = fX (u2 (y)) − fX (u1 (y)) .
dy u1 (y) dy dy

Related to this is an algorithm for sampling at random values from a random variable
X with a desired distribution by using the uniform distribution:

• determine cdf of X, F (x), desired distribution

• find F −1 (y): set F (x) = y, solve for x = F −1 (y)

• generate values y1 , y2 , . . . , yn from Y , assumed to be U (0, 1),

• use x = F −1 (y) to simulate observed x1 , x2 , . . . , xn , from desired distribution.

Exercise 3.5 (Functions of Continuous Random Variables)


98 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

probability, cdf F(x) = P(X < x) = x - x 2/4 probability, cdf F(y) = P(Y < y) = y/2 - y 2/16

1 1

0.75 0.75
0.50 0.50

0.25 0.25

0 1 2 3 4 x 0 1 2 3 4 y = 2x

probability, cdf F(y) = P(Y < y) = (4 + 4y + y )/16


2
probability, cdf F(y) = P(Y < y) = y1/2- y/4

1 1
0.75 0.75
0.50 0.50
0.25 0.25
1/2
-2 -1 0 1 2 y = 2 - 2x 0 1 2 3 4 y=x

Figure 3.12: Distributions of various functions of Y = U (X)

1. Determine pdf of Y = U (X), given pdf of X. Consider density

1 − x2 , 0 ≤ x ≤ 2

fX (x) =
0 elsewhere

(a) Calculate pdf of Y = 2X, version 1.


x t=x
t2 x2
Z   
t
FX (x) = 1− dt = t − = x− ,
0 2 4 t=0 4
 y  y 
FY (y) = P (Y ≤ y) = P (2X ≤ y) = P X ≤ = FX
2 2
y 2

2
y y y
= − 2 = − ,
2 4 2 16
dF
fY (y) = FY0 (y) = =
dy
2y 2y y
(i) 3
2
− 32
(ii) 1
2
− 32
(iii) 1
2
− 32
(iv) 1
2
− y8 ,

where since y = 2x, 0 ≤ x ≤ 2 implies 0 ≤ y2 ≤ 2, or


(i) −1 ≤ y ≤ 2 (ii) −2 ≤ y ≤ 1 (iii) 0 ≤ y ≤ 3 (iv) 0 ≤ y ≤ 4
Section 5. Functions of Continuous Random Variables (LECTURE NOTES 5) 99

(b) Calculate pdf of Y = 2X, version 2.


 y
FY (y) = P (Y ≤ y) = P (2X ≤ y) = P X ≤
2
y  Z y
2
= FX = fX (x) dx
2 0
Z y  
d 2
y  d y  y/2 1
fY (y) = fX (x) dx = fX −0 = 1− · =
dy 0 2 dy 2 2 2

(i) 3
2
− 2y32
(ii) 12 − 2y32
(iii) 12 − 32y
(iv) 21 − y8 ,
Second fundamental theorem of calculus used here, notice cdf of FX (x) is not explicitly calculated.

(c) Calculate pdf for Y = 2 − 2X.


 
2−y
FY (y) = P (Y ≤ y) = P (2 − 2X ≤ y) = P X ≥
2
Z 2
= fX (x) dx,
2−y
2
2    
2−y 2−y
Z
d d
fY (y) = fX (x) dx = 0 − fX
dy 2−y
2
2 dy 2
!
2−y
2 −1
= − 1− · =
2 2

2y 2y y
(i) 3
2
− 32
(ii) 1
2
+ 32
(iii) 1
2
+ 32
(iv) 1
4
+ y8 ,

where since y = 2 − 2x, 0 ≤ x ≤ 2 implies 0 ≤ 2−y


2
≤ 2, or
(i) −1 ≤ y ≤ 2 (ii) −2 ≤ y ≤ 1 (iii) −1 ≤ y ≤ 1 (iv) −2 ≤ y ≤ 2

(d) Calculate pdf for Y = X 2 .



y

Z
2

FY (y) = P (Y ≤ y) = P X ≤ y = P (X ≤ y) = fX (x) dx,
0

Z y  √ 
d √ d √ y 1
fY (y) = fX (x) dx = fX ( y) ( y) − 0 = 1 − · √ =
dy 0 dy 2 2 y

(i) 3
2
− 2y
32
(ii) 2√ 1
y
− 14 (iii) √1y − 14 (iv) 2
1

y
+ 14 ,

Shorten up number of steps, omitted FX ( y).


where y = x2 , 0 ≤ x ≤ 2 implies 0 ≤ y ≤ 2, or
(i) 0 ≤ y ≤ 1 (ii) 0 ≤ y ≤ 2 (iii) 0 ≤ y ≤ 3 (iv) 0 ≤ y ≤ 4

Although there are two roots, ± y, only the positive root fits in the positive interval 0 ≤ y ≤ 2.
100 Chapter 3. Continuous Random Variables (LECTURE NOTES 5)

2. Simulate random sample from desired cdf F (x) using uniform.


Simulate a random sample of size 5 from X with desired cdf (not pdf),

0, x<3
F (x) = 3 3

1− x , x≥3

(a) sample point y


1
chosen at random from
(assume) uniform U(0,1) y
cdf: F(x) = P(X < x)

3 x6 9 12
-1
inverse cdf: F (y) = P(Y < y)
(b) generates sample point x chosen
-1
at random from desired cdf F (y)

Figure 3.13: Distribution y = F (x) and inverse x = F −1 (y)

(a) Determine cdf of X, desired F (x).


3
It is given (does not need to be worked out in this case): F (x) = 1 − x3
3
(b) Find inverse cdf F −1 (y). Since F (x) = 1 − x3 = y, then
3 3 1 1
x = F −1 (y) = (i) x3 (ii) 1− x9 (iii) 1−3(1−y)− 3 (iv) 3(1−y)− 3

(c) Generate values y1 , y2 , y3 , y4 , y5 from Y , assumed to be U (0, 1).


One possible random sample of size 5 from uniform on U (0, 1) gives
0.2189, 0.9913, 0.5775, 0.2846, 0.8068.

y <- runif(5); y
[1] 0.2189 0.9913 0.5775 0.2846 0.8068

(d) Use x = F −1 (y) to generate observed x1 , x2 , . . . , xn , from desired F (x).


Complete the following table.
y sample from U (0, 1) 0.2189 0.9913 0.5775 0.2846 0.8068
1
x = F −1 (y) = 3(1 − y)− 3 3.2575 14.586

y <- c(0.2189,0.9913,0.5775,0.2846,0.8068)
x <- 3*(1-y)^{-1/3}; x

[1] 3.257510 14.586402 3.998027 3.354323 5.189421

You might also like