Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: RE: harrell C statistics comparison


From   Eduardo Nunez <enunezb@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: RE: harrell C statistics comparison
Date   Sun, 6 Sep 2009 13:01:27 -0400

Dear Professor Newson,

Thank you so much for your help.

Best wishes,

Eduardo



On Sun, Sep 6, 2009 at 11:18 AM, Newson, Roger B
<r.newson@imperial.ac.uk> wrote:
>
> In principle, the difference between 2 Harrell's c statistics can be calculated using the -somersd- package (downloadable from SSC), together with -lincom-. As in:
>
> somersd y x1 x2, transf(c) tdist
> lincom x1-x2
>
> which will give you a confidence interval for the difference between Hallell's c of -x1- with respect to -y- and Harrell's c of -x2- with respect to -y-. If -y- is a survival time, with a censorship status variable -cs- which is 1 for right-censored data, -1 for left-censored data and 0 for uncensored data, then you should type
>
> somersd y x1 x2, transf(c) tdist
> lincom x1-x2
>
> This will be a valid comparison if -x1- and -x2- are 2 alternative predictive scores for y, as long as the dataset used to estimate the confidence limits for Harrell's c is not the dataset used to derive the scores in the first place. If you only have 1 dataset, then this should probably be split into a training set (to find the 2 nested models), and a test set (to calculate the 2 Harrell's c-indices and their difference).
>
> I hope this helps.
>
> Best wishes
>
> Roger
>
>
> Roger B Newson BSc MSc DPhil
> Lecturer in Medical Statistics
> Respiratory Epidemiology and Public Health Group
> National Heart and Lung Institute
> Imperial College London
> Royal Brompton Campus
> Room 33, Emmanuel Kaye Building
> 1B Manresa Road
> London SW3 6LR
> UNITED KINGDOM
> Tel: +44 (0)20 7352 8121 ext 3381
> Fax: +44 (0)20 7351 8322
> Email: r.newson@imperial.ac.uk
> Web page: http://www.imperial.ac.uk/nhli/r.newson/
> Departmental Web page:
> http://www1.imperial.ac.uk/medicine/about/divisions/nhli/respiration/popgenetics/reph/
>
> Opinions expressed are those of the author, not of the institution.
>
> -----Original Message-----
> From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Eduardo Nunez
> Sent: 04 September 2009 20:41
> To: statalist@hsphsun2.harvard.edu
> Subject: st: harrell C statistics comparison
>
> Dear Statalister,
>
> A reviewer from a medical journal has asked me to provide 95% CI and
> p-value for two nested Cox models of a Harrell C statistics
> comparison.
> Please help.
>
> For instances,
> Model1= stcox x y, efron
> estat concordance=harrell1
> Model2 stcox x y z, efron
> estat concordance=harrell2
>
> harrell2 - harrell1......xx, 95%CI(xx-yy); p=x.xx
> or maybe harrel2/harrel1 (I am not sure which way)
>
> Maybe boostraping the difference.
>
> I would appreciate any help on this issue.
>
> Eduardo
>
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index