Re: [R] Odd r-squared

2003-11-03 Thread Prof Brian Ripley
On Mon, 3 Nov 2003, Simon Wotherspoon wrote: I would consider the calculation of r-squared in the following to be a bug, but then, I've been wrong before. It seems that R looks to see if the model contains an intercept term, and if it does not, computes r-squared in a way I don't

[R] Problem with batch-file

2003-11-03 Thread Morach Sascha, moracsa1
Hi I have written a graphical user interface using tcltk. Now I would like to run it under rterm.exe with a batch file. I tried the following code written in a .bat file: d:\R\rw1071\bin\Rterm.exe --no-restore --no-save d:\RGui.r d:\RGui.out The problem is that rterm does open but the gui

[R] calling R from Perl

2003-11-03 Thread martin
Hi, I want to call R from Perl to generate plots to be displayed on a webpage. What I found out so far is that there is a package called RSPerl on www.omegahat.org which should do what I need. However in the description it says it has been tested with R 1.3.* the latest. I'm using R 1.8.0

Re: [R] Problem with batch-file

2003-11-03 Thread Prof Brian Ripley
That's exactly what a batch file should do. There is no GUI to run. Please tell us what you actually want to do! On Mon, 3 Nov 2003, Morach Sascha, moracsa1 wrote: Hi I have written a graphical user interface using tcltk. Now I would like to run it under rterm.exe with a batch file. I

Re: [R] calling R from Perl

2003-11-03 Thread Jonathan Baron
On 11/03/03 10:37, [EMAIL PROTECTED] wrote: Hi, I want to call R from Perl to generate plots to be displayed on a webpage. What I found out so far is that there is a package called RSPerl on www.omegahat.org which should do what I need. However in the description it says it has been tested

[R] Questions in R

2003-11-03 Thread vasilis pappas
Hello, I have 4 questions in R and I am looking for 4 answers! the questions are the below: 1) is there a function in R for the Least Significant Difference method for multiple comparisons? 2) Is there a function in R that gives me the studentized deleted residuals and the leverage points? 3) is

Re: [R] Questions in R

2003-11-03 Thread John Fox
At 11:40 AM 11/3/2003 +, =?iso-8859-7?q?vasilis=20pappas?= wrote: Hello, I have 4 questions in R and I am looking for 4 answers! the questions are the below: 1) is there a function in R for the Least Significant Difference method for multiple comparisons? There is a multcomp package that

Re: [R] Questions in R

2003-11-03 Thread Prof Brian Ripley
On Mon, 3 Nov 2003, John Fox wrote: At 11:40 AM 11/3/2003 +, =?iso-8859-7?q?vasilis=20pappas?= wrote: Hello, I have 4 questions in R and I am looking for 4 answers! the questions are the below: 1) is there a function in R for the Least Significant Difference method for multiple

RE: [R] Question about the high dimensional density estimation

2003-11-03 Thread Liaw, Andy
From: Jason Turner [mailto:[EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Hi, I found that the R package KernSmooth can deal with only 1D and 2D data. But now I have a collection of 4-dimensional data (x1,x2,x3,x4) and would like to estimate the mode of the underlying density. What can

[R] second Y axis

2003-11-03 Thread sagar
please send me the manual abnout the [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Visualising Vectors

2003-11-03 Thread Sean O'Riordain
Hi Laura, you should find some useful information in the latest R news Volume 3/2, October 2003 http://cran.r-project.org/doc/Rnews/ refer page 8 where Paul's figure 2 shows some novel symbols showing China Sea Wind Speed, Direction and Temperature... plotted by lat.long. s. Laura Quinn

Re: [R] Odd r-squared

2003-11-03 Thread kjetil
On 3 Nov 2003 at 17:53, Simon Wotherspoon wrote: If you are really interested, you could write your own function to calculate r-squared, deciding from the model matrix if the range space of the model matrix contains a constant vector. If model is the model matrix with n rows, one way is to

[R] lang2(...) with two and more arguments

2003-11-03 Thread Timur Elzhov
Dear R-help, how could I create an R call in C code using lang2 with 2 and more arguments? I tried this code: SEXP f(SEXP fn, SEXP rho) { SEXP R_fcall, x, y; PROTECT(R_fcall = lang2(fn, R_NilValue)); PROTECT(x = allocVector(REALSXP, 1)); PROTECT(y = allocVector(REALSXP, 1));

RE: [R] Weird problem with median on a factor

2003-11-03 Thread RBaskin
Continuing to beat the greasy spot in the road where the dead horse used to be 1) I know that the people building r are working on bigger and better things than this silly question and I appreciate the existence of this complicated package that was dropped in my lap for free. 2) Tony Platt

RE: [R] Calling R from PHP on Win2K?

2003-11-03 Thread Pikounis, Bill
John, I am sorry for the delayed reply...since it sounds like you would need to gain some knowledge on how to integrate components like Apache, MySQL, PHP, and R, I highly recommend one or more of Paul Dubois' books: 1) MySQL 2) MySQL Cookbook 3) MYSQL and Perl for the Web While the primary

[R] FDR in p.adjust

2003-11-03 Thread Arne.Muller
Hello, I've a question about the fdr method in p.adjust: What is the threshold of the FDR, and is it possible to change this threshold? As I understand the FDR (please correct) it adjusts the p-values so that for less than N% (say the cutoff is 25%) of the alternative hypothesis the Null is in

Re: [R] USA map

2003-11-03 Thread kjetil
On 3 Nov 2003 at 9:46, Owen, Jason wrote: library(maps) # since two days available precompiled for windows map('usa') Kjetil Halvorsen R users, In S, there was a function called usa() that would draw the map of the United States, plus it had other options for graphics. I have looked but

Re: [R] USA map

2003-11-03 Thread Martin Wegmann
On Monday 03 November 2003 15:46, Owen, Jason wrote: R users, In S, there was a function called usa() that would draw the map of the United States, plus it had other options for graphics. I have looked but I can't find the equivalent in R. Is there one? Description I don't know the usa()

Re: [R] lang2(...) with two and more arguments

2003-11-03 Thread Prof Brian Ripley
You set the first argument to x and then to y. Why are you surprised? On Mon, 3 Nov 2003, [koi8-r] Timur Elzhov[koi8-r] wrote: Dear R-help, how could I create an R call in C code using lang2 with 2 and more arguments? I tried this code: SEXP f(SEXP fn, SEXP rho) { SEXP R_fcall,

RE: [R] Weird problem with median on a factor

2003-11-03 Thread Prof Brian Ripley
On Mon, 3 Nov 2003 [EMAIL PROTECTED] wrote: Continuing to beat the greasy spot in the road where the dead horse used to be 1) I know that the people building r are working on bigger and better things than this silly question and I appreciate the existence of this complicated package

Re: [R] lang2(...) with two and more arguments

2003-11-03 Thread Peter Dalgaard
Timur Elzhov [EMAIL PROTECTED] writes: Dear R-help, how could I create an R call in C code using lang2 with 2 and more arguments? I tried this code: SEXP f(SEXP fn, SEXP rho) { SEXP R_fcall, x, y; PROTECT(R_fcall = lang2(fn, R_NilValue)); PROTECT(x = allocVector(REALSXP,

[R] comparing characters

2003-11-03 Thread Aurora Torrente
Hi all, I´m having some trouble when trying to compare character values (to check if they are alphabetically ordered). Is it possible to do it in any way? Thanks for your help. Cheers, Aurora __ [EMAIL PROTECTED] mailing list

Re: [R] comparing characters

2003-11-03 Thread Spencer Graves
What sort of trouble are you having? I just got the following from R 1.8.0: ab [1] TRUE ab [1] FALSE aba [1] FALSE aba [1] TRUE hope this helps. spencer graves Aurora Torrente wrote: Hi all, I´m having some trouble when trying to compare character values (to check if they are

[R] problem building MS-Windows package under linux

2003-11-03 Thread Mathieu Ros
hi there, trying to follow the steps of Yan Rossini 2003, I have two problems : first when I 'make CrossCompileBuild', I get : ** make[1]: Entering directory `/home/ros/RWORK/CROSSSOMP/WinR/R-1.8.0/src/gnuwin32' make -f

[R] hclust doesn't return merge details

2003-11-03 Thread Arne Neumann
Dear R-users, I tried to receive the merge details of a clustering by using the summary function of hclust. For illustration I use the Longley data as done by Prof Ripley (Wed 11 Apr 2001) d - dist(longley.y) d - d/max(d) hc - hclust(d, ave) But instead of getting a matrix for $merge I get:

Re: [R] comparing characters

2003-11-03 Thread Marc Schwartz
On Mon, 2003-11-03 at 10:46, Aurora Torrente wrote: Hi all, Im having some trouble when trying to compare character values (to check if they are alphabetically ordered). Is it possible to do it in any way? Thanks for your help. Cheers, Aurora Depending upon what it is you are

[R] Binaries for RORacle and V 1.8. needed

2003-11-03 Thread Claus Gwiggner
Hi, there is no Oracle Client on my machine (Redhat Linux 9.0/686i, R1.8.0). I read the readme.client file and did: I unziped the package ROracle and copied it into the R-library path. I get: library(ROracle) Error in testRversion(descfields) : This package has not been installed properly Who

[R] mva Hclust, heatmap and plotting functions

2003-11-03 Thread Jason Skelton
Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one. When plotting dendrograms created by hclust you can identify clusters by clicking on the graphics and returning a list of what is contained in each cluster.

Re: [R] comparing characters

2003-11-03 Thread Aurora Torrente
That´s what I needed. Thanks! Spencer Graves wrote: What sort of trouble are you having? I just got the following from R 1.8.0: ab [1] TRUE ab [1] FALSE aba [1] FALSE aba [1] TRUE hope this helps. spencer graves __ [EMAIL PROTECTED] mailing

Re: [R] comparing characters

2003-11-03 Thread Giovanni Petris
You could do something like that: x - sample(LETTERS, 10) x [1] K N C F R E L J S Q all.equal(order(x), 1:length(x)) [1] Mean relative difference: 0.5090909 When x is a numeric vector, I usually use any(diff(x) 0) although I don't know whether it's more efficient. HTH, Giovanni

Re: [R] comparing characters

2003-11-03 Thread Uwe Ligges
Aurora Torrente wrote: Hi all, I´m having some trouble when trying to compare character values (to check if they are alphabetically ordered). Is it possible to do it in any way? Thanks for your help. Cheers, a b should work. If that's not your question, you have to be more specific. Uwe Ligges

Re: [R] problem building MS-Windows package under linux

2003-11-03 Thread Paulo Justiniano Ribeiro Jr
I had this same problem before: you need to install the program tidy in your linux machine There is a p[ackage for DEBIAN such that you can use apt-get install tidy don't know about other distros but it is probably available as well P.J. [EMAIL PROTECTED]:~$ apt-cache search tidy perltidy -

RE: [R] hclust doesn't return merge details

2003-11-03 Thread Liaw, Andy
From: Arne Neumann [mailto:[EMAIL PROTECTED] Dear R-users, I tried to receive the merge details of a clustering by using the summary function of hclust. For illustration I use the Longley data as done by Prof Ripley (Wed 11 Apr 2001) d - dist(longley.y) d - d/max(d) hc -

Re: [R] hclust doesn't return merge details

2003-11-03 Thread Thomas W Blackwell
Arne - I have carried out exactly your example below, and I get hc$merge as a matrix with two columns and 15 rows. Do str(hc) to see a useful representation of the contents of the returned list. help(hclust) describes this list in the section Value:. help(Subscript) shows the various

Re: [R] FDR in p.adjust

2003-11-03 Thread James MacDonald
There is no threshold for fdr. The adjusted p-values give the expected proportion of false positives for all comparisons with similar p-values or smaller. In other words, if you choose a p-value of 0.05, you would expect that ~5% of the tests with a p-value of 0.05 or smaller are false positives.

Re: [R]A matrix is full rank is equal to having independent columns?

2003-11-03 Thread Deepayan Sarkar
For any matrix, the following definitions hold: row rank: number of linearly independent rows column rank: number of linearly independent columns There is a theorem stating that these 2 numbers must be the same for any matrix, and (consequently) that number is defined as the 'rank' of the

Re: [R] problem building MS-Windows package under linux

2003-11-03 Thread Prof Brian Ripley
You do need to build R before you build a package, when cross-building. `Small' errors stop make running: try make -k? On Mon, 3 Nov 2003, Mathieu Ros wrote: hi there, trying to follow the steps of Yan Rossini 2003, I have two problems : first when I 'make CrossCompileBuild', I get :

Re: [R]A matrix is full rank is equal to having independent columns?

2003-11-03 Thread kjetil
On 3 Nov 2003 at 13:32, Feng Zhang wrote: This are the same concept, no additional proof is needed. The rank of a matrix is the max number of li columns, or rows (which are the same). Kjetil Halvorsen Dear R listers, Just a simple question. If we say an nxm matrix (nm) is full rank of m,

Re: [R] comparing characters

2003-11-03 Thread Martin Maechler
Giovanni == Giovanni Petris [EMAIL PROTECTED] on Mon, 3 Nov 2003 12:00:41 -0600 (CST) writes: Giovanni You could do something like that: x - sample(LETTERS, 10) x Giovanni [1] K N C F R E L J S Q all.equal(order(x), 1:length(x)) Giovanni [1] Mean relative difference:

[R] svm in e1071 package: polynomial vs linear kernel

2003-11-03 Thread ryszard . czerminski
I am trying to understand what is the difference between linear and polynomial kernel: linear: u'*v polynomial: (gamma*u'*v + coef0)^degree It would seem that polynomial kernel with gamma = 1; coef0 = 0 and degree = 1 should be identical to linear kernel, however it gives

Re: [R] USA map

2003-11-03 Thread kjetil
On 3 Nov 2003 at 16:33, Martin Wegmann wrote: On Monday 03 November 2003 15:46, Owen, Jason wrote: R users, In S, there was a function called usa() that would draw the map of the United States, plus it had other options for graphics. I have looked but I can't find the equivalent in

Re: [R] USA map

2003-11-03 Thread Ray Brownrigg
I don't know the usa() function in S but usa() in R gives a map of the US Where did you find that function? In rw1080 on windows XP I get: Some earlier versions of the maps package (which was then available for Unix only), did have individual functions for each country in the world

Re: [R] USA map

2003-11-03 Thread apjaworski
I believe in R you need library(maps) before using usa(). Of course, the library has to be installed first. Cheers, Andy __ Andy Jaworski Engineering Systems Technology Center 3M Center, 518-1-01 St. Paul, MN 55144-1000 - E-mail: [EMAIL PROTECTED] Tel:

Re: [R] svm in e1071 package: polynomial vs linear kernel

2003-11-03 Thread David Meyer
On Mon, 3 Nov 2003 [EMAIL PROTECTED] wrote: I am trying to understand what is the difference between linear and polynomial kernel: linear: u'*v polynomial: (gamma*u'*v + coef0)^degree It would seem that polynomial kernel with gamma = 1; coef0 = 0 and degree = 1

[R] ROC with GLM?

2003-11-03 Thread itstruei
Hello R-List: Does anybody have code to optimize a logistic regression using ROC curves? I've seen S+ code that does it but never in R. Thanks. __ You Rock! Your E-mail should, too. Visit Rock.com! [[alternative HTML version

[R] Generic Function

2003-11-03 Thread Arend P. van der Veen
Hi, How to I write a generate function that is specific to the class of the argument itself ? For example, I created a data.frame that contains results of an analysis. I would like to write a special plot function with the following usage plot(data.frame) I have been looking through the

[R] write.dta and handling labels

2003-11-03 Thread Arnab mukherji
Hello, I need to write out a data matrix as a STATA 7 file and this happens perfectly with write.dta(), except I cannot seem to export the labelnames to Stata. So far I have tried the following: # X is the data matrix that is to be exported attributes(X)$var.labels - c(apple, banana, cat)

Re: [R] write.dta and handling labels

2003-11-03 Thread Spencer Graves
Have you considered write.table? hope this helps. spencer graves Arnab mukherji wrote: Hello, I need to write out a data matrix as a STATA 7 file and this happens perfectly with write.dta(), except I cannot seem to export the labelnames to Stata. So far I have tried the following: # X is

RE: [R] Generic Function

2003-11-03 Thread Gabor Grothendieck
Here is a simple example using S3 classes. We define a generic function myday. It dispaches, i.e. calls, a function whose name is myday followed by a dot followed by the class of the first argument to myday. After defining myday, define two functions for it to dispatch: myday.POSIXct and

[R] Using _ in names?

2003-11-03 Thread Spencer Graves
I work with data bases that routinely use the underscore, _, as part of names. I've seen discussion on this list of possibly allowing that in a future release of R. How far are we away from that? I got the following from R 1.8.0 under Windows 2000: A_1 Error: syntax error A_B -

Re: [R] write.dta and handling labels

2003-11-03 Thread Thomas Lumley
On Tue, 4 Nov 2003, Arnab mukherji wrote: Hello, I need to write out a data matrix as a STATA 7 file and this happens perfectly with write.dta(), except I cannot seem to export the labelnames to Stata. So far I have tried the following: # X is the data matrix that is to be exported

Re: [R] Using _ in names?

2003-11-03 Thread Thomas Lumley
On Mon, 3 Nov 2003, Spencer Graves wrote: I work with data bases that routinely use the underscore, _, as part of names. I've seen discussion on this list of possibly allowing that in a future release of R. How far are we away from that? About five months. -thomas