* Para utilizar la fórmula, se requiere disponer de la probabilidad del evento A, y la probabilidad condicional del evento B dado el evento A
## Probabilidad de que una pieza se fabricada por el proveedor 1 = 0.60
PA1 <- 0.60
## Probabilidad de que una pieza se fabricada por el proveedor 2 = 0.40
PA2 <- 0.40
## Probabilidad de que una pieza sea buena dada que sea del proveedor 1 es: 0.97
PG.PA1 <- 0.97
## Probabilidad de que una pieza sea mala dada que sea del proveedor 1 es: 0.03
PB.PA1 <- 0.03
## Probabilidad de que una pieza sea buena dada que sea del proveedor 2 es: 0.96
PG.PA2 <- 0.96
## Probabilidad de que una pieza sea mala dada que sea del proveedor 2 es: 0.04
PB.PA2 <- 0.04
##### Cálculo de probabilides condicionales conforme al árbol de decisión para proveedor 1 * Cada una de las variables siguientes se determina en el árbol multiplicando las probabildiades de cada hoja. * Para hallar la probabilidad de cada uno de los resultados experimentales, simplemente se multiplican las probabilidades de las ramas que llevan a ese resultado * Se hace notar que la suma de cada probabilidad en cada paso es 1 o el 100% * Probabilidad de que una pieza sea buena dado el proveedor 1. PA1.I.G = PA1 * PG.PA1 * Probabilidad de que una pieza sea mala dado el proveedor 1. PA1.I.B = PA1 * PB.PA1
PA1.I.G <- PA1 * PG.PA1
PA1.I.B <- PA1 * PB.PA1
cat("La probabilidad de que sea del proveedor 1 y que la pieza sea buena es: ",PA1.I.G)
## La probabilidad de que sea del proveedor 1 y que la pieza sea buena es: 0.582
cat("La probabilidad de que sea del proveedor 1 y que la pieza sea mala es: ",PA1.I.B)
## La probabilidad de que sea del proveedor 1 y que la pieza sea mala es: 0.018
PA2.I.G <- PA2 * PG.PA2
PA2.I.B <- PA2 * PB.PA2
cat("La probabilidad de que sea del proveedor 1 y que la pieza sea buena es: ",PA2.I.G)
## La probabilidad de que sea del proveedor 1 y que la pieza sea buena es: 0.384
cat("La probabilidad de que sea del proveedor 1 y que la pieza sea mala es: ",PA2.I.B)
## La probabilidad de que sea del proveedor 1 y que la pieza sea mala es: 0.016
TB.PA1.B <- (PA1 * PB.PA1) / (PA1 * PB.PA1 + PA2 * PB.PA2)
cat("Conforme al teorema de Bayes (TB), la probabilidad de que sea una pieza mala (Bad) condicionada a que sea primero del proveedor1 es: ", TB.PA1.B)
## Conforme al teorema de Bayes (TB), la probabilidad de que sea una pieza mala (Bad) condicionada a que sea primero del proveedor1 es: 0.5294118
TB.PA2.B <- (PA2 * PB.PA2) / (PA1 * PB.PA1 + PA2 * PB.PA2)
cat("Conforme al teorema de Bayes (TB), la probabilidad de que sea una pieza mala (Bad) condicionada a que sea primero del proveedor2 es: ", TB.PA2.B)
## Conforme al teorema de Bayes (TB), la probabilidad de que sea una pieza mala (Bad) condicionada a que sea primero del proveedor2 es: 0.4705882
* Solo se muestran los valores para las piezas malas (bad) de cada proveedor
tabular <- data.frame('Eventos'=c('A1', 'A2'),
'Prob.Previas'=c(PA1, PA2),
'Prob.Condicionales'=c(PB.PA1, PB.PA2),
'Prob.Conjuntas'=c(PA1.I.B, PA2.I.B),
'Prob.Posteriores'=c(TB.PA1.B, TB.PA2.B))
tabular
## Eventos Prob.Previas Prob.Condicionales Prob.Conjuntas Prob.Posteriores
## 1 A1 0.6 0.03 0.018 0.5294118
## 2 A2 0.4 0.04 0.016 0.4705882
totales <- apply(tabular[-1], 2, sum)
totales <- as.array(c(NA,as.vector(totales)))
tabular <- rbind(tabular, totales)
tabular
## Eventos Prob.Previas Prob.Condicionales Prob.Conjuntas Prob.Posteriores
## 1 A1 0.6 0.03 0.018 0.5294118
## 2 A2 0.4 0.04 0.016 0.4705882
## 3 <NA> 1.0 0.07 0.034 1.0000000
cat("1.- ¿Cuál es la probabilidad de que sea del proveedor 1?", TB.PA1.B)
## 1.- ¿Cuál es la probabilidad de que sea del proveedor 1? 0.5294118
cat("2.- ¿Cuál es la probabilidad de que sea del proveedor 2?", TB.PA2.B)
## 2.- ¿Cuál es la probabilidad de que sea del proveedor 2? 0.4705882