Re: [Rd] k means

2008-05-16 Thread cgenolin
Hi the list I try the flexclust, but I do not manage to see what is wrong in my (very simple) code... Will you have few minutes to check it? Thanks for your help. Christophe --- 8< data <- rbind(c(1,2 ,NA,4 ), c(1,1 ,NA,1 ), c(2,3

Re: [Rd] HoltWinters fitted level parameter not bounded between 0 and 1 (PR#11469)

2008-05-16 Thread John Bodley
Hi, Thanks for the quick response. I upgraded by version of R on Windows to the latest (2.7.0) and re-ran the analysis and get the same result of 48.87989. The original time series was a non-regular zoo() object and I converted it to a time series with daily granularity by setting the values of t

Re: [Rd] basename/dirname produce incorrect results

2008-05-16 Thread ronggui
I have beend realized the problem of reporting with non-English character this way. Just don't know how to do in a better way. Sorry for that. Now, the problem reported has been fixed ( in version (2008-05-15 r45703)). Thanks for your work. Best 2008/5/16 Prof Brian Ripley <[EMAIL PROTECTED]>: >

Re: [Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread Simon Urbanek
Robert, this was discussed before: https://stat.ethz.ch/pipermail/r-devel/2007-December/047594.html and it *is* mentioned in NEWS: o co[rv](use = "complete.obs") now always gives an error if there are no complete cases: they used to give NA if method = "pearson" but an err

Re: [Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread McGehee, Robert
Oops, as I just realized, var does have a 'use' argument in 2.6.2, so I can just use Gabor's suggestion for var. Sorry for that Gabor. -Original Message- From: McGehee, Robert Sent: Friday, May 16, 2008 11:20 AM To: 'Gabor Grothendieck' Cc: R-devel Subject: RE: [Rd] var/sd and NAs in R2.7

Re: [Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread McGehee, Robert
I know I can get around this, I just would prefer that if R is breaking backwards compatibility, then it's intentional (maybe it is, I just don't know). That is, I don't want to require my entire company to upgrade to 2.7.0 just so I can deploy a fix here, and I'd prefer not to check the argument l

Re: [Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread Gerlanc, Daniel
Perhaps _sd_ should take a ... argument. -- Dan -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gabor Grothendieck Sent: Friday, May 16, 2008 11:03 AM To: McGehee, Robert Cc: R-devel Subject: Re: [Rd] var/sd and NAs in R2.7.0 Try var(c(NA, NA, NA), use

Re: [Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread Gabor Grothendieck
Try var(c(NA, NA, NA), use = "pairwise.complete.obs") On Fri, May 16, 2008 at 10:56 AM, McGehee, Robert <[EMAIL PROTECTED]> wrote: > Hello all, > I just upgraded to R 2.7.0 and found that the behavior of 'var' and 'sd' > have changed in the presence NAs (this wasn't explicit in the NEWS file, >

[Rd] var/sd and NAs in R2.7.0

2008-05-16 Thread McGehee, Robert
Hello all, I just upgraded to R 2.7.0 and found that the behavior of 'var' and 'sd' have changed in the presence NAs (this wasn't explicit in the NEWS file, though I see it probably has to do with the change for cor/cov). Anyway, I just want to make sure that it was intentional to produce an error

Re: [Rd] Configure errors say to report here... (PR#11470)

2008-05-16 Thread Simon Urbanek
Adam, this it not a bug in R, this is merely a user error. Please don't file self-inflicted problems as bugs (see posting guide). For Mac-related questions there is R-SIG-Mac. Brian's response was to the point - you break the compilation with the custom flags you use. See details below. O

Re: [Rd] Sweave.sty problems under WinXP (English locale) with default settings

2008-05-16 Thread Duncan Murdoch
On 5/16/2008 9:22 AM, Peter Ruckdeschel wrote: Hi R-devels, I encounter the following problem when I want to build a vignette in the package building process under Windows platform i386-pc-mingw32 arch i386 os mingw32

[Rd] Sweave.sty problems under WinXP (English locale) with default settings

2008-05-16 Thread Peter Ruckdeschel
Hi R-devels, I encounter the following problem when I want to build a vignette in the package building process under Windows platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32

Re: [Rd] interrupting native code

2008-05-16 Thread Luke Tierney
I'm not sure you can make this work as some of the things needed either are or should be private to the core implementation and not available to package code. In any case I would not recommend this approach for two reasons. First, details of what happens in interrupt checking are subject to chan

Re: [Rd] interrupting native code

2008-05-16 Thread Prof Brian Ripley
On Fri, 16 May 2008, Kjell Konis wrote: You mean something like this (I return 1 instead of calling onintr())? Will HAVE_AQUA and Win32 be appropriately defined when building my package (I can't see how to check with R CMD config)? HAVE_AQUA is in Rconfig.h, and WIN32 (not Win32) is defined f

Re: [Rd] HoltWinters fitted level parameter not bounded between 0 (PR#11473)

2008-05-16 Thread david . meyer
I get John's value (48.8789) in 2.7.0 and R-devel (both on Ubuntu). Really seems to be a numeric issue: > HoltWinters(x, beta = 0, gamma = 0)$alpha alpha 48.87989 > HoltWinters(x * 1.01, beta = 0, gamma = 0)$alpha alpha 0.6881547 > HoltWinters(x * 1.001, beta = 0, ga

Re: [Rd] HoltWinters fitted level parameter not bounded between 0 and 1 (PR#11469)

2008-05-16 Thread David Meyer
I get John's value (48.8789) in 2.7.0 and R-devel (both on Ubuntu). Really seems to be a numeric issue: > HoltWinters(x, beta = 0, gamma = 0)$alpha alpha 48.87989 > HoltWinters(x * 1.01, beta = 0, gamma = 0)$alpha alpha 0.6881547 > HoltWinters(x * 1.001, beta = 0, gamma

Re: [Rd] interrupting native code

2008-05-16 Thread Kjell Konis
You mean something like this (I return 1 instead of calling onintr())? Will HAVE_AQUA and Win32 be appropriately defined when building my package (I can't see how to check with R CMD config)? int My_CheckUserInterrupt(void) { R_CheckStack(); #if ( defined(HAVE_AQUA) ) /* R_ProcessEve

Re: [Rd] HoltWinters fitted level parameter not bounded between 0 (PR#11472)

2008-05-16 Thread ripley
It doesn't do it on my system (I get a value of about 0.688 in R 2.7.0 patched on Linux), and 2.5.1 is not current. Does a better starting value help? However, HoltWinters is using optim() in a case it is not designed for (one-dimensional optimization): see the note on its help page. I think

Re: [Rd] HoltWinters fitted level parameter not bounded between 0 and 1 (PR#11469)

2008-05-16 Thread Prof Brian Ripley
It doesn't do it on my system (I get a value of about 0.688 in R 2.7.0 patched on Linux), and 2.5.1 is not current. Does a better starting value help? However, HoltWinters is using optim() in a case it is not designed for (one-dimensional optimization): see the note on its help page. I think

Re: [Rd] interrupting native code

2008-05-16 Thread Prof Brian Ripley
On Fri, 16 May 2008, Kjell Konis wrote: The problem is that my package uses an external pointer to keep track of a structure created by the lp_solve library. If I use R_CheckUserInterrupt in the lp_solve abort function it leaves the structure in a messed-up state after an interrupt occurs. I a

Re: [Rd] interrupting native code

2008-05-16 Thread Kjell Konis
The problem is that my package uses an external pointer to keep track of a structure created by the lp_solve library. If I use R_CheckUserInterrupt in the lp_solve abort function it leaves the structure in a messed-up state after an interrupt occurs. I am not even able to free the memory al

[Rd] HoltWinters fitted level parameter not bounded between 0 and 1 (PR#11469)

2008-05-16 Thread john . bodley
Full_Name: John Bodley Version: 2.5.1 (2007-06-27) OS: Windows XP Submission from: (NULL) (12.144.182.66) I was fitting a number of time series in R using the stats::HoltWinters method to define a single exponential smoothing model, i.e., beta = gamma = 0. I came across an example where the fit

Re: [Rd] Configure errors say to report here... (PR#11470)

2008-05-16 Thread Prof Brian Ripley
What version are your Xcode tools? What X11 do you have installed (I see /sw/include in your path, and that is not normally needed)? Although configure says to report here, it is reporting a bug in your OS, not in R. There will be more information in your config.log, (and perhaps some of it

[Rd] Configure errors say to report here... (PR#11470)

2008-05-16 Thread adik
Hi, I'm attempting to compile R 2.7.0 on my G5, which is running OSX 10.5.2. Yes, I know there is a precompiled binary, but my experience has shown that if I hand-compile ATLAS on my machine and then link R against it when compiling, R runs faster, hence me compiling it myself. Anyway, I

Re: [Rd] ** operator

2008-05-16 Thread Karl Ove Hufthammer
Peter Dalgaard: > Not really, just transcribed during the lexical analysis phase: > > case '*': > if (nextchar('*')) > c='^'; > yytext[0] = c; > yytext[1] = '\0'; > yylval = install(yytext); > return c; > > (There's no "->" function either...) You can also use expression() to see what various e