Re: [R] Rounding problem R vs Excel

2003-06-04 Thread Duncan Murdoch
On Tue, 3 Jun 2003 09:49:44 +0100, you wrote in message [EMAIL PROTECTED]: Duncan If the numbers are not represently exactly how does R resolve problems like the one below? Is there something that needs to be set up in the R environment like the number of significant figures? x-4.145*100+0.5 x

Re: [R] (no subject)

2003-06-04 Thread Roger Koenker
Although not immediately relevant to the present inquiry this might still be an opportune moment to mention again an item that is at the top of my R-wish-listr: The SparseM function slm (for sparse lm) is well suited for problems like this in which the design matrix is quite sparse, so it would

Re: Rtips (was Re: [R] ? building a database with a the great examples

2003-06-04 Thread Paul E. Johnson
I think there is a lot of merit in this idea. I think there is a big question about authentication and protection of Wikis from vandalism. I've set up Wikis for other projects that I started after Rtips. I have not seen your Wiki software before, but it looks pretty nice. I see it does

[R] S+ style implementation of GAM for R?

2003-06-04 Thread Douglas Beare
Hi, I've got the R library mgcv for GAM written by Simon Wood which works well in many instances. However, over the years I got attached to the S+ implementation of GAM which allows loess smoothing in more than 1 dimension as well as spline smoothing. Has anyone ported the S+ GAM library to R?

RE: [R] tseries adf.test

2003-06-04 Thread Pfaff, Bernhard
I have a question regarding the adf.test command in the tseries library. I have a vector of time series observations (2265 daily log prices for the OEX to be exact). I also have this same data in first-differenced form. I want to test both vectors individually for staionarity with an

Re: [R] S+ style implementation of GAM for R?

2003-06-04 Thread Simon Wood
I've got the R library mgcv for GAM written by Simon Wood which works well in many instances. However, over the years I got attached to the S+ implementation of GAM which allows loess smoothing in more than 1 dimension as well as spline smoothing. Has anyone ported the S+ GAM library to R?

[R] libraries in .First

2003-06-04 Thread Laurie Sindlinger
Dear all, I have a question regarding the .First function. I have included help.start() and several libraries in my .First as: .First - function() { help.start(browser = netscape7) library(lattice) library(modreg) library(splines) library(MASS) library(maps)} The libraries maps and splines

[R] X11 not available

2003-06-04 Thread f . grignola
Hi, I just installed R in a Linux box. Unfortunately, I don't have root access and had to install it in my home directory. The software runs fine, but I cannot make it print graphics to the screen. If I type X11 I get: X11 is not available (though it is working for other software). I noticed that

RE: [R] Building an R package under Windows NT

2003-06-04 Thread Benjamin . STABLER
The version of Perl I have is: This is perl, version 5.005_02 built for MSWin32-x86-object Copyright 1987-1998, Larry Wall Binary build 506 provided by ActiveState Tool Corp. http://www.ActiveState.com Built 15:40:37 Oct 27 1998 It was installed by my IS department some time ago. I

Re: [R] libraries in .First

2003-06-04 Thread Prof Brian Ripley
If this is R 1.7.0, the preferred mechanism is options(defaultPackages), as in the following (from the examples in .First) # Example of Rprofile.site local({ old - getOption(defaultPackages) options(defaultPackages = c(old, MASS)) }) .First is run too early to be

RE: [R] Building an R package under Windows NT

2003-06-04 Thread Benjamin . STABLER
The emme2 package is not in the $RHOME/library directory. My sh.exe is the most current one from the tools available at http://www.stats.ox.ac.uk/pub/Rtools/.I searched for sh.exe and that is only one I've got on my system. -Original Message- From: Duncan Murdoch [mailto:[EMAIL

[R] kmeans

2003-06-04 Thread Luis Miguel Almeida da Silva
Dear helpers I was working with kmeans from package mva and found some strange situations. When I run several times the kmeans algorithm with the same dataset I get the same partition. I simulated a little example with 6 observations and run kmeans giving the centers and making just one

[R] Update VR_7.1-6

2003-06-04 Thread Dirk Enzmann
The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, r) : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION' Other packages can be installed without

Re: [R] Update VR_7.1-6

2003-06-04 Thread Thomas Lumley
On Tue, 3 Jun 2003, Dirk Enzmann wrote: The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, r) : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION'

Re: [R] Update VR_7.1-6

2003-06-04 Thread Gavin Simpson
Hi Dirk, dse and VR are bundles or packages. There is a bug in rw1070 in the way it installs bundles. Download the zip file directly and unzip it to the library directory in wherever R is installed on your system. This has been discussed recently (this morning) on the list. Search the

Re: [R] Update VR_7.1-6

2003-06-04 Thread Prof Brian Ripley
I have explained that on R-help earlier today (and when VR_7.1-4 was announced). Please look back a few hours in the R-help postings. On Tue, 3 Jun 2003, Dirk Enzmann wrote: The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the

Re: [R] kmeans

2003-06-04 Thread Prof Brian Ripley
On Tue, 3 Jun 2003, Luis Miguel Almeida da Silva wrote: I was working with kmeans from package mva and found some strange situations. When I run several times the kmeans algorithm with the same dataset I get the same partition. Why does that surprise you? I simulated a little example with

Re: [R] coefficient of logistic regression

2003-06-04 Thread Thomas W Blackwell
Ahmet - In a logistic regression model, fitted probabilities make sense for individual cases (rows in the data set), as well as for future cases (predictions) for which no outcome (success or failure) has been observed yet. Fitted probabilities are calculated from the matrix formula:

Re: [R] rw1070 and package bundles; was: Update VR_7.1-6

2003-06-04 Thread Uwe Ligges
Prof Brian Ripley wrote: I have explained that on R-help earlier today (and when VR_7.1-4 was announced). Please look back a few hours in the R-help postings. On Tue, 3 Jun 2003, Dirk Enzmann wrote: The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files)

RE: [R] Building an R package under Windows NT

2003-06-04 Thread Benjamin . STABLER
I finally figured out the problem. I went ahead and installed Perl 5.8 but that didn't do it. The problem was that an existing program's directory was earlier in the PATH and so one (or more) of the components of the Rcmd INSTALL was an older version. It wasn't zip since I renamed the old zip.

[R] (no subject)

2003-06-04 Thread Gilda Garibotti
Hi, I would like to know if it is possible to get printed output while a loop is taking place. Example: for(i in 1:10){ print(i) some long process } This will print the values of i only after the loop is finished, what I would like is to see them when the process enters the i-th iteration to

Re: [R] (no subject)

2003-06-04 Thread Peter Dalgaard BSA
Gilda Garibotti [EMAIL PROTECTED] writes: Hi, I would like to know if it is possible to get printed output while a loop is taking place. Example: for(i in 1:10){ print(i) some long process } This will print the values of i only after the loop is finished, what I would like is to

Re: [R] coefficient of logistic regression

2003-06-04 Thread John Fox
Dear Ahmet, Sorry for the slow response, but I've been busy all today, coincidentally teaching a workshop on logistic regression. Tom Blackwell sent you a useful suggestion for interpreting coefficients on the odds scale. If you want to trace out the partial relationship of the fitted

[R] Question about looking up names

2003-06-04 Thread Ross Boylan
I think I now understand how R looks up names. Could anyone tell me if I have this right? First it looks up the nested environments created by lexical scoping. Then, if it gets to the top (.GlobalEnv) it also looks through the list of things that have been attached. It never looks in the call

Re: [R] Question about looking up names

2003-06-04 Thread Thomas Lumley
On 3 Jun 2003, Ross Boylan wrote: I think I now understand how R looks up names. Could anyone tell me if I have this right? First it looks up the nested environments created by lexical scoping. Then, if it gets to the top (.GlobalEnv) it also looks through the list of things that have been

[R] Logistic regression problem: propensity score matching

2003-06-04 Thread Paul
Hello all. I am doing one part of an evaluation of a mandatory welfare-to-work programme in the UK. As with all evaluations, the problem is to determine what would have happened if the initiative had not taken place. In our case, we have a number of pilot areas and no possibility of random

Re: [R] Question about looking up names

2003-06-04 Thread Ross Boylan
On Tue, 2003-06-03 at 16:34, Robert Gentleman wrote: Also, note that you can get the effect of lexical scope by doing things like Do you mean you can get the effects of dynamic scope? f- function(x) x+y e1 - new.env() assign(y, 10, env=e1) environment(f) - e1 #now like lexical

Re: [R] Rounding problem R vs Excel

2003-06-04 Thread Marc Schwartz
On Tue, 2003-06-03 at 08:32, Duncan Murdoch wrote: On Tue, 3 Jun 2003 09:49:44 +0100, you wrote in message [EMAIL PROTECTED]: Duncan If the numbers are not represently exactly how does R resolve problems like the one below? Is there something that needs to be set up in the R environment

Re: [R] X11 not available

2003-06-04 Thread Ott Toomet
Hello, | From: [EMAIL PROTECTED] | Date: Tue, 03 Jun 2003 16:23:20 + | | Hi, | | I just installed R in a Linux box. Unfortunately, I don't have root access and | had to install it in my home directory. | The software runs fine, but I cannot make it print graphics to the screen. If I

Re: [R] Your Message Contained a Potential Virus

2003-06-04 Thread Peter Dalgaard BSA
[EMAIL PROTECTED] writes: Result: Virus Detected Virus Name: [EMAIL PROTECTED] File Attachment: screensaver.scr Attachment Status: deleted --- Original message information follows --- From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Date: Tue, 3 Jun 2003 21:38:00 --0400 Subject: Re:

[R] plot rpart tree's from list object

2003-06-04 Thread Christian Schulz
Hello, i want the post plot's from a rpart list object with 18 tree's , getting no error - but getting no files,too? Perhaps i should using assign!? for (i in 1:length(treeList)) { post(treeList[[i]],filename=paste(Tree,i,sep=.ps), title=Arbeitszufriedenheit, digits=getOption(digits) -

Re: [R] (no subject)

2003-06-04 Thread Paul Lemmens
Hoi Peter, --On woensdag 4 juni 2003 0:16 +0200 Peter Dalgaard BSA [EMAIL PROTECTED] wrote: Gilda Garibotti [EMAIL PROTECTED] writes: Hi, I would like to know if it is possible to get printed output while a loop is taking place. Example: for(i in 1:10){ print(i) some long process } This will

Re: [R] Logistic regression problem: propensity score matching

2003-06-04 Thread Prof Brian Ripley
1) Why are you using multinom when this is not a multinomial logistic regression? You could just use a binomial glm. 2) The second argument to predict() is `newdata'. `sample' is an R function, so what did you mean to have there? I think the predictions should be a named vector if `sample'

Re: Rtips (was Re: [R] ? building a database with a the great examples

2003-06-04 Thread Detlef Steuer
On 03-Jun-2003 Paul E. Johnson wrote: I think there is a lot of merit in this idea. I think there is a big question about authentication and protection of Wikis from vandalism. Yes. But I`ll do my daily backups, so major vandalism wouldn't be such a problem. (I hope.) Personally I don't like

[R] Your Message Contained a Potential Virus

2003-06-04 Thread Symantec_AntiVirus_mskavs1
Text: This is an automated message. Please read it carefully. You should know that your recent email message detailed below, to Memorial Sloan-Kettering Cancer Center, was identified as potentially containing a virus. From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] If we have been able to

[R] Mode of MCMC chain

2003-06-04 Thread Patrik Waldmann
Hello, are there any functions in R for estimation of the mode of a MCMC-chain? Best, Patrik Waldmann [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Strip location and grid colour in Lattice

2003-06-04 Thread Mulholland, Tom
I am probably missing something quite obvious, but any help would be appreciated. I am continually getting people misreading the lattice plots because they are expecting the strip (with the factor names in them) to be below the graph. Is there anyway of achieving this. Secondly, from a more

Re: [R] (no subject)

2003-06-04 Thread Christian Hoffmann
Hi everybody, I finally hope to reach the person who started the thread [R] (no subject). The innermost level of text was my original request. Kind regards Christian At 10:50 2003-06-04 +0200, you wrote: Hoi Christian, --On woensdag 4 juni 2003 10:36 +0200 Christian Hoffmann [EMAIL

[R] predict.glm(glm.ob,type=terms)

2003-06-04 Thread orkun
hello pgeo-predict.glm(glm.ob,type=resp) works fine. But I need to get predictions values in terms of each factor variables. pgeo-predict.glm(glm.ob,type=terms) gives Error in rep(1/n,n) %*% model.matrix(object): non conformable arguments Could anyone tell me why ? Ahmet Temiz Turkey

[R] convert factor to numeric

2003-06-04 Thread Philipp Pagel
Hi R-experts! Every once in a while I need to convert a factor to a vector of numeric values. as.numeric(myfactor) of course returns a nice numeric vector of the indexes of the levels which is usually not what I had in mind: v - c(25, 3.78, 16.5, 37, 109) f - factor(v) f [1] 25

Re: [R] Rounding problem R vs Excel

2003-06-04 Thread Duncan Murdoch
On 04 Jun 2003 00:24:08 -0500, you wrote: Excel 2002 (XP): Cell Formula Value = 0.5 - 0.4 - 0.1 0.E+00 =(0.5 - 0.4 - 0.1)-2.7755575615628900E-17 ... What is interesting is the change in the displayed value in Excel when the second formula is surrounded