similar to: Random numbers

Hi R-Developers, I had a question about the random number generator used in the R StandAlone Math Library. The stand-alone library depends on the unif_rand() function for most simulated values, and this function is provided in the sunif.c file in the relevant directory. At present, this program implements the "Marsaglia-Multicarry" algorithm, which is described throughout the R

Dear all, Recently I am generating large random samples (10M) and any duplicated numbers are not desired. We tried several RNGs in R and found Wichmann-Hill did not produce duplications. The duplication problem is the interesting birthday problem. If there are M possible numbers, randomly draw N numbers from them, the average number of dupilcations D = N(N-1)/2/M. For Knuth-TAOCP and

Dear All, I want to use R's random number in a C++ program (I can link libRmath either as shared or static library). I have two questions: 1. If I understand correctly, the underlaying random number generator will be Marsaglia-multicarry, UNLESS I provide my own. In other words (unless I provide it) I cannot use some of the other RNG's available from within R, such as Mersenne-Twister?

Dear Sir/Madam, I am a beginner in R. Here is my questions. 1. Can you give me one test for randomness (a name and descriptive paragraph is sufficient). 2. I have learned a uniform random number generator [e.g. not the algorithms: i)Wichmann-Hill, ii) Marsaglia-Multicarry, iii) Super-Duper (Marsaglia), iv) Mersenne-Twister, v) TAOCP-1997 (Knuth), or vi) TAOCP-2002 (Knuth)] . Is there any other

I have commited fixes for 0.64.1 for (1) From: Mai Zhou <mai@ms.uky.edu> > .Random.seed <- c(1, 1:2) > .Random.seed [1] 1 1 2 > runif(5) Warning: Wrong length .Random.seed; forgot initial RNGkind? set to Wichmann-Hill[1] 0.02253721 0.84832584 ........ Here the length of the seed was being tested before the kind was picked out, so the length of the previous type was used. (2)

I sent this to R-devel early last month, but have received no response, so I guess it really is a bug. This looks like a bug to me, and is a bit hard to describe, but easy to reproduce. ? Basically, if RNGkind is saved as something other than the default, and if the first operation in a session is a set.seed(), the default is reverted to. ?Reproduce by: cafe-rozo> ?R --vanilla R version

This looks like a bug to me, and is a bit hard to describe, but easy to reproduce. Basically, if RNGkind is saved as something other than the default, and if the first operation in a session is a set.seed(), the default is reverted to. Reproduce by: cafe-rozo> R --vanilla R version 2.7.1 (2008-06-23) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. ---559023410-851401618-1218751024=:15885 Content-Type: TEXT/PLAIN; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: QUOTED-PRINTABLE I didn't describe the problem clearly. It's about the number of distinct=20 values. So just

Dear All, With a colleague we are writing a paper where we show how R is a very nice tool to deal with some issues in the analyses of data in evolutionary biology. For the intro, I wonder if 1) Anybody has any rough idea of how many people might be using R or how many people have downloaded R, or similar (I am aware answering this question might require divinatory powers...). 2) Have/are any

I have two questions about using a user-defined random number generator (RNG) with the standalone Rmath library. The default RNG with the standalone Rmath library is the Marsaglia-multicarry generator, which has poor properties. The "R Installation and Administration" manual, in the section "The standalone Rmath library", states that: ``` A little care is needed to use the

Dear All, In the title of a plot, I would like to mix greek letters with numbers, where the numbers are obtained from a particular function to a vector (e.g., max(x)); in each call, the value of this vector can change. Without greek symbols I use something like: title(sub=paste("x1=", deparse(x[1]),"beta = ",deparse(max(x)), "rho = ", deparse(min(x)))) but I'd

Dear All, When generating pdf documentation for a package (using R CMD Rd2dvi --pdf), is there a way to get the date to use the typical international standard of day month year instead of the US one of month day, year? Thanks, Ramon -- Ramón Díaz-Uriarte Dept. Zoology and Statistics University of Wisconsin-Madison Madison, WI 53706-1381 email: rdiazuri at students.wisc.edu (NOTE: starting

Forwarded message: >From BMCCULLO at fcc.gov Fri Nov 6 12:50:07 1998 Message-Id: <s642efa3.006 at fcc.gov> X-Mailer: Novell GroupWise 4.1 Date: Fri, 06 Nov 1998 12:51:24 -0500 From: Bruce McCullough <BMCCULLO at fcc.gov> To: mai at ms.uky.edu Subject: S-Plus RNG Mime-Version: 1.0 Content-Type: text/plain Content-Disposition: inline Mai, Someone forwarded to me your e-mail about

Is it appropriate to use stepAIC (library MASS) with coxph objects (from library survival5) that use "cluster(id)"? It is my understanding that, when using "cluster(id)", we can test for sets of terms by using the methods in Wei et al., (1989; JASA, 84: 1065-1073), or as explained in pp. 53 and ff. of the survival.ps document. But if we use a likelihood ratio test instead

Dear All, The help for prcomp, under "Value" says: sdev: the standard deviation of the principal components (i.e., the eigenvalues of the cov matrix, though the calculation is actually done with the singular values of the data matrix). The way I read it, it implies that the sdev are the eigenvalues, but I think that sdev is actually the square root of the

I am writing a function for the purposes of a simulation. Due to memory problems, the function sometimes crashes. In order to get around this problem, I would like to include to be able to save the "last" seed, so I can pick up with the next run of the simulation after a "crash". I am having trouble understanding what is going on with .Random.seed! For each run of the

I'm running into problems with set.seed--maybe I'm misunderstanding something. I'm running R 1.5.1 on Windows 2000. I'm basically trying to capture the random seed so that I can reproduce a simulation if it's necessary later. Using set.seed, I can certainly get reproducible results, but not the results I get on the first pass. Here's an example: # Generate a random

On Wed, 30-12-2015, at 12:44, Dirk Eddelbuettel <edd at debian.org> wrote: > On 30 December 2015 at 05:00, Dario Strbenac wrote: > | Good day, > | > | In section A.2, the manual advises "Linux users will want the urw-fonts package". However, this package only seems to be available for RedHat Linux and Fedora Linux. What about for Debian or Ubuntu ? There is no

Dear People, I have a vexing problem related to pseudo-random number generation, and would appreciate any help and advice. This problem is not directly related to R, and the only reason I am posting it to this list is that my implementation is using R. Let me describe my problem by giving an example, that is close to what I am trying to do. Suppose we are given a stream of pseudo-random numbers,

I noticed this error in my demo from previous versions as well as R 0.99.0. Is there a way around this one also? Thanks in advance... > demo(dyn.load) demo(dyn.load) ---- ~~~~~~~~ Type <Return> to start : > dyn.load(file.path(R.home(), "demos", "dynload", paste("zero", .Platform$dynlib.ext, sep = ""))) Error in