Re: [R] can we manage memory usage to increase speed?

2005-08-01 Thread Prof Brian Ripley
On Mon, 1 Aug 2005, Haibo Huang wrote: > Please refer to the following post. Which is about Windows only, not Linux. (And on Windows, the answer given is on the help page for memory.size. together with a better one.) > > Ed > > --- Mike Lawrence <[EMAIL PROTECTED]> wrote: > >> Date: Mon, 1 Au

Re: [R] Rgdal windows binary warning message

2005-08-01 Thread Prof Brian Ripley
This is a consequence of how VC++ compiled gdal. R is protecting itself against that. Since R did not crash, there is nothing to worry about. On Tue, 2 Aug 2005, Tony Gill wrote: > Hi all, > > I just downloaded windows binaries of RGDAL (from > http://www.stats.ox.ac.uk/pub/RWin/bin/windows/c

Re: [R] linux compile options (64-bit)

2005-08-01 Thread Prof Brian Ripley
We find it equally strange that you posted this! The advice _is_ in the R-admin manual which the INSTALL file asks you to read if you have any questions. It covers using enhanced BLAS libraries. R builds out of the box on FC3, FC4 and Suse on AMD64. I use Goto's BLAS, but ATLAS can be used (e

Re: [R] can we manage memory usage to increase speed?

2005-08-01 Thread Haibo Huang
Please refer to the following post. Ed --- Mike Lawrence <[EMAIL PROTECTED]> wrote: > Date: Mon, 1 Aug 2005 00:19:06 -0300 > From: Mike Lawrence <[EMAIL PROTECTED]> > To: "Briggs, Meredith M" > <[EMAIL PROTECTED]> > CC: r-help@stat.math.ethz.ch > Subject: Re: [R] How do you increase memeory? >

[R] linux compile options (64-bit)

2005-08-01 Thread Jacob Michaelson
Sorry--sent this before I added a subject line... Begin forwarded message: > From: Jacob Michaelson <[EMAIL PROTECTED]> > Date: August 1, 2005 8:51:58 PM MDT > To: r-help@stat.math.ethz.ch > > > I've been using R for a while under Mac OS X, which thanks to the R > on OS X developers, is probabl

[R] (no subject)

2005-08-01 Thread Jacob Michaelson
I've been using R for a while under Mac OS X, which thanks to the R on OS X developers, is probably the best platform for learning R. I recently built a Linux box with a Pentium D processor, and am running an AMD64 port of Ubuntu with the SMP kernel. After setting up the basics on the box,

Re: [R] Can R do regression with absolute error loss?

2005-08-01 Thread Gabor Grothendieck
Check out the quantreg package. On 8/1/05, Chunyu Yang {msbbb062} <[EMAIL PROTECTED]> wrote: > Dear R users, > > > > I wonder if there are some functions or packages in R that can perform > regression under absolute error loss. I have searched the full manual of > R, but can not find it. I real

[R] Can R do regression with absolute error loss?

2005-08-01 Thread Chunyu Yang {msbbb062}
Dear R users, I wonder if there are some functions or packages in R that can perform regression under absolute error loss. I have searched the full manual of R, but can not find it. I really appreciate your help. Best regards, Ben Yang [[alternative HTML version deleted]] ___

Re: [R] NA when using read.csv

2005-08-01 Thread Wuming Gong
Hi Haibo, Are there any character "#" in the cells? By default, the things at the right side of the character "#" will be masked as comments. Wuming On 8/2/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > On 8/1/05, Haibo Huang <[EMAIL PROTECTED]> wrote: > > Hi, when I used: > > > > Dist=read

Re: [R] nnet package for OS X

2005-08-01 Thread Sean O'Riordain
Hi Lisa, as far as I know, nnet is part of the VR bundle which is a recommended package... refer http://www.stats.bris.ac.uk/R/src/contrib/Descriptions/VR.html so just try library(nnet) and see what happens! cheers! Sean On 02/08/05, Lisa Schweitzer <[EMAIL PROTECTED]> wrote: > Hello every one-

Re: [R] nnet package for OS X

2005-08-01 Thread Lisa Schweitzer
Hello every one-- My sincere apologies if I have missed something obvious, but I need to construct some multinomial logit models. I remember using nnet for this in the past, but I am having trouble locating a OS X version of this (it's been a year or so since I had to use it). I found a

[R] can we manage memory usage to increase speed?

2005-08-01 Thread Zhilin Liu
Hi, Thanks for reading. I am running a process in R for microarray data analysis. RedHat Enterprise Linux 4, dual AMD CPU, 6G memory. However, the R process use only a total of <200M memory. And the CPU usage is total to ~110% for two. The program takes at least 2 weeks to run at the curren

[R] Rgdal windows binary warning message

2005-08-01 Thread Tony Gill
Hi all, I just downloaded windows binaries of RGDAL (from http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.1/) and installed. I ran the example R_HOME\library\rgdal\R-ex\getPixmapGDAL.R and everything seemed to work as expected. However I got the following warning message: DLL attempt

[R] canonical correlations on multiple sets

2005-08-01 Thread Tsjerk Wassenaar
Hi! I was wondering whether anybody could help me with canonical correlation on multiple sets (3 or 4) of variates, e.g. using one of the methods suggested by Kettenring (Biometrika 58(3):433, 1971), starting from covariance or correlation matrices. If someone could point me to an efficient way t

Re: [R] NA when using read.csv

2005-08-01 Thread Gabor Grothendieck
On 8/1/05, Haibo Huang <[EMAIL PROTECTED]> wrote: > Hi, when I used: > > Dist=read.csv("test.csv",header=TRUE) > > to read data from CSV file. For some cells, R > mistakenly put in as NA, while most of the cells still Its not likely that there are errors in this software so its probably not a mi

[R] NA when using read.csv

2005-08-01 Thread Haibo Huang
Hi, when I used: Dist=read.csv("test.csv",header=TRUE) to read data from CSV file. For some cells, R mistakenly put in as NA, while most of the cells still appears to be right, and there is no error message in R. I am pretty sure the csv file is all right, but just can't figure out what went wron

Re: [R] New functions supporting GIF file format in R

2005-08-01 Thread Tuszynski, Jaroslaw W.
Sorry about late response but I was traveling and did not follow R-Help for 2 weeks. I think, I prefer to leave read/write.GIF functions in my package where I have a full control of IO formats. I like to store images as integer or real matrices, and movies as 3D arrays. I believe that 'pixmap' pac

Re: [R] How to hiding code for a package

2005-08-01 Thread Gary Wong
--- Spencer Graves <[EMAIL PROTECTED]> wrote: > Dear Gary Wong: > > Are you trying to defeat the GNU License? If so, > why not make it > available as an S-Plus library? There, you don't > have to provide the > source; binary files will suffice. > > spencer graves > > bogd

Re: [R] Compiling nlme 3.1-62 succeeds if OS X Developer Tools updated; however, compiling mgcv_1.3-4 fails

2005-08-01 Thread Michael Kubovy
Prof. Ripley was right. The problem went away when I updated the Developer Tools. On Aug 1, 2005, at 2:48 AM, Prof Brian Ripley wrote: > This is an internal compiler/assembler error - please check that > your tools are consistent and if so report this as a MacOS X error. > On Sun, 31 Jul 200

Re: [R] converting stata's by syntax to R

2005-08-01 Thread Thomas Lumley
On Mon, 1 Aug 2005, Chris Wallace wrote: > I am struggling with migrating some stata code to R. I have a data > frame containing, sometimes, repeat observations (rows) of the same > family. I want to keep only one observation per family, selecting > that observation according to some other varia

Re: [R] How to hiding code for a package

2005-08-01 Thread Spencer Graves
Dear Gary Wong: Are you trying to defeat the GNU License? If so, why not make it available as an S-Plus library? There, you don't have to provide the source; binary files will suffice. spencer graves bogdan romocea wrote: > There's something else you could try - since y

Re: [R] How to hiding code for a package

2005-08-01 Thread bogdan romocea
There's something else you could try - since you can't hide the code, obfuscate it. Hide the real thing in a large pile of useless, complicated, awfully formatted code that would stop anyone except the most desperate (including yourself, after a couple of weeks/months) from trying to understand it.

Re: [R] Conditional piece-wise dependent regression

2005-08-01 Thread Huntsinger, Reid
In case you want to estimate U and V as well, there's a segmented regression package (called "segmented") on http://cran.r-project.org for this. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Dimitris Rizopoulos Sent: Monday, August 01, 2

Re: [R] converting stata's by syntax to R

2005-08-01 Thread Dimitris Rizopoulos
if you also need to create the `keep' vector, then you could try this approach: fam <- c(1,2,3,3,4,4,4) wt <- c(1,1,0.6,0.4,0.4,0.4,0.2) dat <- data.frame(fam, wt) ### keep <- unlist( lapply(split(wt, fam), function(x){ ind <- rep(FALSE, length(x)) ind[which.max(x)] <- TRU

[R] Clayton-Oakes failure time model

2005-08-01 Thread S. Pohl
Dear R-list, does an R-package exist in which the Clayton-Oakes failure time model is implemented? Thanks for help, Stefan. [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] converting stata's by syntax to R

2005-08-01 Thread Jean Eid
Here is one way this can be done do.call("rbind", by(dat, list(dat$fam) ,function(x) { + if(NROW(x)>1) return(x[which.max(x$wt),]) + else return(x)} + )) and it returns fam wt keep 1 1 1.01 2 2 1.01 3 3 0.61 4 4 0.41 hth, On Mon, 1 Aug 2005, Chris Wallace wrote:

Re: [R] converting stata's by syntax to R

2005-08-01 Thread Peter Dalgaard
Chris Wallace <[EMAIL PROTECTED]> writes: > I am struggling with migrating some stata code to R. I have a data > frame containing, sometimes, repeat observations (rows) of the same > family. I want to keep only one observation per family, selecting > that observation according to some other vari

Re: [R] converting stata's by syntax to R

2005-08-01 Thread ronggui
try > attach(dat) > dat<-dat[order(fam,wt),] #sort the data ,as the stata's byable command does > lis<-by(dat,fam,function(x) x[length(x$fam),]) #equall your stata command ,but return a list. > do.call(rbind,lis) #to make the list to be a matrix-like result. fam wt keep 1 1 1.01 2 2 1.0

Re: [R] reading from posgresql

2005-08-01 Thread Jeff Gentry
> I don't know how to read data from posgresql to R. The RdbiPgSQL package will read from Postgres. Its current incarnation is available via Bioconductor (www.bioconductor.org). __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/list

[R] reading from posgresql

2005-08-01 Thread secretario academico FACEA
Dear all I don't know how to read data from posgresql to R. Does anybody can help me? Thanks a lot Adrián __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/pos

[R] converting stata's by syntax to R

2005-08-01 Thread Chris Wallace
I am struggling with migrating some stata code to R. I have a data frame containing, sometimes, repeat observations (rows) of the same family. I want to keep only one observation per family, selecting that observation according to some other variable. An example data frame is: # construct examp

Re: [R] How to hiding code for a package

2005-08-01 Thread Tuszynski, Jaroslaw W.
I always thought that ability to see and study the code of every package, was a great thing about R and other "open source" environments. So I hope there are no good ways of hiding the code of packages. Jarek \=== Jarek Tuszynski, PhD.

Re: [R] fortran

2005-08-01 Thread Ales Ziberna
See ?dyn.load ?.Fortran and manual "Writing R Extensions"! Then ask again if you have additional questions! Hope this helps, Ales Ziberna - Original Message - From: "javad Ashjari" <[EMAIL PROTECTED]> To: Sent: Monday, August 01, 2005 12:21 PM Subject: [R] fortran > Dear Friends > I

[R] fortran

2005-08-01 Thread javad Ashjari
Dear Friends I am trying to call a fortran subroutine in R, but i can not. How can i use the subroutines in R? Regrads Javad __ [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch

Re: [R] Conditional piece-wise dependent regression

2005-08-01 Thread Dimitris Rizopoulos
Since you want least squares, I think you could use lm() here, i.e., U <- 50 V <- 100 Time <- 1:150 dat <- data.frame(y = rnorm(150), Time, f1 = as.numeric(Time > U), f2 = as.numeric(Time > V)) ### m <- lm(y ~ Time + I(Time - U):f1 + I(Time - V):f2, data = dat) # check also the desig

[R] Conditional piece-wise dependent regression

2005-08-01 Thread Arie
Hi, after reading some R docs, I couldn’t figure out how can I find the solution for the following problem, therefore I would ask this friendly list for an advice. We’re making a least square approximation for an experiment described by the following model: T is the time, Y is some measured val