Re: [R] pca

2003-12-15 Thread Prof Brian Ripley
On Sun, 14 Dec 2003, Douglas Trainor wrote: Somewhere along the line, you have been confused. You're in good company though. Factor analysis and PCA are different entities entirely. Not in SPSS, where the same command is used for both (although we were not told anything like enough about

Re: [R] reverse lexicographic order

2003-12-15 Thread Prof Brian Ripley
On 14 Dec 2003, Peter Dalgaard wrote: Thomas W Blackwell [EMAIL PROTECTED] writes: Murray - If you could guarantee that all of the email addresses have exactly one occurrence of the @ character in them, then something like snip Otherwise, try something like this (I

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread A.J. Rossini
Did you make sure that you loaded the library that contains the C code on the remote hosts? best, -tony Faheem Mitha [EMAIL PROTECTED] writes: On Mon, 1 Dec 2003, A.J. Rossini wrote: Faheem Mitha [EMAIL PROTECTED] writes: So, can this (parallelization at the C level) be done without

Re: [R] pca

2003-12-15 Thread Chuck Cleland
Prof Brian Ripley wrote: On Sun, 14 Dec 2003, Douglas Trainor wrote: Somewhere along the line, you have been confused. You're in good company though. Factor analysis and PCA are different entities entirely. Not in SPSS, where the same command is used for both (although we were not told

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread Faheem Mitha
On Sun, 14 Dec 2003, A.J. Rossini wrote: Did you make sure that you loaded the library that contains the C code on the remote hosts? No, I only loaded the library on the local node (master). However, I'm not sure how I should do so. If I simply start up R on the remote hosts and load up

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread A.J. Rossini
Faheem Mitha [EMAIL PROTECTED] writes: On Sun, 14 Dec 2003, A.J. Rossini wrote: Did you make sure that you loaded the library that contains the C code on the remote hosts? No, I only loaded the library on the local node (master). However, I'm not sure how I should do so. If I simply

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread A.J. Rossini
Whoops. Exhaustion alert. The point is that R gets started up on the remote hosts, but you need to initialize the library on them. So you don't have to manually start up R on the remote hosts, but you do have to tell the remote sessions what to compute (including loading any libraries, etc).

[R] nls arguments

2003-12-15 Thread chenu
Hi all, I've got a problem with the nls function. I have an adjustment which works when I fix one of the argument of my function (Xo=150) : *Xo*=150 f- function (tt*,Xo*,a,b)ifelse(tt*Xo*,a*exp(-b**Xo*),a*exp(-b*tt))

[R] Missing Tcl.h in installation

2003-12-15 Thread Shin Tomoe
Dear Sir, I am a user of R, and I just downloaded rw1081.exe from your site. When I install the said file, it reports that a file is missing. I don’t know the cause. Its the only mising file because when I press ignore, everyting installs just fine. Anyways, Sir, can you send me the lacking

R: [R] nls arguments

2003-12-15 Thread Vito Muggeo
Dear chenu, I am not going to see your code with attention (also I do not understand the `*' symbol you used), however it looks a changepoint-type problem. The package segmented (on CRAN) uses a piecewise linear parameterization to fit regression models with breakpoints. Hope this helps, best,

[R] Week of the Year date conversion

2003-12-15 Thread Wayne Jones
Hello there fellow R-users, I have received some data which comes in the following format: example1-200301 The first 4 digits correspond to the year and the remaining 2 digits correspond to the week of the year. I have tried to convert this to a date by using strptime as follows:

Re: [R] Week of the Year date conversion

2003-12-15 Thread Prof Brian Ripley
On Mon, 15 Dec 2003, Wayne Jones wrote: Hello there fellow R-users, I have received some data which comes in the following format: example1-200301 The first 4 digits correspond to the year and the remaining 2 digits correspond to the week of the year. I have tried to convert this

RE: [R] density plot for very large dataset

2003-12-15 Thread Christophe Pallier
Have you tried the 'sm.density' function from the sm library? I used it for a dataset which 'only' had 13 points. I'm new to R and am trying to perform a simple, yet problematic task. I have two variables for which I would like to measure the correlation and plot versus each other.

Re: [R] help in lme

2003-12-15 Thread Spencer Graves
I can't answer all your questions right now, but I can answer the first: Have you considered gamma: n! = gamma(n+1) gamma(1+1:6) [1] 1 2 6 24 120 720 Regarding the other two, have you consulted Pinhiero and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? Bates

Re: [R] Week of the Year date conversion

2003-12-15 Thread Barry Rowlingson
Prof Brian Ripley wrote: That is not a complete date! Which day of the week is it? example1-200301 strptime(paste(example1, 1, format=%Y%U %d) will work, and you need to do something like that to resolve the ambiguity. Doesn't seem to work: - 6th day of 5th week: strptime(2003 05 06,

Re: [R] Week of the Year date conversion

2003-12-15 Thread Berwin Turlach
BR == Barry Rowlingson [EMAIL PROTECTED] writes: BR Ah ha. Use '%w' for 'day of week': strptime(2003 05 06, format=%Y %U %w) BR [1] 2003-02-02 BR I dont have a calendar to hand to check that 2 Feb is the 6th BR day of the 5th week It's not, it is the first day of the

Re: [R] Week of the Year date conversion

2003-12-15 Thread Prof Brian Ripley
On Mon, 15 Dec 2003, Barry Rowlingson wrote: Prof Brian Ripley wrote: That is not a complete date! Which day of the week is it? example1-200301 strptime(paste(example1, 1, format=%Y%U %d) will work, and you need to do something like that to resolve the ambiguity.

[R] Error with spdep

2003-12-15 Thread Perez Martin, Agustin
Dear useRs: First of all I would like to thank all the responses. I've an error with package spdep. I am working with a Windows XP machine (AMD-2000-XP RAM-256DDR) and 1.8.0. R-version and when I try to load spdep appear the following error: library(spdep) Error in loadNamespace(i, c(lib.loc,

Re: [R] Week of the Year date conversion

2003-12-15 Thread Peter Dalgaard
Berwin Turlach [EMAIL PROTECTED] writes: BR == Barry Rowlingson [EMAIL PROTECTED] writes: BR Ah ha. Use '%w' for 'day of week': strptime(2003 05 06, format=%Y %U %w) BR [1] 2003-02-02 BR I dont have a calendar to hand to check that 2 Feb is the 6th BR day of the

Re: [R] Error with spdep

2003-12-15 Thread Roger Bivand
On Mon, 15 Dec 2003, Perez Martin, Agustin wrote: Dear useRs: First of all I would like to thank all the responses. I've an error with package spdep. I am working with a Windows XP machine (AMD-2000-XP RAM-256DDR) and 1.8.0. R-version and when I try to load spdep appear the following

Re: [R] Error with spdep

2003-12-15 Thread Angel
spdep depends on package 'maptools', install it. - Original Message - From: Perez Martin, Agustin [EMAIL PROTECTED] To: lista R help (E-mail) [EMAIL PROTECTED] Sent: Monday, December 15, 2003 1:17 PM Subject: [R] Error with spdep Dear useRs: First of all I would like to thank all the

[R] Install Design Package

2003-12-15 Thread Marcela XX
Hi! I need to do a Logistic Regresion. I installed the Design Package. I tried lrm(... ) but the program said Object na.delete is not found. What can I do? Thanks!!! Marcela __ [EMAIL PROTECTED] mailing list

Re: [R] help in lme

2003-12-15 Thread Ramon Diaz-Uriarte
Since Spencer Graves already answered the factorial questions, I'll try to answer one of the other two: On Monday 15 December 2003 05:17, [EMAIL PROTECTED] wrote: To anyone who can help, Intelligent question (1) I keep on trying to fit a linear mixed model in R using 'lme(y~fxd.dsgn, data =

Re: [R] Install Design Package

2003-12-15 Thread Frank E Harrell Jr
On Mon, 15 Dec 2003 04:54:12 -0800 (PST) Marcela XX [EMAIL PROTECTED] wrote: Hi! I need to do a Logistic Regresion. I installed the Design Package. I tried lrm(... ) but the program said Object na.delete is not found. What can I do? Thanks!!! Marcela The package has a lot of

Re: [R] packages for ecologists - summary

2003-12-15 Thread Martin Wegmann
Dear R user, I asked for packages which should be mentioned in a ecological related presentation of R and would like to summarize the replies. First of all E. Paradis pointed out, that he would prefer to see R presented in a different way, please read his whole reply

[R] minor documentation typo ?

2003-12-15 Thread Rogers, James A [PGRD Groton]
Hi, I believe the following is a very minor typo in the documentation for merge() : all logical; all=L is shorthand for all.x=L and all.y=L ^^ ^ Looks like 'L' should be 'FALSE'. Someone please alert me if I am wrong. If I am not wrong, is

Re: [R] Basic question on function identical

2003-12-15 Thread Thomas Lumley
On Sat, 13 Dec 2003 [EMAIL PROTECTED] wrote: On 13-Dec-03 Martin Maechler wrote: In general, use == for testing equality of integer numbers (of type integer or not) I hope this is not a suggestion to avoid usage like which(x == max(x)) when x is a vector of

Re: [R] minor documentation typo ?

2003-12-15 Thread Martin Maechler
JARogers == Rogers, James A [PGRD Groton] Rogers on Mon, 15 Dec 2003 10:08:50 -0500 writes: JARogers I believe the following is a very minor typo in JARogers the documentation for merge() : JARogers all logical; all=L is shorthand for all.x=L and all.y=L JARogers

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread Na Li
On 15 Dec 2003, Faheem Mitha stated: So, the question would be how to get the R slaves to load up the library and I don't see any obvious way of doing this (using snow functions). Say you need boot on each slave, call this on the master, clusterEvalQ (cl, library (boot)) where cl is the

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread Faheem Mitha
On Mon, 15 Dec 2003, A.J. Rossini wrote: Yes, and no. Re-read the CPH-statcomp lab, and look at the bootstrap example, which solves the same problem. Look carefully -- it has to initialize the library on each node. If you are just loading the library manually, just do it on each node;

Re: [R] help with random numbers and Rmpi

2003-12-15 Thread Faheem Mitha
On Mon, 15 Dec 2003, Na Li wrote: On 15 Dec 2003, Faheem Mitha stated: So, the question would be how to get the R slaves to load up the library and I don't see any obvious way of doing this (using snow functions). Say you need boot on each slave, call this on the master,

Re: [R] Flury faces.

2003-12-15 Thread Peter Wolf
Eryk Wolski wrote: How to plot if possible flury faces in R? /Eryk __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help Some days ago I was looking for a function to represent multivariate data by faces, too. For

Re: [R] Basic question on function identical

2003-12-15 Thread Ted Harding
On 15-Dec-03 Thomas Lumley wrote: One reason that which.max() exists is that we cannot guarantee which(x==max(x)) to work. It is possible, though rather unlikely, for there to be no x such that x==max(x). One reason is the unpredictable use of 10-byte wide floating point registers on Intel

[R] distribution of second order statistic

2003-12-15 Thread Eugene Salinas (R)
Hi, I am getting some weird results here and I think I am missing something. I am trying to program a function that for a set of random variables drawn from uniform distributions plots that distribution of the second order statistic of the ordered variables. (ie I have n uniform distributions

Re: [R] Week of the Year date conversion

2003-12-15 Thread Gabor Grothendieck
If you don't want to rely on the week calculations in R (which in turn may depend on the OS?) you could try doing it yourself. The formula is pretty simple. Suppose: fw = day of week of first day of each week (0=Sun, 1=Mon, etc.) fy = day of week of first day of year (to be calculated given

Re: [R] distribution of second order statistic

2003-12-15 Thread Prof Brian Ripley
The order statistics have a beta distribution, so pbeta is all you need. On Mon, 15 Dec 2003, Eugene Salinas (R) wrote: I am getting some weird results here and I think I am missing something. I am trying to program a function that for a set of random variables drawn from uniform

[R] mailing list for basic questions

2003-12-15 Thread Martin Wegmann
Dear R-user, I started a few month ago with R and now I think that it is a great program but in the beginning I was close to give up, because I did not manage to read in my files, did not know how the indexing works, where my results are or did not understand how R saves my data and so

Re: [R] Basic question on function identical

2003-12-15 Thread Thomas Lumley
On Mon, 15 Dec 2003 [EMAIL PROTECTED] wrote: On 15-Dec-03 Thomas Lumley wrote: One reason that which.max() exists is that we cannot guarantee which(x==max(x)) to work. It is possible, though rather unlikely, for there to be no x such that x==max(x). One reason is the unpredictable use

[R] Appending intermediate terms in a For Loop

2003-12-15 Thread Jason . L . Higbee
R UseRs: I can't figure out how to append intermediate terms inside a For loop. My use of append(), various indexing, and use of data frames, vectors, matrices has been fruitless. Here's a simplified example of what I'm talking about: i - 1 for(i in 10) { v[i] - i/10 } v [1] 1 NA NA

Re: [R] Appending intermediate terms in a For Loop

2003-12-15 Thread Douglas Bates
[EMAIL PROTECTED] writes: I can't figure out how to append intermediate terms inside a For loop. My use of append(), various indexing, and use of data frames, vectors, matrices has been fruitless. Here's a simplified example of what I'm talking about: i - 1 for(i in 10) { v[i] -

Re: [R] Week of the Year date conversion

2003-12-15 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: On 15 Dec 2003, Peter Dalgaard wrote: Neither of those definitions coincide with ISO, BTW. ISO 8601 week of the year (as distinct from ISO C where these come from) is %V, of course. Doesn't work with RedHat 8 though... -- O__

Re: [R] Appending intermediate terms in a For Loop

2003-12-15 Thread Jason . L . Higbee
Thanks (to those who responded) for the help. Yeah, I forgot to make it a sequence. Must have overlooked that. Thanks again. -Jason [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

[R] Julian Dates

2003-12-15 Thread Ko-Kang Kevin Wang
Hi, I'm a bit confused how julian() works. If I understand right, it returns the number of days since the origin. I have a vector: SLDATX[1:10] [1] 1986-01-06 1986-01-17 1986-02-02 1986-02-04 [5] 1986-02-04 1986-02-21 1986-03-06 1986-03-25 [9] 1986-04-06 1986-04-10 And when I did:

Re: [R] cutree with agnes

2003-12-15 Thread Martin Maechler
[diverted from R-help to R-devel; please follow up on R-devel!] ChrisH == Christian Hennig [EMAIL PROTECTED] on Thu, 11 Dec 2003 15:13:52 +0100 (MET) writes: ChrisH Hi, this is rather a (presumed) bug report than a ChrisH question because I can solve my personal statistical

Re: [R] Basic question on function identical

2003-12-15 Thread Paul Johnson
I hope I am not telling you things you already know. If so, I apologize in advance. There are several C-library addons available to try to deal with the problem that comparisons of floating point numbers can be unpredictable. I think your example with the greater than sign would not be a

[R] Rd files Assignment functions.

2003-12-15 Thread wolski
Hi! Old story again. I put it away for a while because there are always other thinks to do. But I cant deny that I still like to like to comment assignment functions. Under linux I get. Cannot handle Rd file names containing ''. These are not legal file names on all R platforms. Please rename

Re: [R] Julian Dates

2003-12-15 Thread Uwe Ligges
Ko-Kang Kevin Wang wrote: Hi, I'm a bit confused how julian() works. If I understand right, it returns the number of days since the origin. I have a vector: SLDATX[1:10] [1] 1986-01-06 1986-01-17 1986-02-02 1986-02-04 [5] 1986-02-04 1986-02-21 1986-03-06 1986-03-25 [9]

Re: [R] reverse lexicographic order

2003-12-15 Thread Richard A. O'Keefe
Dr Murray Jorgensen ([EMAIL PROTECTED]) wrote: I have some email addresses that I would like to sort in reverse lexicographic order so that addresses from the same domain will be grouped together. How might that be done? Because he wants addresses from the same domain

Re: [R] Julian Dates

2003-12-15 Thread Gabor Grothendieck
What you can do to handle this timezone problem is either to use POSIXt with GMT or use chron (which does not use timezones so can't cause problems like this): Suppose: SLDATX - c( 1986-01-06, 1986-01-17, 1986-02-02, 1986-02-04, ,1986-02-04, 1986-02-21, 1986-03-06, 1986-03-25, ,1986-04-06,

Re: [R] Julian Dates

2003-12-15 Thread Ko-Kang Kevin Wang
Thanks! chron() is very useful indeed. Just out of interest, is it possible to do, say in this case, the number of months (or quarters) after January 1986? i.e. use a different time interval? On Mon, 15 Dec 2003, Gabor Grothendieck wrote: What you can do to handle this timezone problem is

Re: [R] reverse lexicographic order

2003-12-15 Thread Thomas Lumley
On Tue, 16 Dec 2003, Richard A. O'Keefe wrote: Here's a strrev() I knocked together quickly: strrev - function (s) paste(rev(strsplit(s, character(0))[[1]]), collapse=) If anyone can tell me how to vectorise this, I would be glad of the lesson. strrev- function(ss) {

Re: [R] Julian Dates

2003-12-15 Thread Gabor Grothendieck
Assuming you want to find a whole number months and quarters: require(chron) z - chron( SLDATX, format=y-m-d ) months - with( month.day.year(z), 12*(year-1986)+month-1 ) quarters - months %/% 3 --- Date: Tue, 16 Dec 2003 13:58:13 +1300 (NZDT) From: Ko-Kang Kevin Wang [EMAIL

[R] Winedit and R

2003-12-15 Thread Patrick Giraudoux
Hi all, I am trying to install add-on in R (rw1070) to work with WinEdit. Libraries Swinregistry and Rwinedt have been installed via Install Package(s) from local zip files from R. However, when I run library(Rwinedt) I get the following messages: library(Rwinedt) Loading required package:

Re: [R] reverse lexicographic order

2003-12-15 Thread Richard A. O'Keefe
I wrote: If anyone can tell me how to vectorise this, I would be glad of the lesson. where this was strrev - function (s) paste(rev(strsplit(s, character(0))[[1]]), collapse=) Thomas Lumley [EMAIL PROTECTED] suggested strrev- function(ss) { sapply(lapply(

[R] Resampling Stats software

2003-12-15 Thread Brandon Vaughn
Hi, I am new to R (I have most of my experience in SAS and SPSS). I was wondering if anyone has used both Resampling Stats and R, and could comment on strengths/relationships. Also, I have no clue on how to do the various examples from the book Resampling: The New Statistics in R. Can anyone

Re: [R] reverse lexicographic order

2003-12-15 Thread Gabor Grothendieck
Here is a way to do it without using apply. sep must be set to a character not in any of the strings. Below we show its much faster than using apply yet gives the same answer. strRev - function(x, sep = \10) { z - unlist( strsplit( paste( x, sep, sep= ), ) ) z - unlist(

Re: [R] Resampling Stats software

2003-12-15 Thread Jason Turner
Brandon Vaughn wrote: ... I am new to R (I have most of my experience in SAS and SPSS). I was wondering if anyone has used both Resampling Stats and R, and could comment on strengths/relationships. There are a few add-on packages for resampling with R. boot is the one I've used, and can