sNCA.RdThis is the work-horse function for NCA.
sNCA(x, y, dose = 0, adm = "Extravascular", dur = 0, doseUnit = "mg", timeUnit = "h", 
     concUnit = "ug/L", iAUC = "", down = "Linear", R2ADJ = 0.7, MW = 0, Keystring="",
     excludeDelta = 1)usually time
usually concentration
given amount, not amount per body weight
one of "Bolus" or "Infusion" or "Extravascular" to indicate drug administration mode
duration of infusion
unit of dose
unit of time
unit of concentration
interval AUCs to calculate
either of "Linear" or "Log" to indicate the way to calculate AUC and AUMC
Minimum adjusted R-square value to determine terminal slope automatically
molecular weight of the drug
a text string to be shown at the plot in case of manual selection of terminal slope
Improvement of R2ADJ larger than this value could exclude the last point. Default value 1 is for the compatibility with other software.
This replaced previous IndiNCA. Author recommends to use excludeDelta option with about 0.3.
maximum concentration, Cmax
dose normalized Cmax, CMAX / Dose, Cmax / Dose
time of maximum concentration, Tmax
time to observe the first non-zero concentration, for extravascular administration only
last positive concentration observed, Clast
last positive concentration predicted, Clast_pred
time of last positive concentration, Tlast
half-life by lambda z, ln(2)/LAMZ
lambda_z negative of the best-fit terminal slope
earliest time for LAMZ
last time for LAMZ
number of points for LAMZ
correlation of log(concentration) and time
R-squared
R-squared adjusted
back extrapolated concentration at time 0, for intravascular bolus administration only
AUC from 0 to TLST
AUC using all the given points, including trailing zero concentrations
AUC infinity observed
AUCIFO / Dose
AUC infinity predicted using CLSTP instead of CLST
AUCIFP / Dose
AUC % extrapolation observed
AUC % extrapolated for AUCIFP
AUC % back extrapolation observed, for bolus IV administration only
AUC % back extrapolation predicted with AUCIFP, for bolus IV administration only
AUMC to the TLST
AUMC infinity observed using CLST
AUMC infinity determined by CLSTP
AUMC % extrapolated observed
AUMC % extrapolated predicted
mean residence time (MRT) to TLST, for intravascular administration
mean residence time (MRT) infinity using CLST, for intravascular administration
mean residence time (MRT) infinity using CLSTP, for intravascular administration
mean residence time (MRT) to TLST, for extravascular administration
mean residence time (MRT) infinity using CLST, for extravascular administration
mean residence time (MRT) infinity using CLSTP, for extravascular administration
volume of distribution determined by LAMZ and AUCIFO, for intravascular administration
volume of distribution determined by LAMZ and AUCIFP, for intravascular administration
VZO for extravascular administration, VZO/F, F is bioavailability
VZP for extravascular administration, VZP/F, F is bioavailability
clearance using AUCIFO, for intravascular administration
clearance using AUCIFP, for intravascular administration
CLO for extravascular administration, CLO/F, F is bioavailability
CLP for extravascular administration, CLP/F, F is bioavailability
volume of distribution at steady state using CLST, for intravascular administration only
volume of distribution at steady state using CLSTP, for intravascular administration only
Gabrielsson J, Weiner D. Pharmacokinetic and Pharmacodynamic Data Analysis - Concepts and Applications. 5th ed. 2016.
# For one subject
x = Theoph[Theoph$Subject=="1","Time"]
y = Theoph[Theoph$Subject=="1","conc"]
sNCA(x, y, dose=320, doseUnit="mg", concUnit="mg/L", timeUnit="h")
#>           b0         CMAX        CMAXD         TMAX         TLAG         CLST 
#>    2.3687851   10.5000000    0.0328125    1.1200000    0.0000000    3.2800000 
#>        CLSTP         TLST       LAMZHL         LAMZ       LAMZLL       LAMZUL 
#>    3.2801465   24.3700000   14.3043776    0.0484570    9.0500000   24.3700000 
#>      LAMZNPT       CORRXY           R2        R2ADJ       AUCLST       AUCALL 
#>    3.0000000   -0.9999999    0.9999997    0.9999995  148.9230500  148.9230500 
#>       AUCIFO      AUCIFOD       AUCIFP      AUCIFPD       AUCPEO       AUCPEP 
#>  216.6119330    0.6769123  216.6149558    0.6769217   31.2489169   31.2498763 
#>      AUMCLST      AUMCIFO      AUMCIFP      AUMCPEO      AUMCPEP         VZFO 
#> 1459.0711035 4505.5348194 4505.6708646   67.6160287   67.6170065   30.4867482 
#>         VZFP         CLFO         CLFP     MRTEVLST     MRTEVIFO     MRTEVIFP 
#>   30.4863228    1.4772963    1.4772757    9.7974834   20.8000305   20.8003683 
#> attr(,"units")
#>  [1] ""          "mg/L"      "mg/L/mg"   "h"         "h"         "mg/L"     
#>  [7] "mg/L"      "h"         "h"         "/h"        "h"         "h"        
#> [13] ""          ""          ""          ""          "h*mg/L"    "h*mg/L"   
#> [19] "h*mg/L"    "h*mg/L/mg" "h*mg/L"    "h*mg/L/mg" "%"         "%"        
#> [25] "h2*mg/L"   "h2*mg/L"   "h2*mg/L"   "%"         "%"         "L"        
#> [31] "L"         "L/h"       "L/h"       "h"         "h"         "h"        
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y, dose=320, concUnit="mg/L")
#>           b0         CMAX        CMAXD         TMAX         TLAG         CLST 
#>    2.3687851   10.5000000    0.0328125    1.1200000    0.0000000    3.2800000 
#>        CLSTP         TLST       LAMZHL         LAMZ       LAMZLL       LAMZUL 
#>    3.2801465   24.3700000   14.3043776    0.0484570    9.0500000   24.3700000 
#>      LAMZNPT       CORRXY           R2        R2ADJ       AUCLST       AUCALL 
#>    3.0000000   -0.9999999    0.9999997    0.9999995  148.9230500  148.9230500 
#>       AUCIFO      AUCIFOD       AUCIFP      AUCIFPD       AUCPEO       AUCPEP 
#>  216.6119330    0.6769123  216.6149558    0.6769217   31.2489169   31.2498763 
#>      AUMCLST      AUMCIFO      AUMCIFP      AUMCPEO      AUMCPEP         VZFO 
#> 1459.0711035 4505.5348194 4505.6708646   67.6160287   67.6170065   30.4867482 
#>         VZFP         CLFO         CLFP     MRTEVLST     MRTEVIFO     MRTEVIFP 
#>   30.4863228    1.4772963    1.4772757    9.7974834   20.8000305   20.8003683 
#> attr(,"units")
#>  [1] ""          "mg/L"      "mg/L/mg"   "h"         "h"         "mg/L"     
#>  [7] "mg/L"      "h"         "h"         "/h"        "h"         "h"        
#> [13] ""          ""          ""          ""          "h*mg/L"    "h*mg/L"   
#> [19] "h*mg/L"    "h*mg/L/mg" "h*mg/L"    "h*mg/L/mg" "%"         "%"        
#> [25] "h2*mg/L"   "h2*mg/L"   "h2*mg/L"   "%"         "%"         "L"        
#> [31] "L"         "L/h"       "L/h"       "h"         "h"         "h"        
#> attr(,"UsedPoints")
#> [1]  9 10 11
iAUC = data.frame(Name=c("AUC[0-12h]","AUC[0-24h]"), Start=c(0,0), End=c(12,24))
sNCA(x, y, dose=320, doseUnit="mg", concUnit="mg/L", timeUnit="h", iAUC=iAUC)
#>           b0         CMAX        CMAXD         TMAX         TLAG         CLST 
#>    2.3687851   10.5000000    0.0328125    1.1200000    0.0000000    3.2800000 
#>        CLSTP         TLST       LAMZHL         LAMZ       LAMZLL       LAMZUL 
#>    3.2801465   24.3700000   14.3043776    0.0484570    9.0500000   24.3700000 
#>      LAMZNPT       CORRXY           R2        R2ADJ       AUCLST       AUCALL 
#>    3.0000000   -0.9999999    0.9999997    0.9999995  148.9230500  148.9230500 
#>       AUCIFO      AUCIFOD       AUCIFP      AUCIFPD       AUCPEO       AUCPEP 
#>  216.6119330    0.6769123  216.6149558    0.6769217   31.2489169   31.2498763 
#>      AUMCLST      AUMCIFO      AUMCIFP      AUMCPEO      AUMCPEP         VZFO 
#> 1459.0711035 4505.5348194 4505.6708646   67.6160287   67.6170065   30.4867482 
#>         VZFP         CLFO         CLFP     MRTEVLST     MRTEVIFO     MRTEVIFP 
#>   30.4863228    1.4772963    1.4772757    9.7974834   20.8000305   20.8003683 
#>   AUC[0-12h]   AUC[0-24h] 
#>   91.7355220  147.6945866 
#> attr(,"units")
#>  [1] ""          "mg/L"      "mg/L/mg"   "h"         "h"         "mg/L"     
#>  [7] "mg/L"      "h"         "h"         "/h"        "h"         "h"        
#> [13] ""          ""          ""          ""          "h*mg/L"    "h*mg/L"   
#> [19] "h*mg/L"    "h*mg/L/mg" "h*mg/L"    "h*mg/L/mg" "%"         "%"        
#> [25] "h2*mg/L"   "h2*mg/L"   "h2*mg/L"   "%"         "%"         "L"        
#> [31] "L"         "L/h"       "L/h"       "h"         "h"         "h"        
#> [37] "h*mg/L"    "h*mg/L"   
#> attr(,"UsedPoints")
#> [1]  9 10 11
MW = 180.164 # Molecular weight of theophylline
sNCA(x, y/MW, dose=320, doseUnit="mg", concUnit="mmol/L", timeUnit="h")
#> Warning: Molecular weight should be given for more informative results!
#>            b0          CMAX         CMAXD          TMAX          TLAG 
#> -2.8250824530  0.0582802336  0.0001821257  1.1200000000  0.0000000000 
#>          CLST         CLSTP          TLST        LAMZHL          LAMZ 
#>  0.0182056349  0.0182064479 24.3700000000 14.3043775711  0.0484569970 
#>        LAMZLL        LAMZUL       LAMZNPT        CORRXY            R2 
#>  9.0500000000 24.3700000000  3.0000000000 -0.9999998648  0.9999997297 
#>         R2ADJ        AUCLST        AUCALL        AUCIFO       AUCIFOD 
#>  0.9999994593  0.8265971559  0.8265971559  1.2023041953  0.0037572006 
#>        AUCIFP       AUCIFPD        AUCPEO        AUCPEP       AUMCLST 
#>  1.2023209731  0.0037572530 31.2489169405 31.2498763313  8.0985718762 
#>       AUMCIFO       AUMCIFP       AUMCPEO       AUMCPEP          VZFO 
#> 25.0079639629 25.0087190814 67.6160286851 67.6170064935            NA 
#>          VZFP          CLFO          CLFP      MRTEVLST      MRTEVIFO 
#>            NA            NA            NA  9.7974833547 20.8000305256 
#>      MRTEVIFP 
#> 20.8003683211 
#> attr(,"units")
#>  [1] ""            "mmol/L"      "mmol/L/mg"   "h"           "h"          
#>  [6] "mmol/L"      "mmol/L"      "h"           "h"           "/h"         
#> [11] "h"           "h"           ""            ""            ""           
#> [16] ""            "h*mmol/L"    "h*mmol/L"    "h*mmol/L"    "h*mmol/L/mg"
#> [21] "h*mmol/L"    "h*mmol/L/mg" "%"           "%"           "h2*mmol/L"  
#> [26] "h2*mmol/L"   "h2*mmol/L"   "%"           "%"           "L"          
#> [31] "L"           "L/h"         "L/h"         "h"           "h"          
#> [36] "h"          
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y/MW, dose=320, doseUnit="mg", concUnit="mmol/L", timeUnit="h", MW=MW)
#>            b0          CMAX         CMAXD          TMAX          TLAG 
#> -2.8250824530  0.0582802336  0.0001821257  1.1200000000  0.0000000000 
#>          CLST         CLSTP          TLST        LAMZHL          LAMZ 
#>  0.0182056349  0.0182064479 24.3700000000 14.3043775711  0.0484569970 
#>        LAMZLL        LAMZUL       LAMZNPT        CORRXY            R2 
#>  9.0500000000 24.3700000000  3.0000000000 -0.9999998648  0.9999997297 
#>         R2ADJ        AUCLST        AUCALL        AUCIFO       AUCIFOD 
#>  0.9999994593  0.8265971559  0.8265971559  1.2023041953  0.0037572006 
#>        AUCIFP       AUCIFPD        AUCPEO        AUCPEP       AUMCLST 
#>  1.2023209731  0.0037572530 31.2489169405 31.2498763313  8.0985718762 
#>       AUMCIFO       AUMCIFP       AUMCPEO       AUMCPEP          VZFO 
#> 25.0079639629 25.0087190814 67.6160286851 67.6170064935 30.4867482346 
#>          VZFP          CLFO          CLFP      MRTEVLST      MRTEVIFO 
#> 30.4863228055  1.4772962667  1.4772756517  9.7974833547 20.8000305256 
#>      MRTEVIFP 
#> 20.8003683211 
#> attr(,"units")
#>  [1] ""            "mmol/L"      "mmol/L/mg"   "h"           "h"          
#>  [6] "mmol/L"      "mmol/L"      "h"           "h"           "/h"         
#> [11] "h"           "h"           ""            ""            ""           
#> [16] ""            "h*mmol/L"    "h*mmol/L"    "h*mmol/L"    "h*mmol/L/mg"
#> [21] "h*mmol/L"    "h*mmol/L/mg" "%"           "%"           "h2*mmol/L"  
#> [26] "h2*mmol/L"   "h2*mmol/L"   "%"           "%"           "L"          
#> [31] "L"           "L/h"         "L/h"         "h"           "h"          
#> [36] "h"          
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y, dose=320/MW, doseUnit="mmol", concUnit="mg/L", timeUnit="h", MW=MW)
#>           b0         CMAX        CMAXD         TMAX         TLAG         CLST 
#>    2.3687851   10.5000000    5.9116312    1.1200000    0.0000000    3.2800000 
#>        CLSTP         TLST       LAMZHL         LAMZ       LAMZLL       LAMZUL 
#>    3.2801465   24.3700000   14.3043776    0.0484570    9.0500000   24.3700000 
#>      LAMZNPT       CORRXY           R2        R2ADJ       AUCLST       AUCALL 
#>    3.0000000   -0.9999999    0.9999997    0.9999995  148.9230500  148.9230500 
#>       AUCIFO      AUCIFOD       AUCIFP      AUCIFPD       AUCPEO       AUCPEP 
#>  216.6119330  121.9552259  216.6149558  121.9569278   31.2489169   31.2498763 
#>      AUMCLST      AUMCIFO      AUMCIFP      AUMCPEO      AUMCPEP         VZFO 
#> 1459.0711035 4505.5348194 4505.6708646   67.6160287   67.6170065   30.4867482 
#>         VZFP         CLFO         CLFP     MRTEVLST     MRTEVIFO     MRTEVIFP 
#>   30.4863228    1.4772963    1.4772757    9.7974834   20.8000305   20.8003683 
#> attr(,"units")
#>  [1] ""            "mg/L"        "mg/L/mmol"   "h"           "h"          
#>  [6] "mg/L"        "mg/L"        "h"           "h"           "/h"         
#> [11] "h"           "h"           ""            ""            ""           
#> [16] ""            "h*mg/L"      "h*mg/L"      "h*mg/L"      "h*mg/L/mmol"
#> [21] "h*mg/L"      "h*mg/L/mmol" "%"           "%"           "h2*mg/L"    
#> [26] "h2*mg/L"     "h2*mg/L"     "%"           "%"           "L"          
#> [31] "L"           "L/h"         "L/h"         "h"           "h"          
#> [36] "h"          
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y/MW, dose=320/MW, doseUnit="mmol", concUnit="mmol/L", timeUnit="h", MW=MW)
#>          b0        CMAX       CMAXD        TMAX        TLAG        CLST 
#> -2.82508245  0.05828023  0.03281250  1.12000000  0.00000000  0.01820563 
#>       CLSTP        TLST      LAMZHL        LAMZ      LAMZLL      LAMZUL 
#>  0.01820645 24.37000000 14.30437757  0.04845700  9.05000000 24.37000000 
#>     LAMZNPT      CORRXY          R2       R2ADJ      AUCLST      AUCALL 
#>  3.00000000 -0.99999986  0.99999973  0.99999946  0.82659716  0.82659716 
#>      AUCIFO     AUCIFOD      AUCIFP     AUCIFPD      AUCPEO      AUCPEP 
#>  1.20230420  0.67691229  1.20232097  0.67692174 31.24891694 31.24987633 
#>     AUMCLST     AUMCIFO     AUMCIFP     AUMCPEO     AUMCPEP        VZFO 
#>  8.09857188 25.00796396 25.00871908 67.61602869 67.61700649 30.48674823 
#>        VZFP        CLFO        CLFP    MRTEVLST    MRTEVIFO    MRTEVIFP 
#> 30.48632281  1.47729627  1.47727565  9.79748335 20.80003053 20.80036832 
#> attr(,"units")
#>  [1] ""              "mmol/L"        "mmol/L/mmol"   "h"            
#>  [5] "h"             "mmol/L"        "mmol/L"        "h"            
#>  [9] "h"             "/h"            "h"             "h"            
#> [13] ""              ""              ""              ""             
#> [17] "h*mmol/L"      "h*mmol/L"      "h*mmol/L"      "h*mmol/L/mmol"
#> [21] "h*mmol/L"      "h*mmol/L/mmol" "%"             "%"            
#> [25] "h2*mmol/L"     "h2*mmol/L"     "h2*mmol/L"     "%"            
#> [29] "%"             "L"             "L"             "L/h"          
#> [33] "L/h"           "h"             "h"             "h"            
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y/MW, dose=320/MW, doseUnit="mmol", concUnit="mmol/L", timeUnit="h", MW=MW)
#>          b0        CMAX       CMAXD        TMAX        TLAG        CLST 
#> -2.82508245  0.05828023  0.03281250  1.12000000  0.00000000  0.01820563 
#>       CLSTP        TLST      LAMZHL        LAMZ      LAMZLL      LAMZUL 
#>  0.01820645 24.37000000 14.30437757  0.04845700  9.05000000 24.37000000 
#>     LAMZNPT      CORRXY          R2       R2ADJ      AUCLST      AUCALL 
#>  3.00000000 -0.99999986  0.99999973  0.99999946  0.82659716  0.82659716 
#>      AUCIFO     AUCIFOD      AUCIFP     AUCIFPD      AUCPEO      AUCPEP 
#>  1.20230420  0.67691229  1.20232097  0.67692174 31.24891694 31.24987633 
#>     AUMCLST     AUMCIFO     AUMCIFP     AUMCPEO     AUMCPEP        VZFO 
#>  8.09857188 25.00796396 25.00871908 67.61602869 67.61700649 30.48674823 
#>        VZFP        CLFO        CLFP    MRTEVLST    MRTEVIFO    MRTEVIFP 
#> 30.48632281  1.47729627  1.47727565  9.79748335 20.80003053 20.80036832 
#> attr(,"units")
#>  [1] ""              "mmol/L"        "mmol/L/mmol"   "h"            
#>  [5] "h"             "mmol/L"        "mmol/L"        "h"            
#>  [9] "h"             "/h"            "h"             "h"            
#> [13] ""              ""              ""              ""             
#> [17] "h*mmol/L"      "h*mmol/L"      "h*mmol/L"      "h*mmol/L/mmol"
#> [21] "h*mmol/L"      "h*mmol/L/mmol" "%"             "%"            
#> [25] "h2*mmol/L"     "h2*mmol/L"     "h2*mmol/L"     "%"            
#> [29] "%"             "L"             "L"             "L/h"          
#> [33] "L/h"           "h"             "h"             "h"            
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y/MW, doseUnit="mmol", concUnit="mmol/L", timeUnit="h", MW=MW)
#>          b0        CMAX       CMAXD        TMAX        TLAG        CLST 
#> -2.82508245  0.05828023          NA  1.12000000  0.00000000  0.01820563 
#>       CLSTP        TLST      LAMZHL        LAMZ      LAMZLL      LAMZUL 
#>  0.01820645 24.37000000 14.30437757  0.04845700  9.05000000 24.37000000 
#>     LAMZNPT      CORRXY          R2       R2ADJ      AUCLST      AUCALL 
#>  3.00000000 -0.99999986  0.99999973  0.99999946  0.82659716  0.82659716 
#>      AUCIFO     AUCIFOD      AUCIFP     AUCIFPD      AUCPEO      AUCPEP 
#>  1.20230420          NA  1.20232097          NA 31.24891694 31.24987633 
#>     AUMCLST     AUMCIFO     AUMCIFP     AUMCPEO     AUMCPEP        VZFO 
#>  8.09857188 25.00796396 25.00871908 67.61602869 67.61700649  0.00000000 
#>        VZFP        CLFO        CLFP    MRTEVLST    MRTEVIFO    MRTEVIFP 
#>  0.00000000  0.00000000  0.00000000  9.79748335 20.80003053 20.80036832 
#> attr(,"units")
#>  [1] ""              "mmol/L"        "mmol/L/mmol"   "h"            
#>  [5] "h"             "mmol/L"        "mmol/L"        "h"            
#>  [9] "h"             "/h"            "h"             "h"            
#> [13] ""              ""              ""              ""             
#> [17] "h*mmol/L"      "h*mmol/L"      "h*mmol/L"      "h*mmol/L/mmol"
#> [21] "h*mmol/L"      "h*mmol/L/mmol" "%"             "%"            
#> [25] "h2*mmol/L"     "h2*mmol/L"     "h2*mmol/L"     "%"            
#> [29] "%"             "L"             "L"             "L/h"          
#> [33] "L/h"           "h"             "h"             "h"            
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y/MW, dose=as.numeric(NA), doseUnit="mmol", concUnit="mmol/L", timeUnit="h", 
     MW=MW)
#>          b0        CMAX       CMAXD        TMAX        TLAG        CLST 
#> -2.82508245  0.05828023          NA  1.12000000  0.00000000  0.01820563 
#>       CLSTP        TLST      LAMZHL        LAMZ      LAMZLL      LAMZUL 
#>  0.01820645 24.37000000 14.30437757  0.04845700  9.05000000 24.37000000 
#>     LAMZNPT      CORRXY          R2       R2ADJ      AUCLST      AUCALL 
#>  3.00000000 -0.99999986  0.99999973  0.99999946  0.82659716  0.82659716 
#>      AUCIFO     AUCIFOD      AUCIFP     AUCIFPD      AUCPEO      AUCPEP 
#>  1.20230420          NA  1.20232097          NA 31.24891694 31.24987633 
#>     AUMCLST     AUMCIFO     AUMCIFP     AUMCPEO     AUMCPEP        VZFO 
#>  8.09857188 25.00796396 25.00871908 67.61602869 67.61700649          NA 
#>        VZFP        CLFO        CLFP    MRTEVLST    MRTEVIFO    MRTEVIFP 
#>          NA          NA          NA  9.79748335 20.80003053 20.80036832 
#> attr(,"units")
#>  [1] ""              "mmol/L"        "mmol/L/mmol"   "h"            
#>  [5] "h"             "mmol/L"        "mmol/L"        "h"            
#>  [9] "h"             "/h"            "h"             "h"            
#> [13] ""              ""              ""              ""             
#> [17] "h*mmol/L"      "h*mmol/L"      "h*mmol/L"      "h*mmol/L/mmol"
#> [21] "h*mmol/L"      "h*mmol/L/mmol" "%"             "%"            
#> [25] "h2*mmol/L"     "h2*mmol/L"     "h2*mmol/L"     "%"            
#> [29] "%"             "L"             "L"             "L/h"          
#> [33] "L/h"           "h"             "h"             "h"            
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x, y, dose=320, concUnit="mg/L", timeUnit="hr")
#>           b0         CMAX        CMAXD         TMAX         TLAG         CLST 
#>    2.3687851   10.5000000    0.0328125    1.1200000    0.0000000    3.2800000 
#>        CLSTP         TLST       LAMZHL         LAMZ       LAMZLL       LAMZUL 
#>    3.2801465   24.3700000   14.3043776    0.0484570    9.0500000   24.3700000 
#>      LAMZNPT       CORRXY           R2        R2ADJ       AUCLST       AUCALL 
#>    3.0000000   -0.9999999    0.9999997    0.9999995  148.9230500  148.9230500 
#>       AUCIFO      AUCIFOD       AUCIFP      AUCIFPD       AUCPEO       AUCPEP 
#>  216.6119330    0.6769123  216.6149558    0.6769217   31.2489169   31.2498763 
#>      AUMCLST      AUMCIFO      AUMCIFP      AUMCPEO      AUMCPEP         VZFO 
#> 1459.0711035 4505.5348194 4505.6708646   67.6160287   67.6170065   30.4867482 
#>         VZFP         CLFO         CLFP     MRTEVLST     MRTEVIFO     MRTEVIFP 
#>   30.4863228    1.4772963    1.4772757    9.7974834   20.8000305   20.8003683 
#> attr(,"units")
#>  [1] ""           "mg/L"       "mg/L/mg"    "hr"         "hr"        
#>  [6] "mg/L"       "mg/L"       "hr"         "hr"         "/hr"       
#> [11] "hr"         "hr"         ""           ""           ""          
#> [16] ""           "hr*mg/L"    "hr*mg/L"    "hr*mg/L"    "hr*mg/L/mg"
#> [21] "hr*mg/L"    "hr*mg/L/mg" "%"          "%"          "hr2*mg/L"  
#> [26] "hr2*mg/L"   "hr2*mg/L"   "%"          "%"          "L"         
#> [31] "L"          "L/hr"       "L/hr"       "hr"         "hr"        
#> [36] "hr"        
#> attr(,"UsedPoints")
#> [1]  9 10 11
sNCA(x*60, y, dose=320, concUnit="mg/L", timeUnit="min")
#>            b0          CMAX         CMAXD          TMAX          TLAG 
#>  2.368785e+00  1.050000e+01  3.281250e-02  6.720000e+01  0.000000e+00 
#>          CLST         CLSTP          TLST        LAMZHL          LAMZ 
#>  3.280000e+00  3.280146e+00  1.462200e+03  8.582627e+02  8.076166e-04 
#>        LAMZLL        LAMZUL       LAMZNPT        CORRXY            R2 
#>  5.430000e+02  1.462200e+03  3.000000e+00 -9.999999e-01  9.999997e-01 
#>         R2ADJ        AUCLST        AUCALL        AUCIFO       AUCIFOD 
#>  9.999995e-01  8.935383e+03  8.935383e+03  1.299672e+04  4.061474e+01 
#>        AUCIFP       AUCIFPD        AUCPEO        AUCPEP       AUMCLST 
#>  1.299690e+04  4.061530e+01  3.124892e+01  3.124988e+01  5.252656e+06 
#>       AUMCIFO       AUMCIFP       AUMCPEO       AUMCPEP          VZFO 
#>  1.621993e+07  1.622042e+07  6.761603e+01  6.761701e+01  3.048675e+01 
#>          VZFP          CLFO          CLFP      MRTEVLST      MRTEVIFO 
#>  3.048632e+01  2.462160e-02  2.462126e-02  5.878490e+02  1.248002e+03 
#>      MRTEVIFP 
#>  1.248022e+03 
#> attr(,"units")
#>  [1] ""            "mg/L"        "mg/L/mg"     "min"         "min"        
#>  [6] "mg/L"        "mg/L"        "min"         "min"         "/min"       
#> [11] "min"         "min"         ""            ""            ""           
#> [16] ""            "min*mg/L"    "min*mg/L"    "min*mg/L"    "min*mg/L/mg"
#> [21] "min*mg/L"    "min*mg/L/mg" "%"           "%"           "min2*mg/L"  
#> [26] "min2*mg/L"   "min2*mg/L"   "%"           "%"           "L"          
#> [31] "L"           "L/min"       "L/min"       "min"         "min"        
#> [36] "min"        
#> attr(,"UsedPoints")
#> [1]  9 10 11