search for: crr

Hi all I do not understand why I am getting the following error message. Can anybody help me with this? Thanks in advance. install.packages("cmprsk") library(cmprsk) result1 <-crr(ftime, fstatus, cov1, failcode=1, cencode=0 ) one.pout1 = predict(result1,cov1,X=cbind(1,one.z1,one.z2)) predict.crr(result1,cov1,X=cbind(1,one.z1,one.z2)) Error: could not find function "predict.crr" [[alternative HTML version deleted]]

Dear R-help, I'm very sorry to ask 2 questions in a week. I am using the package 'crr' and it does exactly what I need it to when I use the dataset a. However, when I use dataset b I get the following error message: Error in drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) : system is computationally singular: reciprocal condition number = 1....

Hi, I need to fit model using crr, however my covariate is categorical with 3 levels. I use crr(time,status,agesplit,failcode=1,cencode=0) where agesplit is defined as <20,21-29,>30 years, so it takes 0, 1 or 2 for each patient. I hoped to get estimated coefficients for the levels 1 and 2 w.r.t level 0 as in coxph. But, I di...

...go 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 exponentially increasing with increasing number of observations. I therefore wrote a simulation for trying different approaches; check how factors, data frames and matrixes affect the performance so that I could choose the most efficient combination. I have a 1 year old...

Hi everybody, I have a vector of characters and i would like to extract certain parts. My vector is named metr_list: [1] "F:/Naval_Live_Oaks/2005/data//BE.tif" [2] "F:/Naval_Live_Oaks/2005/data//CH.tif" [3] "F:/Naval_Live_Oaks/2005/data//CRR.tif" [4] "F:/Naval_Live_Oaks/2005/data//HOME.tif" And i would like to extract BE, CH, CRR, and HOME in a different vector named "names.id" for example. I read the help files for sub and grep and the likes but i have to recognize that i did not understand it. So i've d...

Hi everybody, I have a vector of characters and i would like to extract certain parts. My vector is named metr_list: [1] "F:/Naval_Live_Oaks/2005/data//BE.tif" [2] "F:/Naval_Live_Oaks/2005/data//CH.tif" [3] "F:/Naval_Live_Oaks/2005/data//CRR.tif" [4] "F:/Naval_Live_Oaks/2005/data//HOME.tif" And i would like to extract BE, CH, CRR, and HOME in a different vector named "names.id" for example. I read the help files for sub and grep and the likes but i have to recognize that i did not understand it. So i've d...

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

Hi all, I am aware that crr lacks the "friendly" command structure of functions such as cph. All is clear to me about including covariates until I want to include a stratification term in the competing risk framework (no nice strat command). I am still a bit of a novice in R - I am looking for an example to help m...

I am trying to use the crr function in the cmprsk package to analyze a large patient dataset (45000 +), The model has 100 + covariates and 5 competing risks. I am finding that R seems to get bogged down and even if I let it run for several hours I don't get anything back. Am I expecting too much, or are there ways to spe...

...ence = 1, absence = 0) of a species using some lidar metrics (i have 4 of these). I am using different models and such .... and when i used gam i got this very weird (for me) result which i thought it is not possible - or i have no idea how to interpret it. > can3.gam <- gam(can>0~s(be)+s(crr)+s(ch)+s(home), family = 'binomial') > summary(can3.gam) Family: binomial Link function: logit Formula: can> 0 ~ s(be) + s(crr) + s(ch) + s(home) Parametric coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 85.39 162.88 0.524 0.6 Approximate s...

...isn't the case after 3 months. ----------------------------------------------- ERROR Errors in running code in vignettes: when running code in ?compete.Rnw? ... > temp$fstat <- as.numeric(event) > temp$msex <- with(temp, 1 * (sex == "M")) > fgfit1 <- with(temp, crr(etime, fstat, cov1 = cbind(age, + msex, mspike), failcode = 2, cencode = 1, variance = TRUE)) When sourcing ?compete.R?: Error: could not find function "crr" Execution halted * checking re-building of vignette outputs ... NOTE Error in re-building vignettes: ... Warning in coxp...

Dear R users, I used crr function in R package 'cmprsk' to fit a competing risks model. There were no any error or warning messages during running the function, but the output was obvious not correct. I saved the model fit as an object called f.crr. I extracted bfitj from the object by doing f.crr$bfitj, and f...

Dear all, Happy New year! I have used the 'crr' function to fit the 'proportional subdistribution hazards' regression model described in Fine and Gray (1999). dat1 is a three column dataset where: - ccr is the time to event variable - Crcens is an indicator variable equal to 0 if the event was achieved, 1 if the event wasn't ac...

...----------------- >> ERROR >> Errors in running code in vignettes: >> when running code in ?compete.Rnw? >> ... >> temp$fstat <- as.numeric(event) >> >> temp$msex <- with(temp, 1 * (sex == "M")) >> >> fgfit1 <- with(temp, crr(etime, fstat, cov1 = cbind(age, >> + msex, mspike), failcode = 2, cencode = 1, variance = TRUE)) >> >> When sourcing ?compete.R?: >> Error: could not find function "crr" >> Execution halted >> >> * checking re-building of vignette outputs ....

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

...uot;treatment" covariate I d like to have it = 0 during time interval 0 to 1 year, and equal to 1 after 1 year. Correct me if I m wrong in trying to do so. First, I have run the following script (R2.7.1 under XPpro) according to previous advices: library(cmprsk) attach(LAMrelapse) fit1<- crr(rel.t, rel.s, treatment, treatment, function(uft) cbind(ifelse(uft<=1,1,0),ifelse(uft>1,1,0)), failcode=1, cencode=0, na.action=na.omit, gtol-06, maxiter) fit1 where: rel.t = time to event (in years) rel.s = status , =1 if disease relapse, =2 if death from non disease related cause (toxicity...

Hello R users, I'm using xyplot (lattice) to plot some data. Unfortunately, the text that goes into each strip is often too big to fit. For instance, one of the factors may be labelled "Satisfaction with waiting time" which is too wide to fit. If I use split.string() to wrap it onto two lines it is then too high to fit in the strip. I'd prefer to apply split.string on the

...ely one editor of a journal where i would like to publish some results insists in giving p-values for the squared deviance i get out from different glm and gam models. I came up with this solution, but sincerely i would like to get yours'all opinion on the matter. p1.glm <- glm(count ~be+ch+crr+home, family = 'poisson') # count - is count of species (vegetation) # be, ch, crr, home - different lidar metrics # calculating d^2 d2.p1 <- round((p1.glm[[12]]-p1.glm[[10]])/p1.glm[[12]],4) d2.p1 0.6705 # calculating f statistics with N = 148 and n=4; f = (N-n-1)/(N-1)(1-d^2) f &l...

Hi, I'm using crrstep package to do stepwise covariate selection for the Fine & Gray competing risks regression model. However, I keep getting an error (please see below). Please help!! > PHstep <- crrstep(years~1+var1+var2+var3+var4+var5,scope.min=~1,censorcmprsk, data=crisk, direction=c("forward&qu...

...eatment" covariate I d like to have it = 0 during time interval 0 to 1 year, and equal to 1 after 1 year. Correct me if I m wrong in trying to do so. First, I have run the following script (R2.7.1 under XPpro) according to previous advices: library(cmprsk) attach(LAMrelapse) fit1<- crr(rel.t, rel.s, treatment, treatment, function(uft) cbind(ifelse(uft<=1,1,0),ifelse(uft>1,1,0)), failcode=1, cencode=0, na.action=na.omit, gtol-06, maxiter) fit1 where: rel.t = time to event (in years) rel.s = status , =1 if disease relapse, =2 if death from non disease related cause (tox...