similar to: discrepancies between stata and r for a cox regression

I would like to produce a complimentary log-log survival plot with only the points appearing on the graph. I am using the code below, taken from the plot.survfit page of help for the the survival package (version 2.35-7). I am running in R 2.10.0 on Windows XP, and the list of packages following the error is loaded. Is there some specific 'type= ' syntax, or an additional parameter that

Hi, I've lost my mind on it... I have to scatterplot two vectors, grouped by a third variable, with two different dimensions according to whether each cell line in the plot is sensitive or resistant to a given drug, and with a different color for each of 9 tissues of origin. Here's what I've done:

Hi everybody, this is a real dummy thing. I sorted a matrix based on a given column, and what I get is right, until it comes to columns of negative and positive values; than, "order" orders everything from max to min in the negative values, and then AGAIN from max to min in the positive values!!! Why isn't everything order from max to min, and that's it? Thank you!!! Attached

Dear friends: I'm interested to make a stacked plot of cumulative incidence. that's, the cuminc model is fitted [fit=cuminc(time, relapse)] and cumulative incidence is in place. I'd like to stack the cuminc plots (relapse of luekemia and death free from leukemia, for example) , then the constituent ratio of leukemia relapse and treatment related mortality is very clear. Can

I am using plot(survfit(Surv(time,status) ~...) and would like to add error bars rather than the confidence intervals. Am I able to do this at specified times? e.g. when time = 20 & 40. leukemia.surv <- survfit(Surv(time, status) ~ x, data = aml) plot(leukemia.surv, lty = 2:3,xlim = c(0,50)) #can i add error bars at times 20 & 40? legend(100, .9, c("Maintenance", "No

dear list, i have a problem using the adaboost function from the package boost. running the example of leukemia works out very well, but so it does not on my own data. i always get the following 22 warnings: 1: cannot compute exact p-value with ties in: wilcox.test.default(x[which(y == 0)], x[which(y == 1)]) can anybody please help me with the meaning of this warning. thnak you stefan

Hi all, I have this microarray large microarray data set (ALL) from which I would like to subset or extract a set of data based on a factor ($mol.biol). I looked up some example of subsetting in, picked up two commands and tried both but I got error messages as follows > testset <- subset(ALL, ALL$mol.biol %in% c("BCR/ABL","ALL1/AF4")) >> Error in

Hi, I am using rpart to do leave one out cross validation, but met some problem, Data is a data frame, the first column is the subject id, the second column is the group id, and the rest columns are numerical variables, > Data[1:5,1:10] sub.id group.id X3262.345 X3277.402 X3369.036 X3439.895 X3886.935 X3939.054 X3953.777 X3970.352 1 32613 HAM_TSP 417.7082 430.4895 619.4776 720.8246

Dear all, I've tried to solve the Es. 12, cap 4 of "Introduction to GLM" by Annette Dobson. It's about the relationship between survival time of leukemia patients and blood cell count. I tried to fit a model with exponential distribution, first by glm (family gamma and then dispersion parameter fixed to 1) and then with survreg. They gave me the same point estimates but the

A postdoctoral fellow position in computational biology is anticipated to open in the fall of 2009 at University of Washington, Seattle, WA. The candidate will join a multi-disciplinary group working on edge-cutting research problems in computational biology. Job Description We seek a highly motivated and diligent researcher at the early stage of academic career with the ability to carry out

Hi, I've been trying unsuccessfully over the past couple of days to get my windows machines to recognize my linux box. I've been working through the troubleshooting guides I can't get the query of the client machine to work using the broadcast address: I have the following 3 machine network: lifebook SUSE 192.168.1.100 renegade WXP 192.168.1.101 maverick W2K

Ref: 14/F/04 Scientific computing systems administrator - IMPPC - Badalona, Spain Job Overview We are seeking a Linux systems administrator who will be responsible for the management of a scientific computing facility that includes a 100 core cluster, several servers, workstations, 96 Tb raw disk storage and backup facilities, to provide support to bioinformatics researchers at the Institute of

Hallo, I just installed all needed packages for my project on my PC. But I cannot load all at one time. I now want to load limma. How can I realize the following plan: I want to install for example limma inclusive all needed other sub packages (add-on). Can anyone tell me the corresponding command? Thanks, Corinna Here is the result of the command library(): Pakete in Library

Dear R users! I want to process a maximum likelihood estimation for a parametric regression survival time model with multiple events per subject. the STATA command for this survival regression is: use survreg stset failure(exercise), id(optionid) local regressors itm posret negret streg `regressors', distribution(weibull) explanation: stset declares data to be survival-time data; exercise

Dear help news reader, I'm trying to draw a Kaplan-Meier curve and would like to ask the news group for some help Supposing I have study comapring two drugs, "A", and "B" and I recorde the time to get to the clinical endpoint (Time), in my case becommming virus free. I have setup the following frame: Time c Drug 1 5 1 A 2 7 1 B 3 2 1 A 4 10 1

Hello, I write to know how can I modify the x axis : when I plot a survival object, R plots a graph with x values = 0, 10, 20, 30 while I want a graph with values 0, 6, 12, 18, 24 in the x axis. How can I do this? In R 2.1.1 version there was "time.inc" in survplot, but in version R 2.2.0 there isn't it! I am sorry for my english and I hope that you understand my problem. Thank you

Hi: I am trying to query pubmed abstracts using the following syntax: url= "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/esearch.fcgi?" search = paste(url, "db=pubmed&term=", queryTerm1, "+AND+", queryTerm2,"+OR+",queryTerm3, "+OR+", queryTerm4, "[abstract]&retmax=100&usehistory=y", sep="") docId <-

Hello, I have a question and need your help urgently. I am new to R but want to learn it. I have several files in a folder which I have imported to R using : temp = list.files(pattern="*.txt") >myfiles = lapply(temp, read.delim) The resulting files are on the workspace stored as List[110]. So they are 110 files in the list. Each file has several different columns and rows. My

Dear Everybody: I'm doing my usual "how does that work in R" thing with some Stata projects. I find a gross gap between the Stata and R in Cox PH models, and I hope you can give me some pointers about what goes wrong. I'm getting signals from R/Survival that the model just can't be estimated, but Stata spits out numbers just fine. I wonder if I should specify initial

Hi, I am conducting a survival analysis and I have been trying, with not success, to do a log-log survival plots. Does anyone know how to do it in R? By the way I am using version 1.0.0 for windows. Many Thanks Marco -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000503/e64aa9ee/attachment.html