
25 Ago
2011
25 Ago
'11
09:08
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