Re: [R] older versions of R (v1.0)

2004-07-15 Thread Uwe Ligges
Kathryn Jones wrote: Hi folks, Does anyone know where I can find downloads for old versions of R, around version 1.0? I found the downloads for unix but not windows. Any help would be greatly appreciated! Thank you, Kathryn Windows binary versions are not archived on CRAN. You may compile it

Re: [R] questions about R

2004-07-15 Thread Torsten Hothorn
On Thu, 15 Jul 2004, [big5] house-ball wrote: Hi, I try to compute critical value for multivariate normal distribution, and I find the crit(fit, const = c(0, 1), d = 1, cov = 0.95, rdf = 0) which seems to compute critical value. However, I can't compute right critical value for

[R] About lattice plots and loops

2004-07-15 Thread Jacques VESLOT
I cannot operate lattice plots, notably those in package gstat, within for() loops. (I can neither read such graphics included in a .R (or .txt) file with the source() function). However same commands executed one by one don't fail. Could you please help me find a means to execute lattice plots

Re: [R] About lattice plots and loops

2004-07-15 Thread Uwe Ligges
Jacques VESLOT wrote: I cannot operate lattice plots, notably those in package gstat, within for() loops. (I can neither read such graphics included in a .R (or .txt) file with the source() function). However same commands executed one by one don't fail. Could you please help me find a means to

Re: [R] Permutations

2004-07-15 Thread Christian T. Steigies
On Wed, Jul 14, 2004 at 07:42:49PM +0200, Jordi Altirriba Gutiérrez wrote: I’ve 12 elements in blocks of 3 elements and I want only to make permutations inter-blocks (no intra-blocks) (sorry if the terminology is not accurate) I am not a mathematician, but this sounds to me a little bit

RE: [R] About lattice plots and loops

2004-07-15 Thread Jacques VESLOT
It works indeed. Thanks a lot. Jacques VESLOT CIRAD -Message d'origine- De : Uwe Ligges [mailto:[EMAIL PROTECTED] Envoyé : jeudi 15 juillet 2004 11:35 À : [EMAIL PROTECTED] Cc : [EMAIL PROTECTED] Objet : Re: [R] About lattice plots and loops Jacques VESLOT wrote: I cannot operate

[R] (Newbie) scope of with()

2004-07-15 Thread Antonio Prioglio
Hi list, As I understand statements within with() are local to what is enclosed within its expression. As some excellent examples given to me previously have illustrated it is nevertheless possible to assign the evaluation of an expression to an external variable like x - with(data, if(..))

RE: [R] Permutations

2004-07-15 Thread Robin Hankin
hi again what a stimulating R discussion! This is R-help at its very best! I think I understand why you don't want pure inter-block permutations. My solution would be to realize that weeding out forbidden permutations is quite difficult and time-consuming (also as several people have pointed out

RE: [R] Permutations

2004-07-15 Thread Jordi Altirriba Gutiérrez
Dear R users, First of all, I want to thank the algorithms , time and suggestions to Rolf, Robert, Marc, Gabor, Adaikalavan, Cliff, Robin, Erich and Fernando. I want to sum up a little bit all the e-mails, the algorithms and the results of a test for 1000 permutations (in my last e-mail is

[R] areg.boot use of inverseTrans and ytype

2004-07-15 Thread Anne
Hi R helpers! I'm still a bit ( alot) confused by the use of inverseTrans and ytype in areg.boot (Hmisc): What I want to do seems very simple, but I do not get the result I want: plot the predicted values in the original scale. (I did not understand the documentation, sorry!) for instance the

[R] formatting tables with long column names via package:xtable within Sweave

2004-07-15 Thread john.gavin
Hi, I use the excellent Sweave tools for writing documents but was wondering how to neatly print a data.frame with long column headings. I cant manage to do this via package:xtable. Typically the labels that I would like use for each column consist of more than one word, but even with just one

[R] Newbie

2004-07-15 Thread Paolo Covelli
Hi, I wish build a R-script (or a R-function) that read a number from the keyboard and then process it. For example: from R I load the function X, that ask me the level of confidence \alpha, by keyboard I write 5 and the function go on. Thanks in advance Paolo [[alternative HTML

Re: [R] Newbie

2004-07-15 Thread Uwe Ligges
Paolo Covelli wrote: Hi, I wish build a R-script (or a R-function) that read a number from the keyboard and then process it. For example: from R I load the function X, that ask me the level of confidence \alpha, by keyboard I write 5 and the function go on. scan() Uwe Ligges Thanks in

Re: [R] formatting tables with long column names via package:xtable within Sweave

2004-07-15 Thread Jonathan Baron
On 07/15/04 11:07, [EMAIL PROTECTED] wrote: Hi, I use the excellent Sweave tools for writing documents but was wondering how to neatly print a data.frame with long column headings. I cant manage to do this via package:xtable. I think a problem with xtable is this:

Re: [R] Newbie

2004-07-15 Thread Adaikalavan Ramasamy
library(sm) x - ask(Enter a number) Enter a number: 5 x [1] 5 On Thu, 2004-07-15 at 11:41, Paolo Covelli wrote: Hi, I wish build a R-script (or a R-function) that read a number from the keyboard and then process it. For example: from R I load the function X, that ask me the level of

[R] indexing a parameter in nls

2004-07-15 Thread John Hendrickx
I'm trying to estimate a diagonal reference model in nls. In its basic form, the model consists of two factors with equal categories and a dependent variable. The model fits the main effects of the two factors such that they are proportional to each other. So the model I want to fit is:

Re: [R] (Newbie) scope of with()

2004-07-15 Thread Gabor Grothendieck
You could copy them back out at the end of your with like this: data(iris) # fetch test data with(iris, { Sepal.Length[Sepal.Length 5] - 5 Sepal.Width - Sepal.Width + 10 for(n in names(iris)) iris[n] - get(n) } ) Antonio Prioglio a.prioglio at city.ac.uk writes: :

Re: [R] (Newbie) scope of with()

2004-07-15 Thread Deepayan Sarkar
On Thursday 15 July 2004 03:01, Antonio Prioglio wrote: Hi list, As I understand statements within with() are local to what is enclosed within its expression. As some excellent examples given to me previously have illustrated it is nevertheless possible to assign the evaluation of an

Re: [R] Permutations

2004-07-15 Thread Ingmar Visser
On 7/14/04 7:42 PM, Jordi Altirriba Gutiérrez [EMAIL PROTECTED] wrote: 4 5 6 | 2 1 3 | 7 8 9 | 10 11 12 NO because it's an intra-block permutation of permutation 3 - - 10 1 7 | 4 8 7 | 5 6 12 | 3 2 9YES---Xth permutation 1 10 7 | 4 8 7 | 5 6 12 | 3 2 9NO because

[R] contriburl argument to install.packages

2004-07-15 Thread John Fox
Dear list members, I can't figure out how to specify the contriburl argument to install.packages() properly when the packages to be installed are in a directory on my local machine. I have in mind a command something like install.packages(missing.packages, contriburl=directory,

Re: [R] Newbie

2004-07-15 Thread Don MacQueen
Also, readline() -Don At 12:41 PM +0200 7/15/04, Paolo Covelli wrote: Hi, I wish build a R-script (or a R-function) that read a number from the keyboard and then process it. For example: from R I load the function X, that ask me the level of confidence \alpha, by keyboard I write 5 and the

Re: [R] contriburl argument to install.packages

2004-07-15 Thread Roger D. Peng
Would this work for you? install.packages(pathtopackagezipfile, CRAN = NULL) Although, you can only do one at a time, I think. -roger John Fox wrote: Dear list members, I can't figure out how to specify the contriburl argument to install.packages() properly when the packages to be installed are in

[R] GHK simulator

2004-07-15 Thread Trenkler, Dietrich
Dear R-community, not to re-invent the wheel I wonder if someone of you has ever written a function to compute the GHK smooth recursive simulator to estimate multivariate normal probabilities. See for instance page 194 of @BOOK{Greene97, author = {William H. Greene}, year = 1997, title =

RE: [R] contriburl argument to install.packages

2004-07-15 Thread John Fox
Dear Roger, -Original Message- From: Roger D. Peng [mailto:[EMAIL PROTECTED] Sent: Thursday, July 15, 2004 9:30 AM To: John Fox Cc: [EMAIL PROTECTED] Subject: Re: [R] contriburl argument to install.packages Would this work for you? install.packages(pathtopackagezipfile, CRAN

RE: [R] contriburl argument to install.packages

2004-07-15 Thread John Fox
Dear Duncan, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Duncan Murdoch Sent: Thursday, July 15, 2004 9:56 AM To: John Fox Cc: [EMAIL PROTECTED] Subject: Re: [R] contriburl argument to install.packages On Thu, 15 Jul 2004 10:09:46 -0400,

Re: [R] contriburl argument to install.packages

2004-07-15 Thread Duncan Murdoch
On Thu, 15 Jul 2004 11:05:43 -0400, John Fox [EMAIL PROTECTED] wrote : CRAN.packages(contriburl=file:c:/temp) Error in file(file, r) : unable to open connection In addition: Warning message: cannot open file `c:/temp/PACKAGES' I *do* have a directory

[R] formula

2004-07-15 Thread solares
Hi, i 'dont understand how to take a general formula, view this: x-1:5 y-c(0,1,1.7,2,2.1.4) dummy-data.frame(x=x,y=y) formula-y~A*log(x)/log(2) formu-as.formula(formula) fm-lm(formu,data=dummy) Error in eval(expr, envir, enclos) : Object A not found but A is the parameter of fitting, why is

[R] color scale to label a plot

2004-07-15 Thread Anne
Hello R-helpers I want to put a color scale in a plot: I've got an xy plot where the values of the response (z=f(x,y)) is symbolically given by colors (like heat or rainbow color scale) I would like to put such a scale with apprpriates labels in the plot, so as to facilitate the

[R] Formula

2004-07-15 Thread Vito Ricci
The '*' operator denotes factor crossing: 'a*b' interpreted as 'a+b+a:b' read ? formula. Best Vito Hi, i 'dont understand how to take a general formula, view this: x-1:5 y-c(0,1,1.7,2,2.1.4) dummy-data.frame(x=x,y=y) formula-y~A*log(x)/log(2) formu-as.formula(formula) fm-lm(formu,data=dummy)

Re: [R] formula

2004-07-15 Thread Luis Torgo
On Thu, 2004-07-15 at 16:28, [EMAIL PROTECTED] wrote: Hi, i 'dont understand how to take a general formula, view this: x-1:5 y-c(0,1,1.7,2,2.1.4) dummy-data.frame(x=x,y=y) formula-y~A*log(x)/log(2) formu-as.formula(formula) fm-lm(formu,data=dummy) You probably meant something similar to

Re: [R] formula

2004-07-15 Thread Adaikalavan Ramasamy
If you want to fit y = a + bx, then you use lm(y ~ x) instead of lm(y ~ A + bx). See the details section of help(formula). x - 1:5 y - c(0, 1.0, 1.7, 2.0, 2.1) lm(x ~ y) Call: lm(formula = x ~ y) Coefficients: (Intercept)y 0.6828 1.7038 PS : I think there is a typo

Re: [R] formula

2004-07-15 Thread Adaikalavan Ramasamy
If you want to fit y = a + bx, then you use lm(y ~ x) instead of lm(y ~ A + bx). See the details section of help(formula). x - 1:5 y - c(0, 1.0, 1.7, 2.0, 2.1) lm(x ~ y) Call: lm(formula = x ~ y) Coefficients: (Intercept)y 0.6828 1.7038 If A was already defined, and

Re: [R] formula

2004-07-15 Thread Adaikalavan Ramasamy
If you want to fit y = a + bx, then you use lm(y ~ x) instead of lm(y ~ A + bx). 'A' is not a parameter but coefficient and you do not need to specify coefficients, which is what the linear model is trying to do anyway ! See the details section of help(formula). x - 1:5 y - c(0, 1.0, 1.7,

Re: [R] color scale to label a plot

2004-07-15 Thread Anne
yes.. legend sould help. I still do not know how to paint my color sacle and put the corresp. labels correctly Thanks Anne - Original Message - From: Berton Gunter [EMAIL PROTECTED] To: Anne [EMAIL PROTECTED] Sent: Thursday, July 15, 2004 5:55 PM Subject: Re: [R] color scale to label a

Re: [R] color scale to label a plot

2004-07-15 Thread Martin Maechler
Anne == Anne Piotet [EMAIL PROTECTED] on Thu, 15 Jul 2004 17:50:17 +0200 writes: Anne Hello R-helpers Anne I want to put a color scale in a plot: Anne I've got an xy plot where the values of the response Anne (z=f(x,y)) is symbolically given by colors (like heat Anne

Re: [R] color scale to label a plot

2004-07-15 Thread Anne
Thanks! That is what I was looking for! Anne - Original Message - From: Martin Maechler [EMAIL PROTECTED] To: Anne [EMAIL PROTECTED] Cc: R list [EMAIL PROTECTED] Sent: Thursday, July 15, 2004 6:42 PM Subject: Re: [R] color scale to label a plot Anne == Anne Piotet [EMAIL

[R] formula and lm

2004-07-15 Thread solares
Hi, don' t understand why the function fomula have this error, i enclose the parameter a with the function I() Thank Ruben x-1:5 y-c( 2 ,4 , 6 , 8 ,11) formu-y~I(a*x) form-formula(formu) dummy-data.frame(x=x,y=y) fm-lm(form,data=dummy) Error in unique(c(AsIs, oldClass(x))) : Object a not found

RE: [R] color scale to label a plot

2004-07-15 Thread Charles Annis, P.E.
Anne: Perhaps this example will help. I have data (several dozen observations at each condition) for different R=stress.ratio. The stress.ratios are: stress.ratio [1] 0.1 0.1 0.1 0.2 0.3 0.5 0.6 0.7 0.8 0.9 0.9 0.9 0.9 0.9 0.9 You will notice multiple datasets for R=0.1 and R=0.9. I

Re: [R] contriburl argument to install.packages

2004-07-15 Thread John Fox
Dear Duncan, On Thu, 15 Jul 2004 11:11:53 -0400 Duncan Murdoch [EMAIL PROTECTED] wrote: On Thu, 15 Jul 2004 11:05:43 -0400, John Fox [EMAIL PROTECTED] wrote : CRAN.packages(contriburl=file:c:/temp) Error in file(file, r) : unable to open connection In addition: Warning

RE: [R] GHK simulator

2004-07-15 Thread Kevin Bartz
Yes. See the mvtnorm package on CRAN. It implements Genz's algorithm, the fastest method in town. (I have a version I converted to S-PLUS as well, if you, or anyone, might need it.) Kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Trenkler, Dietrich

[R] newbie

2004-07-15 Thread Paolo Covelli
Hi List, does it exist a sintax-comand to clear thre console (the screen)? I' ve finded nothing on the standard manuals. (apart from the manual comand CTRL+L) Thank in advance Paolo [[alternative HTML version deleted]] __ [EMAIL PROTECTED]

Re: [R] newbie

2004-07-15 Thread Rolf Turner
Paolo Covelli wrote: does it exist a sintax-comand to clear thre console (the screen)? I've finded nothing on the standard manuals. (apart from the manual comand CTRL+L) You don't say what operating system you're using . Anyhow, as far as I know there is no built-in to clear the screen.

[R] Where is global environment?

2004-07-15 Thread Mike Prager
I am using R 1.9.1 under Windows XP Professional. Is there any R function that returns the filesystem path from which the global environment (workspace) was loaded (and to which it will, by default, be saved)? This seems identical with the current directory, so in any given session it is not

[R] dput and dget

2004-07-15 Thread Mike Prager
I am seeking advice about dput() and dget(). We are using the ascii format supported by these functions as a way to write data from other programs that can be read into R easily. We are able to save complicated results (in the form of an R list) to a single file that can be read trivially into

Re: [R] Where is global environment?

2004-07-15 Thread Mike Prager
I am using R 1.9.1 under Windows XP Professional. Is there any R function that returns the filesystem path from which the global environment (workspace) was loaded (and to which it will, by default, be saved)? This seems identical with the current directory, so in any given session it is not

[R] row naming

2004-07-15 Thread Herman, David (NIH/NIMH)
Hello, If I have a four column data set (with thousands of rows), that doesn't have a header, how do I load in this text file, WITHOUT a row added for naming(i.e. numbering the rows, 1 2 3 4 5..). Also, if a row is added for naming, then will it be actually included in the data?

[R] how to upload [forwarded]

2004-07-15 Thread Federico Calboli
[ This was posted to the R-packages list (which I moderate), but definitely doesn't belong there. Martin Maechler ] --- start of forwarded message --- From: Federico Calboli [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: how to upload Date: 15 Jul 2004 18:30:26 +0100 Dear All, I

Re: [R] Where is global environment?

2004-07-15 Thread Duncan Murdoch
On Thu, 15 Jul 2004 14:15:28 -0400, Mike Prager [EMAIL PROTECTED] wrote : I am using R 1.9.1 under Windows XP Professional. Is there any R function that returns the filesystem path from which the global environment (workspace) was loaded (and to which it will, by default, be saved)? This

[R] More on global environment

2004-07-15 Thread Mike Prager
To follow up on my previous question, suppose a user R session wants to unload one workspace and load another within an R session. Is the following the correct sequence? 1. save.image() to save the current workspace as .Rdata in the current working directory. 2. rm(list=ls()) to remove

Re: [R] row naming

2004-07-15 Thread Mike Prager
At 7/15/2004 02:35 PM Thursday, you wrote: Hello, If I have a four column data set (with thousands of rows), that doesn't have a header, how do I load in this text file, WITHOUT a row added for naming(i.e. numbering the rows, 1 2 3 4 5..). Also, if a row is added for naming, then

Re: [R] More on global environment

2004-07-15 Thread Roger D. Peng
That seems reasonable, although you might use rm(list = ls(all.names = TRUE)) if you are interested in removing objects whose names begin with a period. -roger Mike Prager wrote: To follow up on my previous question, suppose a user R session wants to unload one workspace and load another within

RE: [R] More on global environment

2004-07-15 Thread Liaw, Andy
I've posted the following to R-help before. Hope it helps you. cd - function(dir = tclvalue(tkchooseDirectory()), saveOld=FALSE, loadNew=TRUE) { stopifnot(require(tcltk)) if (saveOld) save.image(compress=TRUE) setwd(dir) rm(list=ls(all=TRUE, envir=.GlobalEnv),

RE: [R] More on global environment

2004-07-15 Thread Mike Prager
At 7/15/2004 03:23 PM Thursday, Andy Liaw wrote: I've posted the following to R-help before. Hope it helps you. cd - function(dir = tclvalue(tkchooseDirectory()), saveOld=FALSE, loadNew=TRUE) { stopifnot(require(tcltk)) flush.console() if (saveOld)

Re: [R] formula and lm

2004-07-15 Thread Spencer Graves
It wants to compute a*x given a and x, and then use ordinarly least squares to estimate b0 and b1 in y = b0 + b1*I(a*x). If that is what you intend, you must supply a. If you want to estimate a, e.g., with no constant, use y~x-1. Does this answer the question? hope this helps. spencer

[R] adding option in the Windows installer for --internet2

2004-07-15 Thread Liaw, Andy
[Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets created. Also, is it possible to make that an option that can be set at

Re: [R] table lookup n R

2004-07-15 Thread Kjetil Halvorsen
?match %in% Kjetil Halvorsen Anne wrote: Hello R helpers! I looked but did not find a table-lookup R-utility. I could use a loop to do the job (old FORTRAN/C habits die hard) but if I have a big table in which I have to search for the values corresponding to a vector, I end up logically with a

[R] Eclipse plugin for R or perhaps S-plus.

2004-07-15 Thread Richard Piper
Does any one know of an eclipse (http://eclipse.org) plugin for R/S. thanks RIchard -- Richard Piper Web: http://icu.rnsh.org Mobile: 0438-120860 Sydney: 9926-8617 Kippax: 6550-5117 Email: [EMAIL PROTECTED] GPG key ID: F6DFFB33 @ wwwkeys.pgp.net Key fingerprint: 3FB8 8355 3AAA 84C3 D548 8FFF EFDB

Re: [R] adding option in the Windows installer for --internet2

2004-07-15 Thread Duncan Murdoch
On Thu, 15 Jul 2004 16:18:43 -0400, Liaw, Andy [EMAIL PROTECTED] wrote : [Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets

Re: [R] adding option in the Windows installer for --internet2

2004-07-15 Thread Adaikalavan Ramasamy
Yes, corporate firewall can be a real pain at times. Duncan does have a point that many users (especially new ones) will not know what this is for. Would it be possible to test during the installation process if --internet2 is suitable and ask for confirmation from the user. A message as follows

RE: [R] adding option in the Windows installer for --internet2

2004-07-15 Thread Liaw, Andy
From: Duncan Murdoch On Thu, 15 Jul 2004 16:18:43 -0400, Liaw, Andy [EMAIL PROTECTED] wrote : [Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the

Re: [R] contriburl argument to install.packages

2004-07-15 Thread Roger D. Peng
Hmm...is this maybe a bug? Check out this section of install.packages(): localcran - length(grep(^file:, contriburl)) 0 if (!localcran) { if (is.null(destdir)) { tmpd - tempfile(Rinstdir) if (!dir.create(tmpd)) stop(Unable to create temp

[R] Functions in a package not visible to a user

2004-07-15 Thread Rui Dantas
Dear all, Sorry if this is obvious, but I couldn't find the answer in any place I looked. I have a package with several R functions. Some are internal auxiliary functions, and should not be available to the user (nor do they, for example, need user documentation). How can I hide them? Best

RE: [R] Functions in a package not visible to a user

2004-07-15 Thread Liaw, Andy
That's one of the purposes of the namespace. See Prof. Tierney's article in the R Newsletter, or the slides of his keynote lecture at the useR! 2004 conference. It's also explained in the `Writing R Extensions' manual. Basically you only export functions that the users should see. The manual

RE: [R] contriburl argument to install.packages

2004-07-15 Thread John Fox
Dear Roger, After Duncan pointed out my error, I was able to get things to work fine. (The context was a discussion on the R-GUI list, in response to which I wrote a tcltk dialog to download missing packages or install them from a local directory. It was the latter that was giving me trouble --

[R] Does AIC() applied to a nls() object use the correct number of estimated parameters?

2004-07-15 Thread Peter.Caley
I'm wondering whether AIC scores extracted from nls() objects using AIC() are based on the correct number of estimated parameters. Using the example under nls() documentation: data( DNase ) DNase1 - DNase[ DNase$Run == 1, ] ## using a selfStart model fm1DNase1 - nls( density ~ SSlogis(

[R] SJava - restart function

2004-07-15 Thread Kathryn Jones
Hi folks, Does anyone know about using event handling in SJava?? I tried examples from here: http://www.omegahat.org/RSJava/examples/ but they are quite old I think and don't work in newer versions of R due to the restart() function being defunct. Does anyone know how to get around this

Re: [R] Does AIC() applied to a nls() object use the correct number of estimated parameters?

2004-07-15 Thread Adaikalavan Ramasamy
I do not know anything about nls(), so apologies if I get it completely wrong. help(AIC) says that AIC is defined to be -2*log-likelihood + k*npar; where k = 2 by default. I think you calculated -2*log-likelihood + k*(npar + 1) instead. Does this help ? On Fri, 2004-07-16 at 03:50, [EMAIL

RE: [R] Does AIC() applied to a nls() object use the correctnumber of estimated parameters?

2004-07-15 Thread Peter.Caley
Thanks Adaikalavan, however the problem remains. Considering AIC() as applied to the linear model in AIC() help documentation: data(swiss) lm1 - lm(Fertility ~ . , data = swiss) AIC(lm1) [1] 326.0716 Clearly this includes the estimation of the residual standard error as an estimated

Re: [R] Random intercept model with time-dependent covariates, results different from SAS [SUMMARY]

2004-07-15 Thread Keith Wong
Thank you Prof Ripley and Dr Bebber for the helpful responses to my post on 4 July 2004, and references to further reading. To close the thread, I summarize the answer to my question. The different results between SAS and R arose from more than one cause. 1. SAS incorrectly assumed that Group