[R] use of glm and cubic splines

2003-10-07 Thread Niel Hens
When I fit a model like glm(gb~bs(age,df=6)), I obtain 7 estimates. Can anybody explain me where they stand for? Thank you Niel Hens Center For Statistics Limburgs Universitair Centrum Universitaire Campus, Building D 3590 Diepenbeek Tel: +32-11-26 82 32 Fax: +32-11-26 82 99 Email:

Re: [R] Installing R in Linux 8.0

2003-10-07 Thread Martyn Plummer
On Tue, 2003-10-07 at 01:00, Peter Dalgaard BSA wrote: Javier Arsuaga [EMAIL PROTECTED] writes: I am trying to install R in Linux 8.0 and I downloaded R-1.7.1-1.i386.rpm and did rpm -hiv R-1.7.1-1.i386.rpm and I am getting the following message: warning: R-1.7.1-1.i386.rpm: V3 DSA

[R] Optimising code

2003-10-07 Thread Crispin Miller
Hi, Does anyone have any advice on speeding up R functions (short of re-implementing them in C :-) )? I have a function that applies a wilcoxon test to 12 sets of about a quarter of a million pairs (and takes about 3 hours). I've replaced the inner loop I had originally with a function call

Re: [R] monotone spline

2003-10-07 Thread alessandro ambrosi
Dear all, Sorry for the bad use of the list I did; I'm very young, only 4 days old to this list :-) Thanks a lot for your help and your hints: At 09.01 06/10/03 +0200, you wrote: mono.con(mgcv) Monotonicity constraints for a cubic regression spline. That is

[R] Installing R with all packages

2003-10-07 Thread Christian Hennig
Hi, I want R to be installed on a UNIX network (Solaris). I am not the system administrator and so I cannot do it myself. The system administrator wants to know which packages I want, and it may be a lot. Is there an easy way to download and install all packages at once? Is it a good idea?

Re: [R] Optimising code

2003-10-07 Thread Uwe Ligges
Crispin Miller wrote: Hi, Does anyone have any advice on speeding up R functions (short of re-implementing them in C :-) )? Some strategies are in every good book on S/R. I have a function that applies a wilcoxon test to 12 sets of about a quarter of a million pairs ... and let me guess:

[R] I need your help....

2003-10-07 Thread Yair Snir
Hello, I have a problem, I can't install the package 'mgu74av2cdf'. I downloaded the zip file, yet when asked the R console to install it from a zip file, I got the answer: Error in file(file, r) : unable to open connection In addition: Warning messages: 1: error -1 in

FW: [R] Optimising code

2003-10-07 Thread Crispin Miller
I have a function that applies a wilcoxon test to 12 sets of about a quarter of a million pairs ... and let me guess: everything is significiant to an almost arbitrary value of \alpha? :-) For each of quarter of a million sets, I do a wilcoxon between two pairs each containing twenty

Re: [R] I need your help....

2003-10-07 Thread Uwe Ligges
Yair Snir wrote: Hello, I have a problem, I can't install the package 'mgu74av2cdf'. I downloaded the zip file, yet when asked the R console to install it from a zip file, I got the answer: Error in file(file, r) : unable to open connection In addition: Warning messages: 1:

Re: FW: [R] Optimising code

2003-10-07 Thread Uwe Ligges
Crispin Miller wrote: I have a function that applies a wilcoxon test to 12 sets of about a quarter of a million pairs ... and let me guess: everything is significiant to an almost arbitrary value of \alpha? :-) For each of quarter of a million sets, I do a wilcoxon between two pairs each

Re: [R] I need your help....

2003-10-07 Thread Peter Dalgaard BSA
Uwe Ligges [EMAIL PROTECTED] writes: Yair Snir wrote: Hello, I have a problem, I can't install the package 'mgu74av2cdf'. I downloaded the zip file, yet when asked the R console to install it from a zip file, I got the answer: Error in file(file, r) : unable to open

[R] If I produce lots of figures at one time, how can I make R graphics keep all the figures already made?

2003-10-07 Thread xinyue ye
If I produce lots of figures at one time, how can I make R graphics keep all the figures already made? thanks, Xinyue Ye, Research Assistant Institute for Geo-spatial Research and Education 125 King Hall Eastern Michigan University Ypsilanti MI 48197 http://ceita.emich.edu/

Re: [R] If I produce lots of figures at one time, how can I make R graphics keep all the figures already made?

2003-10-07 Thread Philippe Glaziou
xinyue ye [EMAIL PROTECTED] wrote: If I produce lots of figures at one time, how can I make R graphics keep all the figures already made? thanks, Open as many graphics devices as needed. See the examples in: ?dev.list -- Philippe __ [EMAIL

RE: FW: [R] Optimising code

2003-10-07 Thread Crispin Miller
Yup, my mistake for not being very clear! Alas the computer is already very fast - looks like it's C for me... -Original Message- From: Patrick Burns [mailto:[EMAIL PROTECTED] Sent: 07 October 2003 12:00 Cc: Crispin Miller Subject: Re: FW: [R] Optimising code I had

Re: [R] If I produce lots of figures at one time, how can I make R graphics keep all the figures already made?

2003-10-07 Thread Thomas Petzoldt
xinyue ye schrieb: If I produce lots of figures at one time, how can I make R graphics keep all the figures already made? thanks, write it to a postscript() or pdf() file or to several bitmaps, see png() for examples. If you are on Windows, there is a Recording function in the menu, if the

Re: [R] Installing R with all packages

2003-10-07 Thread Uwe Ligges
Christian Hennig wrote: Hi, I want R to be installed on a UNIX network (Solaris). I am not the system administrator and so I cannot do it myself. The system administrator wants to know which packages I want, and it may be a lot. Is there an easy way to download and install all packages at

RE: FW: [R] Optimising code

2003-10-07 Thread Adaikalavan RAMASAMY
If you are only interested in calculating the wilcoxon statistic (as one would with calculating permutated p-values for example), the following should suffice. length.na - function(x, ...){ tmp - !(is.na(x) | is.infinite(x)) length(x[tmp], ...) } rank.na - function(x){ y - x[!is.na(x)];

[R] Beginner's query - segmentation fault

2003-10-07 Thread Laura Quinn
I am dealing with a huge matrix in R (20 columns, 54000 rows) and have lots of missing values within the dataset which are currently displayed as the value -999.00 I am trying to create a new matrix (or change the existing one) to display these values as NA so that I can then perform the necessary

Re: [R] If I produce lots of figures at one time, how can I make R graphics keep all the figures already made?

2003-10-07 Thread Henric Nilsson
At 13:31 2003-10-07 +0200, you wrote: If I produce lots of figures at one time, how can I make R graphics keep all the figures already made? thanks, If you are on Windows, there is a Recording function in the menu, if the graphics window is activated. Actually, the recording of plots can be

Re: [R] Installing R with all packages

2003-10-07 Thread Prof Brian Ripley
I once copied all the contrib package sources from CRAN to a local directory (I used an ftp client, but you could use rsync or wget : just don't do it recursively or you will get the Archive too) and ran R CMD INSTALL *.tar.gz The only overhead in installing all packages at once is disc space:

Re: [R] Installing R with all packages

2003-10-07 Thread A.J. Rossini
Christian Hennig [EMAIL PROTECTED] writes: I want R to be installed on a UNIX network (Solaris). I am not the system administrator and so I cannot do it myself. The system administrator wants to know which packages I want, and it may be a lot. Is there an easy way to download and install

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Peter Dalgaard BSA
Laura Quinn [EMAIL PROTECTED] writes: I am dealing with a huge matrix in R (20 columns, 54000 rows) and have lots of missing values within the dataset which are currently displayed as the value -999.00 I am trying to create a new matrix (or change the existing one) to display these values as

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2003, Laura Quinn wrote: I am dealing with a huge matrix in R (20 columns, 54000 rows) and have lots of missing values within the dataset which are currently displayed as the value -999.00 I am trying to create a new matrix (or change the existing one) to display these values as

Re: [R] If I produce lots of figures at one time, how can I make R graphics keep all the figures already made?

2003-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2003, Henric Nilsson wrote: At 13:31 2003-10-07 +0200, you wrote: If I produce lots of figures at one time, how can I make R graphics keep all the figures already made? thanks, If you are on Windows, there is a Recording function in the menu, if the graphics window is

RE: [R] Beginner's query - segmentation fault

2003-10-07 Thread Adaikalavan RAMASAMY
I cannot explain the segmentation fault but try this instead (which works for matrices) temp[which(temp==-999, arr.ind=T)] - NA Are you sure temp is matrix and not a dataframe ? Use class(temp) to find out. Also, if you are getting these -999.00 because you have read files containing them, it

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Uwe Ligges
Laura Quinn wrote: I am dealing with a huge matrix in R (20 columns, 54000 rows) and have lots of missing values within the dataset which are currently displayed as the value -999.00 I am trying to create a new matrix (or change the existing one) to display these values as NA so that I can then

Re: [R] I need your help....

2003-10-07 Thread James MacDonald
This is a question for the Bioconductor listserv, not R-help. In addition, all current versions of affy will automagically download and install any cdfenvs that you may need, so there is no need to do this manually. If anything, you may need to update your version of affy. Jim James W.

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Laura Quinn
thanks, have used temp [temp==0]- NA and this seems to have worked, though it won't let me access individual columns (ie temp$t1 etc) to work on - is there any real advantage in using a matrix, or would i be better advised to deal with dataframes? (I have double checked and temp is currently a

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Uwe Ligges
Adaikalavan RAMASAMY wrote: I cannot explain the segmentation fault but try this instead (which works for matrices) temp[which(temp==-999, arr.ind=T)] - NA No! Please *do* use is.na()- !!! Uwe Ligges Are you sure temp is matrix and not a dataframe ? Use class(temp) to find out. Also, if you

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2003, Laura Quinn wrote: thanks, have used temp [temp==0]- NA and this seems to have worked, though it won't let me access individual columns (ie temp$t1 etc) to work on - is there any real advantage in using a matrix, or would i be better advised to deal with dataframes?

Re: [R] Beginner's query - segmentation fault

2003-10-07 Thread Uwe Ligges
Laura Quinn wrote: thanks, have used temp [temp==0]- NA Please use is.na(temp[temp==0]) - TRUE and this seems to have worked, though it won't let me access individual columns (ie temp$t1 etc) No! temp$t1 is a list element or column of a data.frame, but not a column of a matrix. *PLEASE*,

[R] is.na(x) - TRUE versus x - NA (was: Beginner's query - segmentation fault)

2003-10-07 Thread Duncan Murdoch
On Tue, 07 Oct 2003 14:37:26 +0200, Uwe Ligges [EMAIL PROTECTED] wrote : 2. Use is.na(x) - TRUE instead of x - NA: is.na(temp[temp[ ,t1] == -999.00, t1]) - TRUE I hadn't heard this advice before. The online help ?is.na gives this cryptic advice: Function 'is.na-' may provide a safer

[R] .First.lib doesn't appear to be running after calling library()

2003-10-07 Thread Crispin Miller
Hi - so I've dusted off the C bits of my brain and gotten a library written for my package... It passes R CMD check ok, and I've put a file called '.First.lib.R' in the pacakge's 'R' subdirectory. Its permissions are 644. It says: .First.lib - function(lib,pkg) {

Re: [R] Re: Use of the Foreign package to import Stata files

2003-10-07 Thread Alan Cobo-Lewis
Regarding an old version of foreign allegedly being able to read the file at http://www.ats.ucla.edu/stat/stata/examples/rwg/concord1.dta Thomas Lumley [EMAIL PROTECTED] responded to my call for help: I don't know why it used to work. The file begins with 'h' (0x68), which isn't in my list of

Re: [R] Adjusting for within-cluster correlation: robcov() in Design-package and 'ids' in survey-package

2003-10-07 Thread Thomas Lumley
On Tue, 7 Oct 2003, Bernd Weiss wrote: Dear all, I would like to know if it possible to use the the robcov()-command in the Design- package in order to obtain a robust variance-estimate that adjusts for within-cluster correlation. Does the ids-option in the survey-package the same job?

[R] Re: is.na(x) - TRUE versus x - NA

2003-10-07 Thread Uwe Ligges
Duncan Murdoch wrote: On Tue, 07 Oct 2003 14:37:26 +0200, Uwe Ligges [EMAIL PROTECTED] wrote : 2. Use is.na(x) - TRUE instead of x - NA: is.na(temp[temp[ ,t1] == -999.00, t1]) - TRUE I hadn't heard this advice before. The online help ?is.na gives this cryptic advice: Function 'is.na-'

RE: [R] .First.lib doesn't appear to be running after calling lib rary()

2003-10-07 Thread Liaw, Andy
From the R-exts manual: The R subdirectory contains R code files. The code files to be installed must start with a (lower or upper case) letter and have one of the extensions .R, .S, .q, .r, or .s. We recommend using .R, as this extension seems to be not used by any other software. It should

Re: [R] .First.lib doesn't appear to be running after calling library()

2003-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2003, Crispin Miller wrote: Hi - so I've dusted off the C bits of my brain and gotten a library written for my package... It passes R CMD check ok, and I've put a file called '.First.lib.R' in the pacakge's 'R' subdirectory. Its permissions are 644. It says:

RE: [R] .First.lib doesn't appear to be running after calling lib rary()

2003-10-07 Thread Crispin Miller
Thanks - it is indeed the first '.' that's the problem... Crispin -Original Message- From: Liaw, Andy [mailto:[EMAIL PROTECTED] Sent: 07 October 2003 15:21 To: Crispin Miller Subject: RE: [R] .First.lib doesn't appear to be running after calling lib rary() I put .First.lib in

RE: [R] Problem getting an ifelse statment to work

2003-10-07 Thread Thomas W Blackwell
Greg - I am puzzled that the total counts in table(qs2) and table(qs9) could be different, if these are in fact two columns from the same data frame. I'm guessing that there are NAs in one or both columns, in addition to the digits 1,2,3,4, and that table() by default does not show them.

[R] Sorting matrix or data frame

2003-10-07 Thread Yao, Minghua
Dear all, Could anyone please tell me how to sort a matrix or a data frame against a column/row/component? Many thanks. -MY __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Sorting matrix or data frame

2003-10-07 Thread Uwe Ligges
Yao, Minghua wrote: Dear all, Could anyone please tell me how to sort a matrix or a data frame against a column/row/component? See ?order Uwe Ligges Many thanks. -MY __ [EMAIL PROTECTED] mailing list

[R] R-1.8.0 memory.limit()

2003-10-07 Thread James MacDonald
Using R-1.8.0 (d/l and compiled on 2003-10-01) on WinXP, I seem to be unable to determine the maximum memory allocated to R. The help still says to use memory.limit(size=NA), but this returns the value NA. In addition, I have set --max-mem-size=2G but I run out of memory somewhere around 500Mb

[R] multiple comparisons

2003-10-07 Thread Peter Adler
I'm having trouble finding an R equivalent to the S-Plus multicomp function, which does post-hoc comparisons of treatments means in ANOVAs. Am I missing something obvious? Thanks, Peter Peter Adler, PhD Dept. Ecology, Evolution and Marine Biology University of

Re: [R] Installing R with all packages

2003-10-07 Thread Paulo Justiniano Ribeiro Jr
Hi I use the script in the attached file to do this Cheers P.J. On Tue, 7 Oct 2003, Uwe Ligges wrote: Christian Hennig wrote: Hi, I want R to be installed on a UNIX network (Solaris). I am not the system administrator and so I cannot do it myself. The system administrator wants

Re: [R] R-1.8.0 memory.limit()

2003-10-07 Thread Uwe Ligges
James MacDonald wrote: Using R-1.8.0 (d/l and compiled on 2003-10-01) on WinXP, I seem to be unable to determine the maximum memory allocated to R. The help still says to use memory.limit(size=NA), but this returns the value NA. In addition, I have set --max-mem-size=2G but I run out of memory

Re: [R] multiple comparisons

2003-10-07 Thread Jonathan Baron
On 10/07/03 09:05, Peter Adler wrote: I'm having trouble finding an R equivalent to the S-Plus multicomp function, which does post-hoc comparisons of treatments means in ANOVAs. Am I missing something obvious? The package called multcomp? I don't know if it is the same. -- Jonathan Baron,

Re: [R] Installing R with all packages

2003-10-07 Thread Don MacQueen
If it's really true that you will be the only user of R, then it might be easier to build it in some place where you do have the necessary permissions. You would then follow the suggestion in the provided INSTALL file (from R 1.7.1): quote INSTALLATION You do not need to install R to

[R] plot width in Sweave

2003-10-07 Thread Christoph Lehmann
Hi I didn't find this in the manual: I need to change the width of a plot while I use sweave, so which command/parameters should I insert below, to change the width of a plot \begin{figure}[htbp] \begin{center} echo=TRUE, fig=TRUE= plot(Re(q),ylab =,type=o,col=blue,lwd=1, sub=mystring) @

Re: [R] plot width in Sweave

2003-10-07 Thread Achim Zeileis
On Tuesday 07 October 2003 19:19, Christoph Lehmann wrote: Hi I didn't find this in the manual: I need to change the width of a plot while I use sweave, so which command/parameters should I insert below, to change the width of a plot \begin{figure}[htbp] \begin{center} echo=TRUE,

Re: [R] multiple comparisons

2003-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2003, Jonathan Baron wrote: On 10/07/03 09:05, Peter Adler wrote: I'm having trouble finding an R equivalent to the S-Plus multicomp function, which does post-hoc comparisons of treatments means in ANOVAs. Am I missing something obvious? The package called multcomp? I

[R] Still Cannot Install rimage in R-1.7.1 (RH 9.0) Even With fftw Installed

2003-10-07 Thread Rick Bilonick
I'm still having problems installing rimage - the installation can't find the fftw headers. As suggested, I installed the fftw rpm (for RH 9 from freshrpms). It installed without any errors or warnings. Yet I get exactly the same error message - it can't find the fftw headers. What do I have

Re: [R] Still Cannot Install rimage in R-1.7.1 (RH 9.0) Even With fftw Installed

2003-10-07 Thread Peter Dalgaard BSA
Rick Bilonick [EMAIL PROTECTED] writes: I'm still having problems installing rimage - the installation can't find the fftw headers. As suggested, I installed the fftw rpm (for RH 9 from freshrpms). It installed without any errors or warnings. Yet I get exactly the same error message - it

[R] Compiling/using R on IBM AIX

2003-10-07 Thread Fan
Hi, I plan to use R on IBM AIX machines (4.x.x), and hoping some feedbacks from R users on this OS. I've looked at the R search site, my impression is that there would be very few peoples who are currently using R on AIX, correct me if I'm wrong. To compile R on AIX, I have two little questions:

Re: [R] Still Cannot Install rimage in R-1.7.1 (RH 9.0) Even With fftwInstalled

2003-10-07 Thread James MacDonald
Aren't the headers usually in the devel package (e.g., fftw-devel-2.1.5-0.dag.rh73.i386.rpm)? You might try installing the devel package too. HTH Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI

is.na(v)-b (was: Re: [R] Beginner's query - segmentation fault)

2003-10-07 Thread Richard A. O'Keefe
I am puzzled by the advice to use is.na(x) - TRUE instead of x - NA. ?NA says Function `is.na-' may provide a safer way to set missingness. It behaves differently for factors, for example. However, MAY provide is a bit scary, and it doesn't say WHAT the difference in behaviour is. I

[R] Saving workspace image

2003-10-07 Thread Ted Harding
Hi folks, On quitting R with q(), is it possible to save the workspace to a directory other than the one R was started from? (I sometimes have a project master directory with the major R code and data in that directory, but divisions of the project having their specific stuff in sub-directories.

Re: [R] Saving workspace image

2003-10-07 Thread Marc Schwartz
On Tue, 2003-10-07 at 19:00, [EMAIL PROTECTED] wrote: Hi folks, On quitting R with q(), is it possible to save the workspace to a directory other than the one R was started from? (I sometimes have a project master directory with the major R code and data in that directory, but divisions

[R] Getting alternative editor to stay as a foreground task

2003-10-07 Thread Andrew C. Ward
I apologise in advance if this question ought to be directed elsewhere. I'm trying to use the R sytnax highlighter for JEdit. This was mentioned in R-help a while ago and I'm hoping a few R users will have a (probably obvious) solution. (http://finzi.psych.upenn.edu/R/Rhelp02/archive/8812.html)