TEOREMA-DE-BAYES.R

# UNIVERSIDAD NACIONAL DEL ALTIPLANO
# INGENIERIA ESTADISTICA E INFORMATICA
# CURSO: ESTADISTICA BAYESIANA
# TEMA: TEOREMA DE BAYES

# Forma 1
print("Ingrese el Porcentaje")

## [1] "Ingrese el Porcentaje"

porcentaje <- c(scan())
print("Ingrese la Probabilidad ")

## [1] "Ingrese la Probabilidad "

pro <- c(scan())

TeoremaB <- function(porcentaje,pro){
  suma <- 0.00
  for (i in 1:length(porcentaje)) {
    suma <- suma+porcentaje[i]*pro[i]
  }
  resta <- 0.00
  probabilidad <- c()
  for (i in 1:length(porcentaje)){
    resta <- porcentaje[i]*pro[i]/suma
    probabilidad <- c(probabilidad,resta)
  }
  return(probabilidad)
}
print(TeoremaB(porcentaje, pro))

## numeric(0)

# Forma 2

#Ingreso de las Probabilidades primarias
pro1 <- c(0.6)
pro2 <- c(0.4)
#Ingreso de las Probabilidades secundarias 
pru1 <- c(0.2)
pru2 <- c(0.8)
#Probabilidad de ser culpable
Probabilidad1 <- pru1*pro1/sum(pru1*pro1,pru2*pro2)
#Probabilidad de no ser culpable
Probabilidad2 <- pru2*pro2/sum(pru1*pro1,pru2*pro2)

A <- rbind(Probabilidad1,Probabilidad2)
print(A)

##                    [,1]
## Probabilidad1 0.2727273
## Probabilidad2 0.7272727