Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread ripley
Thank you, will be fixed in the next release of nlme (not yet scheduled, as everythign is code frozen for the release of R 2.5.1 tomorrow). On Tue, 26 Jun 2007, Andrzej Galecki wrote: > > Two attachments: > > 1. getVarCovBugReport.R - Rcode with an example illustrating the problem and > how to

Re: [Rd] slight anomaly in formals<- ? (PR#9758)

2007-06-26 Thread ripley
On Tue, 26 Jun 2007, [EMAIL PROTECTED] wrote: > Hi, > > The R input/output after the following paragraph is from a session with > version.string R version 2.4.0 Patched (2006-11-23 r39958). Which is not recent, but this is unchanged since. The problem is the definition: > `formals<-` function (

Re: [Rd] possible bug in 'scan'

2007-06-26 Thread Prof Brian Ripley
On Tue, 26 Jun 2007, Benjamin Tyner wrote: > R-devel, > > When I run the following code on the attached file, > > tmp <- scan("C:/temp.csv", > what=list("character","numeric"), > sep=",") > > Then tmp[[2]] is a character vector. My impression from the help file > is that it sho

[Rd] trivial typo in ?ks.test (PR#9759)

2007-06-26 Thread bolker
In ?ks.test "if the sample size if less than 100 in" should be "if the sample size is less than 100 in" Ben Bolker --please do not edit the information below-- Version: platform = i486-pc-linux-gnu arch = i486 os = linux-gnu system = i486, linux-gnu status = major = 2 minor = 5.0

Re: [Rd] "Math" group generics for S4, and a bug

2007-06-26 Thread Martin Maechler
> "JMC" == John Chambers <[EMAIL PROTECTED]> > on Tue, 26 Jun 2007 14:47:25 -0400 writes: JMC> Martin Maechler wrote: >>> "MM" == Martin Maechler <[EMAIL PROTECTED]> >>> on Sat, 23 Jun 2007 00:36:43 +0200 writes: >>> >> > {on R-help} >> >> [..

Re: [Rd] "Math" group generics for S4, and a bug

2007-06-26 Thread Martin Maechler
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]> > on Tue, 26 Jun 2007 18:38:27 +0100 (BST) writes: BDR> Remember that you only get S4 dispatch on Math for S4 BDR> generics. log10 is not an S4 generic until you make it BDR> one. From the help page BDR> Note: cu

Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread agalecki
This is a multi-part message in MIME format. --000206090008050103050209 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Two attachments: 1. getVarCovBugReport.R - Rcode with an example illustrating the problem and how to fix it 2. getVarC

[Rd] slight anomaly in formals<- ? (PR#9758)

2007-06-26 Thread georgi . boshnakov
Hi, The R input/output after the following paragraph is from a session with version.string R version 2.4.0 Patched (2006-11-23 r39958). The last element, x$c, of x has no value and after the assignment of x to the formals of f, x$c does not become a formal argument of f. The second assignment doe

Re: [Rd] "Math" group generics for S4, and a bug

2007-06-26 Thread John Chambers
Martin Maechler wrote: >> "MM" == Martin Maechler <[EMAIL PROTECTED]> >> on Sat, 23 Jun 2007 00:36:43 +0200 writes: >> > > {on R-help} > > [.] > [.] > > >> Duncan Murdoch > > DM> You might have better luck with >

Re: [Rd] Behaviour of mle and environments or calling mle inside a function

2007-06-26 Thread Prof Brian Ripley
On Tue, 26 Jun 2007, Iago Mosqueira wrote: > Dear all, > > I would appreciate some help understanding the following behaviour > when stats4::mle is called inside a function. mle seems to look for > its arguments in R_GlobalEnv and not the environment from which it is > called. It is your function

Re: [Rd] "Math" group generics for S4, and a bug

2007-06-26 Thread Prof Brian Ripley
Remember that you only get S4 dispatch on Math for S4 generics. log10 is not an S4 generic until you make it one. From the help page Note: currently those members which are not primitive functions must have been converted to S4 generic functions (preferably _before_ setting an S

[Rd] "Math" group generics for S4, and a bug

2007-06-26 Thread Martin Maechler
> "MM" == Martin Maechler <[EMAIL PROTECTED]> > on Sat, 23 Jun 2007 00:36:43 +0200 writes: {on R-help} [.] [.] >> Duncan Murdoch DM> You might have better luck with DM> log1p(tasa) MM> {very good point, thank you, D

Re: [Rd] conflict between lme4 and RMySQL packages (PR#9753)

2007-06-26 Thread Duncan Murdoch
On 6/26/2007 9:07 AM, Martin Maechler wrote: >> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]> >> on Tue, 26 Jun 2007 10:51:03 +0100 (BST) writes: > > BDR> You have three summary() functions here: > >> library(RMySQL) > BDR> Loading required package: DBI > >> library(lm

Re: [Rd] possible bug in 'scan'

2007-06-26 Thread Dan Davison
Hi Ben, On Tue, Jun 26, 2007 at 11:58:35AM -0400, Benjamin Tyner wrote: > R-devel, > > When I run the following code on the attached file, > > tmp <- scan("C:/temp.csv", >what=list("character","numeric"), >sep=",") > > Then tmp[[2]] is a character vector. My impression f

Re: [Rd] possible bug in 'scan'

2007-06-26 Thread Benjamin Tyner
How embarrassing -- works like a charm. On 6/26/07, Greg Snow <[EMAIL PROTECTED]> wrote: > The what argument looks at what the elements are, not the word they say. > Try this: > > > tmp <- scan("C:/temp.csv", > > what=list("",0), > > sep=",") > > Hope this helps, > > -- > G

[Rd] Behaviour of mle and environments or calling mle inside a function

2007-06-26 Thread Iago Mosqueira
Dear all, I would appreciate some help understanding the following behaviour when stats4::mle is called inside a function. mle seems to look for its arguments in R_GlobalEnv and not the environment from which it is called. library(stats4) lkhd <- function(alpha=1, beta=0.1, sigma=0.1) -

Re: [Rd] possible bug in 'scan'

2007-06-26 Thread Greg Snow
The what argument looks at what the elements are, not the word they say. Try this: > tmp <- scan("C:/temp.csv", > what=list("",0), > sep=",") Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-811

[Rd] possible bug in 'scan'

2007-06-26 Thread Benjamin Tyner
R-devel, When I run the following code on the attached file, tmp <- scan("C:/temp.csv", what=list("character","numeric"), sep=",") Then tmp[[2]] is a character vector. My impression from the help file is that it should be a numeric as specified by 'what' sessionInfo() R

Re: [Rd] R 2.5.0 refuses to print enough digits to recover exact floating point values

2007-06-26 Thread Petr Savicky
I would like to reply to the following email from May in the thread [Rd] R 2.5.0 refuses to print enough digits to recover exact floating point values On Wed, May 23, 2007 at 06:32:36PM +0100, Prof Brian Ripley wrote: > I think this is a bug in the MacOS X runtime. I've checked the C99 > stand

Re: [Rd] conflict between lme4 and RMySQL packages (PR#9753)

2007-06-26 Thread Martin Maechler
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]> > on Tue, 26 Jun 2007 10:51:03 +0100 (BST) writes: BDR> You have three summary() functions here: >> library(RMySQL) BDR> Loading required package: DBI >> library(lme4) BDR> Loading required package: Matrix BDR> Lo

Re: [Rd] boxplot and bxp do not respect xlim by default (PR#9756)

2007-06-26 Thread S . Ellison
What is mystifying is that the issue was not present in previous versions, = so appropriate code already existed. However, I agree that there seem to be a couple of additional issues that = I had missed. =20 I am perfectly happy to look at this again myself, though, and provide = extended code; w

Re: [Rd] boxplot and bxp do not respect xlim by default (PR#9754)

2007-06-26 Thread murdoch
On 6/26/2007 8:16 AM, [EMAIL PROTECTED] wrote: > Full_Name: Steve Ellison > Version: 2.4.1 > OS: Windows, Linux > Submission from: (NULL) (194.73.101.157) > > > bxp() allows specifcation of box locations with at=, but neither adjusts xlim= > to fit at nor does it respect xlim provided explicitly.

Re: [Rd] boxplot and bxp do not respect xlim by default (PR#9754)

2007-06-26 Thread Duncan Murdoch
On 6/26/2007 8:16 AM, [EMAIL PROTECTED] wrote: > Full_Name: Steve Ellison > Version: 2.4.1 > OS: Windows, Linux > Submission from: (NULL) (194.73.101.157) That version is obsolete, but this problem is still present in R-devel. I'll take a look. Duncan Murdoch > > > bxp() allows specifcation

[Rd] boxplot and bxp do not respect xlim by default (PR#9754)

2007-06-26 Thread s . ellison
Full_Name: Steve Ellison Version: 2.4.1 OS: Windows, Linux Submission from: (NULL) (194.73.101.157) bxp() allows specifcation of box locations with at=, but neither adjusts xlim= to fit at nor does it respect xlim provided explicitly. This is because bxp() now includes explicit xlim as c(0.5, n+

Re: [Rd] conflict between lme4 and RMySQL packages (PR#9753)

2007-06-26 Thread Prof Brian Ripley
You have three summary() functions here: > library(RMySQL) Loading required package: DBI > library(lme4) Loading required package: Matrix Loading required package: lattice > find("summary") [1] "package:Matrix" "package:DBI""package:base" and this is simply a namespace issue (note the environ

Re: [Rd] conflict between lme4 and RMySQL packages (PR#9753)

2007-06-26 Thread maechler
Thank you, Dale. I think you've revealed a bug in the method dispatch mechanism. I'm adding here a script which confirms your finding and add slightly more insight: library(RMySQL) library(lme4) data(sleepstudy) fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy) ## MM: First save, then p

Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread ripley
On Mon, 25 Jun 2007, [EMAIL PROTECTED] wrote: > I am using R2.5 under Windows. I presume you mean 2.5.0 (there is no R2.5: see the posting guide). But which version of nlme, which is the relevant fact here? The R posting guide suggests showing the output of sessionInfo() to establish the env

Re: [Rd] R routine standalone

2007-06-26 Thread Prof Brian Ripley
You have not told us your OS nor version of R. The file src/nmath/standalone/README in the R sources will either tell you how to do this or point you to the current documentation (depending on the R version). At a guess (since you use 'project' and don't have the correct case on your file name

[Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread agalecki
Hello, I am using R2.5 under Windows. Looks like the following statement vars <- (obj$sigma^2)*vw in getVarCov.gls method (nlme package) needs to be replaced with: vars <- (obj$sigma*vw)^2 With best regards Andrzej Galecki Douglas Bates wrote: >I'm not sure when the getVarCov.

[Rd] R routine standalone

2007-06-26 Thread Gatsu
I'm writing a C program and I'd like to include in it some R routine (quantile functions, typically). I found the functions I have to use (qnorm.c qbeta.c ecc) e some include (nmath.h rmath.h ecc). But I can't compile the project, I receive linker error. Somebody can help me? I say, if I want

[Rd] conflict between lme4 and RMySQL packages (PR#9753)

2007-06-26 Thread dale . barr
Full_Name: Dale Barr Version: 2.5.1 (patched) OS: Ubuntu linux x86_64 Submission from: (NULL) (138.23.70.108) When RMySQL is loaded in before lme4, the summary() function for lmer objects in the lme4 packages produces the following error: Error in printMer(object) : no slot of name "status" for