[Rd] same test statistic for t-test with and without equal variance assumption

2004-10-01 Thread Adaikalavan Ramasamy
Could some kindly tell me if I am supposed to be getting the same test statistic value with var.equal=TRUE and var.equal=FALSE in t.test ? set.seed(1066) x1 <- rnorm(50) x2 <- rnorm(50) t.test(x1, x2, var.equal=FALSE)$statistic # 0.5989774 t.test(x1, x2, var.equal=TRUE)$statistic # 0.5989774 ???

Re: [Rd] same test statistic for t-test with and without equal variance assumption

2004-10-01 Thread Peter Dalgaard
Adaikalavan Ramasamy <[EMAIL PROTECTED]> writes: > Could some kindly tell me if I am supposed to be getting the same test > statistic value with var.equal=TRUE and var.equal=FALSE in t.test ? > > set.seed(1066) > x1 <- rnorm(50) > x2 <- rnorm(50) > > t.test(x1, x2, var.equal=FALSE)$statistic # 0

[Rd] bzip2 directory won't build on OSF1 due to C99 code and -std1 option (PR#7257)

2004-10-01 Thread hoffman+R
Full_Name: Michael Hoffman Version: 2.0.0-beta-20041001 OS: OSF1 V5.1 Submission from: (NULL) (193.62.199.8) Hello. Building R 2.0.0-beta-20041001 on OSF1 V5.1 failed because the default configure sets R_XTRA_CFLAGS to "-std1 -ieee_with_inexact." The bzip2 directory includes C99

Re: [Rd] same test statistic for t-test with and without equal variance assumption

2004-10-01 Thread Adaikalavan Ramasamy
Peter, thank you! I forgot the to square root in calculating sp. sp <- sqrt( ( (n1-1)*v1 + (n2-1)*v2 )/(n1 + n2 - 2) ) For several simulation runs, the test statistics from both tests are remarkably similar (difference is less than 10e-16). I naively assumed that the statistics value should be sl

(PR#7257) Re: [Rd] bzip2 directory won't build on OSF1 due to C99 code and -std1 option

2004-10-01 Thread Prof Brian Ripley
inlining (and for which there is a configure test). Note that feedback 3 days before release is too late for configuration changes to be made. On Fri, 1 Oct 2004 [EMAIL PROTECTED] wrote: > Full_Name: Michael Hoffman > Version: 2.0.0-beta-20041001 > OS: OSF1 V5.1 > Submission from: (NU

Re: [Rd] same test statistic for t-test with and without equal variance assumption

2004-10-01 Thread Adaikalavan Ramasamy
Apologies, I made an incorrect statement. The statistics from t-test with difference variance assumptions are the same if both groups have the same length. Sorry for troubling everyone again. On Fri, 2004-10-01 at 13:56, Adaikalavan Ramasamy wrote: > Peter, thank you! I forgot the to square root

Re: (PR#7257) Re: [Rd] bzip2 directory won't build on OSF1 due to C99 code and -std1 option

2004-10-01 Thread Michael Hoffman
On Fri, 1 Oct 2004, Prof Brian Ripley wrote: What `C99 code' is that? Specifically it uses the inline keyword. I do not believe that is supported by a strict C89 compiler. Please don't expect us to guess what you have in mind, for we haven't a clue. Please accept my apologies for an inadequate bug

Re: (PR#7257) Re: [Rd] bzip2 directory won't build on OSF1 due to C99 code and -std1 option

2004-10-01 Thread Prof Brian Ripley
I believe this occurs because the configure test for `inline' is not using the R_XTRA_CFLAGS, which is a pervasive design problem (we have noted it occurs with PIC flags as well). If you set -std1 in CFLAGS then it should be used in the configure test which should conclude `inline' is not allowe

Re: (PR#7257) Re: [Rd] bzip2 directory won't build on OSF1 due to C99 code and -std1 option

2004-10-01 Thread Michael Hoffman
If you set -std1 in CFLAGS then it should be used in the configure test which should conclude `inline' is not allowed. OK, I will do as you suggest. I should note that cc (version info below) does support inline although not with -std1 so I'd prefer to have any performance benefits of inline, if po

Re: [Rd] (PR#7254)Documentation: Reference Index (.pdf) -- setOldClass

2004-10-01 Thread Prof Brian Ripley
Have you looked at the 2.0.0 beta version? It is hard for us to know what you are looking at, but I believe this was fixed a while ago (but post 1.9.1). Better PDF cross-refs is one of the changes in 2.0.0-to-be. On Thu, 30 Sep 2004, Pfaff, Bernhard wrote: > Apologies for cross postings or ma

[Rd] Re: [R] R-2.0: roadmap? release statements? plans?

2004-10-01 Thread Duncan Murdoch
(Moved from r-help to r-devel) On Fri, 1 Oct 2004 10:00:05 -0400, "Liaw, Andy" <[EMAIL PROTECTED]> wrote : >Just to keep Valery from thinking that that web page can not be easily >found: Go to www.r-project.org and click on `Developer Page' under the >heading `R Project'. > >A suggestion (for

RE: [Rd] (PR#7254)Documentation: Reference Index (.pdf) -- setOld Class

2004-10-01 Thread Pfaff, Bernhard
Dear Prof. Ripley, thks for your response and no, I have not looked at 2.0.0 beta version; blame on me, but I could not spare time by now for doing so; instead I looked at: http://cran.r-project.org/doc/manuals/fullrefman.pdf as well as on my local one in: 'rw1091/doc/manual/refman.pdf'. Hopeful

Re: [Rd] Problem with _new_ if class "lm" in object representatio n.

2004-10-01 Thread John Chambers
Eric Lecoutre wrote: > > Hi, > > For your lm problem, you could begin to create a virtual class that accepts > either NULL or a lm object: > > > setClassUnion("mylm", c("NULL","lm")) > > setClass("foo", representation(x="numeric", y="mylm")) > > new("foo",x=1,y=NULL) > > Then, you have to ch

[Rd] RCC compatibility patch

2004-10-01 Thread John Garvin
Would you consider the following patch to eval.c to allow compatibility with RCC? (It's in the applyClosure function.) @@ -432,6 +432,14 @@ SEXP f, a, tmp; RCNTXT cntxt; +#ifdef RCC +SEXP comp; +PROTECT(comp = getAttrib(op, install("RCC_CompiledSymbol"))); +if (comp != R_