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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
