#7Vectores aleatorios

Clase 7 - Independencia de variables aleatorias, cambio de variables y distribución Gamma

22 min de lectura

Obs:Obs:

xx e yy absolutamente continuas ̸  ⟹  (x,y)\not\implies (x,y) absolutamente continuo.

Por ejemplo: X∼U[0,1]X \sim \mathcal{U}[0,1] Y=XY=X , (x,y)=(x,x)(x,y)=(x,x), P((x,x)∈s⏟=1)≠∫sf=0P(\underbrace{ (x,x) \in s }_{ =1 }) \neq \int_{s} f =0


Independencia

x1,…,xnx_{1},\dots,x_{n} son independientes.

P(x1∈A1,x2∈A2,…,xn∈An)=∏P(Xi∈Ai)Para todo A1,…,An bolerianosP(x_{1} \in A_{1}, x_{2} \in A_{2},\dots,x_{n} \in A_{n})= \prod P(X_{i} \in A_{i}) \quad \text{Para todo }A_{1},\dots,A_{n}\text{ bolerianos}

Obs:Obs: Si x1,…,xnx_{1},\dots,x_{n}son independientes cualquier subconjunto lo es. Si tengo infinitas variables, son independientes si toda subfamilia lo es.


Prop:Prop: Si x1,…,xnx_{1},\dots,x_{n} son independientes.

Fx1,..,xn(t1,…,tn)=Fx1(t1)…Fxn(tn)F_{x_{1},..,x_{n}}(t_{1},\dots,t_{n}) = F_{x_{1}}(t_{1}) \dots F_{x_{n}} ( t_{n})

Vale que Si FX(t‾)=F1(t1)⋅F2(t2)…Fn(tn)F_{X}(\overline{t})=F_{1}(t_{1})\cdot F_{2}(t_{2})\dots F_{n}(t_{n}) con

lim⁡t→∞Fi(ti)=1,\underset{ t\to \infty }{ \lim } F_{i}(t_{i})=1,

entonces son independientes.


Criterio para discretas

x1,…,xn son independientes  ⟺  Px1,x2,…,xn(x1…xn)=Px1(x1)…Pxn(xn)x_{1},\dots,x_{n} \text{ son independientes} \iff P_{x_{1},x_{2},\dots,x_{n}}(x_{1}\dots x_{n})=P_{x_{1}}(x_{1})\dots P_{x_{n}}(x_{n})

dem:dem:

para n=2n=2   ⟸  )\impliedby) X⊆RX,Y⊆RYX\subseteq R_{X},Y\subseteq R_{Y}

P({x∈A}∩{Y∈B})=∑X∈A,Y∈BPXY(X,Y)=P(\{ x \in A \}\cap \{ Y \in B \})=\sum_{X \in A,Y \in B} P_{XY}(X,Y)= =∑∑PX(X)⋅PY(Y)−∑X∈APX(A)∑Y∈BPY(Y)==\sum \sum P_{X}(X)\cdot P_{Y}(Y)-\sum_{X \in A}P_{X}(A)\sum_{Y \in B}P_{Y}(Y)= =(∑X∈APX(X))⋅(∑Y∈BPY(Y))=\left( \sum_{X \in A}P_{X}(X) \right)\cdot\left( \sum_{Y \in B}P_{Y }(Y) \right) =P(X∈A)⋅P(Y∈B)=P(X \in A)\cdot P(Y \in B)

Tiro una moneda 3 veces. Y=Y= cantidad de caras

X={1Cara la primera0ccX=\begin{cases} 1 & \text{Cara la primera} \\ 0 & cc \end{cases}

Ejemplo:

Y \ X 0 1
1 19\frac{1}{9} 29\frac{2}{9} 13\frac{1}{3}
2 19\frac{1}{9} 29\frac{2}{9} 13\frac{1}{3}
3 19\frac{1}{9} 29\frac{2}{9} 13\frac{1}{3}
13\frac{1}{3} 23\frac{2}{3}

Caso continuo

x1,…xn son independientes  ⟺  fx1,…,xn=fx1⋅fx2⋅⋯⋅fxnx_{1},\dots x_{n} \text{ son independientes}\iff f_{x_{1},\dots,x_{n}} = f_{x_{1}}\cdot f_{x_{2}}\cdot\dots\cdot f_{x_{n}}

Más aún, si

fX(X)=f1(x1)⋅f2(x2)…fn(xn)y∫−∞+∞f1(x) dx=1f_{X}(X)=f_{1}(x_{1})\cdot f_{2}(x_{2})\dots f_{n}(x_{n}) \quad y\quad \int_{-\infty}^{+\infty} f_{1}(x) \, dx = 1

Entonces fx1=fif_{x_{1}}=f_{i} y son independientes.

Ejemplo: (x,y)(x,y) con densidad

fXY(X,Y)=34⋅X2Y⋅1[−1,1](X)⋅1[0,2](Y)f_{XY}(X,Y)= \frac{3}{4}\cdot X^{2}Y\cdot \mathbb{1}_{[-1,1]}(X) \cdot\mathbb{1}_{[0,2]}(Y) =34⋅X2⋅1[−1,1](X)⋅Y⋅1[0,2](Y)=\frac{3}{4}\cdot X^{2}\cdot \mathbb{1}_{[-1,1]}(X) \cdot Y\cdot\mathbb{1}_{[0,2]}(Y) ∫02Y dY=2\int_{0}^{2} Y \, dY = 2 fY(Y)=12⋅Y⋅1[0,2](Y)yfX(X)=32⋅X⋅1[−1,1](X)f_{Y}(Y)= \frac{1}{2} \cdot Y \cdot\mathbb{1}_{[0,2]}(Y)\quad y\quad f_{X}(X)= \frac{3}{2} \cdot X \cdot\mathbb{1}_{[-1,1]}(X)

XX e YY son independientes.


Ejemplo 2:

fXY(X,Y)=e−Y1(0,Y)(X)f_{XY}(X,Y)= e^{-Y} \mathbb{1}_{(0,Y)}(X) fX(X)=e−X1(0,+∞)(X)f_{X}(X)= e^{-X} \mathbb{1}_{(0,+\infty)}(X) fY(Y)=Y⋅e−Y1(0,+∞)(Y)f_{Y}(Y)=Y \cdot e^{-Y} \mathbb{1}_{(0,+\infty)}(Y) fXY≠FX⋅FYf_{XY}\neq F_{X}\cdot F_{Y}

Ejemplo de vector aleatorio.

Uniforme, dado G⊆RnG\subseteq \mathbb{R}^{n}

X∼Usi fX(X)=1(G)(X)vol(G)X \sim \mathcal{U}\quad \text{si } f_{X}(X) = \frac{\mathbb{1}_{(G)(X)}}{vol(G)}

Donde vol(G)=∫Gdx1…dxnvol(G)=\displaystyle\int_{G} dx_{1}\dots dx_{n}


Ejemplo

D={(x,y)∈R2∣x2+y2≤1}D= \{ (x,y) \in \mathbb{R}^{2}\bigm| x^{2}+y^{2}\leq 1 \} X∼U(D)X \sim\mathcal{U}(D) fX(X)=1π⋅1D(X,Y)f_{X}(X)= \frac{1}{\pi}\cdot \mathbb{1}_{D}(X,Y)

XX e YY no son independientes.

soporte rectangular no significa independencia.


Si (X,Y)(X,Y) es un vector discreto y Z=g(X,Y)Z=g(X,Y)

PZ(k)=∑X∈RX,Y∈RY,g(X,Y)=kPXY(X,Y)P_{Z}(k)=\sum_{X \in R_{X},Y \in R_{Y},g(X,Y)=k} P_{XY}(X,Y)

Ejemplo:

X∼P(λ1),Y∼P(λ2),X,Y independientesX \sim P(\lambda_{1}),Y\sim P(\lambda_{2}),\quad X,Y \text{ independientes} Z=X+Y∼P(λ1+λ2)Z=X+Y\sim P(\lambda_{1}+ \lambda_{2}) PZ(k)=∑s∈N0,t∈N0,s+t=kPXY(s,t)=∑s=0kPX,Y(s,k−s)P_{Z}(k)=\sum_{s \in \mathbb{N}_{0},t \in \mathbb{N}_{0},s+t=k} P_{XY}(s,t)=\sum_{s=0}^{k} P_{X,Y}(s,k-s) =∑s=0kλ1ss!⋅e−λ1⋅λ2k−s(k−s)!⋅e−λ2=\sum_{s=0}^{k} \frac{\lambda_{1}^{s} }{s!} \cdot e^{-\lambda_{1}} \cdot \frac{\lambda_{2}^{k-s} }{(k-s)!} \cdot e^{- \lambda_{2}} =(λ1+λ2)kk!⋅e−(λ1+λ2)=\frac{(\lambda_{1}+\lambda_{2})^{k} }{k!} \cdot e^{-(\lambda_{1}+\lambda_{2})}

Ejemplo continuo

X,Y∼D[0,1]Z=XYX,Y \sim\mathcal{D}[0,1] \quad \quad Z = \frac{X}{Y}

X,YX,Y independientes.

P(Z≤t)=P(XY≤t)=P(X≤t⋅Y)P(Z \leq t)= P \left( \frac{X}{Y}\leq t \right) = P(X\leq t \cdot Y) fXY=1[0,1)×[0,1)(X,Y)f_{XY}= \mathbb{1}_{[0,1) \times[0,1)} ( X,Y)

Además Y≥XtY\geq \frac{X}{t}

Si t<1:t<1: Mirando el gráfico de la región vemos que

P(X≤t⋅Y)=t2P(X\leq t \cdot Y)=\frac{t}{2}

Teorema del cambio de variables

U,V∈RnU,V \in \mathbb{R}^{n} g:U⟶V,C1g:U\longrightarrow V, C^{1} con inversa C1C^{1}
A⊆VA\subseteq V

∫g−1(A)f(X) dX=∫Af(g−1(Y))⋅∣Jg−1(Y)∣ dY\int_{g ^{-1}(A)} f(X)\: dX = \int_{A} f(g ^{-1}(Y))\cdot |J_{g ^{-1}} ( Y)|\: dY

Si XX es un vector aleatorio y f=fXf=f_{X}

P(g(X)∈A)=P(X∈g−1(A))=∫AfX(g−1(X))⋅∣Jg−1(X)∣⋅1V(X) dYP(g(X) \in A) = P(X \in g ^{-1}(A)) = \int_{A} f_{X} (g ^{-1}(X)) \cdot |J_{g ^{-1}}(X)| \cdot \mathbb{1}_{V}(X) \: dY fY(Y)=fX(g−1(Y))⋅∣Jg−1(Y)∣⋅1V(Y)⋅1Im(Y)(Y)f_{Y}(Y) = f_{X}(g ^{-1}(Y))\cdot |J_{g ^{-1}}(Y)| \cdot \mathbb{1}_{V}(Y) \cdot \mathbb{1}_{\mathrm{Im}(Y)}(Y)

Lo puedo usar para X,Y∼U[0,1]X,Y \sim \mathcal{U}[0,1] independientes.

Z=XYZ = \frac{X}{Y}

Agrego W=YW=Y así voy de R2\mathbb{R}^{2} a R2\mathbb{R}^{2}

g(X,Y)=(XY,Y)g(X,Y)= \left( \frac{X}{Y},Y \right) g−1(Z,W)=(Z⋅W,W)g ^{-1}(Z,W) = (Z\cdot W,W)

∣Jg−1∣=[WZ01]=W|J_{g ^{-1}}| = \begin{bmatrix}W & Z \\ 0 & 1\end{bmatrix}=W

fZW(Z,W)=1[0,1](Z⋅W)⋅1[0,1](W)⋅Wf_{ZW}(Z,W)= \mathbb{1}_{[0,1]}(Z\cdot W)\cdot \mathbb{1}_{[0,1]}(W) \cdot W

Para rescatar la marginal de ZZ integro:

fZ(Z)=∫−∞+∞fZW(Z,W) dWf_{Z}(Z)=\int_{-\infty}^{+\infty} f_{ZW}(Z,W) \, dW
X,Y∼N(0,1) independientesX,Y \sim N(0,1) \text{ independientes} fX,Y(X,Y)=1(2π)2⋅e−x22⋅e−Y22=12π⋅e−x2−y22=12π⋅e−∣∣X∣∣22f_{X,Y}(X,Y)= \frac{1}{(\sqrt{ 2\pi })^{2}} \cdot e^{ - \frac{x^{2}}{2}} \cdot e^{- \frac{Y^{2}}{2}} = \frac{1}{2\pi} \cdot e^{\frac{-x^{2}-y^{2}}{2}} = \frac{1}{2\pi}\cdot e^{\frac{- || X || ^{2}}{2}}

Obs:Obs: Si A⊆R2×2A \subseteq \mathbb{R}^{2\times2} es ortogonal, A⋅AT=IA\cdot A^{T}=I

Z=A(XY)Z= A\begin{pmatrix} X \\ Y \end{pmatrix} fZ(Z)=12π⋅e−∣∣A−1Z∣∣22det⁡(A−1)=12π⋅12π⋅e−Z122⋅e−Z222f_{Z}(Z)= \frac{1}{2\pi}\cdot e ^{- \frac{||A ^{-1} Z || ^{2}}{2}} \det(A ^{-1})= \frac{1}{\sqrt{ 2\pi }} \cdot \frac{1}{\sqrt{ 2\pi }} \cdot e^{- \frac{Z_{1}^{2}}{2}} \cdot e^{- \frac{Z_{2}^{2}}{2}} Z1,Z2∼N(0,1)independientes.Z_{1},Z_{2} \sim N(0,1) \quad \text{independientes.}

Si ∣∣(a,b)∣∣=1|| (a,b) ||=1 entonces a⋅X+b⋅Y∼N(0,1)a\cdot X+b \cdot Y\sim N(0,1) tomo

A=(abb−a)A=\begin{pmatrix} a & b \\ b & -a \end{pmatrix} Z1=a⋅X+b⋅YZ_{1}= a\cdot X+b\cdot Y

Si a,b∈Ra,b \in \mathbb{R} a⋅X+b⋅Ya \cdot X + b \cdot Y es normal.

Vale que si X∼N(μX,σX2)X \sim N(\mu_X, \sigma_{X}^{2}), Y∼N(μY,σY2)Y \sim N(\mu_{Y}, \sigma_{Y}^{2}) . Independientes.

a⋅X+b⋅Y es normal.a \cdot X + b \cdot Y \text{ es normal.}

Escuchar audio, importante. 00:38:16 parte 2.


Distribución Γ\Gamma

X∼Γ(α,λ)X \sim \Gamma(\alpha,\lambda) fX(X)=λαΓ(α)⋅Xα−1⋅e−λX⋅1(0,+∞)(X)f_{X}(X)=\frac{\lambda ^{\alpha} }{\Gamma(\alpha)} \cdot X^{\alpha-1} \cdot e^{-\lambda X} \cdot \mathbb{1}_{(0,+\infty)}(X) Γ(Y)=∫0+∞XY−1⋅e−Y dX\Gamma(Y)=\int_{0}^{+\infty} X ^{Y-1} \cdot e^{-Y} \, dX Γ(Y+1)=Y⋅P(Y)Γ(n)=(n−1)!Γ(12)=π\begin{array}{c} \Gamma(Y+1)= Y \cdot P(Y) \\ \Gamma(n)=(n-1)! \\ \Gamma\left( \frac{1}{2} \right)=\sqrt{ \pi }\end{array}
X∼Γ(1,λ)X \sim\Gamma(1,\lambda) X∼E(λ)X∼E(λ)Y∼E(λ)\begin{array}{c} X \sim\mathcal{E}(\lambda) \\ X \sim\mathcal{E}(\lambda) \\ Y \sim\mathcal{E}(\lambda) \end{array}

Entonces X+Y∼Γ(2,λ)X+Y \sim \Gamma(2,\lambda) , independientes XX+Y∼U[0,1]\frac{X}{X+Y} \sim \mathcal{U}[0,1]

Temas relacionados