similar to: (Parallel) Random number seed question...

Displaying 20 results from an estimated 2000 matches similar to: "(Parallel) Random number seed question..."

2010 Feb 18
2
sprng 2.0 header file for Rmpi
When I try to install Rmpi, it has a dependency rsprng which, as its description says: Provides interface to SPRNG 2.0 APIs I installed it on Mepis by installing the appropriate debs with minimal drama. In the rpm world of Fedora, no such luxury. Installing rsprng fails because: checking sprng.h usability... no checking sprng.h presence... no checking for sprng.h... no Cannot find sprng 2.0
2007 May 18
1
Cannot install rsprng package (PR#9697)
<<insert bug report here>> I cannot install the R package "rsprng" with the following error message: > install.packages("rsprng") trying URL 'http://www.stathy.com/cran/src/contrib/rsprng_0.3-4.tar.gz' Content type 'application/x-tar' length 35934 bytes opened URL ================================================== downloaded 35Kb * Installing
2011 Mar 14
1
Installing Rmpi on hpc
Hi, I was trying to install the package Rmpi on a hpc cluster running SGE. The command, and the sessionInfo() is as follows: =========================================== > install.packages("Rmpi",dependencies=TRUE) also installing the dependency ‘rsprng’ trying URL ' http://www.ibiblio.org/pub/languages/R/CRAN/src/contrib/rsprng_1.0.tar.gz' Content type
2002 Dec 17
1
rsprng doesn't install on Debian (woody) (PR#2378)
Full_Name: Emmanuel Charpentier Version: 1.6.1 OS: Linux 2.4.20 / Debian Woody Submission from: (NULL) (80.15.78.96) rsprng depends on libraries available in the "genesis" package. However, even when Genesis is installed, the installation of the source package fails for not finding rsprng.h. It turns out that Debian installs it in a subdirectoru of /usr/include, and that the rsprng
2003 Dec 01
2
help with random numbers and Rmpi
Dear People, This may not be the right place to ask a question about Rmpi, but I don't know of a better one. I am trying to get a simple program working using Rmpi with the model of 1 R master and n C slaves. What I am trying to do is have each of the C slaves generate a random number from U[0,1], and then have the master collect all n numbers as a vector and output it. However even doing
2007 Jun 30
1
random numbers
Although this query was inspired by distributed random number generation, one of the questions (#2 below) is a single-machine issue. I call C++ code from R to generate simulated data. I'm doing this on a cluster, and use rmpi and rsprng. While rsprng randomizes R-level random numbers (e.g., from runif), it has no effect on the C code, which is completely SPRNG and MPI ignorant. Currently I
2009 Jul 30
3
user supplied random number generators
?Random.user says (in svn trunk) Optionally, functions \code{user_unif_nseed} and \code{user_unif_seedloc} can be supplied which are called with no arguments and should return pointers to the number of seeds and to an integer array of seeds. Calls to \code{GetRNGstate} and \code{PutRNGstate} will then copy this array to and from \code{.Random.seed}. And it offers as an example void
2011 Dec 12
1
(no subject)
Hi all, I'm trying to get the snow package with mpi working on my linux ubuntu 10.04 machine. When I tried to install it (see below) it already had problems with the dependencies. R version 2.14.0 (2011-10-31) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You
2005 Jun 08
6
Random seed problem in MCMC coupling of chains
Hello! I am performing coupling of chains in MCMC and I need the same value of seed for two chains. I will show demo of what I want: R code, which might show my example is: niter <- 3 nchain <- 2 tmpSeed <- 123 for (i in 1:niter) { # iterations for (j in 1:nchain) { # chains set.seed(tmpSeed) a <- runif(1) cat("iter:", i, "chain:", j,
2009 Aug 24
1
R with MPI
Hello, I plan to use R with my cluster with OpenMPI. I need the packaged 'snow' and 'Rmpi' for that, however, I get an error while downloading and installing them: When I do a: install.packages("Rmpi", dependencies=T) I get this error: checking for mpi.h... no Try to find libmpi.so or libmpich.a checking for main in -lmpi... no libmpi not found.
2011 Jan 07
1
Trouble with installing Rmpi package
Hi, I am having a problem with installing Rmpi package on redhat linux machine. The R I am using is version 2.10.1. Here’s what happens. > install.packages( ''Rmpi'' ) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL ''http://cran.cnr.Berkeley.edu/src/contrib/Rmpi_0.5-9.tar.gz'' Content type
2007 Mar 21
1
Snow Package and SPRNG: Will it solve my problem?
Hello and thanks in advance for your time. I currently have a simulation running on my cluster with the help of snow that relies on global variables being changes regularly to random values. It uses these values, lets call them x1 x2 and x3, in custom functions for logliklyhood and score that gets used in the standard optim function. To get set in the global table on the different
2009 Feb 13
2
Generate random numbers in Fortran
Hi!!! It would like to know if it exists a form to use the functions to generate variates in FORTRAN with the same easiness I use that them in C? Or not? If yes. They would have some example? I would like to use the functions rbeta, rlnorm and others! Sorry my english..rsrsrs Thanks!!!              Fábio Mathias Corrêa    University Federal of the Lavras - Brazil Veja quais são os
2004 Aug 02
2
random seed puzzle
After reading the help page on set.seed, I am unsure about how to solve the following problem. I need to call function f a thousand times. The list of values returned by f, should be as random as possible. f calls g twice: f <- function(){g1 <- g(1); g2 <- g(2); c(g1; g2)} The function g in turn calls sample and returns a number, but also depends on its argument, so, starting from the
2001 Jul 14
1
Calling rand_unif() from fortran
What is the necesary incantation to call rand_unif() from a fortran routine? Just inserting rand_unif() into the fortran code compiles OK, but won't dyn.load(): > dyn.load("forest.so") Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "/home/warnes/R/forest/forest.so": /home/warnes/R/forest/forest.so: undefined symbol:
2010 Oct 24
1
140 packages in R Commander!!
Dear List I just downloaded and installed R 2.12.0 and then installed R Commander . First it got RCmdr and Car, and then suggested for other packages for utilizing the full functionality- I clicked yes! I got 140 packages installed!!! Cran Mirror was UCLA... Here is the list. Is this intentional- I can see some packages like snow and multicore which are desirable but quite optional.(see list
2005 Nov 30
1
RNG stuck via Fortran call
Having not much success with my previous question I try to reformulate it: I'm simulating a Markow chain in Fortran interfaced with R. Each loop of my Fortran calls various functions of the R RNG through the wrapper given below. In a run of 100 iterations of the Markov kernel, after 20 iterations, the RNG seems like frozen. For example, the first call to the RNG in my loop is:
2005 Jun 14
3
Calling C from Fortran
I would like to call C routines from Fortran under linux as suggested in section 5.6 of the "Writing R extensions" documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd() call dblepr("X was", 5, x, 1) call rndend() end but I don't understand
2005 Nov 30
2
R binomial RNG stuck at 1 via Fortran call
Hi, I have some trouble with the result of a fortran function calling the R binomial RNG: The C function rbinom is wrapped as in the file attached. My main fortran program starts as call rndstart() and ends as call rndend() I happen to call the binomial RNG within a loop as b = ggrbinom(1.d0,0.5d0) write(*,*) 'b=',b In certain cases, after a few iterations in the loop, b get
2002 May 08
3
Suggestions for poor man's parallel processing
Almost all of the heavy crunching I do in R is like: > for(i in long.list){ + do.something(i) + } > collect.results() Since all the invocations of do.something are independent of one another, there is no reason that I can't run them in parallel. Since my machine has four processors, a natural way to do this is to divide up long.list into 4 pieces and then start 4 jobs, each of which