Re: [R] missing values and survival analysis

2004-04-12 Thread Frank E Harrell Jr
On Sun, 11 Apr 2004 23:37:22 -0400 [EMAIL PROTECTED] wrote: Hi everyone, I'm analysing a survival analysis data set at the moment with missing values in the covariate and survival vectors (I have about 60 variables). I know there are some functions on the CRAN network to deal with

[R] Windows startup question

2004-04-12 Thread Erich Neuwirth
I would like to be able to double click on an R source file and by that start R and have the file sourced in R. I am using WinXP Pro. So I created an association for file type R and tye following Run action: C:\Program Files\R\rw1090beta\bin\Rgui.exe R_PROFILE=%1 I have an R file with the

[R] R 1.9.0 is release

2004-04-12 Thread Peter Dalgaard
I've rolled up R-1.9.0.tgz a short while ago. This is a new version with a number of new features, most notably a substantial reorganization of the standard packages, a major update of the grid package, and the fact that underscore can now be used as a regular character in variable names. See

[R] question on isoMDS

2004-04-12 Thread Helmut Kudrnovsky
Hello everyone, I have a question on isoMDS. My data set (of vegetation) with 210 samples is in this way: Rotfoehrenau Lavendelweidenau Silberweidenau 067_Breg.7 02 0 071_Dona.4 02

Re: [R] Regression models w/ splines

2004-04-12 Thread Thomas Lumley
On Fri, 9 Apr 2004, Frank E Harrell Jr wrote: On Fri, 09 Apr 2004 15:16:34 -0500 John Fieberg [EMAIL PROTECTED] wrote: Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted

Re: [R] all poss. regression

2004-04-12 Thread Thomas Lumley
On Sat, 10 Apr 2004, Prof Brian Ripley wrote: Look at package leaps. Note that it applies to columns of the design matrix and not terms, so can be nonsensical with categorical predictors (let alone interactions of such). A reasonable compromise is to run regsubsets() (in leaps) returning

Re: [R] (offtopic) I need two sets of 5 different color scales

2004-04-12 Thread Thomas Lumley
On Sat, 10 Apr 2004, Tamas Papp wrote: I need advice in the following: I need two sets of colors for B and F which are easy to distinguish (when printed on a color laser printer), represent cardinality (ie have an intuitive mapping to an interval) or at least ordinality. It's not offtopic.

Re: [R] pasting results into Word/Excel

2004-04-12 Thread Thomas Lumley
On Sun, 11 Apr 2004, Graham M Smith wrote: Anders, This seems a good start, and I will try it get it to work, on of the issues is that the results output is more complex than simply split by comma eg a typical Summary output of my current data set (below) My current approach for model

RE: [R] Only referenc copy when calling C routine?

2004-04-12 Thread Thomas Lumley
On Sun, 11 Apr 2004, Shin, Daehyok wrote: Thanks, Andy. Let me ask a little more about the DUP argument. In the manual, I can see that only reference copy happens when passing an array, if DUP=FALSE. Then, what happens when returning it in a list, if DUP=FALSE? Another reference copy, or

Re: [R] pasting results into Word/Excel

2004-04-12 Thread Graham M Smith
Thomas, My current approach for model summaries is to use the xtable package and write things as HTML tables. It was an HTML approach that Gabor suggested,using R2HTML, and that worked well, but I will try xtable as well. As always with R, there is a solution if you know where to look.

[R] Panel Data Analysis

2004-04-12 Thread Mihai Nica
Greetings: I am trying to find a package/solution for panel (longitudinal) data analysis. Unfortunately it seems I don't know where to start. Could somebody offer a hint? Thanks, Mihai Nica Jackson State University 155 B Parkhurst Dr. Jackson, MS 39202 No good deed will ever remain unpunished

Re: [R] Panel Data Analysis

2004-04-12 Thread Douglas Bates
Mihai Nica [EMAIL PROTECTED] writes: I am trying to find a package/solution for panel (longitudinal) data analysis. Unfortunately it seems I don't know where to start. Could somebody offer a hint? Analysis of longitudinal data is discussed in Pinheiro and Bates (2000), Mixed-effects Models in

[R] R-1.9.0 sources available by rsync

2004-04-12 Thread Douglas Bates
The sources for R-1.9.0 are now available via rsync at rsync.r-project.org. You can, for example, use rsync -avC rsync.r-project.org::r-release ./R-1.9.0 to obtain a copy of the sources. Binary distributions for platforms such as Windows and Mac OS X will be available from the CRAN sites in

Re: [R] Panel Data Analysis

2004-04-12 Thread Mihai Nica
Thanks, this is what I needed! Sincerely, Mihai Nica Jackson State University 155 B Parkhurst Dr. Jackson, MS 39202 No good deed will ever remain unpunished - Original Message - From: Douglas Bates [EMAIL PROTECTED] To: Mihai Nica [EMAIL PROTECTED] Cc: R-Help [EMAIL PROTECTED] Sent:

[R] von Mises Density

2004-04-12 Thread MKG
Dear All, I has posted this question on friday. I am trying to plot von Mises density on the circle. One can use dvm function from the CircStats package, by giving a set of angles, mu and kappa to plot the circular density on the line. Does any one have a macro that does it on the circle? These

[R] [JOB ADV] Computational Statistician

2004-04-12 Thread Warnes, Gregory R
We've just opened a position for a computational statistician in the Non-Clinical Statistics group here at Pfizer. We're looking for someone with strong computational skills in R and SAS, who is interested in creating real applications. The focus will be on statistical genetics, but we

[R] String problems

2004-04-12 Thread MMarques Power
after reading some docs I still can not find anything how to escape the \ string I know that \ in files needs to be escaped in order to be read by R. But in Windows if I try to call any shell util I keep getting errors. example: afile - g:\\simulation\\at works fine with read.table(afile) ...

[R] Random Forest:how to do an automatic rerun using only the important variables

2004-04-12 Thread Hui Han
Hi, I am using the Random Forest Package, and want to do an automatic rerun using only those variables that were most important in the original run. Is there anybody who has experience with this and can give me helpful suggestions? Best regards, Hui Han Department of Computer Science and

Re: [R] String problems

2004-04-12 Thread Uwe Ligges
MMarques Power wrote: after reading some docs I still can not find anything how to escape the \ string I know that \ in files needs to be escaped in order to be read by R. But in Windows if I try to call any shell util I keep getting errors. example: afile - g:\\simulation\\at works fine with

Re: [R] String problems

2004-04-12 Thread James MacDonald
Plus, do you really want your system call to be system(c:\\dir g:\\simulation\\at) not sure what you want, but that is what you will get. Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109

[R] Line numbers in error messages

2004-04-12 Thread Webb Sprague
Hi all, I have searched but to know avail. Is there a way to get a line number when a function crashes? I am doing an edit-source-run cycle. Feel free to cc me as I subscribe to the list in digest. Thanks W __ [EMAIL PROTECTED] mailing list

RE: [R] Random Forest:how to do an automatic rerun using only the important variables

2004-04-12 Thread Liaw, Andy
That's the advantage of having an R interface to RF: you can do such automation rather easily. E.g., twoStageRF - function(x, y, nVar=round(0.1*ncol(x)), ...) { imp - randomForest(x, y, importance=TRUE, ...)$importance[,3] cutoff - sort(imp, decreasing=TRUE)[nVar] randomForest(x[, imp =

[R] Complex sample variances

2004-04-12 Thread Fred Rohde
Hello, Is there a way to get complex sample variances in the survey package on summary statistics other than means? If not, can they be added to a future version? It would be be great to have them on totals, quantiles, ratios, and tables (eg row percent, columns percent, etc). Thanks.

Re: [R] Line numbers in error messages

2004-04-12 Thread Uwe Ligges
Webb Sprague wrote: Hi all, I have searched but to know avail. Is there a way to get a line number when a function crashes? I am doing an edit-source-run cycle. I don't think so, but traceback(), debug() and options(error = recover) are your friends. Also, you might want to check out

Re: [R] Complex sample variances

2004-04-12 Thread Thomas Lumley
On Mon, 12 Apr 2004, Fred Rohde wrote: Hello, Is there a way to get complex sample variances in the survey package on summary statistics other than means? If not, can they be added to a future version? It would be be great to have them on totals, quantiles, ratios, and tables (eg row

[R] [R-pkgs] new version of gregmisc package

2004-04-12 Thread Warnes, Gregory R
gregmisc 0.10.1 --- gregmisc 0.10.1 is now available at ftp://cran.r-project.org/incoming/ and will move to the regular contributed packages location (http://cran.r-project.org/src/contrib/PACKAGES.html#gregmisc) as soon as the CRAN administrators get caught up. Recent Changes

[R] Matrix decomposition

2004-04-12 Thread GIDEON WASSERBERG
I am looking for a manual(s) or any kind of documentation, at the introductory level, regarding matrix algebra (specifically, matrix population models). Any help will be highly appreciated Gideon Gideon Wasserberg (Ph.D.) Wildlife research unit, Department of wildlife ecology, University of

Re: [R] Line numbers in error messages

2004-04-12 Thread Webb Sprague
Hi Patrick, It's very simple using a browser() line in your function somewhere you know your code's OK, then run line by line. The problem is that sometimes you have code of a few hundred lines, to which you have added a strange little line that craps out because of some silly mistake that

Re: [R] Line numbers in error messages

2004-04-12 Thread A.J. Rossini
Uwe Ligges [EMAIL PROTECTED] writes: Webb Sprague wrote: Hi all, I have searched but to know avail. Is there a way to get a line number when a function crashes? I am doing an edit-source-run cycle. I don't think so, but traceback(), debug() and options(error = recover) are your

Re: [R] Line numbers in error messages

2004-04-12 Thread Jason Turner
Webb Sprague wrote: I already use ESS. If it would give line numbers, my life would be perfect! It does, if you don't use source(). Have your R code in one buffer, and the R session running, and (from the code buffer), type C-c C-l. Cheers Jason

Re: [R] Line numbers in error messages

2004-04-12 Thread Thomas Lumley
On Tue, 13 Apr 2004, Jason Turner wrote: Webb Sprague wrote: I already use ESS. If it would give line numbers, my life would be perfect! It does, if you don't use source(). Have your R code in one buffer, and the R session running, and (from the code buffer), type C-c C-l. Yes, yes,

RE: [R] Matrix decomposition

2004-04-12 Thread Liaw, Andy
I don't know anything about matrix population models, but Golub Van Loan's `Matrix Computations' is a standard reference for numerical computations with matrices. There's also David Harville's fairly recent book, `Matrix Algebra from Statisticians' Perspective'. I believe one of James Gentle's

[R] fractal calculation using fdim

2004-04-12 Thread Fred J.
Hello I am getting this Error in slopeopt(AllPoints, Alpha) : Object LineP not found when the dim(data.frame) of the data matrix for which the fdim is being calculated hits a certain values depending on the data set. e.g print(dim(mat)) fd - fdim(mat,q=2) [1] 27432 [1] 27423 [1] 27414

Re: [R] Line numbers in error messages

2004-04-12 Thread Duncan Murdoch
On Mon, 12 Apr 2004 15:20:58 -0700, you wrote: Hi Patrick, It's very simple using a browser() line in your function somewhere you know your code's OK, then run line by line. The problem is that sometimes you have code of a few hundred lines, to which you have added a strange little line that

Re: [R] Line numbers in error messages

2004-04-12 Thread Tony Plate
Isn't source file information often recorded in the source attribute on functions (or calls)? Could either the execution engine or the debugger refer to that information? (Though, in the debugger it might be impossible to uniquely identify expressions that appear multiple times in the

Re: [R] Line numbers in error messages

2004-04-12 Thread Patrick Connolly
On Mon, 12-Apr-2004 at 03:20PM -0700, Webb Sprague wrote: | Hi Patrick, | | It's very simple using a browser() line in your function somewhere you | know your code's OK, then run line by line. | | The problem is that sometimes you have code of a few hundred lines, to | which you have added a

RE: [R] Matrix decomposition

2004-04-12 Thread Bret Collier
Gideon, You might look into Caswell, H. 2001. Matrix population models 2nd edition. Sinauer Associates. Caswell's book covers a majority of matrix population modeling, including some good introductory information and has general matrix manipulation as an appendix. I think most of the work

Re: [R] Line numbers in error messages

2004-04-12 Thread Thomas Lumley
On Mon, 12 Apr 2004, Tony Plate wrote: Isn't source file information often recorded in the source attribute on functions (or calls)? Could either the execution engine or the debugger refer to that information? (Though, in the debugger it might be impossible to uniquely identify expressions

Re: [R] Complex sample variances

2004-04-12 Thread Fred Rohde
Thanks. I'll update the survey package. Sudaan does the standard errors on quantiles using Taylor series. If I can hunt down the formula it uses, could you add that to svyquantile? Fred Thomas Lumley [EMAIL PROTECTED] wrote: On Mon, 12 Apr 2004, Fred Rohde wrote: Hello, Is there a way

[R] Specifying RHOME

2004-04-12 Thread Anar Khan
Hi, I am having trouble installing R-1.8.1 on my Fedora 1 machine. R is already installed in /usr/local on our network, but I wish to be able to work on my own copy. So I installed it locally using: ./configure --prefix=/scratch/akhan/R-1.8.1 make make install The installation looks fine

Re: [R] fractal calculation using fdim

2004-04-12 Thread Robert W. Baer, Ph.D.
Well I must say that I am still struggling to understand exactly how fdim works, and what form the first argument is to take. Presumably the function produces something like 1-slope for a log log plot as a fractal dimension. I guess it is really the makeline() function I don't understand from the

Re: [R] question on isoMDS

2004-04-12 Thread Jari Oksanen
Helmut Kudrnovsky, This issue was discussed on this list recently: Search for messages with Subject: Fwd: MDS problems at and after 26 March. Brief answer to your question is: no, you cannot run isoMDS in MASS if you have zero or negative dissimilarities. You should remove identical cases. After