NAME App::DrugUtils - Utilities related to drugs VERSION This document describes version 0.001 of App::DrugUtils (from Perl distribution App-DrugUtils), released on 2023-10-19. DESCRIPTION This distributions provides the following command-line utilities: * tabulate-drug-concentration FUNCTIONS tabulate_drug_concentration Usage: tabulate_drug_concentration(%args) -> [$status_code, $reason, $payload, \%result_meta] Tabulate drug concentration over time, with one or more dosages of a drug of specified half-life. Examples: * Tabulate concentration of tadalafil (half-life 17.5h) over 1 week period if we give 2.5mg each day for a week: tabulate_drug_concentration( dosages => [ "2.5\@0", "2.5\@24", "2.5\@48", "2.5\@72", "2.5\@96", "2.5\@120", "2.5\@144", ], half_life => 17.5, period => 144 ); Result: [ 200, "OK", [ { time => 0, concentration => "2.500" }, { time => 1, concentration => 2.403 }, { time => 2, concentration => "2.310" }, { time => 3, concentration => "2.220" }, { time => 4, concentration => 2.134 }, { time => 5, concentration => 2.051 }, { time => 6, concentration => 1.971 }, { time => 7, concentration => 1.895 }, { time => 8, concentration => 1.821 }, { time => 9, concentration => "1.750" }, { time => 10, concentration => 1.682 }, { time => 11, concentration => 1.617 }, { time => 12, concentration => 1.554 }, { time => 13, concentration => 1.494 }, { time => 14, concentration => 1.436 }, { time => 15, concentration => "1.380" }, { time => 16, concentration => 1.327 }, { time => 17, concentration => 1.275 }, { time => 18, concentration => 1.225 }, { time => 19, concentration => 1.178 }, { time => 20, concentration => 1.132 }, { time => 21, concentration => 1.088 }, { time => 22, concentration => 1.046 }, { time => 23, concentration => 1.005 }, { time => 24, concentration => 3.466 }, { time => 25, concentration => 3.332 }, { time => 26, concentration => 3.202 }, { time => 27, concentration => 3.078 }, { time => 28, concentration => 2.958 }, { time => 29, concentration => 2.844 }, { time => 30, concentration => 2.733 }, { time => 31, concentration => 2.627 }, { time => 32, concentration => 2.525 }, { time => 33, concentration => 2.427 }, { time => 34, concentration => 2.333 }, { time => 35, concentration => 2.242 }, { time => 36, concentration => 2.155 }, { time => 37, concentration => 2.071 }, { time => 38, concentration => 1.991 }, { time => 39, concentration => 1.914 }, { time => 40, concentration => 1.839 }, { time => 41, concentration => 1.768 }, { time => 42, concentration => 1.699 }, { time => 43, concentration => 1.633 }, { time => 44, concentration => "1.570" }, { time => 45, concentration => 1.509 }, { time => 46, concentration => "1.450" }, { time => 47, concentration => 1.394 }, { time => 48, concentration => "3.840" }, { time => 49, concentration => 3.691 }, { time => 50, concentration => 3.547 }, { time => 51, concentration => "3.410" }, { time => 52, concentration => 3.277 }, { time => 53, concentration => "3.150" }, { time => 54, concentration => 3.028 }, { time => 55, concentration => "2.910" }, { time => 56, concentration => 2.797 }, { time => 57, concentration => 2.688 }, { time => 58, concentration => 2.584 }, { time => 59, concentration => 2.484 }, { time => 60, concentration => 2.387 }, { time => 61, concentration => 2.294 }, { time => 62, concentration => 2.205 }, { time => 63, concentration => "2.120" }, { time => 64, concentration => 2.037 }, { time => 65, concentration => 1.958 }, { time => 66, concentration => 1.882 }, { time => 67, concentration => 1.809 }, { time => 68, concentration => 1.739 }, { time => 69, concentration => 1.671 }, { time => 70, concentration => 1.606 }, { time => 71, concentration => 1.544 }, { time => 72, concentration => 3.984 }, { time => 73, concentration => 3.829 }, { time => 74, concentration => 3.681 }, { time => 75, concentration => 3.538 }, { time => 76, concentration => "3.400" }, { time => 77, concentration => 3.268 }, { time => 78, concentration => 3.141 }, { time => 79, concentration => 3.019 }, { time => 80, concentration => 2.902 }, { time => 81, concentration => 2.789 }, { time => 82, concentration => 2.681 }, { time => 83, concentration => 2.577 }, { time => 84, concentration => 2.477 }, { time => 85, concentration => 2.381 }, { time => 86, concentration => 2.288 }, { time => 87, concentration => 2.199 }, { time => 88, concentration => 2.114 }, { time => 89, concentration => 2.032 }, { time => 90, concentration => 1.953 }, { time => 91, concentration => 1.877 }, { time => 92, concentration => 1.804 }, { time => 93, concentration => 1.734 }, { time => 94, concentration => 1.667 }, { time => 95, concentration => 1.602 }, { time => 96, concentration => "4.040" }, { time => 97, concentration => 3.883 }, { time => 98, concentration => 3.732 }, { time => 99, concentration => 3.587 }, { time => 100, concentration => 3.448 }, { time => 101, concentration => 3.314 }, { time => 102, concentration => 3.185 }, { time => 103, concentration => 3.062 }, { time => 104, concentration => 2.943 }, { time => 105, concentration => 2.828 }, { time => 106, concentration => 2.719 }, { time => 107, concentration => 2.613 }, { time => 108, concentration => 2.512 }, { time => 109, concentration => 2.414 }, { time => 110, concentration => "2.320" }, { time => 111, concentration => "2.230" }, { time => 112, concentration => 2.144 }, { time => 113, concentration => "2.060" }, { time => 114, concentration => "1.980" }, { time => 115, concentration => 1.903 }, { time => 116, concentration => "1.830" }, { time => 117, concentration => 1.758 }, { time => 118, concentration => "1.690" }, { time => 119, concentration => 1.625 }, { time => 120, concentration => 4.061 }, { time => 121, concentration => 3.904 }, { time => 122, concentration => 3.752 }, { time => 123, concentration => 3.606 }, { time => 124, concentration => 3.466 }, { time => 125, concentration => 3.332 }, { time => 126, concentration => 3.202 }, { time => 127, concentration => 3.078 }, { time => 128, concentration => 2.958 }, { time => 129, concentration => 2.844 }, { time => 130, concentration => 2.733 }, { time => 131, concentration => 2.627 }, { time => 132, concentration => 2.525 }, { time => 133, concentration => 2.427 }, { time => 134, concentration => 2.333 }, { time => 135, concentration => 2.242 }, { time => 136, concentration => 2.155 }, { time => 137, concentration => 2.071 }, { time => 138, concentration => 1.991 }, { time => 139, concentration => 1.914 }, { time => 140, concentration => 1.839 }, { time => 141, concentration => 1.768 }, { time => 142, concentration => 1.699 }, { time => 143, concentration => 1.633 }, { time => 144, concentration => "4.070" }, ], { "table.fields" => ["time", "concentration"] }, ] * Same as previous example, but uses -r: tabulate_drug_concentration(half_life => 17.5, period => 144, regular_dosage => "2.5\@24"); Drug concentration over time is calculated using a simple equation: D(t) = D(0) * exp(-t/T) where T is a constant and can be calculated from half-life where D(t)/D(0) = 0.5: 0.5 = 1 * exp(-half_life / T) log(0.5) = -half_life / T T = half_life / log(0.5) When given multiple dosages, each dosage's concentration will be added up. By default, hourly concentration will be calculated for 10 half-lives. This function is not exported. Arguments ('*' denotes required arguments): * dosages => *array[str]* Dosages, each of which in the form of d@t (amount followed by comma and time in hour). * half_life* => *ufloat* Drug half life, in hours. * period => *uint* How many hours to tabulate. If unspecified, defaults to 10 * half-life. * regular_dosage => *str* Regular dosage, in the form of d@t (amount followed by at sign and period in hour, meaning dosage given every t hours). Returns an enveloped result (an array). First element ($status_code) is an integer containing HTTP-like status code (200 means OK, 4xx caller error, 5xx function error). Second element ($reason) is a string containing error message, or something like "OK" if status is 200. Third element ($payload) is the actual result, but usually not present when enveloped result is an error response ($status_code is not 2xx). Fourth element (%result_meta) is called result metadata and is optional, a hash that contains extra information, much like how HTTP response headers provide additional metadata. Return value: (any) HOMEPAGE Please visit the project's homepage at . SOURCE Source repository is at . SEE ALSO Somewhat related utilities: App::VitaminUtils, App::MineralUtils. AUTHOR perlancar CONTRIBUTING To contribute, you can send patches by email/via RT, or send pull requests on GitHub. Most of the time, you don't need to build the distribution yourself. You can simply modify the code, then test via: % prove -l If you want to build the distribution (e.g. to try to install it locally on your system), you can install Dist::Zilla, Dist::Zilla::PluginBundle::Author::PERLANCAR, Pod::Weaver::PluginBundle::Author::PERLANCAR, and sometimes one or two other Dist::Zilla- and/or Pod::Weaver plugins. Any additional steps required beyond that are considered a bug and can be reported to me. COPYRIGHT AND LICENSE This software is copyright (c) 2023 by perlancar . This is free software; you can redistribute it and/or modify it under the same terms as the Perl 5 programming language system itself. BUGS Please report any bugs or feature requests on the bugtracker website When submitting a bug or request, please include a test-file or a patch to an existing test-file that illustrates the bug or desired feature.