Statistics I · Lecture 11
These slides are a free translation and adaptation from the slide deck for Estatística I by Prof. Sandra Custódio and Prof. Teresa Ferreira from the Lisbon Accounting and Business School, Polytechnic University of Lisbon.
A r.v. \(X\) follows a uniform distribution in \([a,b]\subset\mathbb{R}\) with \(-\infty<a<b<\infty\), if its probability density is given by:
\[ f_X(x)=\begin{cases} \frac{1}{b-a} & a\leq x \leq b\\ 0 & otherwise \end{cases} \]
We write \(X\sim U[a,b]\)
The cumulative distribution function is given by:
\[ F_X(x)=\begin{cases} 0 & x< a\\ \frac{x-a}{b-a} & a\leq x < b\\ 1 & x\geq b \end{cases} \]
Note, this distribution is symmetric, and its first two moments are:
The length of small spots in a TV network is a r.v. \(X\) distributed \(U[5,12]\).
Let the r.v. \(X\) be distributed \(U[2,b]\) with \(b>2\). What value must \(b\) take to make \(P(3\leq X\leq 5)=0.4\)?
If \(X\sim U[a,b]\), then \(E[X]\) equals:
A. \(\frac{a+b}{2}\)
B. \(\frac{b-a}{2}\)
C. \(\frac{1}{b-a}\)
D. \(\frac{(b-a)^2}{12}\)
✅ A. The Uniform is symmetric, so its mean sits at the midpoint of the interval. \(\frac{(b-a)^2}{12}\) is the variance.
For \(X\sim U[0,10]\), the probability \(P(X=5)\) is:
A. \(0.1\)
B. \(1\)
C. \(0.5\)
D. \(0\)
✅ D. \(X\) is continuous, so every single point has probability zero, no matter how flat the density is.
Let \(X\sim U[5,15]\).
Compute \(E[X]\), \(V[X]\), \(P(X\leq 8)\) and \(P(7<X\leq 12)\).
\(E[X]=\frac{5+15}{2}=10\) · \(V[X]=\frac{(15-5)^2}{12}=\frac{100}{12}\approx 8.33\)
\(F_X(x)=\frac{x-5}{10}\) on \([5,15]\), so \(P(X\leq 8)=\frac{3}{10}=0.3\)
\(P(7<X\leq 12)=\frac{12-7}{10}=0.5\)
The exponential distribution is rooted in the Poisson distribution, reflecting the waiting time between events originated according to a Poisson process.
Nevertheless, we can apply the exponential distribution to many other phenomena.
A r.v \(X\) is distributed exponentially, with parameter \(\lambda>0\), \(X\sim Exp(\lambda)\), if its probability density function is given by:
\[ f_X(x)=\begin{cases} \lambda e^{-\lambda x} & x\geq 0\\ 0 & x< 0 \end{cases} \]
\(\lambda\) is the rate: the expected number of events per unit of time (or space). The expected waiting time between events is \(\frac{1}{\lambda}\).
The cumulative probability function is:
\[ F_X(x)=\begin{cases} 0 & x< 0\\ 1-e^{-\lambda x} & x\geq 0 \end{cases} \]
And its moments are:
Property
Lack of memory of the exponential distribution:
Let the r.v. \(X\sim Exp(\lambda)\), then:
\[P(X>x+h|X>x)=P(X>h)\]
with \(x,h>0\)
Considering the survival application of this distribution, this property states that, the time left to live is independent of what it already lived.
Grey = time you have already waited with nothing happening. Blue = the original distribution from time zero. Orange dashed = the distribution of what is still to come.
Drag \(x\) anywhere you like. The orange curve slides along with it, but its shape never changes: it is the blue curve translated, not a flatter or steeper one. And the two probabilities in the panel stay equal to four decimals.
Algebraically: \(P(X>x+h\mid X>x)=\frac{e^{-\lambda(x+h)}}{e^{-\lambda x}}=e^{-\lambda h}=P(X>h)\). The \(e^{-\lambda x}\) cancels, and with it everything the past could have told you.
A component that has run 10 years is exactly as likely to survive another year as a brand new one. That is often the wrong model for wear, which is precisely why you must check whether the exponential fits before using it.
In a factory, the execution time of a piece is random variable distributed exponentially with expected value of 5 minutes.
The time it takes until the first consultation, and between consultations, in the clinic of Dr. Shawn are independent and distributed exponentially with \(\lambda=0.1\).
What is the probability that no consultation occurs before the first 10 minutes?
Buses arrive following an Exponential model with \(\lambda=2\) per hour. The expected waiting time for the next bus is:
A. 2 hours
B. 0.5 hours
C. 4 hours
D. 1 hour
✅ B. \(\lambda\) is a rate, here 2 buses per hour. The expected waiting time is \(\frac{1}{\lambda}=0.5\) hours, that is 30 minutes.
Memorylessness says that \(P(X>x+h\,|\,X>x)\) equals:
A. \(P(X>h)\)
B. \(P(X>x)\)
C. \(P(X>x+h)\)
D. \(0\)
✅ A. Having already waited \(x\) tells you nothing. The clock restarts, so only the extra \(h\) matters.
Buses arrive with \(\lambda=0.05\) per minute.
Compute the expected waiting time, \(P(X>10)\), and \(P(X>25\,|\,X>15)\).
\(E[X]=\frac{1}{0.05}=20\) minutes.
\(P(X>10)=1-F(10)=e^{-0.05\times 10}=e^{-0.5}\approx 0.6065\)
By memorylessness, \(P(X>25|X>15)=P(X>10)\approx 0.6065\), the same number. ✅
The Gaussian or Normal Distributions is one of the most used distributions, playing a key role in statistical inference.
The r.v. \(X\) is normally distributed, \(X\sim N(\mu,\sigma^2)\), if its density and cumulative probability distribution functions:
\[ f(x)=\frac{1}{\sigma\sqrt{2\pi}}e^{-\frac{(x-\mu)^2}{2\sigma^2}}\\ P(X\leq a)=F(a)=\int_{-\infty}^a f(t)dt \]
Where \(E[X]=\mu\) and \(Var[X]=\sigma^2\)
If \(\mu=0\) and \(\sigma^2=1\) we call this the standard normal distribution, where \(Z\sim N(0,1)\).
Theorem
Let \(X\sim N(\mu,\sigma^2)\). Let \(Z=\frac{X-\mu}{\sigma}\), then \[Z\sim N(0,1)\]
That is, we can standardize a r.v. distributed Normally. Its probability function is denoted \(\Phi\), and its density is \(\phi\).
\[ \phi(z)= \frac{1}{\sqrt{2\pi}}e^{-\frac{z^2}{2}} \\ \Phi(z)=P(Z\leq z) = \int_{-\infty}^z \frac{1}{\sqrt{2\pi}}e^{-\frac{t^2}{2}} dt \]
How do we know if a r.v. follows a Normal distribution?
Left: your \(X\sim N(\mu,\sigma^2)\). Right: the same problem after standardizing. The two shaded areas are always equal, which is why one table serves every Normal.
Change \(\mu\) and the left curve slides; change \(\sigma\) and it widens or narrows. The right curve never moves: standardizing has already absorbed both.
The orange cut-off tracks the same point in both pictures, at \(x\) on the left and at \(z=\frac{x-\mu}{\sigma}\) on the right, and the shaded areas stay equal to the last decimal.
This is why one page of tabulated \(\Phi(z)\) is enough for every Normal you will ever meet. You never tabulate \(X\), you convert it. Note the readout tells you the exact row and column to look up.
Consider the r.v. \(Z\sim N(0,1)\). Find \(P(Z\leq 1.65)\).
\[P(Z\leq 1.65) = F(1.65)=\Phi(1.65)=0.9505\]
Warning
Consider the r.v. \(X\sim N(6,25)\). Find \(P(X\leq 12)\).
\[P(X\leq 12)=P\left(\frac{X-6}{5}\leq \frac{12-6}{5}\right)\]
Or
\[P(Z\leq 1.2)=\Phi(1.2)=0.8849\]
The exercises go from easy to hard:
Keep the Normal table open.
Let \(X\sim Exp(0.5)\). Find \(P(X>2)\).
\[P(X>2)=1-F(2)=e^{-0.5\times 2}=e^{-1}\approx 0.3679\]
Let \(Z\sim N(0,1)\). Find:
(a) \(P(0<Z\leq 2.05)\) (b) \(P(-1.22\leq Z<1.05)\) (c) \(P(Z\geq -2.05)\)
(d) \(k\) such that \(P(|Z|>k)=0.05\) (e) \(k\) such that \(P(|Z|<k)=0.90\)
(a) \(\Phi(2.05)-\Phi(0)=0.9798-0.5=0.4798\)
(b) \(\Phi(1.05)-\Phi(-1.22)=0.8531-(1-0.8888)=0.7419\)
(c) By symmetry, \(P(Z\geq -2.05)=P(Z\leq 2.05)=0.9798\)
(d) \(0.025\) in each tail, \(\Phi(k)=0.975\Rightarrow k=1.96\) (e) \(0.05\) in each tail, \(\Phi(k)=0.95\Rightarrow k\approx 1.645\)
The number of cars 🚗 checked in an inspection center, per hour, is distributed Poisson with an expected value of 4.
(a) What is the expected time between two cars?
(b) Find the probability of waiting more than 30 minutes for the next car.
If the number of events per hour is \(Poi(4)\), the time between events, in hours, is \(T\sim Exp(4)\).
(a) \(E[T]=\frac{1}{4}\) of an hour, 15 minutes.
(b) \(P(T>0.5)=e^{-4\times 0.5}=e^{-2}\approx 0.1353\). Or, with the Poisson: no cars in 30 minutes, \(P(Y=0)\) with \(Y\sim Poi(2)\). Same number!
The bus 🚌 arrives at your stop at a random moment between 0 and 20 minutes from now, all moments equally likely. Let \(W\) be your waiting time.
(a) Find \(P(W>15)\).
(b) You have waited 10 minutes already. Find \(P(W>15|W>10)\).
(c) Find \(E[W]\) and \(V[W]\).
\(W\sim U(0,20)\)
(a) \(P(W>15)=\frac{20-15}{20}=0.25\) (b) \(P(W>15|W>10)=\frac{P(W>15)}{P(W>10)}=\frac{0.25}{0.5}=0.5\)
(c) \(E[W]=\frac{0+20}{2}=10\) and \(V[W]=\frac{(20-0)^2}{12}\approx 33.33\)
The lifetime of a light bulb 💡, in hours, is distributed exponentially with mean 1000 hours.
(a) Find the probability that a bulb lasts more than 1500 hours.
(b) Find the median lifetime. Is it above or below the mean?
(c) A box has 4 bulbs, independent. Find the probability that at least one of them lasts more than 1500 hours.
\(T\sim Exp\left(\frac{1}{1000}\right)\)
(a) \(P(T>1500)=e^{-1500/1000}=e^{-1.5}\approx 0.2231\)
(b) \(F(m)=1-e^{-m/1000}=0.5\Rightarrow m=1000\ln 2\approx 693\) hours. Below the mean: a few bulbs last very long and pull the mean up.
(c) \(N\sim Bin(4,0.2231)\) counts the bulbs lasting more than 1500 hours: \(P(N\geq 1)=1-0.7769^4\approx 0.6358\)
Go back to the bus 🚌 of Exercise 4. Now suppose the waiting time is instead \(W\sim Exp(0.1)\), so it has the same mean, 10 minutes.
(a) Find \(P(W>15)\) and \(P(W>15|W>10)\).
(b) Compare with Exercise 4. With which bus is it worth to keep waiting?
(a) \(P(W>15)=e^{-1.5}\approx 0.2231\) and, by the lack of memory, \(P(W>15|W>10)=P(W>5)=e^{-0.5}\approx 0.6065\)
(b) With the uniform bus, having waited 10 minutes makes it more likely the bus arrives soon, it must arrive before minute 20. With the exponential bus, the 10 minutes you already waited are worth nothing 😩, you are back at the start.
Problem set 3.2, Questions 1 to 7.
| z | .00 | .01 | .02 | .03 | .04 | .05 | .06 | .07 | .08 | .09 |
|---|---|---|---|---|---|---|---|---|---|---|
| 0.0 | 0.5000 | 0.5040 | 0.5080 | 0.5120 | 0.5160 | 0.5199 | 0.5239 | 0.5279 | 0.5319 | 0.5359 |
| 0.1 | 0.5398 | 0.5438 | 0.5478 | 0.5517 | 0.5557 | 0.5596 | 0.5636 | 0.5675 | 0.5714 | 0.5753 |
| 0.2 | 0.5793 | 0.5832 | 0.5871 | 0.5910 | 0.5948 | 0.5987 | 0.6026 | 0.6064 | 0.6103 | 0.6141 |
| 0.3 | 0.6179 | 0.6217 | 0.6255 | 0.6293 | 0.6331 | 0.6368 | 0.6406 | 0.6443 | 0.6480 | 0.6517 |
| 0.4 | 0.6554 | 0.6591 | 0.6628 | 0.6664 | 0.6700 | 0.6736 | 0.6772 | 0.6808 | 0.6844 | 0.6879 |
| 0.5 | 0.6915 | 0.6950 | 0.6985 | 0.7019 | 0.7054 | 0.7088 | 0.7123 | 0.7157 | 0.7190 | 0.7224 |
| 0.6 | 0.7257 | 0.7291 | 0.7324 | 0.7357 | 0.7389 | 0.7422 | 0.7454 | 0.7486 | 0.7517 | 0.7549 |
| 0.7 | 0.7580 | 0.7611 | 0.7642 | 0.7673 | 0.7704 | 0.7734 | 0.7764 | 0.7794 | 0.7823 | 0.7852 |
| 0.8 | 0.7881 | 0.7910 | 0.7939 | 0.7967 | 0.7995 | 0.8023 | 0.8051 | 0.8078 | 0.8106 | 0.8133 |
| 0.9 | 0.8159 | 0.8186 | 0.8212 | 0.8238 | 0.8264 | 0.8289 | 0.8315 | 0.8340 | 0.8365 | 0.8389 |
| 1.0 | 0.8413 | 0.8438 | 0.8461 | 0.8485 | 0.8508 | 0.8531 | 0.8554 | 0.8577 | 0.8599 | 0.8621 |
| 1.1 | 0.8643 | 0.8665 | 0.8686 | 0.8708 | 0.8729 | 0.8749 | 0.8770 | 0.8790 | 0.8810 | 0.8830 |
| 1.2 | 0.8849 | 0.8869 | 0.8888 | 0.8907 | 0.8925 | 0.8944 | 0.8962 | 0.8980 | 0.8997 | 0.9015 |
| 1.3 | 0.9032 | 0.9049 | 0.9066 | 0.9082 | 0.9099 | 0.9115 | 0.9131 | 0.9147 | 0.9162 | 0.9177 |
| 1.4 | 0.9192 | 0.9207 | 0.9222 | 0.9236 | 0.9251 | 0.9265 | 0.9279 | 0.9292 | 0.9306 | 0.9319 |
| 1.5 | 0.9332 | 0.9345 | 0.9357 | 0.9370 | 0.9382 | 0.9394 | 0.9406 | 0.9418 | 0.9429 | 0.9441 |
| 1.6 | 0.9452 | 0.9463 | 0.9474 | 0.9484 | 0.9495 | 0.9505 | 0.9515 | 0.9525 | 0.9535 | 0.9545 |
| 1.7 | 0.9554 | 0.9564 | 0.9573 | 0.9582 | 0.9591 | 0.9599 | 0.9608 | 0.9616 | 0.9625 | 0.9633 |
| 1.8 | 0.9641 | 0.9649 | 0.9656 | 0.9664 | 0.9671 | 0.9678 | 0.9686 | 0.9693 | 0.9699 | 0.9706 |
| 1.9 | 0.9713 | 0.9719 | 0.9726 | 0.9732 | 0.9738 | 0.9744 | 0.9750 | 0.9756 | 0.9761 | 0.9767 |
| 2.0 | 0.9772 | 0.9778 | 0.9783 | 0.9788 | 0.9793 | 0.9798 | 0.9803 | 0.9808 | 0.9812 | 0.9817 |
| 2.1 | 0.9821 | 0.9826 | 0.9830 | 0.9834 | 0.9838 | 0.9842 | 0.9846 | 0.9850 | 0.9854 | 0.9857 |
| 2.2 | 0.9861 | 0.9864 | 0.9868 | 0.9871 | 0.9875 | 0.9878 | 0.9881 | 0.9884 | 0.9887 | 0.9890 |
| 2.3 | 0.9893 | 0.9896 | 0.9898 | 0.9901 | 0.9904 | 0.9906 | 0.9909 | 0.9911 | 0.9913 | 0.9916 |
| 2.4 | 0.9918 | 0.9920 | 0.9922 | 0.9925 | 0.9927 | 0.9929 | 0.9931 | 0.9932 | 0.9934 | 0.9936 |
| 2.5 | 0.9938 | 0.9940 | 0.9941 | 0.9943 | 0.9945 | 0.9946 | 0.9948 | 0.9949 | 0.9951 | 0.9952 |
| 2.6 | 0.9953 | 0.9955 | 0.9956 | 0.9957 | 0.9959 | 0.9960 | 0.9961 | 0.9962 | 0.9963 | 0.9964 |
| 2.7 | 0.9965 | 0.9966 | 0.9967 | 0.9968 | 0.9969 | 0.9970 | 0.9971 | 0.9972 | 0.9973 | 0.9974 |
| 2.8 | 0.9974 | 0.9975 | 0.9976 | 0.9977 | 0.9977 | 0.9978 | 0.9979 | 0.9979 | 0.9980 | 0.9981 |
| 2.9 | 0.9981 | 0.9982 | 0.9982 | 0.9983 | 0.9984 | 0.9984 | 0.9985 | 0.9985 | 0.9986 | 0.9986 |
| 3.0 | 0.9987 | 0.9987 | 0.9987 | 0.9988 | 0.9988 | 0.9989 | 0.9989 | 0.9989 | 0.9990 | 0.9990 |
In this case the r.v. \(X\) is the time until the first consultation, and \(X\sim Exp(0.1)\).
\[P(X>10)=1-P(X<10)=1-F(10)\]
Or
\[1-(1-e^{-0.1\times 10})=e^{-1}\approx 0.3679\]