2. Simulacion estadistica del Ejercicio 3.7 (distribuciones marginales y condicionadas)

Simulacion:

> y<- c(1,2,3,4)

> x<- c(-1,0,1)

> xy<- expand.grid(X=x,Y=y)

> dpc<- c(0.03,0.01,0,0.02,0.25,0.13,0.05,0.45,0.05,0,0.01,0)

> D <- cbind(xy, dpconjunta = dpc)


> #tabla de distribuciones de probabilidad conjunta de las variables X e Y

> XD <- xtabs(dpconjunta ~ X + Y, data = D)

> XD
    Y
X       1    2    3    4
  -1 0.03 0.02 0.05 0.00
  0  0.01 0.25 0.45 0.01
  1  0.00 0.13 0.05 0.00
  
> #Distribuciones marginales de X

>dmX<-rowSums(XD)

>dmX
  -1    0    1 
0.10 0.72 0.18 

> #Distribuciones marginales de Y

> dmY<-colSums(XD)


> dmY
   1    2    3    4 
0.04 0.40 0.55 0.01 

> #Distribuciones de Y condicionada por X=0, f(y/x=0)

> XD[2,]/dmX[2]
         1          2          3          4 
0.01388889 0.34722222 0.62500000 0.01388889