
Uma pequena alteração para nomear as séries... ### <code r> res <- NULL for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) { id <- paste(n, p, sep="_") sim <- list(rbinom(n, 1,p)); names(sim) <- id res <- c(res, sim) } res ### exemplo de operação com lista - cálculo das proporções sapply(res, function(x) sum(x)/length(x)) ## 10_0.3 10_0.4 10_0.5 20_0.3 20_0.4 20_0.5 30_0.3 30_0.4 30_0.5 ## 0.3000000 0.3000000 0.5000000 0.3000000 0.3000000 0.4000000 0.2333333 0.4000000 0.6000000 ### </code> ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] Em 28 de janeiro de 2016 08:03, Éder Comunello <comunello.eder@gmail.com> escreveu:
Bom dia, Paulo!
Uma forma simples é acumular os resultados em uma lista de vetores. Posteriormente você pode operar a lista para cálculos de interesse.
### <code r> res <- NULL for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) { sim <- rbinom(n, 1,p) res <- c(res, list(sim)) }
res
### exemplo de operação com lista - cálculo das proporções sapply(res, function(x) sum(x)/length(x))
### </code>
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]