similar to: R Kaplan-Meier plotting quirks?

Better would be to use interval censored data. Create your data set so that you have (time1, time2) pairs, each of which describes the interval of time over which the tag was lost. So an animal first captured at time 10 sans tag would be (0,10); with tag at 5 and without at 20 would be (5,20), and last seen with tag at 30 would be (30, NA). Then survit(Surv(time1, time2,

Hello R users I have a question with Kaplan-Meier Curve with respect to my research. We have done a retrospective study on fillings in the tooth and their survival in relation to the many influencing factors. We had a long follow-up time (upto 8yrs for some variables). However, we decided to stop the analysis at the 6year follow up time, so that we can have uniform follow-up time for all the

Hola!! Estoy intentando ejecutar un script com horas, al principio ejecute estos comandos DBx$Date<-strptime(DBx$Date, "%d-%m-%Y") ###Monicap use ; other use Y DBx$Year<-as.POSIXlt(DBx$Date)$year+1900 if(filename!="monicap_50.csv") {DBx$Time<-paste(DBx$Time, ":00", sep="")} Pero me daba el error de que mi base de datos tenia las

When I try to the code from library(survival) of library(ISwR), the following code survfit(Surv(days,status==1)) that could produce Kaplan-Meier estimates shows the following error "Error in survfit(Surv(days, status == 1)) : Survfit requires a formula or a coxph fit as the first argument" How it can be done in R.2.10 -- View this message in context:

Hi All, Please pardon me if I am missing something obvious here. How do I get the Kaplan-Meier estimate function that is created by survfit and plotted by the code. fit <- survfit(Surv(time, status) , data=aml) plot(fit) That is, I need a function that will give me the survival estimate at a given time: \hat{S}(t). Thanks in advance. Ritwik Sinha ritwik.sinha at gmail.com | +12033042111 |

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]]

Dear R-users I am trying to make an adjusted Kaplan-Meier curve (using the Survival package) but I am having difficulty with plotting it so that the plot only shows the curves for the adjusted results. My data come from a randomised controlled trial, and I would like the adjusted Kaplan-Meier curve to only show two curves for the adjusted survival: one for those on treatment (Treatment==1)

Hi all, I am getting an error in my code and I don't know what the problem is. I am using R 2.9 on ubuntu. my code is as follows: ## Libraries ## library(survival) library(foreign) ## reading data ## data<-read.dta("http://psfaculty.ucdavis.edu/bsjjones/cabinet.dta") head(data) attach(data) fit1<-survfit(Surv(durat,censor)) and I get the following error >

Dear all, I have data from 1970 to 1990 for people above age 50. Now I want to calculate survival curves by age starting at age 50 using the Kaplan Meier Estimator. The problem I have is that there are already people in 1970 who are older than 50 years. I guess this is called delayed entry or left truncation (?). I thought the code would be: roland <- survfit(Surv(time=age.enter,

Hola Daniel, si perdona di a responder directamente y no me di cuenta. Ya se donde esta el error, pero queria preguntar si a puede ser que mi R funcione mal o algo porque esta manhana ejecute mi script y funciono perfectamente y ahora volvi a ejecutarlo y me volvio a dar el mismo problema de ayer , despues de reiniciar y demas el tinnR y el R, no se si me vacila o es que tengo algo mal en

Hola a todos, Me gustaria pedir vuestra ayuda a encontrar el error que no consigo encontrar en este archivo. He revisado todo mil veces y probado y no doy con ello.Adjunto el archivo con Google drive porque es muy grande. ? monicap_50.csv <https://docs.google.com/file/d/0B8o2KrPEgG7ATlBMc19lTVk1d3M/edit?usp=drive_web> ? Este es el script, y lo que no entiendo que pasa es que tengo 592044

dear all, I want to plot a kaplan Meier plot with the following functions, but I fail to produce the plot I want: library(survival) tim <- (1:50)/6 ind <- runif(50) ind[ind > 0.5] <- 1; ind[ind < 0.5] <- 0; MS <- runif(50) pred <- vector() pred[MS < 0.3] <- 0; pred[MS >= 0.3] <- 1 df <- as.data.frame(cbind(MS, tim, pred, ind)) names(df) <-

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

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

There are two ways to view weights. One is to treat them as case weights, i.e., a weight of 3 means that there were actually three identical observations in the primary data, which were collapsed to a single observation in the data frame to save space. This is the assumption of survfit. (Most readers of this list will be too young to remember when computer memory was so small that we had to

Dear all, I'm using the "survival" package with R 2.4.0 on Mac OS X 10.4.8. I have two core statistics books (one of which is Altman's medical stats book) which suggest showing the number of individuals at risk at different time intervals on the Kaplan-Meier curve. My plot shows two curves that later cross, because of one significant outlier. I have two queries: Is there an

Is it possible to calculate Kaplan-Meier for left-truncated, right-censored data using survival5? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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

Hello All, ? Am replicating in R an analysis I did earlier using SAS. See this as a test of whether I'm ready to start using R in my day-to-day work. ? Just finished replicating a Kaplan Meier analysis. Everything seems to work out fine except for one thing. The 95% CI around my estimate for the median is substantially larger in R than in SAS. For example, in SAS I have a median of 3.29 with a

Helo: After changing "involuntarily" some of the graphics parameters with the command par() (I did not know that changes with this command are permanent), now when I made a plot of the survival Kaplan-Meier function, the Y axis does not start at 1, and the X axis does starts at 0. The commands that I use are: library(survival) BROWN.SPV = Surv(BROWN$TEMPS, BROWN$DEF)

I have some data which is censored and I want to determine the median. Its actually cost data for a cohort of patients, many of whom are still on treatment and so are censored. I can do the same sort of analysis for a survival curve and get the median survival... ...but can I just use the survival curve functions to plot an X axis that is $ rather than date? If not is there some other way to