similar to: using R's random numbers in another program

Dear All, I am having problems making and using the standalone library, after succesfully installing R from the tar.gz file (I am using Linux, with R-patched, but I have similar problems with R-devel). First, when I do: root at ligarto:/usr/lib/R-patched/src/nmath/standalone > make .... ar: mlutils.o: No such file or directory make[1]: *** [libRmath.a] Error 1 make[1]: Leaving directory

Dear all, I am trying to prepare some figures where I'd like the color of the plot region (i.e., the area bounded by the axes) to be different from the margin area (i.e., the area of the figure not within the axes region). I can use "bg" but that changes the background of everything? How can I accomplish what I want? Thanks, Ramon -- Ramón Díaz-Uriarte Triana 47 28016 Madrid

Dear All, A library (PHYLOGR) that passed the usual tests in R-1.1.1 gives errors with R-devel; my (mis?)understanding of scoping rules is that it should have worked in both. The problems seem related to using the name of the data frame for extracting weights or subsets within a function call. The problems can be reproduced as follows: ********************** datai <- data.frame( y =

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

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 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

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

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

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

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

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 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

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 R experts, Could anyone give some concrete examples of how to call R functions from a program written in C++? Thanks in advance for your help. Nathapong Samlamjiag _________________________________________________________________________ Get Your Private, Free E-mail from MSN Hotmail at http://www.hotmail.com. Share information about yourself, create your own public profile at

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

Hi, I'm seeing something that may be a bug in R's standalone math library, which is packaged by Debian as r-mathlib. I reported it to the Debian BTS as http://bugs.debian.org/657573 I'm using Debian squeeze, and the code was tested with r-mathlib 2.11.1-6 (default on stable) and 2.14.1-1 (from testing/unstable). I summarize this report below. The following code with the R math

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)

Recieving digests. > RNGkind(NULL) [1] "Marsaglia-Multicarry" "Kinderman-Ramage" I would appreciate it if anybody has any comments on the following. Please do not comment on the R functions themselves, since they merely mimic a (bivariate simplification of a) C routine called from S. In particular, I would like to know if anything is available with regard to the

> x <- matrix(1:10, ncol = 2) > xc <- cor(x[, 2], x[, 0]) > cbind(xc, vector()) Process R floating point exception at Tue Feb 11 19:16:51 2003 > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.2 year