strong text#calculate of AUC per each ID
calculate AUC per ID, for CONC>12?
sapply(split(data.f, data.f$ID), function(data.f) trapz(data.f$time, data.f$conc))
data.f=
ID1 time2 conc3
2 5 12
3 7 6
4 12 15