similar to: legend

Hi I have a scatterplot of points with pch=1 and a single point with pch=3, lwd=3. It has a high line width to attract attention to it. The following script plot(rnorm(10),rnorm(10),col="black") points(rnorm(10),rnorm(10),col="red") points(0,0,pch=3,lwd=3) if(TRUE){ legend("bottomleft",c("a","b","Truth"),pch=c(1,1,3),col=c

Dear R-help list, I have a problem with the tick marks of a Kaplan-Meier survival plot. Here is a sample: follow.up<-c(10,20,30,40,50,60,70,80,90,100) #months dead<-c(1,1,1,0,1,1,0,0,0,0) KM <-survfit(Surv(follow.up, dead)) plot(KM) The result is a nice plot. However, our research group thinks it may be a better idea to place the ticks to the years on the time scale, i.e. 0, 12, 24, 36

Hey, I am just starting to learn R now and I typed in this simple survival program: library(survival) t <- c(10,13,18,19,23,30,36,38,54,56,59,75,93,97,104,107,107,107) c <- c(1,0,0,1,0,1,1,0,0,0,1,1,1,1,0,1,0,0) data <- Surv(t,c) km <- survfit(data) summary(km) Call: survfit(formula = data) but everytime I run it I get this error: Error in

Dear Sirs, I want to estimate the survival mean of a few specific teams. I'm trying to calculate it through a Kaplan Meier estimator. For doing so, I load the "survival" package and run the following instructions: "options(survfit.print.mean=TRUE)" allows showing the mean and mean standard error "KM=survfit(Surv(Dias,Censura))"

Hello everybody, does anybody know how the function plot.survfit exactly works? I'd like to plot the log of the cummulative hazard against the log time by using plot.survfit(...fun="cloglog") which does not work correctly. The scales are wrong and there is an error message about infinit numbers. It must have something to do with the censored data, doesn't it? #Example:

Los datos no son de desgaste de cuchilla, sino de consumo de las mismas. Por ello tengo los datos de la siguiente forma: Unidades cambiadas Fecha En unidades cambiadas, suele ser una y en fecha el dia que se hizo el cmabio. Con eso no se muy bien como estructurar los datos para hacer el análisis. Gracias Jesús > Date: Mon, 7 Dec 2015 16:27:18 +0100 > From: griera en yandex.com

Dear all, I am trying to plot KM curves with confidence bands with function survplot under package rms. However, the following codes do not seem to work. The KM curves are produced, but the confidence bands are not there. Any insights? Thanks in advance. library(rms) ########data generation############ n <- 1000 set.seed(731) age <- 50 + 12*rnorm(n) label(age) <- "Age"

Hello I've decided to try and distill an earlier rather ill focused question to try and elicit a response. Any help is greatly appreciated. Why does mod.cox not work with sbrier whilst mod.km does? Can I make it work? > data(DLBCL) > DLBCL.surv<-Surv(DLBCL$time,DLBCL$cens) > > mod.km<-survfit(DLBCL.surv) > mod.cox<-survfit(coxph(DLBCL.surv~IPI, data=DLBCL)) >

Hello, it seems that the main results of survival analysis with package survival are shown only as side effects of the print method. If I compute e.g. a Kaplan-Meier estimate by > km.survdur<-survfit(s.survdur) then I can simply print the results by > km.survdur Call: survfit(formula = s.survdur) n events median 0.95LCL 0.95UCL 100.0 58.0 46.8 41.0 79.3 Is

Dear all, I am doing library(survival) fit <- coxph(Surv(futime,fustat) ~ rx, ovarian) plot(survfit(fit,newdata=ovarian),col=c(1,2)) legend("bottomleft", legend=c("rx = 0", "rx = 1"), lty=c(1,2),col=c(1,2)) Is this correct to compare these two groups? Is the 0.31 the p-value that the median f two groups are equal Why lty does not work here? Many thanks

Hi, I am wondering if anyone can explain to me if cumulative incidence (CI) is just "1 minus kaplan-Meier survival"? Under what circumstance, you should use cumulative incidence vs KM survival? If the relationship is just CI = 1-survival, then what difference it makes to use one vs. the other? And in R how I can draw a cumulative incidence plot. I know I can make a Kaplan-Meier

Hi all, I am trying to draw a Kaplan-Meier curve and I found online that Kaplan - Meier estimates are computed with a function called km in the event package. Is there an update for that because when I choose to download packages in R,. there is no package called event, even though I have selected all the repositories. Thanks in advance, Eleni [[alternative HTML version deleted]]

What is the best way to report the standard error when publishing Kaplan-Meier plots? In my field (Vascular Surgery), practitioners loosely refer to the "10% error" cutoff as the point at which to stop drawing the KM curve. I am interpreting this as the *standard error of the cumulative hazard*, although I'm having a difficult time finding some guidelines about this (perhaps I am

Hi, I am taking a survival class. Recently I need to do the Nelson-Aalen estimtor in R. I searched through the R help manual and internet, but could not find such a R function. I tried another way by calculating the Kaplan-Meier estimator and take -log(S). However, the function only provides the summary of KM estimator but no estimated values. Could you please help me with this? I would

Pero como haría el data frame?? Porque las cuchillas son de la misma referencia. En realidad es para ver cada cuanto se gstan las cuchillas y ver que pedidos hay que hacer de las mismas. La tabla que tengo es: 25 enero-> 1 cuchilla gastada 30 enero -> 1 cuchilla gastada 3 de febrero -> 2 cuchillas gastadas 5 de febrero -> 1 cuchilla gastada Y así.... No tiene necesariamente que ser

Hola a todos, -Estoy estudiando el efecto de dos genotipos (~tratamientos) en la aparición de síndrome metabólico (MetS) con datos longitudinales recogidos a tiempo 0,7,10,15,20 y 25 años. -He hecho un dataframe con las siguientes variables MetS: Síndrome Metabólico (Si=1,No=0) bmi: Indice de masa corporal (IMC) cuando se produce la conversión a MetS+ . Para los que permancen MetS-, esta variable

Hi, Sorry to trouble the list - I would like to ask a question - I can't find the answer in the r-help archives. I am trying to plot 2 survival curves in different colors. What I have tried is this plot(survfit(sim.surv ~ sim.km.smpl),col="red") but both the survival curves turn red... Your advice would be most appreciated! Thank you. Min-Han

Hi, I am attempting to graph a Kaplan Meier estimate for some claims using the survfit function. However, I was wondering if it is possible to plot a cdf of the kaplan meier rather than the survival function. Here is some of my code: library(survival) Surv(claimj,censorj==0) survfit(Surv(claimj,censorj==0)~1) surv.all<-survfit(Surv(claimj,censorj==0)~1) summary(surv.all) plot(surv.all)

Hello, I now receive an error message when obtaining a survival plot, which was not previously received with the same code. I recently updated all my packages. It does not seem to be a peculiarity of my data as I receive the error using data available in R. A plot is produced but I am uncertain regarding the error message. > library(surv2sample) > data(gastric) > fit <-

Buenas, Como pudeo calcular el tiempo de vida? Os cuento, tengo una serie de cuchillas y quiero ver el consumo de las mismas y he pensado en hacer un estudio por tiempo de vida. No se como hacerlo con R Gracias Jesús [[alternative HTML version deleted]]