R Under development (unstable) (2024-09-03 r87093 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > # R CMD BATCH --no-timing --no-restore --no-save fetal_pbtk_testing.R fetal_pbtk_testing.Rout > > # Get rid of anything in the workspace: > rm(list=ls()) > > library(httk) > > head(solve_fetal_pbtk(chem.name = 'bisphenol a', daily.dose = 1, + doses.per.day = 3)) None of the monitored components undergo unit conversions (i.e. conversion factor of 1). AUC is area under the maternal plasma concentration curve in uM*days units. fAUC is area under the fetal plasma concentration curve in uM*days units. The model outputs are provided in the following units: umol: Agutlumen, Atubules, Ametabolized, Aplacenta uM: Cgut, Cliver, Cven, Clung, Cart, Cadipose, Crest, Ckidney, Cplasma, Cplacenta, Cfliver, Cfven, Cfart, Cfgut, Cflung, Cfrest, Cfthyroid, Cfkidney, Cfbrain, Cfplasma unitless: Rblood2plasma, Rfblood2plasma uM*days: AUC, fAUC L/d: Qcardiac, Qthyroid time Agutlumen Cgut Cliver Cven Clung Cart Cadipose Crest Ckidney [1,] 91 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 92 3.111 6.237 8.951 1.647 3.457 1.647 2.707 2.370 9.461 [3,] 93 3.111 12.300 17.730 3.268 6.862 3.268 5.400 4.713 18.780 [4,] 94 3.111 18.190 26.240 4.841 10.170 4.841 8.014 6.987 27.820 [5,] 95 3.111 23.910 34.510 6.366 13.380 6.366 10.550 9.194 36.600 [6,] 96 3.111 29.460 42.540 7.846 16.490 7.846 13.010 11.340 45.130 Cplasma Atubules Ametabolized Rblood2plasma AUC fAUC Aplacenta [1,] 0.0000 0.000 0.000 3.390 0.0000 0.0000 0.0000 [2,] 0.4861 1.496 1.300 3.388 0.2864 0.2226 0.1909 [3,] 0.9649 5.532 4.712 3.387 1.0580 0.8264 0.3872 [4,] 1.4300 12.050 10.170 3.386 2.3010 1.8010 0.5859 [5,] 1.8810 20.980 17.630 3.385 4.0020 3.1350 0.7869 [6,] 2.3190 32.260 27.020 3.384 6.1480 4.8190 0.9902 Cplacenta Cfliver Cfven Cfart Cfgut Cflung Cfrest Cfthyroid Cfkidney [1,] 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 1.750 8.153 1.451 1.451 5.839 3.356 2.264 3.192 8.775 [3,] 3.473 16.190 2.890 2.890 11.590 6.664 4.496 6.339 17.430 [4,] 5.146 23.990 4.296 4.296 17.180 9.876 6.663 9.393 25.820 [5,] 6.770 31.570 5.669 5.669 22.610 12.990 8.767 12.360 33.970 [6,] 8.346 38.920 7.010 7.010 27.870 16.020 10.810 15.240 41.890 Cfbrain Cfplasma Rfblood2plasma Qcardiac Qthyroid [1,] 0.000 0.0000 3.788 8656 118.6 [2,] 2.548 0.3818 3.800 8673 118.7 [3,] 5.060 0.7580 3.812 8690 118.8 [4,] 7.498 1.1230 3.824 8707 118.9 [5,] 9.865 1.4780 3.836 8724 119.0 [6,] 12.160 1.8220 3.848 8741 119.1 > head(solve_fetal_pbtk(chem.cas="80-05-7")) None of the monitored components undergo unit conversions (i.e. conversion factor of 1). AUC is area under the maternal plasma concentration curve in uM*days units. fAUC is area under the fetal plasma concentration curve in uM*days units. The model outputs are provided in the following units: umol: Agutlumen, Atubules, Ametabolized, Aplacenta uM: Cgut, Cliver, Cven, Clung, Cart, Cadipose, Crest, Ckidney, Cplasma, Cplacenta, Cfliver, Cfven, Cfart, Cfgut, Cflung, Cfrest, Cfthyroid, Cfkidney, Cfbrain, Cfplasma unitless: Rblood2plasma, Rfblood2plasma uM*days: AUC, fAUC L/d: Qcardiac, Qthyroid time Agutlumen Cgut Cliver Cven Clung Cart Cadipose Crest Ckidney [1,] 91 1.724e+02 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 92 2.336e-02 6.192 8.954 1.655 3.475 1.655 2.748 2.390 9.509 [3,] 93 3.170e-06 6.016 8.701 1.608 3.376 1.608 2.670 2.323 9.240 [4,] 94 0.000e+00 5.846 8.454 1.562 3.281 1.562 2.595 2.257 8.978 [5,] 95 0.000e+00 5.680 8.215 1.517 3.188 1.517 2.521 2.193 8.724 [6,] 96 0.000e+00 5.519 7.982 1.474 3.097 1.474 2.450 2.131 8.476 Cplasma Atubules Ametabolized Rblood2plasma AUC fAUC Aplacenta [1,] 0.0000 0.000 0.000 3.390 0.0000 0.0000 0.0000 [2,] 0.4885 2.349 2.012 3.388 0.4499 0.3511 0.1919 [3,] 0.4747 4.869 4.107 3.387 0.9315 0.7296 0.1905 [4,] 0.4613 7.321 6.143 3.386 1.4000 1.0970 0.1891 [5,] 0.4482 9.709 8.121 3.385 1.8540 1.4550 0.1875 [6,] 0.4355 12.030 10.040 3.384 2.2960 1.8020 0.1860 Cplacenta Cfliver Cfven Cfart Cfgut Cflung Cfrest Cfthyroid Cfkidney [1,] 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 1.759 8.202 1.459 1.459 5.874 3.376 2.278 3.211 8.828 [3,] 1.709 7.971 1.423 1.423 5.708 3.281 2.214 3.120 8.579 [4,] 1.661 7.745 1.387 1.387 5.546 3.188 2.151 3.032 8.335 [5,] 1.613 7.525 1.351 1.351 5.389 3.097 2.090 2.946 8.099 [6,] 1.568 7.312 1.317 1.317 5.236 3.009 2.031 2.862 7.869 Cfbrain Cfplasma Rfblood2plasma Qcardiac Qthyroid [1,] 0.000 0.0000 3.788 8656 118.6 [2,] 2.564 0.3840 3.800 8673 118.7 [3,] 2.491 0.3731 3.812 8690 118.8 [4,] 2.420 0.3626 3.824 8707 118.9 [5,] 2.352 0.3523 3.836 8724 119.0 [6,] 2.285 0.3423 3.848 8741 119.1 > head(solve_fetal_pbtk(parameters=parameterize_fetal_pbtk(chem.cas="80-05-7"))) None of the monitored components undergo unit conversions (i.e. conversion factor of 1). AUC is area under the maternal plasma concentration curve in uM*days units. fAUC is area under the fetal plasma concentration curve in uM*days units. The model outputs are provided in the following units: umol: Agutlumen, Atubules, Ametabolized, Aplacenta uM: Cgut, Cliver, Cven, Clung, Cart, Cadipose, Crest, Ckidney, Cplasma, Cplacenta, Cfliver, Cfven, Cfart, Cfgut, Cflung, Cfrest, Cfthyroid, Cfkidney, Cfbrain, Cfplasma unitless: Rblood2plasma, Rfblood2plasma uM*days: AUC, fAUC L/d: Qcardiac, Qthyroid time Agutlumen Cgut Cliver Cven Clung Cart Cadipose Crest Ckidney [1,] 91 1.724e+02 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 92 2.336e-02 6.192 8.954 1.655 3.475 1.655 2.748 2.390 9.509 [3,] 93 3.170e-06 6.016 8.701 1.608 3.376 1.608 2.670 2.323 9.240 [4,] 94 0.000e+00 5.846 8.454 1.562 3.281 1.562 2.595 2.257 8.978 [5,] 95 0.000e+00 5.680 8.215 1.517 3.188 1.517 2.521 2.193 8.724 [6,] 96 0.000e+00 5.519 7.982 1.474 3.097 1.474 2.450 2.131 8.476 Cplasma Atubules Ametabolized Rblood2plasma AUC fAUC Aplacenta [1,] 0.0000 0.000 0.000 3.390 0.0000 0.0000 0.0000 [2,] 0.4885 2.349 2.012 3.388 0.4499 0.3511 0.1919 [3,] 0.4747 4.869 4.107 3.387 0.9315 0.7296 0.1905 [4,] 0.4613 7.321 6.143 3.386 1.4000 1.0970 0.1891 [5,] 0.4482 9.709 8.121 3.385 1.8540 1.4550 0.1875 [6,] 0.4355 12.030 10.040 3.384 2.2960 1.8020 0.1860 Cplacenta Cfliver Cfven Cfart Cfgut Cflung Cfrest Cfthyroid Cfkidney [1,] 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 [2,] 1.759 8.202 1.459 1.459 5.874 3.376 2.278 3.211 8.828 [3,] 1.709 7.971 1.423 1.423 5.708 3.281 2.214 3.120 8.579 [4,] 1.661 7.745 1.387 1.387 5.546 3.188 2.151 3.032 8.335 [5,] 1.613 7.525 1.351 1.351 5.389 3.097 2.090 2.946 8.099 [6,] 1.568 7.312 1.317 1.317 5.236 3.009 2.031 2.862 7.869 Cfbrain Cfplasma Rfblood2plasma Qcardiac Qthyroid [1,] 0.000 0.0000 3.788 8656 118.6 [2,] 2.564 0.3840 3.800 8673 118.7 [3,] 2.491 0.3731 3.812 8690 118.8 [4,] 2.420 0.3626 3.824 8707 118.9 [5,] 2.352 0.3523 3.836 8724 119.0 [6,] 2.285 0.3423 3.848 8741 119.1 Warning messages: 1: In solve_model(chem.name = chem.name, chem.cas = chem.cas, dtxsid = dtxsid, : Rblood2plasma not recalculated. Set recalc.blood2plasma to TRUE if desired. 2: In solve_model(chem.name = chem.name, chem.cas = chem.cas, dtxsid = dtxsid, : Clearance not recalculated. Set recalc.clearance to TRUE if desired. > > p <- parameterize_fetal_pbtk(chem.name='bisphenol a')[sort(names(parameterize_fetal_pbtk(chem.name='bisphenol a')))] > # Try to standardize order of variable names > for (this.param in names(p)[order(toupper(names(p)))]) cat(paste(this.param,": ",p[[this.param]],"\r\n",sep="")) adipose_density: 0.95 amnf_density: 1.01 arterial_blood_fraction: 0.16 brain_density: 1.04 BW: 61.1 BW_cubic_theta1: -0.010614 BW_cubic_theta2: 0.029161 BW_cubic_theta3: -0.00050203 Clint: 29.51 Clint.dist: 19.9,10.2,29.8,0 Clmetabolismc: 0.1219 Fabsgut: 0.6441 fblood_weight_ratio: 80 fBW_gompertz_theta0: 0.0018282 fBW_gompertz_theta1: 1.1735 fBW_gompertz_theta2: 0.077577 ffmx_density: 1.1 fhematocrit_cubic_theta1: 4.5061 fhematocrit_cubic_theta2: -0.18487 fhematocrit_cubic_theta3: 0.0026766 Fhep.assay.correction: 0.6744 Fraction_unbound_plasma_fetus: 0.03111 Funbound.plasma: 0.02434 Funbound.plasma.adjustment: 0.6321 Funbound.plasma.dist: 0.0385,0.0334,0.0456 gut_density: 1.04 hematocrit_quadratic_theta0: 39.192 hematocrit_quadratic_theta1: -0.10562 hematocrit_quadratic_theta2: -0.00071045 Kadipose2pu: 230.8 Kfbrain2pu: 214.6 Kfgut2pu: 491.7 Kfkidney2pu: 739 Kfliver2pu: 686.6 Kflung2pu: 282.6 Kfplacenta2pu: 147.2 Kfrbc2pu: 301.9 Kfrest2pu: 190.7 Kfthyroid2pu: 268.8 Kgut2pu: 520.6 kgutabs: 0.3711 kidney_density: 1.05 Kkidney2pu: 800.3 Kliver2pu: 753.2 Klung2pu: 292.2 Kplacenta2pu: 147.9 Krbc2pu: 301.5 Krest2pu: 200.9 Kthyroid2pu: 268.9 liver_density: 1.05 lung_density: 1.05 MA: 28840 million.cells.per.gliver: 110 MW: 228.3 pH_Plasma_fet: 7.28 pH_Plasma_mat: 7.38 pKa_Accept: NA pKa_Donor: 9.46 placenta_density: 1.02 Pow: 2089 pre_pregnant_BW: 61.1 Qadipose_percent_initial: 8.5 Qadipose_percent_terminal: 7.8 Qbrain_percent: 12 Qcardiac_cubic_theta0: 301.78 Qcardiac_cubic_theta1: 3.2512 Qcardiac_cubic_theta2: 0.15947 Qcardiac_cubic_theta3: -0.0047059 Qfbrain_percent: 14.3 Qfda_logistic_theta0: 1125.3 Qfda_logistic_theta1: 0.18031 Qfda_logistic_theta2: 35.939 Qfdv_gompertz_theta0: 1.892 Qfdv_gompertz_theta1: 0.098249 Qfdv_gompertz_theta2: 0.0064374 Qfgut_percent: 6.8 Qfkidney_percent: 5.4 Qfliver_percent: 6.5 Qflvtl_logistic_theta0: 506.3 Qflvtl_logistic_theta1: 0.21916 Qflvtl_logistic_theta2: 30.231 Qfnonplacental_percent: 75 Qfplacenta_logistic_theta0: 262.2 Qfplacenta_logistic_theta1: 0.22183 Qfplacenta_logistic_theta2: 28.784 Qfrvtl_logistic_theta0: 2466.5 Qfrvtl_logistic_theta1: 0.14837 Qfrvtl_logistic_theta2: 43.108 Qfthyroid_percent: 1.5 Qgfr_quadratic_theta0: 113.73 Qgfr_quadratic_theta1: 3.5784 Qgfr_quadratic_theta2: -0.067272 Qgut_percent: 16 Qgut_percent_initial: 17 Qgut_percent_terminal: 12.5 Qkidney_cubic_theta0: 53.248 Qkidney_cubic_theta1: 3.6447 Qkidney_cubic_theta2: -0.15357 Qkidney_cubic_theta3: 0.0016968 Qkidney_percent: 18 Qliver_percent_initial: 27 Qliver_percent_terminal: 20 Qplacenta_linear_theta1: 0.059176 Qthyroid_percent_initial: 1.5 Qthyroid_percent_terminal: 1.1 term: 40 thyroid_density: 1.05 Vamnf_logistic_theta0: 822.34 Vamnf_logistic_theta1: 0.26988 Vamnf_logistic_theta2: 20.15 venous_blood_fraction: 0.595 Vgutc: 0.019 Vkidneyc: 0.00458 Vliverc: 0.0233 Vlungc: 0.0158 Vplacenta_cubic_theta1: -1.7646 Vplacenta_cubic_theta2: 0.91775 Vplacenta_cubic_theta3: -0.011543 Vplasma_mod_logistic_theta0: 1.2406 Vplasma_mod_logistic_theta1: 0.31338 Vplasma_mod_logistic_theta2: 17.813 Vplasma_mod_logistic_theta3: 2.4958 Vthyroidc: 0.000283 Wadipose_linear_theta0: 17.067 Wadipose_linear_theta1: 0.14937 Wfbrain_gompertz_theta0: 0.01574 Wfbrain_gompertz_theta1: 0.70707 Wfbrain_gompertz_theta2: 0.064827 Wfgut_gompertz_theta0: 0.00081828 Wfgut_gompertz_theta1: 0.65028 Wfgut_gompertz_theta2: 0.047724 Wfkidney_gompertz_theta0: 6.3327e-05 Wfkidney_gompertz_theta1: 1.0409 Wfkidney_gompertz_theta2: 0.076435 Wfliver_gompertz_theta0: 0.0074774 Wfliver_gompertz_theta1: 0.65856 Wfliver_gompertz_theta2: 0.061662 Wflung_gompertz_theta0: 0.00030454 Wflung_gompertz_theta1: 1.0667 Wflung_gompertz_theta2: 0.084604 Wfthyroid_gompertz_theta0: 0.0038483 Wfthyroid_gompertz_theta1: 0.30799 Wfthyroid_gompertz_theta2: 0.0398 > > quit("no") > proc.time() user system elapsed 8.23 0.40 8.62