[R] Comparing differences in AUC from 2 different models

From: Monica Pisica <pisicandru_at_hotmail.com>
Date: Thu, 17 Jul 2008 14:06:10 +0000


I would like to compare differences in AUC from 2 different models, glm and gam for predicting presence / absence. I know that in theory the model with a higher AUC is better, but what I am interested in is if statistically the increase in AUC from the glm model to the gam model is significant. I also read quite extensive discussions on the list about ROC and AUC but I still didn't find my answer.

To calculate the AUC and plot the ROC I used the package PresenceAbsence. The help file for auc() says: " The standard errors from auc are only valid for comparing an individual model to random assignment (i.e. AUC=.5). To compare two models to each other it is necessary to account for correlation due to the fact that they use the same test set. If you are interested in pair wise model comparisons see the Splus ROC library from Mayo clinic. auc is a much simpler function than what is available from the Splus ROC library from Mayo clinic."

I did download this library but I don't have access to S-PLUS and even if supposedly the code is very similar between S-PLUS and R I still don't quite understand what is going on because I am a little bit confused what some parameters represent . For example "markers" and "status", although I think "status" represent my original data (all coded 0 and 1) and "markers" might be the probabilities obtained from my 2 models. The confusion may also steam from the fact that I don't have a medical or biological training and maybe "markers" and "status" do have a special meaning for these 2 disciplines.

I will really appreciate if you can help in finding a way to compare differences in AUC.




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 17 Jul 2008 - 14:12:01 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 17 Jul 2008 - 18:31:43 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