similar to: Randomness not due to seed

Hi, How do I use a function inside a regression formula? When I tried to use min() function inside a regression formula, I ran into errors. Here are the commands to reproduce the error: mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/binary.csv")) attach(mydata) names(mydata) mylm2 = lm(admit~gre+I(max(gpa,rank))) *Error in model.frame.default(formula = admit ~

I am looking for a way to save a plot (graphics contents) to a file after the plot has been calculated but before it has been rendered. More specifically, assume that I made a plot that took a very long time to produce, I would like to save this plot to a generic file that I can later, on a different machine, render to either PDF, PNG, SVG using the usual R graphics devices, without recalculating

I am looking for a way to save a plot (graphics contents) to a file after the plot has been calculated but before it has been rendered. More specifically, assume that I made a plot that took a very long time to produce, I would like to save this plot to a generic file that I can later, on a different machine, render to either PDF, PNG, SVG using the usual R graphics devices, without recalculating

Dear all, I was doing an experiment to disprove some theory therefore performing lot of random simulation. Goal is to show the audience that although something has very rare chance to occur but it doesn't mean that event would be impossible. In this case after getting that rare event I need to show that same scenario for multiple times to explain other audience. Hence I need to somehow

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

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

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