9. Simulacion estadistica del Ejercicio 2.16(Teorema de Bayes)

Simulacion con R-Project del Teorema de Bayes.

Simulacion:

> poblac <- c("Maquina1","Maquina2","Maquina3");

> piezas<- c(200,100,50);

> ppiezas <- piezas/sum(piezas)

> DefectuosaSegunMaquina <- c(0.3,0.2,0.1)

> #Teorema de Bayes

> #Probabilidad maquina siendo pieza defectuosa 

> MaquinaSiendoDefectuosa<- (ppiezas*DefectuosaSegunMaquina)/sum(ppiezas
*DefectuosaSegunMaquina)

> A<-rbind(MaquinaSiendoDefectuosa)

> colnames(A)<-poblac

> A 
                         Maquina1  Maquina2   Maquina3
MaquinaSiendoDefectuosa 0.7058824 0.2352941 0.05882353