Re: [R] Building package on Windows: No rule to make target '-llapack'

2004-06-10 Thread Prof Brian Ripley
On Wed, 9 Jun 2004, Zhu Wang wrote: I have a problem to build a package on Windows XP while there is no problem on Linux. The Makefile is something like: There would be a problem on Linux, if that Makefile were used. I suspect it is not used. ### LIBNAME=cts PKG_LIBS =

Re: [R] Re: R equivalent of Splus rowVars function

2004-06-10 Thread Diethelm Wuertz
Rmetrics availalble from www.rmetrics.org has implemented several R functions which deal with column statistics, row statistics and rolling analysis. Diethelm Wuertz David Brahm wrote: Mark Leeds [EMAIL PROTECTED] wrote (to S-News): does anyone know the R equivalent of the SPlus rowVars

[R] Failure to compile on Itanium

2004-06-10 Thread F. Tusell
I am trying to compile R on an Itanium machine running Red Hat 7.2 and gcc version 2.96. The build fails with the following symptoms: g77 -fPIC -g -O2 -c cmplx.f -o cmplx.lo g77 -fPIC -g -O2 -c cmplxblas.f -o cmplxblas.lo gcc -shared -L/usr/local/lib -o libRlapack.so dlapack0.lo dlapack1.lo

[R] question about hierclust {multiv}

2004-06-10 Thread Xinan Yang
my major is bioinformatics, and i'm trying to cluster ( agglomerate the closest pari of observations ) in R. i have already got my own similarities metric, but do not know how to clust it based on similarities instead of dissimilarities. since the help document of hierclust mentions the parameter

Re: [R] Help with a Lattice plot that fails with an empty unique combination

2004-06-10 Thread Martin Maechler
Deepayan == Deepayan Sarkar [EMAIL PROTECTED] on Wed, 9 Jun 2004 18:04:38 -0500 writes: Deepayan On Wednesday 09 June 2004 01:58, Tom Mulholland wrote: While using Lattice I received the following error. Error in if (xx != 0) xx/10 else z/10 : argument is of length zero

Re: [R] Failure to compile on Itanium

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, F. Tusell wrote: I am trying to compile R on an Itanium machine running Red Hat 7.2 and gcc version 2.96. The build fails with the following symptoms: g77 -fPIC -g -O2 -c cmplx.f -o cmplx.lo g77 -fPIC -g -O2 -c cmplxblas.f -o cmplxblas.lo gcc -shared

Re: [R] Getting Pr from Summary(lm)

2004-06-10 Thread Martin Maechler
David == David J Netherway [EMAIL PROTECTED] on Thu, 10 Jun 2004 11:19:05 +0930 writes: David Thanks for all the relies. David I recently discovered names and applied it to lm David objects but did not think to apply it to the David summary object. ok, and now you should

Re: [R] question about hierclust {multiv}

2004-06-10 Thread Wolski
Hi! I would convert the simmilarities into dissimilarities by myself. Its quite easy. For example to make a dissimilarity from correlation which is a simmilarity measures you can: dcorr = 1 - cor and dcorr will be a dissimilarity. More general. Make the smallest value the largest and the

[R] nls and R scoping rules

2004-06-10 Thread joerg van den hoff
I apologize for posting this in essence the second time (no light at the end of the tunnel yet..): is there a way to enforce that nls takes both, the data *and* the model definition from the parent environment? the following fragment shows the problem. # cut here== wrapper -

Re: [R] question about hierclust {multiv}

2004-06-10 Thread Martin Maechler
As I've already said on the R-devel list, MM why on earth are you using hierclust() from the ORPHANED package MM 'multiv', when there's hclust() in the core 'stats' package MM and 'agnes' in the recommended 'cluster' package ? {and your question is not about hierclust but about

Re: [R] nls and R scoping rules

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, joerg van den hoff wrote: I apologize for posting this in essence the second time (no light at the end of the tunnel yet..): is there a way to enforce that nls takes both, the data *and* the model definition from the parent environment? the following fragment shows

RE: [R] how to initialize random seed properly ?

2004-06-10 Thread Lutz Prechelt
I want to start R processes on multiple processors from single shell script and I want all of them to have different random seeds. The usual way would be using the process ID for the initialization. Or the shell script could just hand out consecutive numbers 1, 2, 3, ... Lutz

[R] Manova and specifying the model

2004-06-10 Thread Stefanie von Felten
Hi, I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory

Re: [R] nls and R scoping rules

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, Prof Brian Ripley wrote: Around R 1.2.x the notion was introduced that variables should be looked for in the environment of a formula. Functions using model.frame got converted to do that, but nls did not. I guess that the best way forward is to ensure that nls (and

[R] Table to Data Frame

2004-06-10 Thread Ko-Kang Kevin Wang
[Forwarding on behalf of a colleague] She's got a list with several tables: tab - list() for(i in 1:6) { + tab[[i]] - table(freq[i]) + } tab [[1]] 0 0.17 0.3 0.51 2.534 1962515252 [[2]] 0 0.17 0.31 2.534 199136

Re: [R] Manova and specifying the model

2004-06-10 Thread Petr Pikal
Hallo On 23 May 2004 at 13:44, Stefanie von Felten wrote: Hi, I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call.

Re: [R] nls and R scoping rules

2004-06-10 Thread joerg van den hoff
Prof Brian Ripley wrote: On Thu, 10 Jun 2004, Prof Brian Ripley wrote: Around R 1.2.x the notion was introduced that variables should be looked for in the environment of a formula. Functions using model.frame got converted to do that, but nls did not. I guess that the best way forward is

Re: [R] nls and R scoping rules

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, joerg van den hoff wrote: Prof Brian Ripley wrote: On Thu, 10 Jun 2004, Prof Brian Ripley wrote: Around R 1.2.x the notion was introduced that variables should be looked for in the environment of a formula. Functions using model.frame got converted to do

Re: [R] Table to Data Frame

2004-06-10 Thread Uwe Ligges
Ko-Kang Kevin Wang wrote: [Forwarding on behalf of a colleague] She's got a list with several tables: tab - list() for(i in 1:6) { Most easy solution: Make freq[i] a factor with all the levels that may appear. Then all tables have same dimension, e.g.: tab - lapply(freq, function(x) table(

Re: [R] Table to Data Frame

2004-06-10 Thread Wolski
Hi! First as Uwe Ligges explained and finally. do.call(rbind,tab) Sincerely Eryk *** REPLY SEPARATOR *** On 6/10/2004 at 10:17 PM Ko-Kang Kevin Wang wrote: [Forwarding on behalf of a colleague] She's got a list with several tables: tab - list() for(i in 1:6) { +

[R] Clustering Categorial and Continuous Variables

2004-06-10 Thread Wayne Jones
Hi there fellow R users, R has many different clustering packages (e.g. mclust,cluster,e1071). However, can anyone recommend a method to deal with data sets that contain categorial and continuous variables? Regards Wayne KSS Ltd Seventh Floor St James's Buildings 79 Oxford Street

Re: [R] Table to Data Frame

2004-06-10 Thread Adaikalavan Ramasamy
# Generate data universe - c(0, 0.17, 0.5, 1, 2.5, 3, 4) prob - c(94, 1, 1, 1, 1, 1, 1)/1000 tab - list(NULL) for(i in 1:6){ tab[[i]] - table(sample(universe, 218, prob=prob, rep=TRUE)) } # Re-table all - unique(unlist(sapply(tab, function(x) names(x retab - t(sapply( tab, function(x)

Re: [R] Clustering Categorial and Continuous Variables

2004-06-10 Thread Wolski
Hi! You need a apropriate dissimilarity measure. look for daisy in package cluster help(daisy,package=cluster) x: numeric matrix or data frame. Dissimilarities will be computed between the rows of 'x'. Columns of mode 'numeric' (i.e. all columns when 'x' is a matrix) will

[R] Help with plotmath

2004-06-10 Thread Steve Roberts
There must be a simple answer. I want to plot an expression, where the expression is held in a string variable. The obvious solution along the lines of ex-x^2 plot( c(0,1), c(0,1), main=as.expression(ex) ) gives me the a title x^2 - ie doesn't treat it like an expression. I suspect I don't

Re: [R] Help with plotmath

2004-06-10 Thread Prof Brian Ripley
You *parse* strings to form expressions. plot( c(0,1), c(0,1), main=parse(text=ex) ) On Thu, 10 Jun 2004, Steve Roberts wrote: There must be a simple answer. I want to plot an expression, where the expression is held in a string variable. The obvious solution along the lines of ex-x^2

Re: [R] moving data and output?

2004-06-10 Thread Martin Klaffenboeck
Am 09.06.2004 22:23:51 schrieb(en) Peter Alspach: 2. How can I make outputs of the grafics (plot, hist, ...) into a file? what platform (Windows or Unix)? Unix. (Gentoo Linux in special) Martin __ [EMAIL PROTECTED] mailing list

Re: [R] Help with plotmath

2004-06-10 Thread Petr Pikal
Hi On 10 Jun 2004 at 12:28, Steve Roberts wrote: There must be a simple answer. I want to plot an expression, where the expression is held in a string variable. The obvious solution along the lines of ex-x^2 plot( c(0,1), c(0,1), main=as.expression(ex) ) Plotmath example led me to

RE: [R] question related to S-Plus

2004-06-10 Thread Pikounis, Bill
Hi Rick, learning R. A question to this end, though. Would all of my existing S-Plus arrays, functions, and so on have to be re-created from scratch in R, or is there a way to copy them into the .RData directory? The answer to this question has major I am not sure if you have had this

[R] ordered probit or logit / recursive regression

2004-06-10 Thread VUILLEUMIER Mathieu
I make a study in health econometrics and have a categorical dependent variable (take value 1-5). I would like to fit an ordered probit or ordered logit but i didn't find a command or package who make that. Does anyone know if it's exists ? other question : i would like to make a least square

RE: [R] question related to S-Plus

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, Pikounis, Bill wrote: Hi Rick, learning R. A question to this end, though. Would all of my existing S-Plus arrays, functions, and so on have to be re-created from scratch in R, or is there a way to copy them into the .RData directory? The answer to this

Re: [R] ordered probit or logit / recursive regression

2004-06-10 Thread Prof Brian Ripley
On Thu, 10 Jun 2004, VUILLEUMIER Mathieu wrote: I make a study in health econometrics and have a categorical dependent variable (take value 1-5). I would like to fit an ordered probit or ordered logit but i didn't find a command or package who make that. Does anyone know if it's exists ?

Re: [R] nls and R scoping rules

2004-06-10 Thread Luke Tierney
On Thu, 10 Jun 2004, Prof Brian Ripley wrote: On Thu, 10 Jun 2004, Prof Brian Ripley wrote: Around R 1.2.x the notion was introduced that variables should be looked for in the environment of a formula. Functions using model.frame got converted to do that, but nls did not. I guess

Re: [R] how to initialize random seed properly ?

2004-06-10 Thread ryszard . czerminski
Here is my best solution so far using $RANDOM in bash or maybe somebody has pure R solution ? #!/bin/bash for i in 1 2; do p=tmp$i cat $p.R EOF set.seed($RANDOM) cat('rnorm(3) =', rnorm(3), '\n') EOF nice R CMD BATCH --no-save --no-restore $p.R $p.log done sleep 3 grep rnorm tmp[12].log Best

Re: [R] Help with plotmath

2004-06-10 Thread Steve Roberts
No that isn't it I'm afraid - my expression is pasted together from various bits and pieces and is at some stage a string variable. I needed parse as Brian pointed out Thanks anyway, Steve. From: Petr Pikal [EMAIL PROTECTED] To: Steve Roberts [EMAIL

Re: [R] lsoda with arbitrary zero thresholds (with psuedo-solution)

2004-06-10 Thread Setzer . Woodrow
Dear Hank, Last question first: really, only you can say for sure if 4e-281 and 5e-11 are small enough; it depends on the units you measure your state variables in. However, this strategy cannot get the state variables to exactly 0. Obviously, you could get closer to 0.0 faster by setting

Re: [R] how to initialize random seed properly ?

2004-06-10 Thread Itay Furman
perldoc -f srand has a quick and good advice on initializing a random sequence. I second the recommendation to prepare the list of random numbers in advance for sake of reproducibility. Itay -- [EMAIL PROTECTED]

Re: [R] how to initialize random seed properly ?

2004-06-10 Thread A.J. Rossini
[EMAIL PROTECTED] writes: Here is my best solution so far using $RANDOM in bash or maybe somebody has pure R solution ? Sure -- use the rsprng or rlecuyer packages from CRAN, which support (independent?) streams of random numbers from a common set of seeds. best, -tony -- [EMAIL PROTECTED]

[R] Response to questions raised in Mar 17 reply

2004-06-10 Thread Jeff Yanosky
Hi, I would like to repsond to a few questions raised in a reply to a question posted by another user on March 17, 2004. The entire message is copied at the end of this email. The relevant questions and statements are as follows: What did you not understand about help(memory.size)? This is

[R] Lattice::qqmath -- groups option question

2004-06-10 Thread richard . kittler
Does the 'groups' option on qqmath just color the points differently in the main distribution or does it actually overlay separate quantile plots for each subset? I would like to be able to do the latter. --Rich Richard Kittler AMD TDG 408-749-4099

Re: [R] Lattice::qqmath -- groups option question

2004-06-10 Thread Deepayan Sarkar
On Thursday 10 June 2004 09:21, [EMAIL PROTECTED] wrote: Does the 'groups' option on qqmath just color the points differently in the main distribution or does it actually overlay separate quantile plots for each subset? I would like to be able to do the latter. The 'groups' option doesn't

[R] Informal discussion group about R

2004-06-10 Thread Baize, Harold
I've started a tribe for discussing R and sharing scripts. Tribe.net is one of the popular on-line social communities, like Friendster. Visit and see if it is a forum that you find useful. To join the tribe you will need to register with Tribe.net. I hope it will be of help to newbies,

Re: [R] Building package on Windows: No rule to make target '-llapack'

2004-06-10 Thread Zhu Wang
On Thu, 2004-06-10 at 06:40, Prof Brian Ripley wrote: On Wed, 9 Jun 2004, Zhu Wang wrote: I have a problem to build a package on Windows XP while there is no problem on Linux. The Makefile is something like: There would be a problem on Linux, if that Makefile were used. I suspect it is

[R] Everything was ok, but I failed to get .o file in g77

2004-06-10 Thread Rui
Hi folks, I tried to compile a FORTRAN routine using the command “g77 –o2 –c all.f”, it seems everything is fine, there is no error message at all, but finally I can not get the file all.o. Therefore, I can not move on to the next step using “g77 –shared –c all.dll all.o” to get the file

RE: [R] Lattice::qqmath -- groups option question

2004-06-10 Thread richard . kittler
Thanks for the info. Given that I rule out qqmath would the best method be to make repeated calls to qqnorm without plotting and then overlay the results or is there a more elegant method? --Rich Richard Kittler AMD TDG 408-749-4099 -Original Message- From: Deepayan Sarkar

Re: [R] Everything was ok, but I failed to get .o file in g77

2004-06-10 Thread Prof Brian Ripley
I am afraid your message is in an unreadable character set and none was specified by your mailer. So I don't really know what you used, but I have ^V where - should be, and so on. What you would need is g77 -O2 -c all.f ^ note You appear to have asked for a file '2' here via -o2. g77

[R] odesolve: lsoda vs rk4

2004-06-10 Thread Chris Knight
I'm trying to use odesolve for integrating various series of coupled 1st order differential equations (derived from a system of enzymatic catalysis and copied below, apologies for the excessively long set of parameters). The thing that confuses me is that, whilst I can run the function rk4: out

[R] Questions about Preserving registers

2004-06-10 Thread Rui
Hi folks, I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file. However, when I used the command dyn.load(“test.dll”), I got the following message: NULL Warning message: DLL attempted to change FPU control word from 9001f to 90003 I read the instruction on Duncan Murdoch’s

Re: [R] Questions about Preserving registers

2004-06-10 Thread Uwe Ligges
Rui wrote: Hi folks, I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file. And what is the reason not to use the recommended gcc? Uwe Ligges However, when I used the command dyn.load(test.dll), I got the following message: NULL Warning message: DLL attempted to change FPU

[R] displaying a table in full vs. typing variable name

2004-06-10 Thread Rishi Ganti
I have a data frame called totaldata that is 10,000 rows by about 9 columns. The data frame contains many zero entries which are important. If I type totaldata it only prints out the first two columns. I expect (and want) it to print out all 9 columns. (I am actually sinking this to a text

Re: [R] Questions about Preserving registers

2004-06-10 Thread Thomas Lumley
On Thu, 10 Jun 2004, Rui wrote: Hi folks, I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file. However, when I used the command dyn.load(“test.dll”), I got the following message: NULL Warning message: DLL attempted to change FPU control word from 9001f to 90003 I read

Re: [R] displaying a table in full vs. typing variable name

2004-06-10 Thread Uwe Ligges
Rishi Ganti wrote: I have a data frame called totaldata that is 10,000 rows by about 9 columns. If about 9 equals 2, the behaviour reported below is expected. The data frame contains many zero entries which are important. If I type totaldata it only prints out the first two columns. I expect

Re: [R] displaying a table in full vs. typing variable name

2004-06-10 Thread Marc Schwartz
On Thu, 2004-06-10 at 11:26, Uwe Ligges wrote: Rishi Ganti wrote: I have a data frame called totaldata that is 10,000 rows by about 9 columns. If about 9 equals 2, the behaviour reported below is expected. That is, of course, for sufficiently large values of about... ;-) Marc

Re: [R] Lattice::qqmath -- groups option question

2004-06-10 Thread Deepayan Sarkar
On Thursday 10 June 2004 10:38, [EMAIL PROTECTED] wrote: Thanks for the info. Given that I rule out qqmath would the best method be to make repeated calls to qqnorm without plotting and then overlay the results or is there a more elegant method? I don't think there's any good way currently.

Re: [R] odesolve: lsoda vs rk4

2004-06-10 Thread Peter Dalgaard
Chris Knight [EMAIL PROTECTED] writes: I'm trying to use odesolve for integrating various series of coupled 1st order differential equations (derived from a system of enzymatic catalysis and copied below, apologies for the excessively long set of parameters). The thing that confuses me is

Re: [R] Questions about Preserving registers

2004-06-10 Thread Duncan Murdoch
On Thu, 10 Jun 2004 09:12:05 -0700 (PDT), Thomas Lumley [EMAIL PROTECTED] wrote : DLL attempted to change FPU control word from 9001f to 90003 I read the instruction on Duncan Murdoch’s website about preserving registers, but I still don’t understand it. For example, This code is for Delphi.

Re: [R] off topic publication question

2004-06-10 Thread Roger D. Peng
I always thought that the use of we was a reference to the author(s) and the reader. -roger Erin Hodgess wrote: Dear R People: Please excuse the off topic question, but I know that I'll get a good answer here. If a single author is writing a journal article, should she use We performed a test or

Re: [R] How to Describe R to Finance People

2004-06-10 Thread Ko-Kang Kevin Wang
Thank you very much to those who contributed to this rather interesting discussion/debate. I was very surprised (and almost overwhelmed) by the volume of replies based on this topic! I have prepared some slides and put the draft version on www.stat.auckland.ac.nz/~kwan022/tmp/RFin.ppt. Feel

[R] a scope problem

2004-06-10 Thread Rajarshi Guha
Hi, I have some code that looks like: dftc - df[sets$tcset,] pt - numeric(nrow(dftc)) sub - 1:nrow(dftc) for (i in 1:nrow(dftc)) { n - nnet( fmla, data=dftc, weights=wts, subset=sub[-i], size=4, decay=0.01) pt[i] - predict( n, dftc[ i, ], type='class' ) }

[R] tryCatch() and preventing interrupts in 'finally'

2004-06-10 Thread Henrik Bengtsson
With tryCatch() it is possible to catch interrupts with tryCatch(). Then you can use a 'finally' statement to clean up, release resources etc. However, how can I protect against additional interrupts? This is a concern when the hold down Ctrl+C and generates a sequence of interrupts. Example:

Re: [R] a scope problem

2004-06-10 Thread Tony Plate
This looks like it probably is a scope problem with non-standard evaluation rules for the argument subset= of nnet. Instead of subset=sub[-i], try data=dftc[-i,] (I've not tested this since I don't have the data objects you used.) hope this helps, Tony Plate At Thursday 04:38 PM 6/10/2004,

[R] running R UNIX in a mac computer

2004-06-10 Thread Tiago R Magalhaes
Hi to you all My question is: there is a package written in UNIX for which there is no Mac version. I would like to know if it's possible to install the R UNIX version on the MacOSX and run that UNIX package on my Mac (through this UNIX R Vresion on a Mac) I have seen a porfile for r version

Re: [R] running R UNIX in a mac computer

2004-06-10 Thread Ulises Mora Alvarez
Hi: Please tell us: which Mac OS are you using? If you are using any of the 10.x, then you should download R-1.9.0 at: http://cran.us.r-project.org/bin/macosx If you have anything older (MacOS 9 or lower), then you should visit: http://cran.us.r-project.org/bin/macos Regards... On Thu, 10 Jun

Re: [R] tryCatch() and preventing interrupts in 'finally'

2004-06-10 Thread Luke Tierney
More complete support for interrupt management is still on the to do list, unfortunately. onexit actions have similar issues (and have for a long time). What I was hoping to do is provide a user level mechanism for turning off interrupts around a piece of code and then optionally enabling them

Re: [R] running R UNIX in a mac computer

2004-06-10 Thread Tiago R Magalhaes
Hola ulisses thank you very much for replying to my posting I am using the last Mac Version: OS X 10.3 Panther I use R in the Mac and it works fine The thing is that I want to use a package that was made in UNIX and I wanted to know how easy it would be to run the UNIX version of R in the Mac

Re: [R] running R UNIX in a mac computer

2004-06-10 Thread Ulises Mora Alvarez
Hola! If you are already using Panther with R.dmg, then is pretty straight forward to use a Unix package (panther it's a kind of Unix). If the package is already on CRAN, all you have to do is to search the package (either from the command line or from the menus) and then install. If what