[R] AUC steady state calculations

From: Bart Joosen <bartjoosen_at_hotmail.com>
Date: Sun, 08 Jun 2008 14:22:03 +0200


a collegue has send me an excel sheet with some plasma values, and now he wants to know the AUC steady state. I took a look at the CRAN taskviews, and came up with PK, PKtools, ... The AUC calculation, no problem with that, but how do I calculate the steady state? One way of thinking was with the aid of multiple t-test, to find where we couldn't find difference between the different measuring points, and take this as steady state, but I'm not sure that this is the right way.

As I'm not home in the pharmacokinetic world, I was hoping someone with some more experience can shed some light on this (altough for me) dark material.

Kind Regards


PS: Here some example data:
dat<- structure(list(Sample = c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L), tijd = c(12L, 36L, 60L, 84L, 108L, 132L, 12L, 36L, 60L, 84L, 108L, 132L), conc = c(0.621518061431366, 0.87531564366726,

0.916311538568891, 0.880947260781843, 0.852202744098934, 0.218909173985895, 
1.22305496551836, 1.30075841227452, 0.995918019674464, 1.33214099618361, 
1.42613784296527, 0.290928921761672)), .Names = c("Sample", "Time", 
"conc"), row.names = c(NA, -12L), class = "data.frame")
	[[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 Sun 08 Jun 2008 - 12:26:47 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 Sun 08 Jun 2008 - 12:30:38 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