## Calculation of AUCt, AUCall, AUCinf... [NCA / SHAM]

Hi Yu!

Thanks, nice to read...

This topic comes up from time to time – not only on this forum, but also at the PKPD-list.

Let’s start with an example – a picture may tell more than a thousand words:

I generated data for a simple 1-compartment model, where k

Since in BE we

F

AUC based on the parameters (calculating the definite integral) is 100/k

Unfortunately we don’t have theoretical curves, but discrete sampling intervals, a lot of noise (biological and – to a much (!) lesser extent – analytical), and an artificial cut-off called LLOQ…

I used this algorithm (t

Obviously the concentration at the 24h sampling point was just above LLOQ for the reference (0.82), and BLQ (0.78) for the test.

We can calculate the AUC by numerical integration up to 24h (t

An enlarged view of the the last section of the profile is given below:

Now let’s start with the calculations (based on the linear trapezoidal rule for simplicity). We will assess the methods based on the bias (deviation from the theoretical value of T/R = 95%).

Calculation up to the last measured concentration for each formulation (AUC

Calculation up to the last time point where both formulations show concentrations above LLOQ (AUC

Like Method #1, but a concentration of zero is

Extrapolation from t

Extrapolation from t

(erroneous; see this post for corrections)

Now let’s look at the Pros and Cons.

Method #1 showed the largest bias, but has its merits. The method is available in all standard software (both commercial and open-source); the bias seen in the example most likely diminishes in ‘real world’ data sets. Pragmatically speaking: this method has already been applied in tenthousands of BE studies…

Although in this simple example Method #2 showed the smallest bias,

Extrapolation based on Method #4 is recommended by the FDA, but considered suboptimal (nicely speaking) by PK experts. Using this method we rely on the concentration with the largest variability (both inaccuracy and precision).

It’s definitely better to go with Method #5,

Conclusion:

It’s like already discussed in this thread on sampling times; whatever method you choose, lay it down in the protocol and calculate all AUCs in the same manner. AUC is a quite robust metric, unless you have a relatively high LLOQ and large differences between formulations deviations simply mean out.

My personal recommendation would be AUC

Oh, there are just a few native speakers here; we do all our best to communicate – somehow.

❝ Hi, its very nice to meet you.

❝ The informations from this forum are very helpful. I very appreciate it.

Thanks, nice to read...

❝ how do you calculate the AUCt and do statistical analysis if there are concentrations below LOQ (at early or late sampling point)? As I know, the't' time in AUCt should be the same time. Someone recommend to make it zero in the early time. And they also recommend to delete it in the late time to calculate the AUCinf. Is it meaningful to do statistical analysis of AUCt with different 't' value? Please let me know how to proceed the study in such a case.

This topic comes up from time to time – not only on this forum, but also at the PKPD-list.

Let’s start with an example – a picture may tell more than a thousand words:

I generated data for a simple 1-compartment model, where k

_{a}= 1, k_{el}= 0.2, and F×V/CL = 100 (reference) / 95 (test).Since in BE we

*assume*V and CL to be constant, and*set*D_{test}= D_{reference}(in most regulations)F

_{rel}= AUC_{test}/ AUC_{reference}.AUC based on the parameters (calculating the definite integral) is 100/k

_{el}-100/k_{a}for the reference and 95/k_{el}-95/k_{a}or AUC_{reference}= 400 and AUC_{test}= 380. Therefore F_{rel}= 380/400 = 95%. Fine.Unfortunately we don’t have theoretical curves, but discrete sampling intervals, a lot of noise (biological and – to a much (!) lesser extent – analytical), and an artificial cut-off called LLOQ…

I used this algorithm (t

_{1}0.5, t_{n}24, 10 sampling points in-between) to establish a sampling schedule. Sampling times were rounded to 2 decimal places, as well as the simulated concentrations; LLOQ was set to 0.8.Obviously the concentration at the 24h sampling point was just above LLOQ for the reference (0.82), and BLQ (0.78) for the test.

`┌──────┬──────┬──────┐`

│ t │ ref │ test |

├──────┼──────┼──────┤

│ 0.00 │ 0.00 │ 0.00 |

│ 0.50 │29.83 │28.34 |

│ 0.71 │37.60 │35.72 |

│ 1.01 │45.29 │43.02 |

│ 1.44 │51.28 │48.72 |

│ 2.05 │53.49 │50.82 |

│ 2.91 │50.43 │47.91 |

│ 4.14 │42.10 │39.99 |

│ 5.89 │30.51 │28.99 |

│ 8.37 │18.73 │17.79 |

│11.90 │ 9.25 │ 8.79 |

│16.92 │ 3.39 │ 3.22 |

│24.00 │ 0.82 │(0.78)|

└──────┴──────┴──────┘

We can calculate the AUC by numerical integration up to 24h (t

_{last}) for the reference, and only up to 16.92h (t_{last}) for the test…An enlarged view of the the last section of the profile is given below:

Now let’s start with the calculations (based on the linear trapezoidal rule for simplicity). We will assess the methods based on the bias (deviation from the theoretical value of T/R = 95%).

**Method 1**:Calculation up to the last measured concentration for each formulation (AUC

_{t})`┌───────────┬───────────┬────────┬────────┐`

│ ref │ test │ T/R │ % Bias │

├───────────┼───────────┼────────┼────────┤

│ 401.88595 │ 367.64200 │ 91.48% │ -3.71% │

└───────────┴───────────┴────────┴────────┘

**Method 2**Calculation up to the last time point where both formulations show concentrations above LLOQ (AUC

_{p})`┌───────────┬───────────┬────────┬────────┐`

│ ref │ test │ T/R │ % Bias │

├───────────┼───────────┼────────┼────────┤

│ 386.98255 │ 367.64200 │ 95.00% │ 0.00% │

└───────────┴───────────┴────────┴────────┘

**Method 3**Like Method #1, but a concentration of zero is

*assumed*at the time point after t_{last},*i.e.*, a small triangle is added to AUC_{t}. The method is available in WinNonlin; this metric was named by Pharsight ‘AUC_{all}’ (you don’t find it in any textbook on PK – it’s Pharsight’s ‘invention’).`┌───────────┬───────────┬────────┬────────┐`

│ ref │ test │ T/R │ % Bias │

├───────────┼───────────┼────────┼────────┤

│ 401.88595 │ 379.04080 │ 94.32% │ -0.72% │

└───────────┴───────────┴────────┴────────┘

**Method 4**Extrapolation from t

_{last}to t=∞ based on the*observed*concentration C_{last}(AUC_{inf}, AUC_{inf-obs})`┌───────────┬───────────┬────────┬────────┐`

│ ref │ test │ T/R │ % Bias │

├───────────┼───────────┼────────┼────────┤

│ 405.98047 │ 383.74829 │ 94.52% │ -0.50% │

└───────────┴───────────┴────────┴────────┘

**Method 5**Extrapolation from t

_{last}to t=∞ based on the*predicted*concentration C_{last}(AUC_{inf}, AUC_{inf-pred})`┌───────────┬───────────┬────────┬────────┐`

│ ref │ test │ T/R │ % Bias │

├───────────┼───────────┼────────┼────────┤

│ 405.93053 │ 383.74829 │ 94.54% │ -0.49% |

└───────────┴───────────┴────────┴────────┘

(erroneous; see this post for corrections)

Now let’s look at the Pros and Cons.

Method #1 showed the largest bias, but has its merits. The method is available in all standard software (both commercial and open-source); the bias seen in the example most likely diminishes in ‘real world’ data sets. Pragmatically speaking: this method has already been applied in tenthousands of BE studies…

Although in this simple example Method #2 showed the smallest bias,

^{1}I would not recommend it because:- You would just throw away valuable information (all concentrations of the ‘longer’ profile after t
_{last}of the ‘shorter’ profile will not be used).

- The method is not standard in any software I know. It would be rather tricky to get the right numbers out from a list of partial areas.

- You would make the life of inspectors miserable.

*before*the dose (the first triangle which is used in all methods), in some situations AUC_{all}> AUC_{inf}– which lead to a lot of questions in the past. This happens if a next sampling interval is quite large, and the estimated half life is comparable small. See also this post.Extrapolation based on Method #4 is recommended by the FDA, but considered suboptimal (nicely speaking) by PK experts. Using this method we rely on the concentration with the largest variability (both inaccuracy and precision).

It’s definitely better to go with Method #5,

^{2}which is standard according to the WHO, European, and many other countries’ guidelines.Conclusion:

It’s like already discussed in this thread on sampling times; whatever method you choose, lay it down in the protocol and calculate all AUCs in the same manner. AUC is a quite robust metric, unless you have a relatively high LLOQ and large differences between formulations deviations simply mean out.

My personal recommendation would be AUC

_{t}calculated by the lin-up/log-down method followed by extrapolation by Method #5.❝ Sorry for the poor english.

Oh, there are just a few native speakers here; we do all our best to communicate – somehow.

- Edit (years later). Why not?

Fisher D, Kramer W, Burmeister Getz E.*Evaluation of a Scenario in Which Estimates of Bioequivalence Are Biased and a Proposed Solution: t*J Clin Pharm. 2016;56(7):794–800. doi:10.1002/jcph.663. free resource._{last}(Common).

- Method #5 is the default in EquivTest/PK (though #4 is available as well).

—

Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮

Science Quotes

*Dif-tor heh smusma*🖖🏼 Довге життя Україна!_{}Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮

Science Quotes

### Complete thread:

- Calculation of AUCt yckim 2007-11-27 02:46 [NCA / SHAM]
- Calculation of AUCt, AUCall, AUCinf...Helmut 2007-11-27 22:02
- Calculation of AUCt, AUCall, AUCinf... yckim 2007-11-28 02:44
- Calculation of AUCt, AUCall, AUCinf... Helmut 2007-11-28 03:49
- Calculation of AUCt, AUCall, AUCinf... jmlee 2008-10-10 08:29
- Calculation of AUCt, AUCall, AUCinf... Helmut 2008-10-10 17:48
- Calculation of AUCt, AUCall, AUCinf... mittyri 2014-03-19 19:15
- Errors happen d_labes 2014-03-20 08:26
- Errors happen fno 2014-03-20 11:03

- Errors happen d_labes 2014-03-20 08:26

- Calculation of AUCt, AUCall, AUCinf... jmlee 2008-10-10 08:29

- Calculation of AUCt, AUCall, AUCinf... Helmut 2007-11-28 03:49
- Calculation of AUCt, AUCall, AUCinf... Ohlbe 2007-11-28 10:37
- Calculation of AUCt, AUCall, AUCinf... Helmut 2007-11-28 12:28
- References Helmut 2007-11-28 17:58

- Calculation of AUCt, AUCall, AUCinf... Helmut 2007-11-28 12:28

- Calculation of AUCt, AUCall, AUCinf... yckim 2007-11-28 02:44

- Calculation of AUCt, AUCall, AUCinf...Helmut 2007-11-27 22:02