similar to: nearest neighbors

Hi, I am having trouble with the survival library, particualrily the coxph function. the following works coxph(jtree9$cph.call,z,rep(1,dim(z)[1])) Call: coxph(formula = jtree9$cph.call, data = z, weights = rep(1, dim(z)[1])) coef exp(coef) se(coef) z p SM 0.2574 1.294 0.0786 3.274 1.1e-03 Sex -0.1283 0.880 0.1809 -0.709

I downloaded the latest version of R (1.5.0) and issued a read.table() command that used to work with (1.3.1). Unfortunately I erased 1.3.1 from my hard drive after I downloaded R 1.5.0. Here's the command and the error message I get under 1.5.0: > golub1<-read.table("C:/Microarrays/code/data_set_ALL_AML_train.txt",sep="\t",q

At a snail's pace I keep on translating an introduction to R into italian; I have reached the section describing the glm() function, in which some example code is presented. The very last line of code, before the beginning of the section on Poisson models is: ldp <- ld50(coef(fmp)); ldl <- ld50(coef(fmp)); c(ldp, ldl) which of course gives results 43.663 and 43.663; the correct code

At a snail's pace I keep on translating an introduction to R into italian; I have reached the section describing the glm() function, in which some example code is presented. The very last line of code, before the beginning of the section on Poisson models is: ldp <- ld50(coef(fmp)); ldl <- ld50(coef(fmp)); c(ldp, ldl) which of course gives results 43.663 and 43.663; the correct code

Hello, I thought some of you might like to know about a GNU project that is complementary to R in some ways, namely gretl http://ricardo.ecn.wfu.edu/gretl (GNU Regression, Econometrics and Time-series Library). gretl (a library with cli and gui clients, the gui using GTK) is designed to be very user-friendly, and suitable for teaching econometrics. It has a fairly wide variety of least-squares

I'm trying to do a Windows build of Ray Brownrigg's port of the maps library, and running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c

I'm trying to do a Windows build of Ray Brownrigg's port of the maps library, and running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c

I tend to use table() quite a bit for quick "diagnostics", summary, etc. I have wished for a more automatic way of plotting these. One possibility would be something like the following function; The question is if (something like) the following is worth providing (and then maintaining...) at all : plot.table <- function(x, type = "h", ylim = c(0, max(x)), lwd = 2,

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample

I'm trying to build the 1.4.0 development version of R for Windows. I'll be posting problems as they arise, until I give up (or succeed), or someone asks me to stop. 1. I don't have zlib, the compression library. The source file gnuwin32/unzip/WHERE says to get it from ftp://ftp.info-zip.org/pub/infozip/zlib/zlib.html; that web site says it has been closed and redirects me to

Just thought I would mention that there is a new O'Reilly book out, "Developing Bioinformatics Computer Skills", which is ok, pretty superficial about some things but at least it lets you know what is happening and where. R gets a few pages in there, 394-396, mostly nice press and accurate except that Bill Venables (hi Bill) may be surprised to find out that he is a member of the

Hi Duncan, I've been looking around at GUI toolkits to use with Python (and by extension R). I initially tried Tk, but was disappointed because there is no consistent set of high level widgets. I've subsequently done some looking at QT, wxWindows, and GTK. QT is out because I am a *commercial* user, and I don't want to have a license hassle. GTK seemed to have a lower level of

Hi, Kurt Hornik and I will start to publish a quarterly electronic ``R Newsletter'' beginning with January 2001 (if all goes well). The idea is to have a mixture of articles describing 1) new features in R itself 2) contributed add-on packages 3) nice applications 4) more general statistical computing issues related to R 5) or ... OK, now comes my christmas wish to you: If you are

Dear R-help, I think I have kinda isolated the problem I had to the following: i <- 0 while( {i <- i + 1} < 5) { if(i < 3) next print(i) } This seems to go into an infinite loop. After I break the execution, i has the value 1. At the R prompt, if I start from i <- 0 and keep typing {i <- i + 1} < 5, it eventually evaluate to TRUE. So why does the while loop not work?

Prof. Tierney, Thanks very much for the info. Why does the loop work if I move the assignment out of the condition? E.g., the following works: i <- 0 while(i < 5) { i <- i + 1 if(i < 3) next print(i) } Regards, Andy > -----Original Message----- > From: Luke Tierney [mailto:luke at nokomis.stat.umn.edu] > Sent: Tuesday, October 02, 2001 3:36 PM > To: Liaw, Andy

Byron Ellis has been making some progress on a hdf5 library for microarray data (and any other thing you want to put in there). In doing so some issues have arisen that are of more general interest. 1) hdf5 supports annotation (through comments) so it would be nice if the comment function in R became generic. I think this is backward compatible and basically not really an issue of any

Hi, The following function does work optimize.W<-function(W,k,G,Groups,cph.call,z){ n<-length(Groups) grp.wt<-rep(0,n) for(i in 1:(length(G))){ ind<-Groups == G[i] if(G[i]!=k){ grp.wt[ind]<-W[i] } elsegrp.wt[ind]<-1 } z<-data.frame(cbind(z,grp.wt=grp.wt)) #needed to make the case weights #part of the data

It was my fault, as you have seen. Y matrix must contain successes and failures, and not successes and trials. Sorry for my ignorance and thanks a lot for your help. Oscar -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20001220/341c2519/attachment.html

> From: Robert Gentleman [mailto:rgentlem@jimmy.harvard.edu] > > Greg, > you might be better off writing a wrapper, > > R->foo where foo has less than 65 args and does nothing more than > unpack them and calls bar (the original) with as many args as you > want, not much code change and not much cost... Hi Rob, [BTW, sorry about the crosspost

Prof. Gentleman (and R-help), Here's an example of what didn't work. I still don't understand why. Function: trycon <- function(file, n) { f.con <- file(file, open="rt") on.exit(close(f.con)) i <- 0 while( length(readln <- readLines(f.con, 1)) > 0 ) { x <- unlist(strsplit(readln, " ")) if(length(x) <= 6 && x[3] ==