Boa Tarde pessoal !

Peguei do manual do Bioestat o exemplo da Regra de Bayes para fazer no R e consegui na mão grande, logo abaixo. Tentei então colocar em loop para expandir para além de 3 situações e aí deu erro no cálculo do denominador (o loop está no final). Como posso resolver ?
Obrigado,
Mauricio Cardeal
UFBA

# Regra de Bayes
# Dados populacionais
PB1 <- 0.15 # asma
PB2 <- 0.01 # enfisema
PB3 <- 0.009 # ca

# Dados clínicos
PAdadoB1 <- 0.26 # dispneia dado asma
PAdadoB2 <- 0.75 # dispneia dado enfisema
PAdadoB3 <- 0.43 # dispneia dado ca

PB1dadoA <- (PAdadoB1 * PB1) / (PAdadoB1*PB1 + PAdadoB2*PB2 + PAdadoB3*PB3)
PB2dadoA <- (PAdadoB2 * PB2) / (PAdadoB1*PB1 + PAdadoB2*PB2 + PAdadoB3*PB3)
PB3dadoA <- (PAdadoB3 * PB3) / (PAdadoB1*PB1 + PAdadoB2*PB2 + PAdadoB3*PB3)

rbind(Asma_dado_dispneia=PB1dadoA,Asma_dado_enfisema=PB2dadoA,Asma_dado_ca=PB3dadoA)

                       [,1]
Asma_dado_dispneia 0,774270
Asma_dado_enfisema 0,148898
Asma_dado_ca       0,076831


# LOOP

PB <- c(0.15,0.01,0.009)
PAdadoB <- c(0.26,0.75,0.43)
dados <- data.frame(PB,PAdadoB)
n = dim(dados)[1]
dados$PBdadoA <- 0
dados$den <- 0
dados
for(i in 1:n)
{
dados$num <- dados$PAdadoB[i] * dados$PB[i]
    for(j in 1:n)
      {
        dados$den <- dados$PAdadoB[j]*dados$PB[j] + dados$den
      }
dados$PBdadoA <- dados$num[i] / dados$den[i]
}

dados

     PB PAdadoB  PBdadoA     den     num
1 0,150    0,26 0,025610 0,15111 0,00387
2 0,010    0,75 0,025610 0,15111 0,00387
3 0,009    0,43 0,025610 0,15111 0,00387