Function for AUC?

RSiteSearch("AUC")

would lead you to

http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46416.html

Hello

Is there an easy way, i.e. a function in a package, to calculate the
area under the curve (AUC) for drug serum levels?

Thanks for any advice
--
Armin Goralczyk, M.D.
--
Universit?tsmedizin G?ttingen
Abteilung Allgemein- und Viszeralchirurgie
Rudolf-Koch-Str. 40
39099 G?ttingen
--
Dept. of General Surgery
University of G?ttingen
G?ttingen, Germany
--
http://www.chirurgie-goettingen.de

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?...
Nom : non disponible
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20071217/1e517255/attachment.pl

Hi all

Hi Armin,
Do you know the rocr package ? This is very easy to draw ROC curves and to
calculate AUC with it.
http://rocr.bioinf.mpi-sb.mpg.de/
Hope this will help.

I know ROCR although I am not familiar with ROC. I have seen the AUC
function of ROCR, but are you sure that this is applicable to the
stated problem, i.e. calculate AUC of serum levels of a drug over
time?! ROC seems to handle completely different problems.

I have been reading the following post (thanks for the hint Stephen):

On Dec 17, 2007 2:58 AM, Stephen Weigand <weigand.stephen at gmail.com> wrote:

RSiteSearch("AUC")

would lead you to

http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46416.html

I tried it:

y<-c(1,2,3,4,5);x<-c(10,15,10,5,0)
trap.rule <- function(x,y) sum(diff(x)*(y[-1]+y[-length(y)]))/2
trap.rule(x,y)

[1] -45

It is not the correct value, but the formula seems applicable and I
changed it to

auc <- function(x,y) sum((x[-length(x)] + x[-1]) * (y[-1]-y[-length(y)]))/2
auc(x,y)

[1] 35

which seems to be correct. I hope everyone agees.
I didn't know it's that simple. I guess I don't need another function.
Thank's for all the help and all the suggestions.

Armin Goralczyk, M.D.
--
Universit?tsmedizin G?ttingen
Abteilung Allgemein- und Viszeralchirurgie
Rudolf-Koch-Str. 40
39099 G?ttingen
--
Dept. of General Surgery
University of G?ttingen
G?ttingen, Germany
--
http://www.chirurgie-goettingen.de

Armin Goralczyk <agoralczyk at gmail.com> [Mon, Dec 17, 2007 at 11:07:25AM CET]:
[AUC]

I tried it:

y<-c(1,2,3,4,5);x<-c(10,15,10,5,0)

Are you sure you don't have x and y wrong? Normally the x values
should be monotonically increasing.

Johannes H?sing               There is something fascinating about science. 
                              One gets such wholesale returns of conjecture 
mailto:johannes at huesing.name  from such a trifling investment of fact.                
http://derwisch.wikidot.com         (Mark Twain, "Life on the Mississippi")