3. Simulacion estadistica del Ejercicio 2.7 (Probabilidad cartas enviadas)

Tendremos que utilizar el package prob y combinat del R-Project

Simulacion:

> library(prob)

> library(combinat)

> L <- c("C1","C2","C3");

> M <- urnsamples(L, size=3, replace = FALSE, ordered = TRUE);

> N <- probspace(M);

> N
  X1 X2 X3     probs
1 C1 C2 C3 0.1666667
2 C1 C3 C2 0.1666667
3 C3 C1 C2 0.1666667
4 C3 C2 C1 0.1666667
5 C2 C3 C1 0.1666667
6 C2 C1 C3 0.1666667

> #La carta 1 llegue a su destinatario correspondiente.

> A <- subset(N,X1 == "C1")

> #La carta 2 llegue a su destinatario correspondiente.

> B <- subset(N,X2 == "C2");

> #La carta 3 llegue a su destinatario correspondiente.

> C <- subset(N,X3 == "C3");

> #Probabilidad de que las cartas lleguen correctamente a sus destinatarios.

> prob(union(A,union(B,C)))
[1] 0.6666667