From: Dukka k.c. <dukkakc_at_gmail.com>

Date: Thu, 12 Jun 2008 10:10:25 -0400

Date: Thu, 12 Jun 2008 10:10:25 -0400

*> 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
*

Thank You

-- Dukka KC UNCC [[alternative HTML version deleted]] ______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help 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.
*