Re: [R] Updating to nlme 3.1-62 failing from source (OS X)

2005-08-01 Thread Prof Brian Ripley
This is an internal compiler/assembler error - please check that your tools are consistent and if so report this as a MacOS X error. (It is the sort of thing which happens if gcc targetted for one assembler is used with another, for example on Solaris with GNU vs Sun assemblers.) Given that

[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

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 design

[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] 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: R-help@stat.math.ethz.ch Sent: Monday, August 01, 2005 12:21 PM Subject: [R] fortran

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.

[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

[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!

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

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.01

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 variable.

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: I

[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

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)] - TRUE

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,

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

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

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

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 2005,

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 bogdan romocea

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'

[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

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 mistake but

[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

[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

[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

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

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] 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:

[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] 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 really

[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,

[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 probably the

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?