Displaying 20 results from an estimated 4000 matches similar to: "Survival-Analysis: How to get numerical values from survfit (and not just a plot)?"
2006 Jul 11
2
new object
Um texto embutido e sem conjunto de caracteres especificado associado...
Nome: n?o dispon?vel
Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060711/1da40382/attachment.pl
2006 Jun 10
1
Calculating survival for set time intervals
Hello friends and fellow R users,
I have successfully tabulated and entered my survival data into R and
have generated survival curves. But I would like to be able to determine
what the survival rates are now at one month, three months, six months
and one year.
I have a data set, via.wall, which I have entered into R, and which
generates the following Surv object:
Surv(Days,Status==1)
[1]
2008 Dec 06
1
Kaplan-Meier function from survfit
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 |
2010 Sep 10
2
survfit question
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)
2007 Oct 15
2
some question about partial prediction in survival
Hi there:
i got a problem to get the prediction from a model recently. for
example if i use a survival analysis to predict the risk. i use the code
like below: i found the the prediction is not equal to (coef * x + coef
* sex) , could someone help me with why this happened? and can someone
explain to me how this command "predict(f, type="terms")" works? is
every partial
2006 May 24
2
median of a survfit object
Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?...
Nom : non disponible
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060524/312404be/attachment.pl
2008 Mar 13
1
survival curve for only certain values of a factor
Hello:
Using the built-in dataset aml as an example:
data(aml)
If I use instead dummy variables:
aml$x1 = (aml$x=="maintained")aml$x2 = (aml$x=="unmaintained")
and I want to plot the survival curve using x1, x2, and I just want the 2 levels, rather than 4 curves from:
fit <- survfit(Surv(time, status) ~ x1+x2, data=aml)
plot(fit)
I guess because there are 2 levels
2004 Oct 05
2
Nelson-Aalen estimator in R
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
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
2006 Jan 20
3
command in survival package
Hi there,
I have a question about one command sentence when I follow the example
in the book of "Survival analysis in S":
> aml1<-aml[aml$group==1]
but I got the error warning: NULL data frame with 23 rows
Thus, I couldn't keep going on the next command:
esf.fit<-survfit(Surv(aml1,status)~1).
and also when I try
> aml1<-aml[aml$group==1,]
2015 Dec 07
3
Tiempo de vida
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]]
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)
2010 Dec 27
1
Problem using pkg "survival"
Hello all.
I've been attempting to utilize the "survival" pkg (
http://cran.r-project.org/web/packages/survival/index.html), while reading
through this guide (http://www.ms.uky.edu/~mai/Rsurv.pdf). I figured working
through the guide would be the best way to go, before attempting my own
data.
I tried to utilize the Kaplain-Meier estimator as shown in the guide:
2015 Dec 07
2
Tiempo de vida
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
2004 Apr 21
1
Boot package
Dear mailing list,
I tried to run the example for the conditional bootstap written in the help file
of censboot. I got the following result:
STRATIFIED CONDITIONAL BOOTSTRAP FOR CENSORED DATA
Call:
censboot(data = aml, statistic = aml.fun, R = 499, F.surv = aml.s1,
G.surv = aml.s2, strata = aml$group, sim = "cond")
Bootstrap Statistics :
original bias std. error
t1*
2006 May 05
2
How to access results of survival analysis
Hi List,
A friend of mine recently asked the same question as Heinz T?chler. Since
I've already written the code I'd like to share with the list.
# x is an object returned by "survfit";
# "smed" returns a matrix of 5 columns of
# n, events, median, 0.95LCL, 0.95UCL.
# The matrix returned has rownames as the
# group labels (eg., treatment arms) if any.
smed <-
2015 Dec 08
2
Tiempo de vida
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
2009 Feb 20
1
log-minus-log plot
Dear experts,
I would like to know how to plot the log-minus-log plot for survival
analysis (to check the proportional assumption) in R.
Using the AML example.
fit <- survfit(Surv(time, status) ~ x, data=aml)
length(fit$surv) #20
as the length of fit$surv is shorter than aml$x and aml$time. I don't
know how to plot.
Thank you.
Regards,
CH
--
CH Chan
Research Assistant - KWH
2020 Sep 29
5
2 KM curves on the same plot
Hello,
Can anyone suggest a simple way to generate a Kaplan-Meier plot with 2 survfit objects, just like this one:?
https://drive.google.com/file/d/1fEcpdIdE2xYtA6LBQN9ck3JkL6-goabX/view?usp=sharing
Suppose I have 2 survfit objects: fit1 is for the curve on the left (survtime has been truncated to the cutoff line: year 5), fit2 is for the curve on the right (minimum survival time is at the
2011 Apr 18
1
time dependent hazard ratios
Hi, I am new to time-dependent Cox model to estimate time dependent hazard
ratios. Let me use aml dataset from survival package:
> aml3<-survSplit(aml2,cut=c(5,10,20),end="time",start="start",
event="status",episode="i")
If I want to esimate hazard ratio for each of the time intervals 0-5, 5-10,
10-20 and >=20, would the following calculate