[R] C-statistic comparison with partially paired datasets

Hanneke Wijnhoven hanneke.wijnhoven at falw.vu.nl
Wed Aug 12 14:09:19 CEST 2009


Thank you for your quick response!

I want to compare the discriminative capacity of different 
anthropometric measures in predicting mortality, focussing on the "thin" 
site of these measures.
Since these associations are not linear (U shaped for BMI and inversily 
J-shaped for mid-upper arm circumference) and I do not want to include 
the prediction by "obesity", I am using all values below the median of 
each separate measure to calculate a C-statistic (below the median, the 
association is approximately linear).
As a result, some different and some overlapping cases are included.
I understand your point though.

Any suggestion is welcome.


Frank E Harrell Jr schreef:
> Hanneke Wijnhoven wrote:
>> 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 for any help,
>> Hanneke Wijnhoven
> Hanneke,
> I'm having trouble seeing how the unpaired observations can contribute 
> information in general.  If for example all of the observations were 
> unpaired, one C-statistic might be larger because it came from a 
> dataset with more extreme observations that were easier to discriminate.
> Frank

Hanneke A.H. Wijnhoven (PhD) 
Institute of Health Sciences 
Vrije Universiteit Amsterdam 
De Boelelaan 1085 
1081 HV Amsterdam 
The Netherlands 
Tel. +31 (0) 20 5989951 
Fax. +31 (0) 20 5986940 
hanneke.wijnhoven at falw.vu.nl

More information about the R-help mailing list