similar to: case weights-coxph (solved)

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

Is there an implementation of a reasonable k-nearest neighbor finder already in one of the packages? -- +---------------------------------------------------------------------------+ | Robert Gentleman phone : (617) 632-5250 | | Associate Professor fax: (617) 632-2444 | | Department of Biostatistics office: not yet

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

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

Hi, I am trying to automate the creation of tables for some simply analyses. There are lots and lots of tables, thus the creation of a user-defined function to make and output them to excel. My problem is that some of the analyses have convergence issues, which I want captured and included in the output so the folks looking at them know how to view those estimates. I am successfully able to do

1. I did not attempt to sort through your voluminous code. But I suspect you are trying to reinvent wheels. 2. I don't understand this: "I've failed to find a solution after much searching of various R related forums." A web search on "error handling in R" **immediately** brought up ?tryCatch, which I think is what you want. If not, you should probably explain why it

You can capture warnings by using withCallingHandlers. Here is an example, its help file has more information. dataList <- list( A = data.frame(y=c(TRUE,TRUE,TRUE,FALSE,FALSE), x=1:5), B = data.frame(y=c(TRUE,TRUE,FALSE,TRUE,FALSE), x=1:5), C = data.frame(y=c(FALSE,FALSE,TRUE,TRUE,TRUE), x=1:5)) withWarnings <- function(expr) { .warnings <- NULL # warning handler will

tryCatch() is good for catching errors but not so good for warnings, as it does not let you resume evaluating the expression that emitted the warning. withCallingHandlers(), with its companion invokeRestart(), lets you collect the warnings while letting the evaluation run to completion. Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, Mar 6, 2018 at 2:45 PM, Bert Gunter <bgunter.4567 at

Hi, Is there a command that will give the total size of an R object in bytes? thanks Nicholas -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at

Dear R community, I need to plot the results of some simulations I did using QTL Cartographer. I am plotting LOD scores over three chromosomes. The three plot have to be one next to the other. The procedure I am using is: par(mfrow=c(1,3)) plot(x$x, x$y, ylim=c(0,35), type="l", col="blue", las=1, xaxs="i", yaxs="i", xlab="X Chromosome",

After update from centos 6.6 to centos 6.7 and reboot it, I have get a lot of this error into /var/log/messages: > May??3 11:27:20 s-virt kernel: EDAC MC0: CE row 2, channel 1, label > "": (Branch=0 DRAM-Bank=2 RDWR=Read RAS=6093 CAS=896, CE Err=0x10000 > (Correctable Patrol Data ECC)) > May??3 11:27:21 s-virt kernel: EDAC MC0: CE row 2, channel 1, label > "":