similar to: Packages/functions for competing risk analysis

Hello again, I m trying to use timereg package as you suggested (R2.7.1 on XP Pro). here is my script based on the example from timereg for a fine & gray model in which relt = time to event, rels = status 0/1/2 2=competing, 1=event of interest, 0=censored random = covariate I want to test library(timereg) rel<-read.csv("relapse2.csv", header = TRUE, sep = ",",

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

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 ?).

Hi all, I've been reading and using the information from the list for some time but this is my first question here. English is not my primary language, so sorry in advance for any language mistakes. :) I'm working with the "timereg" package to analize survival data. I want to perform a multivariate analisis of clinical information similar to the Cox regression but taking

useR's, I am trying to do a power calculation for a survival analysis using a logrank test and I need some help properly doing this in R. Here is the information that I know: - I have 2 groups, namely HG and LG - Retrospective analysis with subjects gathered from archival data over 20 years. No new recruitment of subjects and no estimated time to target accrual and accrual rate. - Survival

I'm sorry for the really vague subject line but I am not sure how to succinctly describe what I am doing and what the problem is. But, here goes: 1. I have data with two-way data with frequencies. Below is an example, though in reality I am looking at about 10 different variables that I am crossing so the values of X1 and X2 change. X1 and X2 are place holders. Here's the dataset

Hi all, anyone has an idea how I could fix this? \donttest{ ## Set colors from colorspace package with a fallback col <- try(colorspace::rainbow_hcl(5), silent = TRUE) %||% rainbow(5) } The problem is that this makes R CMD check freak out (http://www.r-project.org/nosvn/R.check/r-devel-linux-x86_64-debian-clang/falsy-00check.html) if the colorspace package is not declared as a dependency.

Hi R-Helpers, I don't think I need to post a dataset for this question but if I do, I can. Anyway, I am having a lot of trouble with the ifelse command. Here is my code: vn$PM.DIST_flag <- ifelse( (vn$PM.EXP > 0.0) & (vn$PM.DIST.TOT != 1.0), 1, 0 ) And here is my output that doesn't make ANY sense: PM.EXP PM.DIST.TOT PM.DIST_flag 0 0 0 0 0 0 0 0 0 177502 1 0 31403

Hi All, I have some data where I am doing fairly simple calculations, nothing more than adding, subtracting, multiplying and dividing. I’m running into a problem when I divide one variable by another and when they’re both 0 I get NaN. I realize that if you divide a non-zero by 0 then you get Inf, which is, of course, correct. But in my case I never get Inf, just NaN because of the structure

Dear R Users, Do you have a sample code for developing a nomogram with competing-risks? Any help is appreciated. Kind regards, ND Nguyen

Hi, I posted this question on stats.stackexchange.com 3 days ago but the answer didn't really address my question concerning the speed in competing risk regression. I hope you don't mind me asking it in this forum: I?m doing a registry based study with almost 200 000 observations and I want to perform a competing risk analysis. My problem is that the crr() in the cmprsk package is

Dear all, I'm interested in fitting survival trees with competing risk analysis. The tree should show the cumulative incidence function for each terminal node . I read several paper illustrating this possibility, but to the best of my knowledge no R code are reported. There is any R package for this fit? Thank you very much in advance. Mario Petretta

Dear R users, is there a way, I mean a package, to perform a competing risk model which can handle time dependent covariates ? my main covariate (additional treatment to patients) appears not to follow the proportional hazards assumption, its effect being observed after one year of treatment but not before (this is expected / makes sense on a clinical point of view). SO I was planning to use a

I have a competing risk data where a patient may die from either AIDS or Cancer. I want to compare the cox model for each of the event of interest with a competing risk model. In the competing risk model the cumulative incidence function is used directly. I used the jackknife (pseudovalue) of the cumulative incidence function for each cause (AIDS or Cancer) in a generalized estimating equation. I

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

This message was also sent to the MEDSTATS mailing list, so here is the reply I posted to that: Philippe, The machinery to use is to split follow-up time so finely that you can safely assume that rates are constant in each interval, and then just stuff it all into a Poisson model. This allows you to use any kind of time-dependent variables as well as accommodating competing risks. In the Epi

Dear R-helpers: I am trying to use aalen{timereg} to run a model and found that it generates the following message. the data set I used is > dim(indf) [1] 13743 10 #R code: fit_aalen <- aalen(Surv(durdxdeathp, die) ~ 1+factor(gender), indf) #start of error message below Non unique survival times: break ties ! Break ties yourself Error in aalenBase(times, ldata, X, status, id, clusters,

Hi, I am trying to automate the creation of tables for some simply analyses. There are lots and lots of tables, thus the creation of a user-defined function to make and output them to excel. My problem is that some of the analyses have convergence issues, which I want captured and included in the output so the folks looking at them know how to view those estimates. I am successfully able to do

Dear all, I am trying to use the crr function in the cmprsk package version 2.2 to analyse 198 observations.I have receive the error in solve.default. Can anyone give me some insights into where the problem is? Thanks here is my script : cov=cbind(x1,x2) z<-crr(ftime,fstatus,cov)) and data file: x1 x2 fstatus ftime 0 .02 1 263 0 .03 1 113 0 .03 1 523

Hello, I have a problem regarding calculation of Cumulative Incidence Function. The event of interest is failure of bone-marrow transplantation, which may occur due to relapse or death in remission. The data set that I have consists of- lifetime variable, two indicator variables-one for relapse and one for death in remission, and the other variables are donor type (having 3 categories), disease