Re: [R] factor always have type integer

2004-09-09 Thread Prof Brian Ripley
On Wed, 8 Sep 2004, Erich Neuwirth wrote: The function I need is valtype-function(x) typeof(ifelse(is.factor(x),levels(x),x)) It is easy enough to write. Are there any other special cases where the values and the storage mode differ? All classed objects are more than the internal

Re: [R] factor always have type integer

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Erich Neuwirth wrote: The simple answer to my problem (are the values of a vector numeric or not) is is.numeric, and that is enough for what I need right now. But this way I do not get an answer discriminating between integers and and doubles. What is the canonical way of

Re: [R] Cox regression for prevalence estimates

2004-09-09 Thread The Michaelson Institute
Thnks to Prof. Frank E Harrell Jr , Prof. Bernardo Rangel Tura, and Prof. Thomas Lumley for their response. Tomas Karpati (MD) The Michaelson Institute for the prevention of Blindness Hadassah Medical Org. email: [EMAIL PROTECTED] Phone: +972-2-6256458, fax: +972-2-6232895 On Tue, 07 Sep 2004

Re: [R] Problem installing R

2004-09-09 Thread Prof Brian Ripley
On Wed, 8 Sep 2004, Niggley, Jason wrote: As I am installing R for windows, I get the following error: C:\ProgramFiles\R\rw1091\library\stats\chtml\stats.chm Is that really not Program Files? An error occurred while trying to copy a file: The source file is corrupted. Can you help

[R] Evaluating R function in C code

2004-09-09 Thread Manoj - Hachibushu Capital
Dear All, For the following C code: SEXP testR(SEXP rho) { SEXP evalStr,evalRes; Char cmdstr[1024]; SEXP val; /* Get the value of variable 'val' defined already in the R environ */ PROTECT(val =

RE: [R] isoMDS

2004-09-09 Thread Jari Oksanen
On Wed, 2004-09-08 at 21:31, Doran, Harold wrote: Thank you. Quick clarification. isoMDS only works with dissimilarities. Converting my similarity matrix into the dissimilarity matrix is done as (from an email I found on the archives) d- max(tt)-tt Where tt is the similarity matrix. With

Re: [R] isoMDS

2004-09-09 Thread Jari Oksanen
On Thu, 2004-09-09 at 04:53, Kjetil Brinchmann Halvorsen wrote: Mardia, kent Bibby defines the standard transformation from a similarity matrix to a dissimilarity (distance) matrix by d_rs - sqrt( c_rr -2*c_rs + c_ss) where c_rs are the similarities. This assures the diagonal of the

[R] Handling the windows clipboard/32KB limit

2004-09-09 Thread Ritter, Christian C MCIL-CTANL/S
(R 1.9.1; Windows 2000;) I'm just comparing ease of use, speed, etc for methods of transferring data frames in the Excel, MySQL, R triangle. It turns out that going from Excel to R (when doing this carefully). Using the clipboard is actually quite fast and efficient (2 seconds for

Re: [R] Installing packages on OS X

2004-09-09 Thread Jari Oksanen
On Wed, 2004-09-08 at 21:25, hadley wickham wrote: On my computer, it seems that (binary?) packages installed through the GUI in RAqua are not used available to the command line version of R, while (source) packages installed with R CMD INSTALL are available to both. This is a problem when I

Re: [R] rodbc windows doesn't find dsn

2004-09-09 Thread Jack Tanner
Prof Brian Ripley wrote: 'odbcConnect' establishes a connection to the dsn, and 'odbcDriverConnect' allows a more flexible specification via a connection string. 'odbcConnect' uses the connection string 'DSN=dsn;UID=uid;PWD=pwd'. The point here seems to be that odbcConnect is

[R] filling patterns in image and examples

2004-09-09 Thread Marta Rufino
Dear R-users, I would like to do a image plot with patterns filling the squares instead of colors. How can I do it? Is there manuals with examples of image.plot? I search but could not find many. Thank you All the best Marta [[alternative HTML version deleted]]

[R] Indexing dataframe

2004-09-09 Thread Jacques VESLOT
I am sorry to ask such question, but I can't find a solution... I have a dataframe 'd2004' and I want to remove two columns: 'd2004$concentration' and 'd2004$stade. I could do it just as follows: names(d2004) [1] Localite Date parcelle maille presence.plant

Re: [R] Handling the windows clipboard/32KB limit

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Ritter, Christian C MCIL-CTANL/S wrote: (R 1.9.1; Windows 2000;) I'm just comparing ease of use, speed, etc for methods of transferring data frames in the Excel, MySQL, R triangle. It turns out that going from Excel to R (when doing this carefully). Using the clipboard

Re: [R] Indexing dataframe

2004-09-09 Thread Thomas Petzoldt
Jacques VESLOT wrote: I am sorry to ask such question, but I can't find a solution... I have a dataframe 'd2004' and I want to remove two columns: 'd2004$concentration' and 'd2004$stade. d2004$concentration - NULL d2004$stade - NULL Hope it helps! Thomas P.

AW: [R] Indexing dataframe

2004-09-09 Thread TEMPL Matthias
Hi, Use Newdata - subset(d2004, select=-c(concentration,stade)) See ?subset for details Best, Matthias -Ursprüngliche Nachricht- Von: Jacques VESLOT [mailto:[EMAIL PROTECTED] Gesendet: Donnerstag, 09. September 2004 12:38 An: [EMAIL PROTECTED] Betreff: [R] Indexing dataframe

Re: [R] Installing packages on OS X

2004-09-09 Thread Prof Brian Ripley
Here's my guess: R CMD INSTALL installs to .Library, and does not look at R_LIBS unless set in its environment. install.packages() installs at .libPaths()[1], and uses all the library trees available. So if you have R_LIBS set in ~/.Renviron (as I have), or .libPaths() set in .Rprofile or

RE: [R] isoMDS

2004-09-09 Thread Doran, Harold
Thank you. I use the same matrix on cmdscale as I did with isoMDS. I have reproduced my steps below for clarification if this happens to shed any light. Here is the original total matrix (see opening thread if you care how this is created) a b c d e f g h a 4 4 2 4 1 2 0 0 b 4 4 2 4 1 2 0 0

Re: [R] degrees of freedom (lme4 and nlme)

2004-09-09 Thread Ramon Diaz-Uriarte
Dear Elizabeth, When I looked for this a couple of years ago, I found DF's to be discussed in the book by Pinheiro Bates Mixed effects models for S and S-Plus, as well as the documentation for SAS's PROC MIXED (I believe that the discussion on df's on the SAS manual was more complete than on

RE: [R] isoMDS

2004-09-09 Thread Jari Oksanen
On Thu, 2004-09-09 at 14:25, Doran, Harold wrote: Thank you. I use the same matrix on cmdscale as I did with isoMDS. I have reproduced my steps below for clarification if this happens to shed any light. --- snip --- Doran, Your data clarified things. It seems to me now, that your data are not

RE: [R] isoMDS

2004-09-09 Thread Hanke, Alex
I get the following message: Error in isoMDS(tt) : zero or negative distance between objects 1 and 2 This makes sense since a and b are identical in their relationship to c to h. Drop row 1 and col 1 and you get isoMDS(tt[2:8,2:8]) initial value 14.971992 iter 5 value 8.027815 iter 10 value

[R] Adding GSL library path to SHLIB

2004-09-09 Thread Vicky Landsman
Dear R-list people, I asked a similar question a few hours before. I will try to be more specific. We like to add the GSL library to the file SHLIB in order to make it possible to run the C code using GSL functions from R. We read that the path to the libgsl.a should be added to the line

RE: [R] Handling the windows clipboard/32KB limit

2004-09-09 Thread Ritter, Christian C MCIL-CTANL/S
Thanks Brian for the hint. Indeed, writeClipboard does not suffer from the 32KB limit. Here is a line which can put large quantities of data on the clipboard: writeClipboard(gsub(NA,\#N/A,apply(Alldata,1,paste,collapse=\t))) this can then be pasted either manually or via a VBA macro into a

Re: [R] locator() in a multi-figure setting using mfrow(): SOLVED

2004-09-09 Thread Christoph Lehmann
thanks to some great hints by Paul Murrel I could solve it: here we are with one solution (code needs to be cleaned and simplified, but maybe one can understand it) ### ## create a multifigure setting nr - 4 nc - 2 opar - par(mfrow = c(nr, nc)) slices - 8 m - matrix(runif(100),10,10)

Re: [R] Adding GSL library path to SHLIB

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Vicky Landsman wrote: Dear R-list people, I asked a similar question a few hours before. I will try to be more specific. We like to add the GSL library to the file SHLIB in order to make it possible to run the C code using GSL functions from R. We read Read where?

[R] tree architecture in Random Forest

2004-09-09 Thread Ute
Hello, I'm doing classification with random Forest using numeric predictor variables. I would like to know a detail in the tree architecture: to which side of a node data are sent if , for the split variable, their value is higer than the split point? best regards, Ute Müller

RE: [R] tree architecture in Random Forest

2004-09-09 Thread Liaw, Andy
From: Ute Hello, I'm doing classification with random Forest using numeric predictor variables. I would like to know a detail in the tree architecture: to which side of a node data are sent if , for the split variable, their value is higer than the split point? The right, as one

RE: [R] confidence intervals

2004-09-09 Thread Rogers, James A [PGRD Groton]
Robert, I have this quick hack to obtain approximate Shewhart prediction intervals for variance component models fit with lme (to nitpick slightly, confidence intervals have the interpretation of containing parameters, while prediction and tolerance intervals have the interpretation of

[R] kolmogorov-smirnov for discrete ordinal scale data

2004-09-09 Thread Gila Lithwick
Hi, I was wondering whether there is an implementation of the Kolmogorov-Smirnov goodness of fit test for discrete, ordinal scale data in R - I've only managed to find the test for continuous data. Thanks! Gila __ [EMAIL PROTECTED] mailing list

[R] Copulas

2004-09-09 Thread J W
Is anyone already in the process of developing copula estimation methods or porting Prof. Rene Carmona's EVANESCE (http://www.orfe.princeton.edu/~rcarmona/SVbook/svbook.html)? JW __ [EMAIL PROTECTED] mailing list

RE: [R] kolmogorov-smirnov for discrete ordinal scale data

2004-09-09 Thread Hanke, Alex
Hi I think the answer is no. However, I have written a script that implements the test described in Testing for shifts in the Vertical Distribution of Plankton using a robust Kolmogorov-Smirnov like Statistic by Smith, Beet and Solow (1998). The test has the properties you are looking for. If this

[R] discriminant analysis

2004-09-09 Thread Rogério Rosa da Silva
Dear all, I'm having difficulty getting the procedure for discriminant analysis. How can I specify the forward or backward stepwise analyses of variables into the model? I'm using lda in MASS and discrimin (ade4). I run R version 1.9.1. Thanks, Rogério

Re: [R] Adding GSL library path to SHLIB

2004-09-09 Thread Vicky Landsman
Dear Prof. Ripley, We read the archive thread http://maths.newcastle.edu.au/~rking/R/help/02b/0547.html Thank you for your help, we will try to create the Makevars file. Vicky. - Original Message - From: Prof Brian Ripley [EMAIL PROTECTED] To: Vicky Landsman [EMAIL PROTECTED] Cc:

[R] SJava, Client X11

2004-09-09 Thread Laurent Buffat
Hi, I have installed the SJava package (see http://www.omegahat.org/RSJava/ ) without any problem. I was able to run the examples (see http://www.omegahat.org/RSJava/examples/index.html ) from R, on my “R Linux server” and directly on the display of the server, But, if I run theses

[R] scoping rules

2004-09-09 Thread Whit Armstrong
Can someone help me with this simple example? sq - function() { y - x^2 y } myfunc - function() { x - 10 sq() } myfunc() executing the above in R yields: myfunc() Error in sq() : Object x not found I understand that R's scoping rules cause it to look for x in

Re: [R] Adding GSL library path to SHLIB

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Vicky Landsman wrote: Dear Prof. Ripley, We read the archive thread http://maths.newcastle.edu.au/~rking/R/help/02b/0547.html Things have changed in 2 years in R -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics,

[R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Vicky Landsman
Following the recommendation of Prof. Ripley, I have created the Makevars file with the line: PKG_LIBS=-L/usr/lib/libm -lm -L/usr/local/lib/libgsl -lgsl -L/usr/local/lib /libgslcblas -lgslcblas in the working directory. Now I have the code file Example3.c which computes the Bessel function value

Re: [R] Installing packages on OS X

2004-09-09 Thread hadley wickham
Thanks to you both. Setting R_LIBS='~/Library/R/library/' in ~/.Renviron did the trick. Thanks, Hadley __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] Proposal for New R List: Criticism? Comments?

2004-09-09 Thread Berton Gunter
Folks: I would like to propose a new R list, tentatively labeled r-contents. I wish to briefly explain the purpose and format here and solicit public comments, pro or con, so feel free to criticize or suggest a better name and other improvements or alternatives. R presently consists of a suite

Re: [R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Vicky Landsman wrote: Following the recommendation of Prof. Ripley, I have created the Makevars file with the line: PKG_LIBS=-L/usr/lib/libm -lm -L/usr/local/lib/libgsl -lgsl -L/usr/local/lib /libgslcblas -lgslcblas in the working directory. No, that's not what I

RE: [R] Skipping panels in Lattice

2004-09-09 Thread Hanke, Alex
I have the same problem. As far as I can see, the only thing you can do is : attach(df2) group=paste(facb,facc,sep= ) bwplot( dv ~ faca | factor(group)) Alex -Original Message- From: Leon Barmuta [mailto:[EMAIL PROTECTED] Sent: September 9, 2004 1:19 AM To: [EMAIL PROTECTED] Subject:

Re: [R] Proposal for New R List: Criticism? Comments?

2004-09-09 Thread Jonathan Baron
I think that a lot of posts on r-help are exactly of the form you suggest: How do I do X? Answer: Use Y. (Or maybe, Use Y. And next time RTFM. But so what. The answer is still there.) Often, when the answer is not of that form, the question is unclear. In other cases, the questioner is

[R] R course in Princeton and Boston

2004-09-09 Thread eugene dalt
Anyone saved the upcoming R course announcement for Princeton and Boston? Please email me... Thanks Eugene. __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Vicky Landsman
Dear Prof. Ripley, The GSL reference book explains that lm should be added as well (as far as I understand, and this is what I saw when I run the file in C in a standard way). No matter, with lm and without, I get the same picture. Here is what I get after R CMD SHLIB Example3.c make: Warning:

[R] Confused about loading other packages from a package

2004-09-09 Thread hadley wickham
In my package, I create a new method for plot with the following signature: setMethod(plot, signature(x=marrayNorm, y=formula), plot.ma) where marrayNorm is a class defined in the marray package. After building and installing my package, I get the following warnings when I load my package

Re: [R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Paul Roebuck
On Thu, 9 Sep 2004, Vicky Landsman wrote: Following the recommendation of Prof. Ripley, I have created the Makevars file with the line: PKG_LIBS=-L/usr/lib/libm -lm -L/usr/local/lib/libgsl -lgsl -L/usr/local/lib /libgslcblas -lgslcblas in the working directory. That should probably look

Re: [R] Confused about loading other packages from a package

2004-09-09 Thread Roger D. Peng
In R 2.0.0-to-be I think putting marray in the Depends: field will solve your problem. Otherwise, you might want to import the marrayNorm class in your NAMESPACE file. -roger hadley wickham wrote: In my package, I create a new method for plot with the following signature: setMethod(plot,

Re: [R] scoping rules

2004-09-09 Thread Duncan Murdoch
On Thu, 9 Sep 2004 12:01:56 -0400 , Whit Armstrong [EMAIL PROTECTED] wrote : Can someone help me with this simple example? sq - function() { y - x^2 y } myfunc - function() { x - 10 sq() } myfunc() executing the above in R yields: myfunc() Error in sq() : Object

Re: [R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Vicky Landsman wrote: Dear Prof. Ripley, The GSL reference book explains that lm should be added as well (as far as I understand, and this is what I saw when I run the file in C in a standard way). No matter, with lm and without, I get the same picture. Here is what I get

Re: [R] R/S-plus Course***In Princeton Boston, ***R/Splus Fundamentals and Programming Techniques, september - October, 2004

2004-09-09 Thread Spencer Graves
[EMAIL PROTECTED] wrote: R/S-plus Course***In Princeton Boston,***R/Splus Fundamentals and Programming Techniques, september - October, 2004 XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce 2-day R/S-plus Fundamentals and Programming Techniques. Princeton,NJ

Re: [R] Confused about loading other packages from a package

2004-09-09 Thread hadley wickham
Unfortunately, marrayClasses doesn't seem to have a namespace: library(maVis) Error in loadNamespace(imp[[1]], c(lib.loc, .libPaths()), keep.source) : package 'marrayClasses' does not have a name space Error in library(maVis) : package/namespace load failed Any other ideas? Hadley

[R] Blom's approximation to rankits?

2004-09-09 Thread Lisa Wang
Hello, My name is Lisa and I'm a statistician at Princess Margare Hospital. I wonder if there is any function in R that calculate the Normal rankits based on Blom's approximation? Thank you very much Lisa Wang Msc. Princess Margaret hospital Toronto, Ca

Re: [R] Confused about loading other packages from a package

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, hadley wickham wrote: Unfortunately, marrayClasses doesn't seem to have a namespace: library(maVis) Error in loadNamespace(imp[[1]], c(lib.loc, .libPaths()), keep.source) : package 'marrayClasses' does not have a name space Error in library(maVis) :

[R] GEE model

2004-09-09 Thread Lisa Wang
Hello there, My name is Lisa and I'm doing analysis using genelized linear model-GEE model because I have some repeated measurements on the same patients (in a bilateral eye disease case). It would be very appreciated if you could shed some light on what function to use in R. Thank you very

Re: [R] Blom's approximation to rankits?

2004-09-09 Thread Bob Wheeler
You need not use Blom's approximations. Exact values to several decimal places are given by normOrder() in SuppDists. Lisa Wang wrote: Hello, My name is Lisa and I'm a statistician at Princess Margare Hospital. I wonder if there is any function in R that calculate the Normal rankits based on

[R] Rslides.sty

2004-09-09 Thread Doran, Harold
I recently came across Rslides.sty during a web search and was not aware of it. However, I can't seem to get it to actually work (I use miktex and the texniccenter distribution for Windows). Does anyone have a sample tex file with this style to share? Second, I looked around the R web site to

RE: [R] Dyn.load of sharing object with GSL library

2004-09-09 Thread Liaw, Andy
As Prof. Ripley said, twice already: seek local expert help. If you don't know who that is, start with the person that setup your user account on that machine. Such system-specific problem is much easier to resolve that way, instead of having others on the list guessing. Andy From: Vicky

[R] RDA file

2004-09-09 Thread T. Murlidharan Nair
What is an RDA file ? I was trying to work through a turorial on cluster analysis using R and bioconductor. The data files in the library are stored as RDA files. Its part of the golubEsets library. Cheers ../Murli __ [EMAIL PROTECTED] mailing list

[R] Efficient Cartesian product of data.frames

2004-09-09 Thread Timothy W. Victor
Hello List, I am looking for efficient code to produce the Cartesian product of two or more data.frames. I'd like to be able to do this without resorting to looping. I have searched the FAQ, web, etc without luck. That being said, the help page for merge says that the function can produce what

RE: [R] R conversion

2004-09-09 Thread John Fox
Dear Mark, See ?factanal for a factor-analysis function, in the standard stats package (though factanal does ML factor analysis and not principal axes). Note that help.search(factor analysis) turns this up. With respect to your last question, it's not possible to know the source of the

RE: [R] Proposal for New R List: Criticism? Comments?

2004-09-09 Thread John Fox
Dear Bert, I believe that you've identified an important issue -- and one that's occasionally been discussed on this list previously -- but I'm not sure that another email list is a good solution. Some method of indexing functions in packages that would allow people to more easily locate them

RE: [R] RDA file

2004-09-09 Thread Liaw, Andy
See ?save and ?load. Andy From: T. Murlidharan Nair What is an RDA file ? I was trying to work through a turorial on cluster analysis using R and bioconductor. The data files in the library are stored as RDA files. Its part of the golubEsets library. Cheers ../Murli

Re: [R] R conversion

2004-09-09 Thread Mark Strivens
Thanks John for the post - I had found factanal However factanal it seems you have specify the number of factors to be fitted up front, whereas with the SAS procedure you don't - this is apparently important to the analysis! Perhaps I could emulate this function by breaking it down using factor

Re: [R] Efficient Cartesian product of data.frames

2004-09-09 Thread Gabor Grothendieck
Timothy W. Victor tvictor at dolphin.upenn.edu writes: I am looking for efficient code to produce the Cartesian product of two or more data.frames. First create some test data consisting of a list of n=2 data frames. data(iris) L - list(iris1 = iris[1:3,1:2], iris2 = iris[1:3,3:4])

Re: [R] converting to data frame

2004-09-09 Thread Prof Brian Ripley
On Thu, 9 Sep 2004, Murli Nair wrote: Is there a method that converts a csv file to a data frame ? There are R _functions_, if that is what you are asking. (method is a technical term in R.) Try help.search(csv), which gives me one entry. -- Brian D. Ripley, [EMAIL