薬剤師のプログラミング学習日記

プログラミングやコンピュータに関する記事を書いていきます

eGFRとクレアチニンクリアランスの計算式をグラフ描画する

腎機能評価の指標としてeGFR(推算糸球体濾過量)やCCr(クレアチニンクリアランス)、シスタチンC(Cys-C)などがありますが、病院等で働く薬剤師であれば、仕事中に薬の投与量のチェックをするためにこれらの値を計算することも多いかと思います。 ネットで…

薬の蓄積率について-投与間隔と半減期からわかること-

ネットで蓄積率を検索すると『投与間隔/消失半減期が3以下なら定常状態のある薬、4以上なら定常状態のない薬』などの記事が見つかります。今回はこの3や4という数字の意味について見ていきたいと思います。 蓄積率とは 蓄積率の式をグラフでみる 具体的な例…

ニュートン法で方程式の近似解を求める

前回「トラフ値からkeを推定する」で定常状態のCss_minを求める式からkeを求めました。このときscipy.optimizeモジュールのnewton()を使いましたが、これはニュートン法(ニュートン・ラフソン法)というアルゴリズムを実装した関数です(引数に導関数を与え…

トラフ値からkeを推定する

今回の記事は参考文献『4ステップ 臨床力UPエクササイズ4 TDM領域』の『Step3 薬物治療の個別化(薬物投与設計)の実践』の項目を基に書いています。 この項目の概要としては、喘息症状が悪化した患者のスクリーニングを通して薬物クリアランスの変動を考…

経口繰り返し投与のグラフを描く

n回繰り返し投与後の血中濃度 経口投与の定常状態における最高・最低血中濃度とAUC Sympyによる微分と定積分 経口繰り返し投与のクラス 血中濃度推移の描画 参考 n回繰り返し投与後の血中濃度 投与間隔(タウ)でn回繰り返し投与後のt時間経過時の血中濃…

経口投与のフリップ・フロップをグラフで見る

Flip-Flop(フリップ・フロップ) 各条件で用いる数式 ke<kaの場合 ke>kaの場合 ke=kaの場合 TmaxとCmax グラフと文字の更新処理 フリップフロップを確認する 参考 Flip-Flop(フリップ・フロップ) 前回記事の最後に、消失相の傾きはkaとkeの大小関係によってどちらが反映さ</kaの場合>…

血中濃度データに経口投与モデルを当てはめる

以前1-コンパートメントと2-コンパートメントの静注モデルでもやりましたが、今回は前回記事の1次速度の吸収過程のある1-コンパートメントモデルでやってみようと思います。 血中濃度データの片対数プロット プロットから消失速度定数(ke)と吸収速度定数(ka…

経口投与1-コンパートメントモデルのグラフを描く

1次速度の吸収過程のある1-コンパートメントモデル 吸収過程と消失過程 参考 1次速度の吸収過程のある1-コンパートメントモデル 経口投与時の血中濃度推移はこのモデルで表すことができます。※ka: 吸収速度定数, ke: 消失速度定数 この連立微分方程式を解…

テイコプラニンの負荷投与のグラフを描く(2-コンパートメントモデル)

半減期の長い薬物の場合は定常状態への到達が遅れます。このような薬物の場合、開始時期に投与量や回数を増やすことで目標とする血中濃度への到達を早める投与設計を行うことがあります。 今回はこれまで使ってきた2-コンパートメント点滴静注モデルのプログ…

Pythonによるシミュレーション高速化のためのコードの改善

はじめに 定常状態の血中濃度を得る部分を最小限に 演算に関わる部分はNumpy配列のみを使う コード改善後の処理時間 はじめに 前回まで書いていたモンテカルロシミュレーションのプログラムを改善して、どれだけ処理時間が短縮できるかを試みよう、という記…