similar to: How to obtain seed after generating random number?

I am working on a reproducible computing platform for which I would like to be able to _exactly_ reproduce an R object. However, I am experiencing unexpected randomness in some calculations. I have a hard time finding out exactly how it occurs. The code below illustrates the issue. mylm1 <- lm(dist~speed, data=cars); mylm2 <- lm(dist~speed, data=cars); identical(mylm1, mylm2); #TRUE

Hello Timo, I'm running dovecot-1.2.11/Maildir and plan to migrate to single UID mailboxes some day, but for now, I've got "system" users and I'm testing permissions handling in order to set up shared mailboxes. The private namespace mailboxes location is location = maildir:/courriel/boites/%u:CONTROL=/courriel/meta/%u:INDEX=/var/dovecot-test/indexes/%1u/%u I've

I run a decent amount of 3ware hardware, all under centos-4. There seems to be some sort of fundamental disagreement between ext3 and 3ware's hardware RAID5 mode that trashes write performance. As a representative example, one current setup is 2 9550SX-12 boards in hardware RAID5 mode (256KB stripe size) with a software RAID0 stripe on top (also 256KB chunks). bonnie++ results look

Hello, I'm trying to set up, as explained in the wiki, shared mailboxes with symlink. I'm running dovecot-1.1.8 with LDAP pass/userdb. The Maildir layout goes like this : /courriel/boites/foobar : Maildir of user foobar (nfs) /courriel/meta/foobar : control files of user foobar (nfs) /var/dovecot/indexes : indexes (locally stored) I created : /courriel/boites/doveshared/.box/ :

Is there an easy way to compile a packages's fortran so that it has the information for ddd or gdb to follow the lines of the code? I seem to remember seeing something about this, but I cannot find it in Writing R Extensions. (I see about Valgrind, debugging a segfault, and some other tricks, but what I really want is to follow the fortran logic.) Thanks for any pointers, Paul

Hello all, I have downloaded samba-2.2.8a-1-sol8-suncc-32bit.pkg.gz and I am trying to gunzip it but I get gunzip: samba-2.2.8a-1-sol8-gcc295.pkg.gz: not in gzip format what could be the problem ? ********************************************************************************* Ce courriel ainsi que ses pieces jointes sont strictement reserves a l'usage de la ou du destinataire et peut

Is CRAN having trouble or is it just me? (The web interface is very slow and install.packages() is timing out and giving me "service not available"). Paul ==================================================================================== La version française suit le texte anglais. ------------------------------------------------------------------------------------ This email

I am trying to extend an S3 method to work with an S4 object as well as the S3 objects it works with, but UseMethod does not seem to recognize the S4 class and dispatches to the default method. Is this to be expected or should I be looking for an error in my code? If it is not an error in my code, is there an easy way to do this, or do I have to convert the generic to S4 and then make those

Hi, R friends. I had some difficulty with the following code: qqnorm(freq, log='y') qqline(freq) as the line drawn was seemingly random. The exact data I used appears below. After wandering a bit within the source code for "abline", I figured out I should rather write: qqnorm(freq, log='y') par(ylog=FALSE) qqline(log10(freq)) par(ylog=TRUE)

I am trying to check my packages with R-rc and latex is failing to find utf8.def. I presume my latex installation is defective or too old, or is this file distributed with R and I just don't have a path set correctly? Paul ______________ * checking PDF version of manual ... WARNING LaTeX errors when creating PDF version. This typically indicates Rd problems. LaTeX errors found: !

I was trying to subscribe to an RSS feed from <http://developer.r-project.org/RSSfeeds.html> but IE chokes on DTDs and Firefox also has trouble with this feed. (I'm not sure if it has trouble with DTDs in general.) I went to one of our web developers to see what I should use to read the feed, and he pointed out that the mix of RSS 2.00 and rss-0.91 looks like a mistake: <!--

Hi, I could try latest 7zip (4.57) on a triple boot machine (Debian SID i386, x86_64 and native Windows XP SP2), and found something surprising in the results : amd64 at 2400MHz linux32 -> 1814 linux64 -> 1946 win32 -> 2004 So performance is much better in Linux 64 bit than in 32 bit, while both are slower than native windows system. Is it normal, I mean does anyone else obtain

Could anyone please explain how to set up contrasts between means in R. I want to know if "before I conduct an experiment and believe the mean for 1 and 2 will be different from means 3 and 4, Is this true?" That is what I have to prove or disprove, I thought that contrasts would be the way to go. Thanks for the help. Igor [[alternative HTML version deleted]]

I am having the following problem building R-beta_2009-10-12_r50045.tar.gz on Ubuntu 9.04 - the Jaunty Jackalope Paul ________ ... make[2]: Entering directory `/home/paul/toolchain/R/R-beta/src/library/Recommended' begin installing recommended package boot gzip: invalid option -- 'x' Try `gzip --help' for more information. tar: This does not look like a tar archive

Is it possible in R to call a fortran routine that sets variables in a common block and expect the values to persist when a call is made from R to a second routine that uses the common block? If not (as I suspect), is it possible to use a common block in a group of routines that are used together, for as long the routines do not return to R? Paul

Hi everyone, Can anyone give me more details about the 'lp', 'risk', 'terms', and 'expected' options from the predict.coxph function? I actually found that risk is obtained from exp(lp). I can't find information on 'expected' type however. Is anybody can expand the information on this function? Thank you in advance, I will appreciate it.

Unfortunately, RNGkind is buggy. It will not generate warnings except the full name "Buggy Kinderman-Ramage" is supplied for normal.kind. match.arg is supposed to be called before "==" comparison. ======================================== Shengqiao Li Research Associate The Department of Statistics PO Box 6330 West Virginia University Morgantown, WV 26506-6330

Hi, How to install Windows 2003 on Xen 3 without vt-enabled hardware? Thanks to your help. Jean-Marc Renson ********************************************************************************************* Ce courriel et les documents qui y sont attachés peuvent contenir des informations confidentielles. Si vous n''etes pas le destinataire escompté, merci d''en informer

I'm working on the next release of coxme (ready to start some testing), and have written major chunks of it using noweb -- very similar to Sweave except I'm generating code and documentation for the code rather than vingetes. The question: I have directory with .Rnw objects and a Makefile therein that generates most of the .R files. I don't see any guidance in the documentation

Dear users, I'm working on binary classification problem using Support Vector Machines (SVM). My objective is to train a series of SVM models on a grid of hyperparameters and then select those that maximize the AUC based on an independent validation sample. My attempted code is shown below. It runs well on "small" data sets but when I use it on a slightly larger sample (e.g., my