[R-br] Otimizar tempo de processamento

[Benilton Carvalho]

Tiago,

isso e' resultado de y sendo igual a NA ou integer(0).

Vc vai precisar trocar o %dopar% por %do% para debugar. Uma sugestao
e' comecar por adiciionar umas condicoes para teste dentro do for()

for (k in 1:n) y <- rbinom(1, y*f, p)

usar algo como

for (k in 1:n){
  y <- rbinom(1, y*f, p)
  if (length(y) == 0){
     browser()
  }else{
     if (!is.finite(y))
        browser()
  }
}

b