[R] About Mcneil Hanley test for a portion of AUC!

From: Dukka k.c. <dukkakc_at_gmail.com>
Date: Thu, 12 Jun 2008 10:10:25 -0400

Dear all

 I am trying to compare the performances of several methods using the AUC0.1 and
not the whole AUC. (meaning I wanted to compare to AUC's whose x axis only goes to
0.1 not 1)

I came to know about the Mcneil Hanley test from Bernardo Rangel Tura and I referred to the original paper for the calculation of "r" which is an argument of the function
cROC. I can only find the value of "r" for the whole AUC's .

> seROC<-function(AUC,na,nn){
> a<-AUC
> q1<-a/(2-a)
> q2<-(2*a^2)/(1+a)
> se<-sqrt((a*(1-a)+(na-1)*(q1-a^2)+(nn-1)*(q2-a^2))/(nn*na))
> se
> }
> cROC<-function(AUC1,na1,nn1,AUC2,na2,nn2,r){
> se1<-seROC(AUC1,na1,nn1)
> se2<-seROC(AUC2,na2,nn2)
> sed<-sqrt(se1^2+se2^2-2*r*se1*se2)
> zad<-(AUC1-AUC2)/sed
> p<-dnorm(zad)
> a<-list(zad,p)
> a

Could somebody kindly suggest me how to calculate the value of "r" or some ways to calculate the statistical significance measure for the differences of auc for a part of
the curve like AUC0.1.

Thank You

Dukka KC

	[[alternative HTML version deleted]]

R-help_at_r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Thu 12 Jun 2008 - 15:51:12 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 12 Jun 2008 - 17:30:42 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive