17 de Setembro de 2015
Distribuição de probabilidade de uma \(Bin(10, p)\), com \(p=0.2, 0.5\) e \(0.8\).
Definição: a função de densidade de probabilidade de uma v.a. \(X\) é uma função \(f\) que verifica:
Toda v.a. \(X\) à qual seja possÃvel associar uma função de densidade de probabilidade será chamada de v.a. contÃnua.
A probabilidade de que uma v.a. \(X\) contÃnua pertença a um intervalo da reta (a,b], \(a<b\) é dada por:
\[P(a<X\leq b)=\int_{a}^{b}f(x)dx\]
Obs: Quando \(X\) é v.a. contÃnua, \(P(a<X\leq b)=P(a\leq X\leq b)=P(a<X< b)=P(a\leq X <b).\)
Notação: se \(X\) v.a. contÃnua com função de densidade de probabilidade \(f\), no lugar de \(f\) denotaremos \(f_{X}\).
A função de distribuição acumulada (f.d.a.) de uma v.a. contÃnua \(X\), que denotaremos por \(F_{X}\), é a probabilidade de que \(X\) pertença ao intervalo da reta (\(-\infty,x\)], ou seja, \[F_{X}(x)=P(X\leq x)=\int_{-\infty}^{x}f_X(u)du\]
Exemplo: Seja \(X\) uma v.a. contÃnua com a seguinte função de densidade:
\[ f_{X}(x) = \left\{ \begin{array}{ll} x, & 0\leq x < 1 \\ 2-x, & 1 \leq x \leq 2 \\ 0, & \mbox{caso contrário} \\ \end{array} \right.\]
Pela definição de função de densidade:
Propriedade: toda v.a. \(X\) contÃnua (ou seja, que possui \(f_{X}\) como função de densidade) tem probabilidade pontual nula: \(P(X=x)=0.\)
Resumindo: \(F(x)=P(X\leq x)\)
Dada a função
\[ f(x) = \left\{ \begin{array}{rcl} 0 & \mbox{se} & x < 0\\ 2e^{-2x} & \mbox{se} & x \geq 0 \\ \end{array} \right. \]
Mostre que esta é uma função de densidade de probabilidade.
Calcule a probabilidade de \(X>10\).
Fonte: Morettin & Bussab, EstatÃstica Básica \(5^a\) edição, pág 166.
Uma função de densidade de probabilidade deve satisfazer as seguintes propriedades:
\(f(x)\geq0,\) para todo \(x \in \mathbb{R}.\)
\(\displaystyle \int_{-\infty}^{\infty} f(x) dx = 1.\)
Note que \(e^{-x}\) é positiva para qualquer \(x\) e, consequentemente, \(2e^{-2x}\) também.
Resta mostrar que sua integral é 1:
\[ \int 2e^{-2x} dx = -e^{-2x}\]
Note que a função está definida nesta forma para \(x \geq 0\); para \(x<0\), ela é 0.
Então a integral é:
\[ \int_{-\infty}^{\infty} f(x) dx = \int_{-\infty}^0 0 dx + \int_0^{\infty} 2e^{-2x} dx = \]
\[ = \left [ -e^{-2x} \right ] _{0}^{\infty} = \lim_{x\rightarrow \infty} -e^{-2x} - \left ( -e^{-0} \right) = 1 \]
A probabilidade é dada por: \[P(X>10) = \int_{10}^{\infty} 2e^{-2x} dx = \lim_{x\rightarrow \infty} -e^{-2x} - \left ( -e^{(-2)10} \right) = \frac{1}{e^{20}}\]
Definição: Seja \(X\) uma v.a. contÃnua com densidade \(f_{X}\), a esperança de \(X\) é dada por: \[\mathbb E(X)= \int_{-\infty}^{\infty}x\,f_{X}(x)dx \]
Definição: Seja \(X\) uma v.a. contÃnua com densidade \(f_{X}\), o \(k\)-ésimo momento de \(X\) é dado por:
\[\mathbb E(X^k)= \int_{-\infty}^{\infty}x^k\,f_{X}(x)dx \]
Definição: Seja \(X\) uma v.a. com valor esperado \(\mathbb E(X)\), definimos por variância, a quantidade:
\[Var(X)= \mathbb E[(X - \mathbb E(X))^2] = \mathbb E(X^2) - [\mathbb E(X)]^2\]
E definimos como desvio padrão: \[\sigma(X) = \sqrt{Var(X)}\]
Note que assim como no caso discreto, ambas as quantidades oferecem medidas de dispersão da variável \(X\) em relação ao valor esperado \(\mathbb E(X)\).
Considere a seguinte função de densidade: \[ f_{X}(x) = \begin{cases} x, & 0\leq x < 1 \\ 2-x, & 1 \leq x \leq 2 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
Encontre a \(\mathbb E(X)\) e \(Var(X).\)
\[\mathbb E(X)= \int_{-\infty}^{\infty}xf_{X}(x)dx= \int_{0}^{1}x^{2}dx+ \int_{1}^{2}x(2-x)dx= 1\]
\[Var(X)= \int_{-\infty}^{\infty}[x-\mathbb E(X)]^{2}f_{X}(x)dx= \int_{-\infty}^{\infty}[x-1]^{2}f_{X}(x)dx =\frac{1}{6}\]
Uma v.a. \(X\) tem distribuição triangular no intervalo \([0,1]\) se sua função densidade for dada por:
\[ f(x) = \begin{cases} 0 & \mbox{se} & x < 0\\ Cx & \mbox{se} & 0 \leq x \leq 1/2 \\ C(1-x) & \mbox{se} & 1/2 \leq x \leq 1 \\ 0 & \mbox{se} & x > 1\\ \end{cases} \]
Fonte: Morettin & Bussab, EstatÃstica Básica \(5^a\) edição, pág 166.
Item 1 - Devemos escolher \(C\) de modo que \(f(x)\) satisfaça:
\(f(x)\geq0,\) para todo \(x \in \mathbb{R}\).
\(\int_{-\infty}^{\infty} f(x) dx = 1\).
Pela primeira condição, temos que \(C>0\). Agora, para que \(C\) satisfaça a segunda condição, devemos integrar \(f(x)\):
\[ \begin{aligned} \int_{-\infty}^{\infty} f(x) dx &= \int_{-\infty}^{0}\!\! 0 dx + \int_0^{1/2}\!\! Cx dx + \int_{1/2}^1\!\! C(1-x) dx + \int_{1}^{\infty}\!\! 0 dx\\ & = C \int_0^{1/2}\!\! x dx + C \int_{1/2}^1\!\! (1-x) dx = C \left ( \left[ \frac{x^2}{2} \right]_0^{1/2} + \left[ x- \frac{x^2}{2} \right]_{1/2}^1 \right )\\ & = C \left( \frac{1}{8} +1 - \frac{1}{2} - \frac{1}{2} + \frac{1}{8} \right) = \frac{C}{4} \quad \Rightarrow \quad C \mbox{ deve ser igual a 4.} \end{aligned} \]
Item 2 - Densidade da distribuição triangular no intervalo \([0,1]\)
Item 3 - Para encontrarmos as probabilidades dos eventos, basta integrar nas regiões correspondentes: \[P\left(X \leq \frac{1}{2} \right) = \int_0 ^{1/2} f(x)dx = \int_0 ^{1/2} 4x dx = \frac{1}{2}\]
Note que: \(P\left(X > \frac{1}{2} \right) = 1-P\left(X \leq \frac{1}{2} \right) = 1-\frac{1}{2} = \frac{1}{2}\)
\[\begin{aligned} P\left(\frac{1}{4} \leq X \leq \frac{3}{4}\right) &= \int_{1/4} ^{3/4} f(x)dx \\ &= \int_{1/4} ^{1/2} 4xdx + \int_{1/2} ^{3/4} 4(1-x)dx = \frac{3}{4} \end{aligned} \]
Calcule a esperança, a variância e a f.d.a. da v.a. \(X\) com a densidade triangular em \([0,1]\).
Fonte: Morettin & Bussab, EstatÃstica Básica \(5^a\) edição, pág 171.
Esperança \[ \begin{aligned} \mathbb E(X) &= \int_{-\infty}^{\infty} \!\! x f(x) dx = \int_0^{1/2} \! 4x^2dx + \int_{1/2}^{1} \! 4x(1-x)dx\\ & = \left[ \frac{4x^3}{3} \right]_0^{1/2} + \left[ \frac{2}{3} x^2(3-2x) \right]_{1/2}^{1} = \frac{1}{2} \end{aligned} \]
Variância \[ \begin{aligned} \mbox{Var}(X)& = \int_{-\infty}^{\infty} \!\! (x-\mathbb{E}(X))^2 f(x) dx \\ &= \int_0^{1/2} \! 4 \left(x-\frac{1}{2}\right)^2 \! xdx + \int_{1/2}^{1} \! 4 \left(x-\frac{1}{2}\right)^2 \! (1-x)dx\\ &= \left[ x^4-\frac{4}{3}x^3 + \frac{1}{2}x^2 \right]_0^{1/2} + \left[ -x^4+\frac{8}{3}^3-\frac{5}{2}x^2 +x \right]_{1/2}^{1}\\ &= \frac{1}{24} \end{aligned} \]
A função de distribuição acumulada de uma variável aleatória contÃnua é dada por \[F(x) = \int_{-\infty}^x f(t) dt\]
Temos que para \(x \in [0,1/2)\), \(F(x)\) é dada por \[F(x) = \int_0^x 4t dt = 2x^2\]
Para \(x \in [1/2,1]\), a acumulada é dada por \[F(x) = \int_0^{1/2} 4t dt + \int_{1/2}^x 4(1-t)dt = -2x^2 + 4x -1\]
Para valores de \(x \geq 1\), a acumulada assume valor 1.
O gráfico de \(F(x)\) é dado por:
\[ f_{X}(x) = \begin{cases} \frac{1}{b-a}, & a \leq x \leq b \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
Notação: \(X \sim U[a,b]\) ou \(X \sim U(a,b)\).
Cálculo da função de distribuição acumulada:
\[ F_{X}(x) = \left\{ \begin{array}{ll} 0, & x < a \\ \int_{a}^{x}\frac{1}{b-a}dt=\frac{x-a}{b-a}, & a \leq x \leq b \\ 1, & x > b \\ \end{array} \right.\]
Gráficos da função de densidade de probabilidade e da função de distribuição acumulada de uma v.a. \(X \sim U[a,b]\):
Se \(X \sim U[a,b]\), então a esperança e variância são dadas por:
Esperança: \[\mathbb E(X)= \int_{a}^{b}x\frac{1}{(b-a)}dx=\frac{(b+a)}{2}\]
Variância \[\mathbb E(X^{2})= \int_{a}^{b}x^{2}\frac{1}{(b-a)}dx= \frac{(a^{2}+ab+b^{2})}{3}\]
\[Var(X)=\mathbb E(X^{2})-[\mathbb E(X)]^{2}= \frac{(a^{2}+ab+b^{2})}{3}-\frac{(b+a)^{2}}{4}= \frac{(b-a)^{2}}{12}\]
A dureza \(H\) de uma peça de aço pode ser pensada como uma v.a. com distribuição uniforme no intervalo \([50,70]\) da escala de Rockwel. Então:
\[ f_{H}(h) = \begin{cases} \displaystyle \frac{1}{70-50}=\frac{1}{20}, & 50 \leq x \leq 70 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
Calcule a probabilidade de que uma peça tenha dureza entre 55 e 60.
\[\begin{aligned} P(55<H<60) &= \int_{55}^{60}\frac{1}{20}dh \\ &= \frac{1}{20}(60-55)= \frac{1}{4} \end{aligned} \]
Seja \(X\) uma v.a. distribuÃda uniformemente, com média 15 e variância 25/3.
Encontre a função de densidade de \(X\).
Qual é a probabilidade de que \(X > 14\)?
Lembre-se que a esperança e variância de uma v.a. \(X \sim U[a,b]\) são dadas por \[\mathbb E(X) = \frac{a+b}{2} \qquad \mbox{e} \qquad Var(X) = \frac{(b-a)^2}{12}\]
Temos o seguinte sistema, portanto:
\[ \begin{cases} \displaystyle \frac{a+b}{2} &=& 15 \\ \displaystyle \frac{(b-a)^2}{12} &=& \displaystyle \frac{25}{3} \\ \end{cases} \qquad \Longleftrightarrow \qquad \begin{cases} a+b &=& 30 \\ (b-a)^2 &=& 100 \\ \end{cases} \]
Ou simplesmente (você é capaz de dizer por que tomamos a raiz positiva apenas, neste sistema não-linear?)
\[ \left\{ \begin{array}{ccr} a + b & = & 30 \\ b - a & = & 10 \\ \end{array} \right. \]
O sistema tem solução \(a=10\), \(b=20\), o que nos mostra que \(X \sim U[10,20]\) e
\[ f_{X}(x) = \begin{cases} \frac{1}{10}, & 10 \leq x \leq 20 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
A probabilidade de que \(X>14\) é dada por \[P(X>14) = \int_{14}^{20}\frac{1}{10}dx = \frac{(20-14)}{10} = 0.6\]
\[ f_{X}(x) = \begin{cases} \lambda e^{-\lambda x}, & x \geq 0 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
\[ F_{X}(x) = \begin{cases} \int_{0}^{x}\lambda e^{-\lambda t}dt=1-e^{-\lambda x}, & x \geq 0 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
Gráficos da função de densidade de probabilidade (esquerda) e da função de distribuição acumulada (direita) de \(X \sim Exp(0.5)\):
Se \(X \sim Exp(\lambda)\), a esperança e variância são dadas por:
Esperança: \[\mathbb E(X)=\int_{0}^{\infty}x\lambda e^{-\lambda x}dx=\frac{1}{\lambda}\]
Variância:
\[\mathbb E(X^{2})=\int_{0}^{\infty}x^{2}\lambda e^{-\lambda x}dx=\frac{2}{\lambda^{2}}\]
\[Var(X)=\mathbb E(X^{2})-[\mathbb E(X)]^{2}=\frac{2}{\lambda^{2}}-\frac{1}{\lambda^{2}}=\frac{1}{\lambda^{2}}\]
O tempo de vida, \(X\), em horas, de um componente eletrônico segue uma distribuição exponencial de tal forma que \(P(X\leq 1000) = 0.75\).
Qual é o tempo médio de vida do componente?
Sabemos que se \(X \sim exp(\lambda)\), então \[F(x) = P(X \leq x) = 1-e^{-\lambda x} \quad \mbox{ e } \quad \mathbb E(X) = \lambda^{-1}\]
Basta então observarmos que
\[P(X \leq 1000) = 1 - e^{-\lambda 1000} = 0.75 \Leftrightarrow \lambda = \frac{\ln(4)}{1000} = 0.0013863\]
Concluimos então que o tempo médio de vida (em horas) é \[\mathbb E(X) = \frac{1}{0.0013863} = 721.3475,\] e que \(75\%\) dos componentes duram 1000 horas ou menos.
Um antiga fábrica de tubos de TV determinou que a vida média dos tubos de sua fabricação é de \(800\) horas de uso contÃnuo e segue uma distribuição exponencial.
Qual a probabilidade de que a fábrica tenha que substituir um tubo gratuitamente, se oferece uma garantia de \(300\) horas de uso?
\(X\): vida útil do tubo de TV.
\(\mathbb E[X]= 800\). Como \(X \sim Exp(\lambda)\) e \(\mathbb E[X]= \frac{1}{\lambda}=800\), portanto \(\lambda= \frac{1}{800}.\)
\[ f_{X}(x) = \begin{cases} \lambda e^{-\lambda x}, & x \geq 0 \\ 0, & \mbox{caso contrário} \\ \end{cases} \quad = \quad \begin{cases} \frac{1}{800} e^{-\frac{x}{800}}, & x \geq 0 \\ 0, & \mbox{caso contrário} \\ \end{cases} \]
Se \(X<300\), a fábrica tem que substituir gratuitamente.
\[P(X<300)= \int_0^{300}\frac{1}{800}e^{-\frac{x}{800}}dx =\left[-e^{-\frac{x}{800}} \right]_0^{300}= 0.3127\]
A seguinte função de distribuição de probabilidade \[ f_{X}(x) =
\begin{cases}
2 e^{-2 x}, & x \geq 0 \\
0, & \mbox{caso contrário} \\
\end{cases}
\]
representa a distribuição do Ãndice de acidez (\(X\)) de um determinado produto alimentÃcio.
O produto é consumÃvel se este Ãndice for menor do que \(2\).
O setor de fiscalização apreendeu \(30\) unidades do produto. Qual a probabilidade de que pelo menos \(10\%\) da amostra seja imprópria para consumo?
Produto liberado para consumo se: \(P(X<2)= \int_0^22 e^{-2 x}dx= \left[-e^{-2x} \right]_0^{2}=0.98\)
Produto não consumÃvel com probabilidade \(p=P(X \geq 2) = 1-P(X<2)= 0.02\)
\(Y\): número de unidades impróprias para consumo na amostra de 30
Então \(Y\sim \mbox{Bin}(30,0.02)\)
Probabilidade de que pelo menos 10% de uma amostra de 30 unidades seja imprópria:
\[ \begin{aligned} P(Y\geq3) &= 1-P(Y<3)=1-[P(Y=0)+P(Y=1)+P(Y=2)]\\ &= 1-\left[\binom{30}{0}(0.02)^0(0.98)^{30} + \binom{30}{1}(0.02)^1(0.98)^{29}+\binom{30}{2}(0.02)^2(0.98)^{28} \right]\\ &= 0.022 \end{aligned} \]