<style type="text/css">p { margin-bottom: 0.21cm; }</style>

<p class="western" style="margin-bottom: 0cm;">Boa Tarde pessoal
!<br><br>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 ? <br>Obrigado,<br>Mauricio Cardeal<br>UFBA <br><br># Regra
de Bayes<br># Dados populacionais<br>PB1 <- 0.15 # asma<br>PB2 <-
0.01 # enfisema<br>PB3 <- 0.009 # ca<br><br># Dados
clínicos<br>PAdadoB1 <- 0.26 # dispneia dado asma <br>PAdadoB2 <-
0.75 # dispneia dado enfisema<br>PAdadoB3 <- 0.43 # dispneia dado
ca<br><br>PB1dadoA <- (PAdadoB1 * PB1) / (PAdadoB1*PB1 +
PAdadoB2*PB2 + PAdadoB3*PB3)<br>PB2dadoA <- (PAdadoB2 * PB2) /
(PAdadoB1*PB1 + PAdadoB2*PB2 + PAdadoB3*PB3)<br>PB3dadoA <-
(PAdadoB3 * PB3) / (PAdadoB1*PB1 + PAdadoB2*PB2 +
PAdadoB3*PB3)<br><br>rbind(Asma_dado_dispneia=PB1dadoA,Asma_dado_enfisema=PB2dadoA,Asma_dado_ca=PB3dadoA)<br><br>                      
[,1]<br>Asma_dado_dispneia 0,774270<br>Asma_dado_enfisema
0,148898<br>Asma_dado_ca      
0,076831<br><br><br># LOOP<br><br>PB <- c(0.15,0.01,0.009)<br>PAdadoB
<- c(0.26,0.75,0.43)<br>dados <- data.frame(PB,PAdadoB)<br>n =
dim(dados)[1]<br>dados$PBdadoA <- 0<br>dados$den <-
0<br>dados<br>for(i in 1:n)<br>{<br>dados$num <- dados$PAdadoB[i]
* dados$PB[i]<br>    for(j in 1:n)<br>     
{<br>        dados$den <-
dados$PAdadoB[j]*dados$PB[j] + dados$den<br>     
}<br>dados$PBdadoA <- dados$num[i] / dados$den[i]<br>}<br><br>dados<br><br>    
PB PAdadoB  PBdadoA     den    
num<br>1 0,150    0,26 0,025610 0,15111 0,00387<br>2
0,010    0,75 0,025610 0,15111 0,00387<br>3 0,009   
0,43 0,025610 0,15111 0,00387<br><br>
</p>