Re: [Rd] Recursion error after upgrade to R_2.11.1 [Sec=Unclassified] [Sec=Unclassified]

2010-10-07 Thread Martin Maechler
TR == Troy Robertson troy.robert...@aad.gov.au on Thu, 7 Oct 2010 13:50:49 +1100 writes: On 10/06/2010 06:12 PM, Troy Robertson wrote: Hi all, After no replies to my previous message I thought I might show some code to demonstrate the change and again seek

[Rd] long double, C, fortran

2010-10-07 Thread Richard D. Morey
I'm using .Call() to call C code from R under Windows (on an Intel Core 2 duo). The C code involves some very small numbers, and I think I'm losing precision using doubles. I thought I might use long doubles to see if I can get that precision back. I have a few questions: 1. Does this affect

Re: [Rd] long double, C, fortran

2010-10-07 Thread Simon Urbanek
On Oct 7, 2010, at 2:28 PM, Richard D. Morey wrote: I'm using .Call() to call C code from R under Windows (on an Intel Core 2 duo). The C code involves some very small numbers, and I think I'm losing precision using doubles. The magnitude of numbers has really nothing to do with precision

Re: [Rd] cleanup_pkg() in tools:::.build_packages() is broken

2010-10-07 Thread Hervé Pagès
This has been fixed in R-beta. Thanks! H. On 10/01/2010 06:31 PM, Hervé Pagès wrote: Hi, The cleanup_pkg() function defined the big tools:::.build_packages() function in tools/R/build.R is currently broken. When Makefiles are used cleanup_pkg() doesn't clean anything because of the way

Re: [Rd] tar problem when using 'R CMD build' on Windows

2010-10-07 Thread Hervé Pagès
Hi, On 09/30/2010 10:45 AM, Hervé Pagès wrote: [...] If that is not desirable, could the tools:::.build_packages() function not put quotes around the value of the TAR environment variable (the code does shQuote(TAR)), at least on Windows, so I can set this variable to tar --no-same-owner

Re: [Rd] long double, C, fortran

2010-10-07 Thread Matt Shotwell
In working with small numbers, I usually look for a solution involving a log transformation. For instance, when a or b are very small, we can compute the ratio a / ( a + b ) more safely by 1 / ( 1 + exp( log(b) - log(a) ) ) Alternative data types are sure to cause issues with portability. It's

Re: [Rd] Recursion error after upgrade to R_2.11.1 [Sec=Unclassified] [Sec=Unclassified]

2010-10-07 Thread Troy Robertson
-Original Message- From: Martin Maechler [mailto:maech...@stat.math.ethz.ch] Sent: Thursday, 7 October 2010 9:50 PM To: Troy Robertson Cc: 'Martin Morgan'; 'r-devel@R-project.org'; 'John Chambers' Subject: Re: [Rd] Recursion error after upgrade to R_2.11.1 [Sec=Unclassified]

[Rd] model.frame deficiency

2010-10-07 Thread Terry Therneau
The model.frame function has trouble with a certain type of really long formula. Here is a test: tname - paste('var', 1:50, sep='') tmat - matrix(rnorm(500), ncol=50, dimnames=list(NULL, tname)) tdata - data.frame(tmat) temp1 - paste( paste(tname, tname, sep='='), collapse=', ') temp2 -

Re: [Rd] Recursion error after upgrade to R_2.11.1 [Sec=Unclassified] [Sec=Unclassified]

2010-10-07 Thread Peter Dalgaard
On 10/07/2010 12:58 PM, Duncan Murdoch wrote: R 2.12.2 (currently beta) has in its NEWS : Just pointing out a typo: the version of the beta is 2.12.0. Duncan Murdoch o A facility for defining reference-based S4 classes (in the OOP style of Java, C++, etc.) has been added

Re: [Rd] Recursion error after upgrade to R_2.11.1 [Sec=Unclassified] [Sec=Unclassified]

2010-10-07 Thread Kasper Daniel Hansen
Troy For what it's worth, the ExpressionSet / eSet class in Biobase from Bioconductor has for a long time implemented a pass-by-reference system, using amongst other stuff locked environments to make sure things are read only. You might find this useful to look at. In general it is hard to