Re: [Rd] model.matrix evaluation challenges

2009-08-10 Thread Felix Andrews
how about... nrow(with(cc, model.matrix(params, data=environment( cheers -Felix 2009/8/10 Ben Bolker : >  I am having difficulty with evaluation/environment construction > for a formula to be evaluated by model.matrix().  Basically, I > want to construct a model matrix that first looks in "n

Re: [Rd] Bug in nlm, found using sem; failure in several flavors (PR#13883)

2009-08-10 Thread Peter Dalgaard
William Dunlap wrote: From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-project.org] On Behalf Of Ben Bolker Sent: Monday, August 10, 2009 12:52 PM To: r-devel@r-project.org Subject: Re: [Rd] Bug in nlm, found using sem; failure in several flavors (PR#13883) A wild guess, bas

Re: [Rd] Bug in nlm, found using sem; failure in several flavors (PR#13883)

2009-08-10 Thread William Dunlap
> From: r-devel-boun...@r-project.org > [mailto:r-devel-boun...@r-project.org] On Behalf Of Ben Bolker > Sent: Monday, August 10, 2009 12:52 PM > To: r-devel@r-project.org > Subject: Re: [Rd] Bug in nlm, found using sem; failure in > several flavors (PR#13883) > > > >A wild guess, based on

Re: [Rd] Bug in nlm, found using sem; failure in several flavors (PR#13883)

2009-08-10 Thread Ben Bolker
A wild guess, based on some recent "non-deterministic" behavior in lme4: could this be driven by a bug in optimized linear algebra libraries (BLAS/LAPACK)? (Unfortunately, I'm not clear on how to check the versions of BLAS/LAPACK being used and whether they are subject to bugs, but perhaps so

Re: [Rd] identical(0, -0)

2009-08-10 Thread Petr Savicky
On Mon, Aug 10, 2009 at 05:47:57AM -0400, Duncan Murdoch wrote: > I wouldn't mind a "strict" option. It would compare bit patterns, so > would distinguish +0 from -0, and different NaN values. I think that a logical option "strict" in the above meaning could be useful for debugging. The default

Re: [Rd] identical(0, -0)

2009-08-10 Thread Duncan Murdoch
For people who want to play with these, here are some functions that let you get or set the "payload" value in a NaN. NaN and NA, Inf and -Inf are stored quite similarly; these functions don't distinguish which of those you're working with. Regular finite values give NA for the payload value,

[Rd] Rcpp and variable protection

2009-08-10 Thread David Scherrer
Hi all, is there already a variable protection (as use else when linking to C/C++ e.g. PROTECT(a = AS_NUMERIC(a)); ) built in when using RcppExport from Rcpp? Many thanks, David [[alternative HTML version deleted]] __ R-devel@r-project.org mai

Re: [Rd] identical(0, -0)

2009-08-10 Thread Duncan Murdoch
On 8/10/2009 9:55 AM, Simon Urbanek wrote: On Aug 10, 2009, at 5:47 , Duncan Murdoch wrote: Petr Savicky wrote: On Sat, Aug 08, 2009 at 10:39:04AM -0400, Prof. John C Nash wrote: I'll save space and not include previous messages. My 2 cents: At the very least the documentation needs a fix.

Re: [Rd] multicore mclapply hangs

2009-08-10 Thread Simon Urbanek
On Aug 10, 2009, at 6:17 , Markus Schmidberger wrote: Hello Rune, please use the mailinglist 'R-sig-hpc' (https://stat.ethz.ch/mailman/listinfo/r-sig-hpc , do not forget to register!) for this topic. .. or even better ask the package maintainer (with more precise details such as the exa

Re: [Rd] identical(0, -0)

2009-08-10 Thread Simon Urbanek
On Aug 10, 2009, at 5:47 , Duncan Murdoch wrote: Petr Savicky wrote: On Sat, Aug 08, 2009 at 10:39:04AM -0400, Prof. John C Nash wrote: I'll save space and not include previous messages. My 2 cents: At the very least the documentation needs a fix. If it is easy to do, then Ted Harding's s

Re: [Rd] Error Building R + X11 on SuSe 10S

2009-08-10 Thread rory . winston
Thanks Duncan. I have since found that building cairo and adding '--with-cairo' to the configure command solves this problem, and also gives (at least to my eyes) much nicer looking 2D png() plots. Cheers -- Rory On Aug 10, 2009 12:54pm, Duncan Murdoch wrote: > rory.wins...@gmail.com wrote:

Re: [Rd] Error Building R + X11 on SuSe 10S

2009-08-10 Thread Duncan Murdoch
rory.wins...@gmail.com wrote: Following on... Sorry folks, it looks like I have misdiagnosed the issue. When I connect to the server using an X client and start up the newly built R instance, I see capabilities() shows X11 support as expected. So this changes the character of my query somewh

Re: [Rd] Error Building R + X11 on SuSe 10S

2009-08-10 Thread rory . winston
Following on... Sorry folks, it looks like I have misdiagnosed the issue. When I connect to the server using an X client and start up the newly built R instance, I see capabilities() shows X11 support as expected. So this changes the character of my query somewhat: As I was hoping to run R in

Re: [Rd] multicore mclapply hangs

2009-08-10 Thread Markus Schmidberger
Hello Rune, please use the mailinglist 'R-sig-hpc' (https://stat.ethz.ch/mailman/listinfo/r-sig-hpc, do not forget to register!) for this topic. And please provide some more information: * in R: sessionInfo() * and your used operation system Best Markus Rune Schjellerup Philosof schrieb:

[Rd] Error Building R + X11 on SuSe 10

2009-08-10 Thread rory . winston
Hi all I am having problems building R + x11 support on the following system: # uname -a Linux 2.6.16.60-0.27-bigsmp #1 SMP Mon Jul 28 13:06:32 UTC 2008 i686 i686 i386 GNU/Linux # cat /etc/SuSE-release SUSE Linux Enterprise Server 10 (i586) VERSION = 10 PATCHLEVEL = 2 I am using the source do

[Rd] multicore mclapply hangs

2009-08-10 Thread Rune Schjellerup Philosof
When I execute mclapply it creates the needed processes, but these processes never begin computing anything, they just wait indefinitely. I recently upgraded to version 2.9.1, which might have caused the problem. -- Med venlig hilsen Rune Schjellerup Philosof Ph.d.-studerende, Statistik, IST, SDU

Re: [Rd] identical(0, -0)

2009-08-10 Thread Duncan Murdoch
Petr Savicky wrote: On Sat, Aug 08, 2009 at 10:39:04AM -0400, Prof. John C Nash wrote: I'll save space and not include previous messages. My 2 cents: At the very least the documentation needs a fix. If it is easy to do, then Ted Harding's suggestion of a switch (default OFF) to check for s

[Rd] Duplicated lines in configure

2009-08-10 Thread rory . winston
Hi There seems to be two functionally equivalent lines in the configure script https://svn.r-project.org/R/trunk/configure: ... --with-gnu-ld assume the C compiler uses GNU ld [default=no] --with-x use the X Window System --with-gnu-ld assume the C compiler uses GNU ld default=no ... These see

Re: [Rd] identical(0, -0)

2009-08-10 Thread Petr Savicky
On Sat, Aug 08, 2009 at 10:39:04AM -0400, Prof. John C Nash wrote: > I'll save space and not include previous messages. > > My 2 cents: At the very least the documentation needs a fix. If it is > easy to do, then Ted Harding's suggestion of a switch (default OFF) to > check for sign difference w

Re: [Rd] Why is strptime always returning a vector of length 9 ?

2009-08-10 Thread Martin Maechler
{Correcting thinko below .. } > "MM" == Martin Maechler > on Mon, 10 Aug 2009 08:55:52 +0200 writes: > "l" == laurent > on Sun, 09 Aug 2009 21:45:07 +0200 writes: l> Thanks. It seems that the source of my confusion comes l> from using first using str() (and