similar to: Need help interpreting output from rcorrp.cens with Cox regression

I am reposting my message from April 8th because I never received a response to the original post: Dear R-listers, I am an MD and clinical epidemiologist developing a measure of comorbidity severity for patients with liver disease. Having developed my comorbidity score as the linear predictor from a Cox regression model I want to compare the discriminative ability of my comorbidity measure with

I have two alternative Cox models with C-statistics 0.72 and 0.78. My question is if 0.78 is significantly greater than 0.72. I'm using rcorrp.cens. I cannot find the U statistics in the output of the function. This is the output of the help example: > x1 <- rnorm(400) > x2 <- x1 + rnorm(400) > d.time <- rexp(400) + (x1 - min(x1)) > cens <- runif(400,.5,2) > death

Hi R-users, I'm having some problems in using the Hmisc package. I'm estimating a cox ph model and want to test whether the drop in concordance index due to omitting one covariate is significant. I think (but I'm not sure) here are two ways to do that: 1) predict two cox model (the full model and model without the covariate of interest) and estimate the concordance index (i.e. area

Dear, I wanna to compare AUC generated by two distribution models using the same sample. I tried improveProb function's example code below. set.seed(1) library(survival) x1 <- rnorm(400) x2 <- x1 + rnorm(400) d.time <- rexp(400) + (x1 - min(x1)) cens <- runif(400,.5,2) death <- d.time <= cens d.time <- pmin(d.time, cens) rcorrp.cens(x1, x2, Surv(d.time, death))

Dear R-help, Would it be appropriate to do the following to calculate a p-value for the difference between c-ind of x1 and c-inx of x2 using the output from rcorrp.senc() > r<-rcorrp.senc(x1,x1,y) > pValue<-1-pnorm((r[11]-r[12])/(r[2]/r[5])*1.96) Osman O. Al-Radi, MD, MSc, FRCSC Chief Resident, Cardiac Surgery University of Toronto, Canada

If there are any other users who use AHRQ's SAS code comoanaly2010 and comformat2010 to create comorbidity variables, I thought you might be interested in the following PRELIM code we wrote to mimic its functionality in R. It seems to yield similar results, but may contain errors. Please feel free to comment (kindly) or enhance. I'm sure there are better ways to skin this cat, but we at

Hi, I'm trying to calculate 95% confidence interval of C statistic of logistic regression model using rcorr.cens in rms package. I wrote a brief function for this purpose as the followings; CstatisticCI <- function(x) # x is object of rcorr.cens. { se <- x["S.D."]/sqrt(x["n"]) Low95 <- x["C Index"] - 1.96*se Upper95 <- x["C

Dear R-help, This is an example in the {Hmisc} manual under rcorr.cens function: > set.seed(1) > x <- round(rnorm(200)) > y <- rnorm(200) > round(rcorr.cens(x, y, outx=F),4) C Index Dxy S.D. n missing uncensored Relevant Pairs Concordant Uncertain 0.4831 -0.0338 0.0462 200.0000

Dear all, I am interested to evaluate reclassification using net reclassification improvement and Integrated Discrimination Index IDI after survival analysis (Cox proportional hazards using stcox). I search a R package or a R code that specifically addresses the categorical NRI for time-to-event data in the presence of censored observation and, if possible, at different follow-up time points. I

Dear R Group: Based on the following data, how to do a great Spider (Radar) Plot? Any advice is greatly appreciated. HospID Rate Age Charlson NIHSS 1 0.2 49 3.5 0 2 0.1 48 1.8 12 3 0.4 56 2.1 5 4 0.3 77 0 7 5 0.2 67 6.5 3 6 0.1 62 4.8 4.6 7 0.1 64 12 5.2 8 0.3 61 3 2.8 9 0.15 69 4.5 1.9 10 0.22 80 0 6.7 11 0.34 61 6 4.2 12 0.18 63 3 6.1 13 0.09 64 8 15 12 0 56 10 11 15 0.1 70 11 7 16

Hi all, I have a difficulty to calculate the PCA scores. The PCA scores I calculated doesn't match with the scores generated by R, mypca<-princomp(mymatrix, cor=T) myscore<-as.matrix(mymatrix)%*%as.matrix(mypca$loadings) Does anybody know how the mypca$scores were calculated? Is my formula not correct? Thanks a lot! Phoebe [[alternative HTML version deleted]]

Dear R users, For the purpose of validating a prediction model using validate() from the rms package, I am running into some trouble with using the fastbw() function breaking up natural groups of variables. Is there any way I can specify to keep certain variable together? In particular, if interactions are included I would also like to keep the main effects in the model. Another example is a

Hi, I am trying to calculate net reclassification improvement (NRI) and Inegrated Discrimination Improvement (IDI) for a survival dataset to compare 2 risk models. It seems that the improveProb() in Hmisc package does this only for binary outcome, while rcorrp.cens() does take survival object, but doesn't output NRI or IDI. Can anyone suggest any other packages that can calculate NRI and IDI

Does anyone know of an R-function or method to compare two C-statistics (Harrells's C - rcorr.cens) obtained from 2 different models in partially paired datasets (i.e. some similar and some different cases), with one continuous independent variable in each separate model? (in a survival analysis context)? I have noticed that the rcorrp.cens function can be used for paired data. Thanks

Greetings, I'm learning R and I'm stuck on a basic concept: how to specify a logical condition once and then perform multiple transformations under that condition. The program below is simplified to demonstrate the goal. Its results are exactly what I want, but I would like to check the logical state of gender only once and create both (or any number of) scores at once.

I downloaded data from http://www.hcp.med.harvard.edu/ncs/ Which provides data in DTA (STATA), XPT (SAS), and POR (SPSS) formats all of which I have tried to read with the foreign package but I am not able to load any of them. I have 2 gb of RAM, but R crashes when the memory gets just over 1 GB. I am using Windows version 2.1.1. The size of the DTA file is 48 MB; the xpt file is 188

Hello R Community, I'm continuing to work through logistic regression (thanks for all the help on score test) and have come up against a new opposition. I'm trying to compute Somers Dyx as some suggest this is the preferred method to Somers Dxy (Demaris, 1992). I have searchered the [R] archieves to no avail for a function or code to compute Dyx (not Dxy). The overview of Hmisc has

Hi everyone, I've been bashing my head against this for days now, and can't figure out what to do. I have the following dataframe header appetitive stimulus aversive stimulus chaining contingency discriminative stimulus extinction intermittent reinforcement negative reinforcer operant response place learning positive reinforcer punishment reinforcement

Dear R-helpers, I recently switched from SAS to R, in order to model the occurrence of rare events through logistic regression. Is there a package available in R to calculate the Goodman-Kruskal Gamma? After searching a bit I found a function "rcorr.cens" which should do the job, but it is not clear to me how to define the input vectors? Is "x" a vector with the fitted

I am using the rcorr.cens function from the Hmisc package and the time-dependent ROC curve obtained using tdrocc in the survcomp package. I understand that the C statistic from rcorr.cens has to be subtracted from 1 if high values of the risk variable lower survival. Given that I wonder what the connection is between that C statistic and the AUC from the tdrocc object. If they are substantially