have you tried
res <- rnorm(10^6)-rnorm(10^6)
?
it take 0.5 sec for me...
b
On Feb 28, 2008, at 3:30 PM, Minimax wrote:
Dear useRs,
Suppose we have loop:
res <- c()
for (i in 1:100000) {
x <- rnorm(2)
res <- c(res,x[2]-x[1])
}
and this loop for 10^5 cases runs about - for example 5 minutes.
When I add one zero (10^6) this loop will not end overnight but
probably
hangs. This occurs regardless of calculated statistics in such
simulation, always above 10^5 times. Nested loops do not help.
Any suggestions for collecting larger amount of Monte Carlo data ?
Regards
Minimax
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
