similar to: RNG.c: unif_rand, MARSAGLIA_MULTICARRY (PR#2437)

I would like to suggest a modification of initial scrambling of the seed in RNG_Init (called from do_setseed). The modified code is equivalent, but faster. Patch against R-devel_2007-09-22 follows --- R-devel-orig/src/main/RNG.c 2007-09-02 07:49:35.000000000 +0200 +++ R-devel-modif/src/main/RNG.c 2007-09-23 10:51:59.234566440 +0200 @@ -216,8 +216,8 @@ BM_norm_keep = 0.0; /* zap Box-Muller

Just out of curiosity, here is a question that might be embarrasingly stupid. Hmm ... What I need is to call runif() (more precisely, `runif(1, 0, 1)') from inside a C function. Is this possible using call_S? If so, how? -k =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

hello, everybody, I met a problem that I want to generate a random uniform number by using function with c interface. What I found is only "unif_rand()", but its range is [0,1]. How can I create a uniform random variable by using a function with c interface. Who can tell me what function it is and how to use? I will appreciate for it very much. Thanks in advance! xiaotong wang

While I untar a large archive on xfs , ext3 (ver 1.3 and ver 1.4) file systems , on ppc processor and kernel ver 2.6.21 , I get an error. Also sometimes, on ext3 (1.3 and 1.4) the file system goes read-only while untarring. The same tar file when untarred on a i386 machine works properly. ERROR: -------------- tar: Skipping to next header gzip: stdin: invalid compressed data--crc error tar:

Dear R-users: Recently I found my simulation run into an apparently infinite loop. After a few days of tracing and chasing, I believe it is caused by the built-in unif_rand() and exp_rand() functions: unif_rand() can produce a value of 0 which causes the following part of exp_rand() running into an infinity loop u = unif_rand(); for (;;) { u += u; if (u > 1.0) break;

https://bugzilla.mindrot.org/show_bug.cgi?id=2437 Bug ID: 2437 Summary: ssh with ControlMaster and ControlPath hangs on 2nd session in same terminal Product: Portable OpenSSH Version: 6.7p1 Hardware: amd64 OS: Linux Status: NEW Severity: normal Priority: P5

>I am having a little problem with executing logon scripts on Win95 >clients, >logging into Samba Netlogon service. >When a user logges in, a MS-DOS prompt shows up and quicky disappears. >I have done the logon script in MS-DOS editor. I have the identical copy >of >the batch file on the Win95 machine. If, I run the logon.bat script in >ms-dos prompt manually, it all

Hi: I was wondering why Marsaglia's new ziggurat method for generating deviates from the standard normal distribution has not been implemented in the R base package. I know that it is available in SuppDists pacakage of Bob Wheeler, as "rziggurat". According my timing tests, it is about 6 to 7 times faster (on a Pentium 2.4 MHz) machine than the default Inversion method used in

Hi, Quick question. With an old phone system a receptionist receiving a call has 1 button to push to transfer calls to a specific extension, with Asterisk, a receptionist would actually put the caller on hold, pick up another line, call the extension, ask if the person is available, hang up pick up the caller again and transfer. To me it's seems a long way to simply do a receptionist

On r-help Bill Simpson wrote: >Marsaglia's multiply-with-carry generator is very easy to implement--he >has posted some C code to newsgroups--and seems to test out OK. I wonder >if it would make sense to use one of his generators instead. I would like a system where the random number generator can be specified (perhaps in options or a .RandomGenerator variable). Most often everyone

Hi, I am working on a population of an invasive clam. The data are the size of each clam per station (2mm on average). Each station is found at a different distance from a power nuclear station, so at different water temperatures. The fist step I want to do is to identify cohort size at each station or (zone of water temperature). The second step will be to see whether the size or number of

I've been beating my head against this one and just can't figure it out. I hope someone here may have an answer. The employees using the workstations on this network are getting increasingly upset with this problem. The problem is wildly varying logon and logoff times over the network. This is definitely not a matter of long profile transfers. An individual can log onto a workstation

Hi All. I have R code whose functionality is being replicated within a C+ program. The outputs are to be compared to validate the conversion somewhat - however (as is always the case) I have stuffed my code with random number calls. Random uniform numbers in C+ are being produced using the (Boost) mersenne-twister generators (mt11213b & mt19937) - which is the default type of generator

I hope someone can enlighten me on this. Situation: NT network, Samba PDC, about 20 NT 4.0 workstations. log(s).smbd are created per machine for easier analysis (as log.smbd.<NetBios Name>). Not often, but often enough to be of concern, are errors in request_oplock_break that seem to indicate that another smbd process that should be listening for break requests on a UDP socket

Hi, I would like to extract the fitted values from a model using LMER but only for the fix portion of the model and not for the fix and random portion (e.g it is the procedure outpm or outp in SAS). I am aware of the procedure fitted() but I not sure it give the fitted values both for the fixed and random or only the fixed. I looked in the r help and the r list and I haven?t not found much

Dear Samba and Linux list users - Several people have commented on the volume of spam reaching the list. Andrew (the list maintainer, I am but a lowly assistant) has been advised of your comments and will take action if he thinks it's necessary. People wishing to make further comment should do so directly to the list admins and NOT, repeat NOT to the list itself! My thanks to those list

Thanks for telling me that you could not get my message, I hope this work better... so my question was: I built a population matrix to which I applied the fonction eigen in order to find the main parameters about my population. I know that the first eigen value correspond to lambda or exponential growth rate of my population. My problem is that I want to have the 95% confidence interval of the

Tom Dickson (bombcar@bombcar.com) wrote: > >-----BEGIN PGP SIGNED MESSAGE----- >Hash: SHA1 > >Here is the issue, and why the Samba list is basically unable to do >anything: > >If you send a message to the Samba list, it gets resent to everyone on >the list. This message will contain a FROM: line, showing an email address. Which should be <samba@lists.samba.org>

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

Hi, I have a few questions about the RNG in R; apologies if these are dumb questions: 1. It is my understanding that, among the three types of random number generators available in R now, the best one is the Marsaglia Multicarry. Is this correct? 2. How does the best RNG in R compare (in terms of quality) to the RNG in SPlus? (based on Marsaglia's Super Duper)? Does the Super-Duper in