Displaying 20 results from an estimated 10000 matches similar to: "survfit is too slow! Looking for an alternative"
2011 Nov 04
1
survfit function?
Hi, I
am working on fitting a proportional hazard model to predict the
probability of default for mortgage loans. I have a question regarding
survfit function. My
historical data set is a pool of loans with monthly observed default
status for the next 24 months. The data is left truncated (delayed entry
to observation window after the loan is opened) and right censored. I
would like to
2006 Oct 25
1
Incorrect 'n' returned by survfit()
I've a data set with 60000 rows of data representing 6000+ distinct loans. I did a coxph() regression on it (see call below), but a subsequent survfit() call on the coxph object is almost certainly wrong. It gives n=6 when it should be
more like 6000+ (I think)
> survfit(resultag)
Call: survfit.coxph(object = resultag)
n events median 0.95LCL 0.95UCL
6 489 Inf
2009 Feb 25
3
survival::survfit,plot.survfit
I am confused when trying the function survfit.
my question is: what does the survival curve given by plot.survfit mean?
is it the survival curve with different covariates at different points?
or just the baseline survival curve?
for example, I run the following code and get the survival curve
####
library(survival)
fit<-coxph(Surv(futime,fustat)~resid.ds+rx+ecog.ps,data=ovarian)
2011 Oct 01
4
Is the output of survfit.coxph survival or baseline survival?
Dear all,
I am confused with the output of survfit.coxph.
Someone said that the survival given by summary(survfit.coxph) is the
baseline survival S_0, but some said that is the survival S=S_0^exp{beta*x}.
Which one is correct?
By the way, if I use "newdata=" in the survfit, does that mean the survival
is estimated by the value of covariates in the new data frame?
Thank you very much!
2007 Sep 27
2
center option of basehaz in survfit
I have a very general question about what the centering option in basehaz does to factors. (basehaz computes the baseline cumulative hazard for a coxph object using the Breslow estimator).
Lets say I'm interested in a survival model with two (dichotomous) factors and a continuous covariate.
Variable Possible Values
Factor1 0 or 1
Factor2 0 or 1
2007 Nov 13
2
plotting coxph results using survfit() function
i want to make survival plots for a coxph object using survfit
function. mod.phm is an object of coxph class which calculated results
using columns X and Y from the DataFrame. Both X and Y are
categorical. I want survival plots which shows a single line for each
of the categories of X i.e. '4' and 'C'. I am getting the following
error:
> attach(DataFrame)
>
2006 Dec 29
2
Survfit with a coxph object
I am fitting a coxph model on a large dataset (approx 100,000 patients), and
then trying to estimate the survival curves for several new patients based
on the coxph object using survfit. When I run coxph I get the coxph object
back fairly quickly however when I try to run survfit it does not come
back. I am wondering if their is a more efficient way to get predicted
survival curves from a coxph
2010 Nov 12
3
predict.coxph
Since I read the list in digest form (and was out ill yesterday) I'm
late to the discussion.
There are 3 steps for predicting survival, using a Cox model:
1. Fit the data
fit <- coxph(Surv(time, status) ~ age + ph.ecog, data=lung)
The biggest question to answer here is what covariates you wish to base
the prediction on. There is the usual tradeoff between too few (leave
out something
2012 Oct 13
4
Problems with coxph and survfit in a stratified model with interactions
I?m trying to set up proportional hazard model that is stratified with
respect to covariate 1 and has an interaction between covariate 1 and
another variable, covariate 2. Both variables are categorical. In the
following, I try to illustrate the two problems that I?ve encountered, using
the lung dataset.
The first problem is the warning:
To me, it seems that there are too many dummies
2013 Jan 31
1
obtainl survival curves for single strata
Dear useRs,
What is the syntax to obtain survival curves for single strata on many subjects?
I have a model based on Surv(time,response) object, so there is a single row per subject and no start,stop and no switching of strata.
The newdata has many subjects and each subject has a strata and the survival based on the subject risk and the subject strata is needed.
If I do
newpred <-
2012 Nov 27
4
Fitting and plotting a coxph with survfit, package(surv)
Hi Dear R-users
I have a database with 18000 observations and 20 variables. I am running
cox regression on five variables and trying to use survfit to plot the
survival based on a specific variable without success.
Lets say I have the following coxph:
>library(survival)
>fit <- coxph(Surv(futime, fustat) ~ age + rx, data = ovarian)
>fit
what I am trying to do is plot a survival
2005 Aug 17
1
GLM/GAM and unobserved heterogeneity
Hello,
I'm interested in correcting for and measuring unobserved
heterogeneity ("missing variables") using R. In particular, I'm
searching for a simple way to measure the amount of unobserved
heterogeneity remaining in a series of increasingly complex models
(adding additional variables to each new model) on the same data.
I have a static database of 400,000 or
2010 Jun 23
1
Probabilities from survfit.coxph:
Hello:
In the example below (or for a censored data) using survfit.coxph, can
anyone point me to a link or a pdf as to how the probabilities appearing in
bold under "summary(pred$surv)" are calculated? Do these represent
acumulative probability distribution in time (not including censored time)?
Thanks very much,
parmee
*fit <- coxph(Surv(futime, fustat) ~ age, data = ovarian)*
2010 Sep 23
2
extending survival curves past the last event using plot.survfit
Hello,
I'm using plot.survfit to plot cumulative incidence of an event.
Essentially, my code boils down to:
cox <-coxph(Surv(EVINF,STATUS) ~ strata(TREAT) + covariates, data=dat)
surv <- survfit(cox)
plot(surv,mark.time=F,fun="event")
Follow-up time extends to 54 weeks, but the last event occurs at week
30, and no more people are censored in between. Is there a
2009 Dec 22
1
Slow survfit -- is there a faster alternative?
Using R 2.10 on Windows:
I have a filtered database of 650k event observations in a data frame
with 20+ variables.
I'd like to be able to quickly generate estimate and plot survival
curves. However the survfit and cph() functions are extremely slow.
As an example: I tried
results.cox<-coxph(Surv(duration, success) ~ start_time + factor1+
factor2+ variable3, data=filteredData) #(took a
2009 May 10
2
plot(survfit(fitCox)) graph shows one line - should show two
R 2.8.1
Windows XP
I am trying to plot the results of a coxph using plot(survfit()). The plot should, I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong?
My code is reproduced below, my figure is attached to this EMail message.
John
> #Create simple survival object
>
2006 May 24
1
multiple destinations in duration (survival) analysis
Hi,
I'm trying to estimate a (parametric) competing risks model in the context
of duration (or survival, if you wish) analysis. That is, instead of
studying the transition of subjects to "death", I wish to study the
transitions to multiple *destinations* (which is different from studying
multiple *durations*, or recurrent events). I am more interested in the
hazard function rather
2012 Jan 24
1
Plotting coxph survival curves
Hi,
I am attempting to plot survival curves estimated by cox proportional
hazards regression model. The formula for the model is this:
F.cox.weight <- coxph(Surv(Lifespan, Status) ~ MS + Weight + Laid + MS:Laid
+ Weight:Laid, data = LongF)
MS = Mating status (mated/virgin)
Weight = adult female weight, continuous covariate
Laid = number of eggs laid by each female, continuous covariate
I
2011 Feb 03
3
coxph fails to survfit
I have a model with quant vars only and the error message does not make sense:
(mod1 <- coxph(Surv(time=strt,time2=stp,event=(resp==1))~ +incpost+I(amt/1e5)+rate+strata(termfac),
subset=dt<"2010-08-30", data=inc,method="efron"))
Call:
coxph(formula = Surv(time = strt, time2 = stp, event = (resp ==
1)) ~ +incpost + I(amt/1e+05) + rate + strata(termfac),
2009 Feb 25
3
survival::predict.coxph
Hi,
if I got it right then the survival-time we expect for a subject is the
integral over the specific survival-function of the subject from 0 to t_max.
If I have a trained cox-model and want to make a prediction of the
survival-time for a new subject I could use
survfit(coxmodel, newdata=newSubject) to estimate a new
survival-function which I have to integrate thereafter.
Actually I thought