Re: [R] 64 bit R slower than 32 bit R on Sun Sparc Solaris?

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Roger D. Peng wrote: > Are you using an optimized BLAS for both builds? That's one > possibility. Also, 64-bit builds use up more memory initially since the > pointers are bigger. I've tried both 64-bit and 32-bit builds on > Sparc/Solaris and haven't seen any slowdown. It

Re: [R] C function name garbled

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Sigal Blay wrote: > On Tue, Sep 07, 2004 at 04:57:57PM -0400, Duncan Murdoch wrote: > > On Tue, 7 Sep 2004 12:46:55 -0700, S Blay <[EMAIL PROTECTED]> wrote : > > > > >I wrote an R wrapper function (phylpro) around a C function > > >(Rphylpro). > > >The first time I'm running

RE: [R] Thanks

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Liaw, Andy wrote: > A small drawback to Gabor's proposal: > > > "aa" %:% "dd" > [1] "b" "c" "d" > > "00" %:% "99" > character(0) > > Defining it as a seq() method for characters is probably not a terribly good > idea, as one should expect it to work for any character vectors

RE: [R] Thanks

2004-09-07 Thread Liaw, Andy
A small drawback to Gabor's proposal: (B (B> "aa" %:% "dd" (B[1] "b" "c" "d" (B> "00" %:% "99" (Bcharacter(0) (B (BDefining it as a seq() method for characters is probably not a terribly good (Bidea, as one should expect it to work for any character vectors as input. (B(One can argue that

Re: [R] Thanks

2004-09-07 Thread Gabor Grothendieck
Gabor Grothendieck myway.com> writes: : : åæ capitalbio.com> writes: : : : Thanks a lot for your timely rely. : : I still wonder whether I can use "a":"d" instead of 1:4. : : I remember I fulfill it successfully according to the guidance of some : materials on R,but fail to find it now. : :

Re: [R] Thanks

2004-09-07 Thread Gabor Grothendieck
åæ capitalbio.com> writes: : Thanks a lot for your timely rely. : I still wonder whether I can use "a":"d" instead of 1:4. : I remember I fulfill it successfully according to the guidance of some materials on R,but fail to find it now. You could define your own seq function and operator: R> "%

Re: [R] download of 1.9.1

2004-09-07 Thread Kevin Wang
Hi, On Tue, 7 Sep 2004, Kathy DeRiemer wrote: > Hi there, > > I downloaded version 1.9.1 WinZip file, 9.367 bytes ( so I can update my R Uit doesn't sound right, 9 bytes > and use the vcd package) but cannot get the download from your Web site to > unzip and install. So silly! > >

[R] download of 1.9.1

2004-09-07 Thread Kathy DeRiemer
Hi there, I downloaded version 1.9.1 WinZip file, 9.367 bytes ( so I can update my R and use the vcd package) but cannot get the download from your Web site to unzip and install. So silly! The error message is :" Error reading header after processing 0 entries." I tried downloading from 3 dif

Re: [R] 64 bit R slower than 32 bit R on Sun Sparc Solaris?

2004-09-07 Thread Roger D. Peng
Are you using an optimized BLAS for both builds? That's one possibility. Also, 64-bit builds use up more memory initially since the pointers are bigger. I've tried both 64-bit and 32-bit builds on Sparc/Solaris and haven't seen any slowdown. -roger Jason Liao wrote: Hello, everyone! I guess n

[R] 64 bit R slower than 32 bit R on Sun Sparc Solaris?

2004-09-07 Thread Jason Liao
Hello, everyone! I guess no one is still using R on Sun Sparc these days. But our department has a (pretty new) two-CPU Sun server. We recently compiled R as a 64 bit application and expected a performance boost. But it runs 25-30% slower than the 32 bit version of R. Anyone knows why this is so? T

Re: [R] C function name garbled

2004-09-07 Thread Sigal Blay
On Tue, Sep 07, 2004 at 04:57:57PM -0400, Duncan Murdoch wrote: > On Tue, 7 Sep 2004 12:46:55 -0700, S Blay <[EMAIL PROTECTED]> wrote : > > >I wrote an R wrapper function (phylpro) around a C function > >(Rphylpro). > >The first time I'm running my function, it runs with no errors. > >The second

[R] Thanks

2004-09-07 Thread 孟欣
Thanks a lot for your timely rely. I still wonder whether I can use "a":"d" instead of 1:4. I remember I fulfill it successfully according to the guidance of some materials on R,but fail to find it now. Thanks again! __ [EMAIL PROTECTED] mailing list ht

Re: [R] a little question about R

2004-09-07 Thread Kevin Wang
Hi, On Wed, 8 Sep 2004, [gb2312] ÃÏÐÀ wrote: > Hello,sir: Here's a little question about R which needs your help.Thanks in advance. > If I wanna make a sequence just like a,b,c,d (In other words,a vector consists of 4 > characters :a,b,c,d ).How can I do it in a shortcut manner? Yes,I can do it

Re: [R] a little question about R

2004-09-07 Thread Alec Stephenson
letters[1:4] Alec Stephenson Department of Statistics Macquarie University NSW 2109, Australia >>> " $(AÃÏÐÀ" <[EMAIL PROTECTED]> 09/08/04 10:41am >>> Hello,sir: Here's a little question about R which needs your help.Thanks in advance. If I wanna

RE: [R] a little question about R

2004-09-07 Thread Wiener, Matthew
letters[1;4] LETTERS[1:4] for capitals. Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Tuesday, September 07, 2004 8:42 PM To: [EMAIL PROTECTED] Subject: [R] a little question about R Hello,sir: Her

Re: [R] a little question about R

2004-09-07 Thread Roger D. Peng
Try letters[1:4] -roger åæ wrote: Hello,sir: Here's a little question about R which needs your help.Thanks in advance. If I wanna make a sequence just like a,b,c,d (In other words,a vector consists of 4 characters :a,b,c,d ).How can I do it in a shortcut manner? Yes,I can do it as following: c("a"

RE: [R] a little question about R

2004-09-07 Thread John Fox
You can use letters[1:4]. (B (BI hope that this helps, (B John (B (B> -Original Message- (B> From: [EMAIL PROTECTED] (B> [mailto:[EMAIL PROTECTED] On Behalf Of ?? (B> Sent: Tuesday, September 07, 2004 7:42 PM (B> To: [EMAIL PROTECTED] (B> Subject: [R] a little question about R 

Re: [R] pairwise comparisons

2004-09-07 Thread Roger D. Peng
Are you talking about repeated rows? Does unique() do what you want? -roger Rebecca Young wrote: Hello, I am a new R user. I am trying to calculate vector correlations for all pairwise comparisons in my data frame without repeats. I am familiar with the expand.grid function, but this includes re

[R] a little question about R

2004-09-07 Thread 孟欣
Hello,sir: Here's a little question about R which needs your help.Thanks in advance. If I wanna make a sequence just like a,b,c,d (In other words,a vector consists of 4 characters :a,b,c,d ).How can I do it in a shortcut manner? Yes,I can do it as following: c("a","b","c","d") and the result is:

[R] pairwise comparisons

2004-09-07 Thread Rebecca Young
Hello, I am a new R user. I am trying to calculate vector correlations for all pairwise comparisons in my data frame without repeats. I am familiar with the expand.grid function, but this includes repeats. Is there a way to use expand.grid and eliminate repeats? Or is there another function th

RE: [R] heatmap help

2004-09-07 Thread Liaw, Andy
Sorry. I think you need the argument distfun=as.dist. E.g., heatmap(x, distfun=as.dist, symm=TRUE, scale="none") Andy > From: Liaw, Andy > > I believe you want to pass the (symmetric) matrix as is, > rather than wrapped > in as.dist(). E.g., > > x <- as.matrix(dist(matrix(rnorm(100), 20, 5

RE: [R] heatmap help

2004-09-07 Thread Liaw, Andy
I believe you want to pass the (symmetric) matrix as is, rather than wrapped in as.dist(). E.g., x <- as.matrix(dist(matrix(rnorm(100), 20, 5))) heatmap(x, symm=TRUE, scale="none") HTH, Andy > From: Paul Lepp > > Dear R wizards, > Hopeful someone can help me with what I believe is a > p

[R] heatmap help

2004-09-07 Thread Paul Lepp
Dear R wizards, Hopeful someone can help me with what I believe is a pretty simple task. I pretty new to R so some (much) of the obvious escapes me. How do I get a distance matrix into heatmap? What do I tell distfun if what I'm trying to map is already an ordered distance matrix? I trie

Re: [R] Multiple comparisons in a non parametric case

2004-09-07 Thread Spencer Graves
Great summary, Rolf. Just one minor issue that recently bit me: In a data mining application with hundred of p-values, people want to make subtle distinctions based on extremely small p-values. In such applications, even a modest amount of skewness (to say nothing of outliers) migh

Re: [R] Multiple comparisons in a non parametric case

2004-09-07 Thread Rolf Turner
It looks to me like what you are doing is trying to judge significance of differences by non-overlap of single-sample confidence intervals. While this is appealing, it's not quite right. I just looked into my copy of Applied Nonparametric Statistics (second ed.) by Wayne W. Daniel (Duxbury, 1990)

Re: [R] using text on the x axis ticks rather than numbers

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Adaikalavan Ramasamy wrote: > Out of curiosity, why is srt not respected in mtext ? See examples below It doesn't make a lot of sense. mtext is about putting text on >m It would be nice to have 'srt' as an argument in mtext in future. You haven't been listening to the discus

Re: [R] using text on the x axis ticks rather than numbers

2004-09-07 Thread Adaikalavan Ramasamy
Out of curiosity, why is srt not respected in mtext ? See examples below plot(1:10, xaxt="n"); par(srt=45) mtext( paste("Point", LETTERS[1:10]), side=1, at=1:10 ) plot(1:10, xaxt="n"); mtext( paste("Point", LETTERS[1:10]), side=1, at=1:10, srt=45 ) It would be nice to have 'srt' as an argumen

Re: [R] C function name garbled

2004-09-07 Thread Duncan Murdoch
On Tue, 7 Sep 2004 12:46:55 -0700, S Blay <[EMAIL PROTECTED]> wrote : >I wrote an R wrapper function (phylpro) around a C function >(Rphylpro). >The first time I'm running my function, it runs with no errors. >The second time I'm trying to run it, I get an error message >with the first argument

[R] Multiple comparisons in a non parametric case

2004-09-07 Thread Marco Chiarandini
Dear all, I am conducting a full factorial analysis. I have one factor consisting in algorithms, which I consider my treatments, and another factor made of the problems I want to solve. For each problem I obtain a response variable which is stochastic. I replicate the measure of this response valu

[R] C function name garbled

2004-09-07 Thread S Blay
I wrote an R wrapper function (phylpro) around a C function (Rphylpro). The first time I'm running my function, it runs with no errors. The second time I'm trying to run it, I get an error message with the first argument to .Call garbled. Set up: > dyn.load("Phylpro.so") > source("phylpro.R") >

Re: [R] using text on the x axis ticks rather than numbers

2004-09-07 Thread Marc Schwartz
On Tue, 2004-09-07 at 12:52, Rajarshi Guha wrote: > Hello, > is there a way in which I can use text labels rather than numbers on > the x axis ticks? I basically have a vector of (say) 8 points and I want > to plot these sequentially. Thus the x axis would have ticks at 1 .. 8. > Rather than havi

Re: [R] using text on the x axis ticks rather than numbers

2004-09-07 Thread Eryk Wolski
Hi! first supress axis drawing with param axes=FALSE in your plot function. The use axis function to add your axes. ?axis /E On Tue, 7 Sep 2004, Rajarshi Guha wrote: > Hello, > is there a way in which I can use text labels rather than numbers on > the x axis ticks? I basically have a vector of

Re: [R] Cox regression for prevalence estimates

2004-09-07 Thread Thomas Lumley
robust standard errors are available from coxph with the option robust=TRUE. -thomas On 6 xxx -1, The Michaelson Institute wrote: > Hello, I'm an MD working in an eye clinic. I'm learning by myself to use R > for use in my research works and for implementation in a software project. > T

[R] using text on the x axis ticks rather than numbers

2004-09-07 Thread Rajarshi Guha
Hello, is there a way in which I can use text labels rather than numbers on the x axis ticks? I basically have a vector of (say) 8 points and I want to plot these sequentially. Thus the x axis would have ticks at 1 .. 8. Rather than having the labels 1 .. 8 I would like to have some arbitrary tex

Re: [R] gridBase and heatmap

2004-09-07 Thread Sean Davis
Kevin, Thanks for the advice. However, I was pretty careful to include a par(new=T) [see below]. Sean On Sep 7, 2004, at 1:14 PM, Kevin Bartz wrote: The problem is that most base plotting functions first wipe the graphics device clean. To do what you want, you need to use par(new = T) libe

RE: [R] gridBase and heatmap

2004-09-07 Thread Kevin Bartz
The problem is that most base plotting functions first wipe the graphics device clean. To do what you want, you need to use par(new = T) liberally between plots. Does that work for you? Kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Sean Davis Sent

[R] contrast/test statements in polynomial mixed model

2004-09-07 Thread Dede Greenstein
Hello I am having some trouble figuring out how to write contrast/test syntax for the following mixed model regression y ~ group*age + group*age^2 random =~1|ID (i have centered the age variable & there are only 2 groups) the contrast i am trying to construct should compare the overall shape o

[R] gridBase and heatmap

2004-09-07 Thread Sean Davis
I would like to use gridBase to place four separate heatmaps (actually, a stripped-down heatmap.2 from ght gregmisc package that contains only the "image" part) into four different viewports. I can get the placement correct, but I keep 'losing' the previous plot. Any suggestions? Here is

Re: [R] rodbc windows doesn't find dsn

2004-09-07 Thread Prof Brian Ripley
On Sat, 4 Sep 2004, Jack Tanner wrote: > Under WinXP, I have a system DSN called foo. It has stored > username/password information. (I can click Configure in ODBC > Administrator and then Test Data Source connects to MySQL without me > re-entering data.) > > Using RODBC, I can connect to the

Re: [R] Further png() question

2004-09-07 Thread Laura Quinn
Thanks all! It appears I just had the par in the wrong place in my loop. Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax: +44 113 343 6716 mail: [EMAIL PROTECTED] On Tue, 7 Sep 2004, Petr Pikal wrote: > >

Re: [R] Further png() question

2004-09-07 Thread Petr Pikal
On 7 Sep 2004 at 15:37, Laura Quinn wrote: > Ok, I have reinstalled R-1.9.0 and this appears to have fixed the > problems I was having with png(). However, I have a further question > regarding png() > > Is it possible to pass a par() argument to the png() command? I am > wanting to produce 4 p

RE: [R] Further png() question

2004-09-07 Thread Liaw, Andy
Insert the par() call after png(). par() is device-specific. Andy > From: Laura Quinn > > Ok, I have reinstalled R-1.9.0 and this appears to have fixed the > problems I was having with png(). However, I have a further question > regarding png() > > Is it possible to pass a par() argument to th

Re: [R] Further png() question

2004-09-07 Thread partha_bagchi
The following is not a problem in R 1.9.1: > png(file = "test.png") > par(mfrow = c(2,2)) > plot(c(1:10)) > plot(c(1:10)) > plot(c(1:10)) > plot(c(1:10)) > dev.off() null device 1 > Have you tried that? HTH, Partha Laura Quinn <[EMAIL PROTECTED]> Sent by: [EMAIL PROTECTED] 09/07

Re: [R] Further png() question

2004-09-07 Thread Uwe Ligges
Laura Quinn wrote: Ok, I have reinstalled R-1.9.0 and this appears to have fixed the problems I was having with png(). However, I have a further question regarding png() Is it possible to pass a par() argument to the png() command? I am wanting to produce 4 plots per object, which I normally acheiv

Re: [R] debugging an S4 method

2004-09-07 Thread John Chambers
"Bickel, David" wrote: > > Does anyone know how to use the equivalent of debug() on an S4 method? I would like > R to enter the browser not for the generic function, but for the method of the class > that I specify. This is what the trace() function does, if you supply the signature= argument s

Re: [R] Further png() question

2004-09-07 Thread Sean Davis
Are you calling png before your par command--I think you should. > png (file='filename') > par(mfrow=c(1,4)) > plot(1:10) > plot(1:10) > plot(1:10) > plot(1:10) > dev.off() See if that works for you. Sean On Sep 7, 2004, at 10:37 AM, Laura Quinn wrote: Ok, I have reinstalled R-1.9.0 and this appear

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Bernardo Rangel Tura
At 06:17 07/09/2004, you wrote: Hi guys. Check this out: NaN +NA [1] NaN NA + NaN [1] NA I thought "+" was commutative by definition. What's going on? In my version, both cases is NA: > NaN +NA [1] NA > NA + NaN [1] NA > R.version _ platform i386-pc-mingw32 arch i386 os min

[R] Further png() question

2004-09-07 Thread Laura Quinn
Ok, I have reinstalled R-1.9.0 and this appears to have fixed the problems I was having with png(). However, I have a further question regarding png() Is it possible to pass a par() argument to the png() command? I am wanting to produce 4 plots per object, which I normally acheive on an X window b

Re: [R] Cox regression for prevalence estimates

2004-09-07 Thread Bernardo Rangel Tura
At 09:27 07/09/2004, you wrote: Bernardo Rangel Tura wrote: At , The Michaelson Institute wrote: How can R be used to calculate the prevalence ratios using Cox regression + robust variance estimates ? Well, In Design package have a command: cph This command have a option "robsut" with default=FALSE

Re: [R] Error when running configure on xeon with r-1-9-patched

2004-09-07 Thread Gavin Simpson
Peter Dalgaard wrote: Gavin Simpson <[EMAIL PROTECTED]> writes: configure:25282: gcc -o conftest -D__NO_MATH_INLINES -mieee-fp -O2 -g -pipe -march=i386 -mcpu=i686 -I/usr/local/include -L/usr/local/lib conftest.c -lreadline -ldl -lncurses -lm -L/usr/local/lib -L/usr/lib/gcc-lib/i386-redhat-linux/

Re: [R] Cox regression for prevalence estimates

2004-09-07 Thread Frank E Harrell Jr
Bernardo Rangel Tura wrote: At , The Michaelson Institute wrote: How can R be used to calculate the prevalence ratios using Cox regression + robust variance estimates ? Well, In Design package have a command: cph This command have a option "robsut" with default=FALSE, but in help is write: " ...

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Prof Brian Ripley
On 7 Sep 2004, Peter Dalgaard wrote: > Robin Hankin <[EMAIL PROTECTED]> writes: > > > Dear Professor Ripley > > > > thank you for your reply. > > > > > > > >> > NaN +NA > > >> [1] NaN > > >> > NA + NaN > > >> [1] NA > > >> > > >> I thought "+" was commutative by definition. What's goin

Re: [R] Cox regression for prevalence estimates

2004-09-07 Thread Bernardo Rangel Tura
At , The Michaelson Institute wrote: How can R be used to calculate the prevalence ratios using Cox regression + robust variance estimates ? Well, In Design package have a command: cph This command have a option "robsut" with default=FALSE, but in help is write: " ... robust if TRUE a robust varian

Re: [R] Error when running configure on xeon with r-1-9-patched

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Gavin Simpson wrote: > Prof Brian Ripley wrote: > > Is libg2c.so in your library path? That's one common problem. > > Many thanks for your reply Prof. Ripley. How can I check what my library > path is? > > cat /etc/ld.so.conf gives > > include ld.so.conf.d/*.conf > /usr/X1

Re: [R] Error when running configure on xeon with r-1-9-patched

2004-09-07 Thread Peter Dalgaard
Gavin Simpson <[EMAIL PROTECTED]> writes: > configure:25282: gcc -o conftest -D__NO_MATH_INLINES -mieee-fp -O2 -g > -pipe -march=i386 -mcpu=i686 -I/usr/local/include -L/usr/local/lib > conftest.c -lreadline -ldl -lncurses -lm -L/usr/local/lib > -L/usr/lib/gcc-lib/i386-redhat-linux/3.3.3 > -L/usr

RE: [R] how to add error bar to the data in R?

2004-09-07 Thread Warnes, Gregory R
Alternatively, use the plotCI() function from the gregmisc package. -G -Original Message- From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Sent: 9/7/04 5:08 AM Subject: Re: [R] how to add error bar to the data in R? Hello, Use the function xYplot from hmisc package

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Peter Dalgaard
Robin Hankin <[EMAIL PROTECTED]> writes: > Dear Professor Ripley > > thank you for your reply. > > > > >> > NaN +NA > >> [1] NaN > >> > NA + NaN > >> [1] NA > >> > >> I thought "+" was commutative by definition. What's going on? > > > >It is clearly not under your compiler/OS. We could

Re: [R] Error when running configure on xeon with r-1-9-patched

2004-09-07 Thread Gavin Simpson
Prof Brian Ripley wrote: Is libg2c.so in your library path? That's one common problem. Many thanks for your reply Prof. Ripley. How can I check what my library path is? cat /etc/ld.so.conf gives include ld.so.conf.d/*.conf /usr/X11R6/lib /usr/lib/qt-3.3/lib /usr/local/lib and I don't have LD_LIB

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Robin Hankin
Dear Professor Ripley thank you for your reply. > NaN +NA [1] NaN > NA + NaN [1] NA I thought "+" was commutative by definition. What's going on? It is clearly not under your compiler/OS. We could add a configure test for broken systems and fix it in arithmetic.c but it hardly seems wort

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Jari Oksanen
On Tue, 2004-09-07 at 12:47, Prof Brian Ripley wrote: > On Tue, 7 Sep 2004, Robin Hankin wrote: > > > Check this out: > > I am unable to reproduce it on any of the 7 different systems I checked > (Solaris, Linux, Windows with various compilers). > > > > NaN +NA > > [1] NaN > > > NA + NaN > > [

Re: [R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Prof Brian Ripley
On Tue, 7 Sep 2004, Robin Hankin wrote: > Check this out: I am unable to reproduce it on any of the 7 different systems I checked (Solaris, Linux, Windows with various compilers). > > NaN +NA > [1] NaN > > NA + NaN > [1] NA > > I thought "+" was commutative by definition. What's going on?

[R] noncommutative addition: NA+NaN != NaN+NA

2004-09-07 Thread Robin Hankin
Hi guys. Check this out: NaN +NA [1] NaN NA + NaN [1] NA I thought "+" was commutative by definition. What's going on? R.version _ platform powerpc-apple-darwin6.8 arch powerpc os darwin6.8 system powerpc, darwin6.8 status major1 minor9.0 year 2004 month04

Re: [R] how to add error bar to the data in R?

2004-09-07 Thread Marta Rufino
Hello, Use the function xYplot from hmisc package cheers, Marta __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Error when running configure on xeon with r-1-9-patched

2004-09-07 Thread Martyn Plummer
On Mon, 2004-09-06 at 20:18, Gavin Simpson wrote: > Dear list, > > I'm having a problem configuring R-1-9-patched on a Intel Xeon machine > running Fedora Core 2. I used svn to checkout the source, altered > config.site to add some optimisation flags to match the Fedora rpm on > CRAN (the flags