similar to: Suppressing a warning from library()

Dear R users I would like to calculate the Cumulative incidence for an event adjusting for competing risks and adjusting for covariates. One way to do this in R is to use the cmprsk package, function crr. This uses the Fine & Gray regression model. However, a simpler and more classical approach would be to implement the Kalbfleisch & Prentice method (1980, p 169), where one fits cause

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

Dear all, I was wondering if anyone can help? I am an R user but recently I have resorted to SAS to calculate the probability of the event (and the associated confidence interval) for the Cox model with combinations of risk factors. For example, suppose I have a Cox model with two binary variables, one for gender and one for treatment, I wish to calculate the probability of survival for the

Hi All, I'm trying to figure out the cumulative incidence curve in R in some limited time. I found in package "cmprsk", the command "plot.cuminc" can get this curve. But I noticed that there is no mark for the censored time there, comparing with the KM curve by "plot.survfit". Here are my codes (attached is the data): ----------------

------------------------------ On Sat, May 16, 2015 8:04 AM BST Uwe Ligges wrote: >Not sure why this goes to R-devel. You just could have asked the >maintainer. Terry Therneau is aware of it and promised he will fix it. > The quickest fix is to add cmprsk to the recommended list , and that's is an R-devel issue. >On 16.05.2015 07:22, Hin-Tak Leung wrote: >> 'make

'make check-all' for current R has been showing this error in the middle for a few months now - any thought on fixing this? I think cmprsk should be either included in the recommended bundle, or the survival vignette to not depend on it. Having 'make check-all' showing glaring ERROR's for a few months seems to defeat the purpose of doing any checking at all via 'make

Dear R-users Has anybody implemented a function/package that will compute an individual's risk of an event in the presence of competing risks, adjusted for the individual's covariates? The only thing that seems to come close is the cuminc function from cmprsk package, but I would like to adjust for more than one covariate (it allows you to stratify by a single grouping vector). Any

Dear Prof Gray and everyone, As our package developers discussed about incompatibility between Design and survival packages, I faced another problem with cmprsk- a survival dependent packacge. The problem is exactly similar to what happened to the Design package that when I just started running cuminc function, R was suddenly closed. These incidents suggest that maybe many other survival

Dear everyone, My question involves the use of the survival object. We can have Surv(time,time2,event, type=, origin = 0) (1) As detailed on p.65 of: http://cran.r-project.org/web/packages/survival/survival.pdf My data (used in my study) is 'right censored' i.e. my variable corresponding to 'event' indicates whether a person is alive (0) or dead (1) at date last seen

> Dear All, > > I have built a survival cox-model, which includes a covariate * time interaction. (non-proportionality detected) > I am now wondering how could I most easily get survival predictions from my model. > > My model was specified: > coxph(formula = Surv(event_time_mod, event_indicator_mod) ~ Sex + > ageC + HHcat_alt + Main_Branch + Acute_seizure +

> On Jun 1, 2019, at 12:59 PM, Peter Langfelder <peter.langfelder at gmail.com> wrote: > > On Sat, Jun 1, 2019 at 3:22 AM Therneau, Terry M., Ph.D. via R-devel > <r-devel at r-project.org> wrote: >> >> In the next version of the survival package I intend to make a non-upwardly compatable >> change to the survfit object. With over 600 dependent packages

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

Dear R users, I d like to assess the effect of "treatment" covariate on a disease relapse risk with the package cmprsk. However, the effect of this covariate on survival is time-dependent (assessed with cox.zph): no significant effect during the first year of follow-up, then after 1 year a favorable effect is observed on survival (step function might be the correct way to say that ?).

On the local machine the command R11 CMD Rd2pdf survfit.Rd works fine. R12 CMD Rd2pdf survfit.Rd fails with the message below. Converting Rd files to LaTeX ... survfit.Rd Creating pdf output from LaTeX ... Error in texi2dvi("Rd2.tex", pdf = (out_ext == "pdf"), quiet = FALSE, : Running 'texi2dvi' on 'Rd2.tex' failed. Messages: sh: texi2dvi: command not

This works ok: > cox = coxph(surv ~ bucket*(today + accor + both) + activity, data = data) > fit = survfit(cox, newdata=data[1:100,]) but using strata leads to problems: > cox.s = coxph(surv ~ bucket*(today + accor + both) + strata(activity), > data = data) > fit.s = survfit(cox.s, newdata=data[1:100,]) Error in model.frame.default(data = data[1:100, ], formula = ~bucket + :

Hello, I have a question regarding competing risk regression using cmprsk package (function crr()). I am using R2.9.1. How can I do to assess the effect of qualitative predictor (gg) with more than two categories (a,b,c) categorie c is the reference category. See above results, gg is considered like a ordered predictor ! Thank you for your help Jan > # simulated data to test > set.seed(10)

Is the phrase "index <- match(x, sort(unique(x)))" reliable, in the sense that it will never return NA? Context: Calculation of survival curves involves the concept of unique death times. I've had reported cases in the past where survfit failed, and it was due to the fact that two "differ by machine precision" values would sometimes match and sometimes not,

A problem has been pointed out by a French user of the survival package and I'm looking for a pointer. > options(OutDec= ",") > fit <- survfit(Surv(1:6 /2) ~ 1) > fit$time [1] NA 1 NA 2 NA 3 A year or two ago some test cases that broke survfit were presented to me. The heart of the problem was numbers that were almost identical, where table(x) and unique(x) gave

First a statistical issue: The survfit routine will produce predicted survival curves for any requested combination of the covariates in the original model. This is not the same thing as an "adjusted" survival curve. Confusion on this is prevalent, however. True adjustment requires a population average over the confounding factors and is closely related to the standardized

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),