Displaying 20 results from an estimated 1000 matches similar to: "c-statiscs 95% CI for cox regression model"
2010 Mar 12
3
How to format dates (with no century)?
Hi,
I have dates in this kind of format (day, month, year):
> dput(head(dates, 10))
c("6.4.7", "29.12.98", "19.10.91", "20.6.92", "2.9.3", "23.6.3",
"13.7.93", "23.3.7", "26.6.95", "15.2.10")
So, as you can see, there is no century. How can I change this
character data into dates? Any help
2010 Aug 15
1
calibration plot labels
Dear all,
when i do the calibration plot, i put the x label y label , there is
some labels are i did not put it , like "resample optimism added ..." i
want to get rid of the these label , is any body know how can i get rid of these label.
these are the following command i used
cal <- calibrate(f, u=12, method=c("boot"), B=100,m=70, data=a1)
2010 Apr 08
2
C-index and Cox model
Dear all R users,
I am building a Cox PH model on a small dataset. I am wondering how to
measure the predictive power of my cox model? Normally the ROC curve or Gini
value are used in logistic regression model. Is there any similar
measurement suitable for Cox model?
Also if I use C-index statistic to measure the predictive power, is it a
time-dependent value (i.e. do I need to calculate it for
2011 May 22
1
How to calculate confidence interval of C statistic by rcorr.cens
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
2008 Mar 03
1
Cox model+ROCR
Dear list,
I am trying to build a cox model and then perform ROC analysis in order to
retrieve some genes that are correlated with breast cancer. When I calculate
the hazard score taking into account different numbers of genes and their
coefficients ( I am trying to find the pest predictor number of genes), I
retrieve from around 1 values (for few genes included ) to size of e+80
values (for many
2011 Mar 01
1
which does the "S.D." returned by {Hmisc} rcorr.cens measure?
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
2009 Mar 09
1
rcorr.cens Goodman-Kruskal gamma
Dear r-helpers!
I want to classify my vegetation data with hierachical cluster analysis.
My Dataset consist of Abundance-Values (Braun-Blanquet ordinal scale; ranked) for each plant species and relev?.
I found a lot of r-packages dealing with cluster analysis, but none of them is able to calculate a distance measure for ranked data.
Podani recommends the use of Goodman and Kruskals' Gamma for
2010 Jul 26
1
After writing data in MMF using SEXP structure, can i reference in R?
Hi all,
After writing data in MMF(Memory Map File) using SEXP structure, can i
reference in R?
If input data is larger than 2GB, Can i reference MMF Data in R?
my work environment :
R version : 2.11.1
OS : WinXP Pro sp3
Thanks and best regards.
Park, Young-Ju
from Korea.
---------[ ???????? ???????? ???????? ]----------
???????? : R-help Digest, Vol 89,
2007 Dec 19
1
using rcorr.cens for Goodman Kruskal gamma
Dear List,
I would like to calculate the Goodman-Kruskal gamma for the predicted
classes obtained from an ordinal regression model using lrm in the Design
package. I couldn't find a way to get gamma for predicted values in Design
so have found previous positings suggesting to use :
Rcorr.cens(x, S outx = TRUE) in the Hmisc package
My question is, will this work for predicted vs
2004 Jun 04
1
use of "rcorr.cens" with binary response?
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
2013 Jan 24
4
Difference between R and SAS in Corcordance index in ordinal logistic regression
lrm does some binning to make the calculations faster. The exact calculation
is obtained by running
f <- lrm(...)
rcorr.cens(predict(f), DA), which results in:
C Index Dxy S.D. n missing
0.96814404 0.93628809 0.03808336 32.00000000 0.00000000
uncensored Relevant Pairs Concordant Uncertain
32.00000000
2013 Sep 27
1
Problems when moving to Imports from Depends
Hi all,
one of my packages uses the rcorr.cens function from the Hmisc
package. Until now I have simply put the Hmisc package into Depends:,
but prodded on by new CRAN requirements, I tried to moving it into
Imports:. However, this fails because rcorr.cens calls the function
is.Surv from survival, which does not seem to be on the search path
when Hmisc is "imported from" rather then
2011 Jun 13
1
Somers Dyx
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
2004 Jul 19
2
Evaluating the Yield of Medical Tests
Hello,
I'm a biostatistician in Toronto. I would like to know if there is
anything in survival analysis developed in R for the method "Evaluating
the Yield of Medical Test" (JAMA. May 14,1982--Vol 247, No.18 Frank E.
Harrell, Jr,PhD; Robert M. Califf, MD; David B. Pryor, MD;Kerry L.Lee,
PhD; Robert A. Rosait,MD.)
Hope to hear from you and thanks
Lisa Wang, MSc
Project Organiser
2003 Mar 11
1
Goodman / Kruskal gamma
The Goodman/Kruskal gamma is a nice descriptive rank-order
correlation statistic, often used in psychology. It is nice
because it is easy to understand. It takes all pairs of values
of each variable and asks whether they are congruent (S+ is the
number in the same order for both variables) or discordant (S-,
opposite ranking). The statistic is (S+ - S-)/(S+ + S-). It is
like tau except for the
2005 Jul 11
1
validation, calibration and Design
Hi R experts,
I am trying to do a prognostic model validation study, using cancer
survival data. There are 2 data sets - 1500 cases used to develop a
nomogram, and another of 800 cases used as an independent validation
cohort. I have validated the nomogram in the original data (easy with
the Design tools), and then want to show that it also has good results
with the independent data using 60
2012 Sep 12
2
How to resolve the following error: (list) object cannot be coerced to type 'double'
Hello,
I'd like to test a significance of two variables in their correlation using
rcorr, which gave me an error of format incompatibility.
Below are the lines that I typed in the R window and let me know if anyone
knows how to resolve this.
Var=read.csv("03apr10ab_corr_matrix_in_overlaps.csv",header=F)
colnames(Var)=c("D Prime","T statistics")
D Prime
2012 Sep 12
1
digit precision in p value of rcorr
Hi all,
Sorry about posting a really novice question.
I was able to run rcorr after converting the list to a matrix by your help.
I'm though wondering if there is any way to find out an exact p value as
the output only gave me 0 for P value as shown below.
I've added options(digits=10), which doesn't seem to help at all. Any help
would be appreciated.
P
D Prime T
2011 Aug 19
1
Hmisc::rcorr on a 'data.frame'?
Dear all
?Hmisc::rcorr states that it takes as main argument "a numeric
matrix". But is it normal that it fails in such an ugly way on a data
frame? (See below.) If the function didn't attempt any conversion to a
matrix, I would have expected it to state that in the error message
that it didn't accept 'data.frame' objects in its input. Also, I
vaguely remember having used
2008 Mar 26
2
pseudo R square and/or C statistic in R logistic regression
Dear all,
I am now doing the logistic regression using R. (glm, family=binomial). Besides the standardize summary statistics generated from R, I am also interested in some more informations concerning the model fitting / prediction etc; Particularly I am interested in "pseudo R squar" and "C statistic". I searched the R- help and could only get very limited information. (Post