search for: futime

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung, dist='weibull') plot (lung.wbs) Returns the error msg: Error in xy.coords(x, y, xlabel, ylabel, log) : x and y length...

...s if I just do some regular graphing tasks. But as long as I run lines(fit) with "fit" object is a survfit object, this would reset to default plot region. See the ovarian example below: > > > > library(survival) > > ovarian1<-ovarian > > ovarian1$fustat[ovarian$futime>450]<-0 > > ovarian1$futime[ovarian$futime>450]<-450 > > ovarian2<-subset(ovarian,futime>450) > > > > fit1 <- survfit(Surv(futime, fustat) ~ rx, data = ovarian1) > > fit2 <- survfit(Surv(futime, fustat) ~ rx, data = ovarian2) > > > >...

[code]>require("survival") > coxph(Surv(futime,fustat)~age + strata(rx),ovarian) Call: coxph(formula = Surv(futime, fustat) ~ age + strata(rx), data = ovarian) coef exp(coef) se(coef) z p age 0.137 1.15 0.0474 2.9 0.0038 Likelihood ratio test=12.7 on 1 df, p=0.000368 n= 26, number of events= 12 > coxph(Surv(futime,fu...

Hi all, Using: R version 2.8.1 Patched (2009-03-07 r48068) on OSX (10.5.6) with survival version: Version: 2.35-3 Date: 2009-02-10 I get the following using the first example in ?summary.survfit: > summary( survfit( Surv(futime, fustat)~1, data=ovarian)) Call: survfit(formula = Surv(futime, fustat) ~ 1, data = ovarian) time n.risk n.event survival std.err lower 95% CI upper 95% CI 59 26 1 0.962 0.0377 0.890 1.000 115 25 1 0.923 0.0523 0.826 1.000 156...

...Executing it, gives the error "Error in floor(temp) : Non-numeric argument to mathematical function" best regards, Heinz library(survival) ## Example from help page of survdiff ## Expected survival for heart transplant patients based on ## US mortality tables expect <- survexp(futime ~ ratetable(age=(accept.dt - birth.dt), sex=1,year=accept.dt,race="white"), jasa, cohort=FALSE, ratetable=survexp.usr) Error in floor(temp) : Non-numeric argument to mathematical function sessionInfo('survival') R version 2.12...

...has to do with what happens when I specify a time point t for the upper limit of the time range under consideration (we want to avoid using the right-end tail of the KM curve). Copying from the example in the help file: TR <- ovarian[1:16,] TE <- ovarian[17:26,] train.fit <- coxph(Surv(futime, fustat) ~ age,x=TRUE, y=TRUE, method="breslow", data=TR) lpnew <- predict(train.fit, newdata=TE) Surv.rsp <- Surv(TR$futime, TR$fustat) Surv.rsp.new <- Surv(TE$futime, TE$fustat) If time is left NULL, the maximum futime in the TE data is assumed: UnoC(Surv.rsp, Su...

...he principle should be the same with left truncation I hope) library(survival) library(eha) # COMPARE coefs between survreg ('survival' pkg) and aftreg ('eha' pkg) #Fitting NULL models (no covariates) results in (approximately) the same coefs (which is good!) m1_NULL=survreg(Surv(futime/365, status==1) ~ 1, data=pbcseq) m2_NULL=aftreg(Surv(futime/365, status==1) ~ 1, data=pbcseq) c(m1_NULL$coef, 1/m1_NULL$scale) #--> intercept= 3.878656 , shape = 1.478177 c(m2_NULL$coef[1], exp(m2_NULL$coef[2])) #--> intercept= 3.878859 , shape=1.478150 # NOW I adjust for covariates m...

...nf)==0 are well defined. Below is a short example to reproduce the problem. Does anyone know why survreg() must behave that way? Is there an alternate solution to this problem? Sincerely, Jerome Asselin library(survival) data(ovarian) newovarian <- ovarian newovarian$lower59 <- newovarian$futime-59 newovarian$time59 <- Surv(newovarian$lower59,newovarian$futime, event=rep(3,nrow(newovarian)),type="interval") survreg(time59~ecog.ps+rx,data=newovarian,dist="lognormal") #THIS DOES NOT WORK BECAUSE ONE OF THE LOWER BOUNDS IS ZERO #Error in survreg(time59 ~ ecog.ps + r...

Hola Estoy tratando de correr un survival analysis usando un Cox regression model. Tengo una duda respecto a la organizacion del script. Tengo una variable que es -tamano del individuo- y quiero ver si hay diferencia en sobrevivencia respecto a tamano. Como diseno de campo los tamanos fueron ubicados de forma aleatoria en bloques al azar. Cuado planteo el script tengo algo como:

Hi, I've been trying to get LLVM working as a JIT compiler on Android for a while. It works now, except that the setLastModificationAndAccessTime function won't build because the Bionic C library lacks both futimes and futimens. There doesn't appear to be any "nice" workaround for this problem: the only ways I can think of to make setLastModificationAndAccessTime work are either to go by /proc/self/fd or invoke the system call directly (the second option isn't helped by the fact that t...

...p.c > cc -g -D__i386__ -o makedep makedep.o > ../tools/makedep -C. -S.. -T.. -I/usr/include/freetype2 fnt2bdf.c fnt2fon.c make_ctests.c makedep.c relpath.c sfnt2fnt.c > `makedep' is up to date. > ../../tools/makedep -C. -S../.. -T../.. ffs.c fstatvfs.c futimes.c getopt.c getopt1.c getpagesize.c gettid.c interlocked.c isinf.c isnan.c lstat.c memcpy_unaligned.c memmove.c mkstemps.c poll.c pread.c pwrite.c readlink.c spawn.c statvfs.c strcasecmp.c strerror.c strncasecmp.c symlink.c usleep.c > cc -c -I. -I. -I....

...in the list of covariates in the output of survreg function with standard error, z, p.value etc? Does it mean that intercept was fitted with the covariates? Does Value column represent coefficients or some thing else? Regards, ------------------------------------------------- tmp = survreg(Surv(futime, fustat) ~ ecog.ps + rx, ovarian, dist='weibull',scale=1) > summary(tmp) Call: survreg(formula = Surv(futime, fustat) ~ ecog.ps + rx, data = ovarian, dist = "weibull", scale = 1) Value Std. Error z p (Intercept) 6.962 1.322 5.267 1.39e-07 e...

...sorry). I can't seem to find a package or function that can handle 2 competing risks for the main event. I realize I may be using the well-known ones, timereg and cmprisk, wrong. And if I am, then I'd appreciate a correction. Here's the model I tried (and it crashes) comp.risk(Surv(futime, retention==2) ~ sex, data=dta, dta$retention, causeS=c(1,3), resample.iid=1, n.sim=100) I apologize, but I am cannot create the example dataset using code. I am not able. The data, though is one to one. ID retention sex futime ------------------------------------ 1 0 1...

Dear All: I have some questions about the output of coxph. Below is the input and output: ---------------------------------------- > coxph(formula = Surv(futime, fustat) ~ age + rx + ecog.ps, data = + ovarian, x = TRUE) Call: coxph(formula = Surv(futime, fustat) ~ age + rx + ecog.ps, data = ovarian, x = TRUE) coef exp(coef) se(coef) z p age 0.147 1.158 0.0463 3.17 0.0015 rx -0.815 0.443 0.6342 -1.28 0.2000...

...1 1 1 ... $ chemo : Factor w/ 2 levels "0","1": 1 1 2 1 1 2 2 1 1 2 ... $ mets : Factor w/ 2 levels "0","1": 1 2 2 1 2 2 2 1 1 2 ... $ status : Factor w/ 4 levels "1","2","3","4": 1 3 2 1 3 3 3 1 1 3 ... $ futime : num 10.46 1.15 2.43 2.83 6.82 ... $ smk : Factor w/ 2 levels "0","1": 2 2 2 1 2 1 2 1 2 2 ... $ surg.yr : int 88 94 92 93 86 85 95 98 91 85 ... $ nodes : Factor w/ 2 levels "0","1": 1 1 2 1 1 2 1 1 1 2 ... $ os : num 0 1 0 0 1 1...

..._SETRESUID + /* Define if you have the setrlimit function. */ #undef HAVE_SETRLIMIT --- configure.orig Thu Jun 7 14:12:53 2001 +++ configure Thu Jun 7 14:13:40 2001 @@ -3184,7 +3184,7 @@ fi -for ac_func in arc4random atexit b64_ntop bcopy bindresvport_sa clock fchown fchmod freeaddrinfo futimes gai_strerror getcwd getaddrinfo getgrouplist getnameinfo getrlimit getrusage getttyent getusershell glob inet_aton inet_ntoa inet_ntop innetgr login_getcapbool md5_crypt memmove mkdtemp on_exit openpty realpath rresvport_af setdtablesize setenv setegid seteuid setlogin setproctitle setresgid setre...

futimens() is available in Android KitKat. We work extensively with LLVM and this allowed us to remove our ugly workaround for not having this functionality. Steve On Mon, Nov 25, 2013 at 8:30 AM, James Lyon <jameslyon0 at gmail.com> wrote: > Hi, > > I've been trying to get LLVM wo...

Dear all, I am using R 3.4.3 on Windows 10. I am preparing some teaching materials and I'm having trouble matching the by-hand version with the R code. I have fitted a Cox model - let's use the ovarian data as an example: library(survival) data(ovarian) ova_mod <- coxph(Surv(futime,fustat)~age+rx,data=ovarian) If I want to make predict survival for a new set of individuals at 100 days then that is trivial using predict.coxph e.g.: newdata <- data.frame(futime=rep(100,5),fustat=rep(1,5),age=c(45,50,55,60,65),rx=c(1,2,1,2,1)) preds <- predict(ova_mod,newdata,type="e...

On 25/11/2013 18:27, Stephen Hines wrote: > futimens() is available in Android KitKat. We work extensively with > LLVM and this allowed us to remove our ugly workaround for not having > this functionality. Hi Steve, Is your work on a branch somewhere? James is putting a lot of time into this, and the LLVM community has been making best...

...mate and plot hazard funciton. However function "muhaz" always gives error message "Error in Surv(times, delta) : object "times" not found". I could not even run their sample codes in the user's manual as follows: data(ovarian) attach(ovarian) fit1 <- muhaz(futime, fustat) it gave the same error message. By the way, the "survival" package is installed and functions well on my computer. Does anyone have the same problem? Thanks! Deming