f ( t ) = 1 e − b e − c t {\displaystyle f(t)=a\mathrm {e} ^{-b\mathrm {e} ^{-ct}}} どこ
a は漸近線である。リム t → ∞ 1 e − b e − c t = 1 e 0 = 1 {\textstyle \lim _{t\to \infty }a\mathrm {e} ^{-b\mathrm {e} ^{-ct}}=a\mathrm {e} ^{0}=a} bは x 軸方向の変位を設定します(グラフを左または右に平行移動します)。cは 成長率(y スケーリング)を設定しますeはオイラー数 (e = 2.71828...)です。
導出 関数曲線は、絶対死亡率(減衰率)が現在のサイズとともに指数関数的に減少するというゴンペルツ・マケハム死亡率法則から導き出すことができる。数学的には、
k r ∝ 1 y ( t ) {\displaystyle k^{r}\propto {\frac {1}{y(t)}}}
どこ
r = y ′ ( t ) y ( t ) {\textstyle r={\frac {y'(t)}{y(t)}}} 成長率k は任意の定数である。
使用例 ゴンペルツ曲線の使用例としては、以下のようなものがあります。
携帯電話の 普及は、当初はコストが高かったため普及が遅く、その後急速な成長期を経て、飽和状態に達したため普及が鈍化した[ 6 ]。 限られた空間での人口増加は、出生率が最初は上昇し、資源の限界に達すると低下する[ 7 ]。 腫瘍の成長のモデリング[ 8 ] 宇宙 における水素の再電離のモデリング[ 9 ] 金融における市場影響のモデリング[ 10 ] および集約された地方融資の動態[ 11 ] 捕食者と被食者の関係に着目し、被食動物の個体数増加を詳細に分析する。 集団内の細菌細胞のモデリング 疾病の蔓延を調査する 英語版ウィキペディアのサイズは、ある程度ゴンペルツ関数と修正関数でモデル化できる[ 12 ]。
アプリケーション
ゴンペルツ曲線 個体群生物学で は、特にゴンペルツ関数が重要視されています。この関数は、特定の生物集団の急速な増加を記述するのに非常に有用であり、環境収容力 (細胞数/個体数)が決定されれば、最終的な水平漸近線も説明することができます。
以下のようにモデル化されます。
N ( t ) = N 0 exp ( ln ( N 私 / N 0 ) ( 1 − exp ( − b t ) ) ) {\displaystyle N(t)=N_{0}\exp(\ln(N_{I}/N_{0})(1-\exp(-bt)))}
どこ:
t {\textstyle t} 時間が来たN 0 {\textstyle N_{0}} は細胞の初期密度ですN 私 {\textstyle N_{I}} プラトー細胞/個体群密度はb {\textstyle b} 腫瘍の初期増殖速度この関数はプラトー細胞数を考慮することで、実際の集団動態を 正確に模倣するのに役立ちます。この関数は、集団の成長を一般的に詳細に記述する最も広く受け入れられている慣例であるシグモイド関数にも準拠しています。さらに、この関数は、 対数期 を経て急速に増殖する細菌や癌細胞の集団でよく見られる初期成長率を利用します。人気があるにもかかわらず、腫瘍成長の初期成長率は、患者に存在するさまざまな微小環境、または集団生物学の場合のさまざまな環境要因を考慮すると、事前に決定することが困難です。癌患者の場合、年齢、食事、民族、遺伝的素因、代謝、ライフスタイル、 転移 の起源などの要因が、腫瘍成長率の決定に役割を果たします。収容力もこれらの要因に基づいて変化すると予想されるため、このような現象を記述することは困難です。
代謝機能は、生物体内の代謝速度を説明することに特に重点を置いています。この機能は腫瘍細胞のモニタリングに適用できます。代謝速度は動的で柔軟性が高く、がんの増殖をより正確に詳細に把握できます。代謝曲線は、組織の維持と生成に体が供給するエネルギーを考慮に入れています。このエネルギーは代謝とみなすことができ、細胞分裂において特定のパターンに従います。エネルギー保存則は 、質量や発生時間の違いに関わらず、このような成長をモデル化するために使用できます。すべての分類群は 同様の成長パターンを共有しており、このモデルは、結果として、腫瘍発生の基盤となる細胞分裂を考慮しています。
B = ∑ C ( N C B C ) + ( E C d N C d t ) {\displaystyle B=\sum _{C}(N_{C}B_{C})+\left(E_{C}{\operatorname {d} \!N_{C} \over \operatorname {d} \!t}\right)}
B {\textstyle B} = 生物が安静時に使用するエネルギーN C {\textstyle N_{C}} =当該生物の細胞数B C {\textstyle B_{C}} = 個々の細胞の代謝率N C B C {\textstyle N_{C}B_{C}} =既存の組織 を維持するために必要なエネルギーE C {\textstyle E_{C}} =個々の細胞から新しい組織を作り出すのに必要なエネルギー安静時のエネルギー消費量と代謝率によるエネルギー消費量を区別することで、モデルは成長速度をより正確に決定できます。安静時のエネルギー消費量は組織の維持に必要なエネルギーよりも少なく、両者を合わせると既存組織の維持に必要なエネルギーを表します。これら2つの要素と、新しい組織の生成に必要なエネルギーを組み合わせることで、成長速度を包括的に把握でき、さらに、遅延期 を正確に表現することができます。
腫瘍の増殖 1960年代にAK Laird [ 13 ] は、腫瘍の成長データに適合させるために初めてゴンペルツ曲線をうまく使用しました。実際、腫瘍は栄養素の利用可能性が制限された限られた空間で成長する細胞集団です。腫瘍のサイズをX(t)とすると、ゴンペルツ曲線は次のように表すことができます。
X ( t ) = K exp ( ログ ( X ( 0 ) K ) exp ( − α t ) ) {\displaystyle X(t)=K\exp \left(\log \left({\frac {X(0)}{K}}\right)\exp \left(-\alpha t\right)\right)} どこ:
X ( 0 ) {\textstyle X(0)} は、観察開始時点での腫瘍の大きさです。K {\textstyle K} これは収容力、つまり利用可能な栄養素で到達できる最大サイズです。実際には次のようになります。リム t → + ∞ X ( t ) = K {\displaystyle \lim _{t\rightarrow +\infty }X(t)=K} X(0)>0 とは独立して。治療などがない場合、通常は X(0) < K であるのに対し、治療がある場合は X(0) > K になる可能性があることに注意してください。
α {\textstyle \alpha } これは細胞の増殖能力に関連する定数である。ログ ( ) {\textstyle \log()} 自然対数 を指します。X(t)のダイナミクスは、ゴンペルツ微分方程式によって支配されることが示される。
X ′ ( t ) = α ログ ( K X ( t ) ) X ( t ) {\displaystyle X^{\prime }(t)=\alpha \log \left({\frac {K}{X(t)}}\right)X(t)}
つまり、分解すると次の形式になります。
X ′ ( t ) = F ( X ( t ) ) X ( t ) 、 と F ′ ( X ) ≤ 0 、 {\displaystyle X^{\prime }(t)=F\left(X(t)\right)X(t),\quad {\mbox{ただし}}\quad F^{\prime }(X)\leq 0,}
F(X) は細胞集団の瞬間的な増殖率であり、その減少傾向は、ロジスティック成長率と同様に、細胞集団の増加に伴う栄養素の競合によるものです。しかし、根本的な違いがあります。ロジスティックの場合、細胞集団が小さいときの増殖率は有限です。
F ( X ) = α ( 1 − ( X K ) ν ) ⇒ F ( 0 ) = α < + ∞ {\displaystyle F(X)=\alpha \left(1-\left({\frac {X}{K}}\right)^{\nu }\right)\Rightarrow F(0)=\alpha <+\infty }
一方、ゴンペルツの場合、増殖率は無制限である。
リム X → 0 + F ( X ) = リム X → 0 + α ログ ( K X ) = + ∞ \lim_{X\rightarrow 0^{+}}F(X)=\lim_{X\rightarrow 0^{+}}\alpha \log \left({\frac {K}{X}}\right)=+\infty }
Steel [ 14 ] や Wheldon [ 15 ] が指摘したように、細胞集団の増殖速度は最終的に細胞分裂 時間によって制限されます。したがって、これは Gompertz 方程式が小さな腫瘍の成長をモデル化するのに適していないという証拠となる可能性があります。さらに最近では[ 16 ] 、免疫系 との相互作用を含めると、Gompertz や F(0) が制限されない他の法則では免疫監視の可能性が排除されることが指摘されています。
Fornalski らによる理論的研究[ 17 ] は、放物線関数の方が適切なごく初期の段階を除き、癌の成長に対する Gompertz 曲線の生物物理学的基礎を示した。彼らはまた、Gompertz 曲線が、癌の動態関数の幅広いファミリーの中で最も典型的なケースを記述していることを発見した。
ゴンペルツ成長とロジスティック成長 ゴンペルツ微分方程式
X ′ ( t ) = α ログ ( K X ( t ) ) X ( t ) {\displaystyle X^{\prime }(t)=\alpha \log \left({\frac {K}{X(t)}}\right)X(t)}
これは一般化ロジスティック微分方程式 の極限ケースである。
X ′ ( t ) = α ν ( 1 − ( X ( t ) K ) 1 ν ) X ( t ) {\displaystyle X^{\prime }(t)=\alpha \nu \left(1-\left({\frac {X(t)}{K}}\right)^{\frac {1}{\nu }}\right)X(t)}
(どこν > 0 {\displaystyle \nu >0} は正の実数であるため
リム ν → + ∞ ν ( 1 − x 1 / ν ) = − ログ ( x ) {\displaystyle \lim _{\nu \rightarrow +\infty }\nu \left(1-x^{1/\nu }\right)=-\log \left(x\right)} 。
さらに、一般化ロジスティック関数の グラフ には、
X ( t ) = ( ν ν + 1 ) ν K {\displaystyle X(t)=\left({\frac {\nu }{\nu +1}}\right)^{\nu }K}
そして、ゴンペルツ関数のグラフでは、
X ( t ) = K e = K ⋅ リム ν → + ∞ ( ν ν + 1 ) ν {\displaystyle X(t)={\frac {K}{e}}=K\cdot \lim _{\nu \rightarrow +\infty }\left({\frac {\nu }{\nu +1}}\right)^{\nu }} 。
ゴンプエックスの成長法則 上記の考察に基づき、ウェルドン[ 15 ] は、ゴンペルツの法則をわずかに修正した腫瘍増殖の数理モデルであるGomp-Exモデルを提案した。Gomp-Exモデルでは、初期段階では資源の競合がないため、細胞集団は指数関数的に増加すると仮定されている。しかし、ある臨界サイズ閾値が存在する。X C {\displaystyle X_{C}} そのため、X > X C {\displaystyle X>X_{C}} 資源をめぐる競争がないという前提は、ほとんどのシナリオにおいて当てはまります。しかし、制約要因 によって影響を受ける可能性があり、そのためにはサブファクター変数を作成する必要があります。
成長はゴンペルツの法則に従う。
F ( X ) = 最大 ( 1 、 α ログ ( K X ) ) {\displaystyle F(X)=\max \left(a,\alpha \log \left({\frac {K}{X}}\right)\right)}
となることによって:
X C = K exp ( − 1 α ) 。 {\displaystyle X_{C}=K\exp \left(-{\frac {a}{\alpha }}\right).}
ここでは、いくつかの数値推定値[ 15 ] がある。X C {\displaystyle X_{C}} :
X C ≈ 10 9 {\textstyle X_{C}\approx 10^{9}} ヒト腫瘍の場合X C ≈ 10 6 {\textstyle X_{C}\approx 10^{6}} マウス 腫瘍の場合
逆ゴンペルツ関数 ゴンペルツ関数は1対1対応(全単射関数 とも呼ばれる)であるため、その逆関数は従来の関数表記法で単一の 連続関数 として明示的に表現できます。次の形式のゴンペルツ関数が与えられた場合:
f ( t ) = 1 e − e b − c t + d {\displaystyle f(t)=a\mathrm {e} ^{-\mathrm {e} ^{b-ct}}+d} どこ
d は水平漸近線の底辺です。リム t → − ∞ 1 e − e b − c t + d = 1 e − ∞ + d = d {\textstyle \lim _{t\to -\infty }a\mathrm {e} ^{-\mathrm {e} ^{b-ct}}+d=a\mathrm {e} ^{-\infty }+d=d} a は基線から第 2 漸近線までの距離です。リム t → ∞ 1 e − e b − c t + d = 1 e 0 + d = 1 + d {\textstyle \lim _{t\to \infty }a\mathrm {e} ^{-\mathrm {e} ^{b-ct}}+d=a\mathrm {e} ^{0}+d=a+d} bは x 軸方向の変位を設定します(グラフを左または右に平行移動します)。cは 成長率(y スケーリング)を設定しますeは オイラー数 (e = 2.71828...)です。対応する逆関数は次のように表すことができます。
f − 1 ( t ) = 1 c [ b − ln ( ln ( 1 t − d ) ) ] {\displaystyle f^{-1}(t)={\frac {1}{c}}\left[b-\ln \left(\ln \left({\frac {a}{t-d}}\right)\right)\right]}
逆関数は、順方向ゴンペルツ関数のように水平ではなく垂直な2つの漸近線の間の実数 の範囲でのみ数値を生成します。垂直漸近線によって定義される範囲外では、逆関数は負の数の対数を計算する必要があります。このような理由などから、特にフィッティングに使用できるデータ点が比較的少ない場合は、逆ゴンペルツ関数をデータに直接当てはめようとすると、多くの場合非現実的です。代わりに、データの転置関係を順方向ゴンペルツ関数に当てはめ、上記で示した2つの関数の関係を使用して、それを同等の逆関数に変換することができます。
このように、逆関数には多くの用途があります。例えば、特定のELISAアッセイで は、濃度をゴンペルツ関数で吸光度に非常に良く適合させることができる 標準曲線 が用いられます。このように標準値をゴンペルツ関数に適合させた後、アッセイ中のサンプルの未知濃度を測定された吸光度から計算するには、標準曲線を適合させた際に得られたゴンペルツ関数の逆関数を用います。
参考文献 ↑ カークウッド、TBL (2015) 「死の解読:ゴンペルツ(1825)『人間の死亡法則を表現する機能の性質、および生命の偶発性の価値を決定する新しい方法』の解説」 「 . Philosophical Transactions of the Royal Society of London B . 370 (1666). doi : 10.1098/rstb.2014.0379 . PMC 4360127 . PMID 25750242 . ↑ ゴンペルツ、ベンヤミン( 1825)。 「人間の死亡法則を表現する機能の性質と、生命の偶発性の価値を決定する新しい方法について」 。 ロンドン 王立 協会の哲学的トランザクション 。115 : 513–585。doi : 10.1098/rstl.1825.0026。S2CID 145157003。 2021年1 月 20日にオリジナルから アーカイブ 。 2020年9月4日 に 取得 。 ↑ ド・モアブル、アブラハム (1725)。 年金保険… 。ロンドン、イングランド:フランシス・フェイラム、ベンジャミン・モット、W・ピアソン。 第2版は1743年に、第3版は1750年に、第4版は1752年に発行された。↑ Greenwood, M. (1928). "生物学的観点から見た死亡率の法則" . Journal of Hygiene . 28 (3): 267– 294. doi : 10.1017/S002217240000961X . PMC 2167778 . PMID 20475000 . ↑ Makeham, William Matthew (1860). 「死亡率の法則と年金表の作成について」 。The Assurance Magazine, and Journal of the Institute of Actuaries . 8 (6): 301–310 . doi : 10.1017/S204616580000126X 。 ↑ Islam T、Fiebig DG、Meade N ( 2002)。「限られたデータによる多国籍通信需要のモデリング」。International Journal of Forecasting。18 ( 4): 605–624。doi : 10.1016/S0169-2070(02) 00073-0 。 ↑ Zwietering MH、Jongenburger I、Rombouts FM、van 't Riet K (1990 年 6 月)。 「 細菌 増殖 曲線 のモデリング」 。 応用 環境微生物学 。56 ( 6 ) : 1875–81。Bibcode : 1990ApEnM..56.1875Z。doi : 10.1128 / AEM.56.6.1875-1881.1990。PMC 184525。PMID 16348228 。 。↑ Sottoriva A、Verhoeff JJ 、 Borovski T、McWeeney SK、Naumov L、Medema JP、et al. (2010 年 1 月)。 「 癌幹細胞腫瘍モデルは、浸潤性形態と表現型の異質性の増加を明らかにする」 。Cancer Research。70 ( 1 ) : 46– 56。doi : 10.1158/0008-5472.CAN-09-3663。PMID 20048071 。 ↑ Montero-Camacho, Paulo; Li, Yin; Cranmer, Miles (2024). "Five parameters are all you need (in ΛCDM)". arXiv : 2405.13680 [ astro-ph.CO ]. ↑ Caravelli F、Sindoni L 、Caccioli F 、Ududec C (2016 年8 月)。「有限 の 収容力を持つ最適な成長軌道 」 。Physical Review E。94 ( 2–1 ) 022315。arXiv : 1510.05123。Bibcode : 2016PhRvE..94b2315C。doi : 10.1103 / PhysRevE.94.022315。PMID 27627325。S2CID 35578084 。 。↑ Rocha LS、Rocha FS、Souza TT (2017-10-05)。 「 あなたの国の公共部門は拡散借入者ですか?ブラジルから の 実証 的 証拠 」 。PLOS ONE。12 ( 10 ) e0185257。arXiv : 1604.07782。Bibcode : 2017PLoSO..1285257R。doi : 10.1371 / journal.pone.0185257。PMC 5628819。PMID 28981532 。 ↑ 「Wikipedia:Wikipediaの成長モデル化」 、 Wikipedia 、2023年3月18日、 2025年3月6日のオリジナルから アーカイブ、 2023年3月23日 取得 ↑ Laird AK (1964年9月) 「腫瘍増殖のダイナミクス」 British Journal of Cancer 13 ( 3): 490– 502. doi : 10.1038/bjc.1964.55 . PMC 2071101 . PMID 14219541 . ↑ Steel GG (1977). 腫瘍の成長速度論 . オックスフォード: クラレンドン・プレス. ISBN 0-19-857388-X 。1 2 3 Wheldon TE (1988). Mathematical Models in Cancer Research . Bristol: Adam Hilger. ISBN 0-85274-291-6 。↑ d'Onofrio A (2005). "腫瘍免疫系競合と免疫療法をモデル化するための一般的な枠組み:数学的解析と生物医学的推論". Physica D . 208 ( 3– 4): 220– 235. arXiv : 1309.3337 . Bibcode : 2005PhyD..208..220D . doi : 10.1016/j.physd.2005.06.032 . S2CID 15031322 . ↑ フォルナルスキー KW、レシュチェンスカ J、ドブジンスキー L、ヴィソッキ P、ヤニアク MK (2020)。 「ガンペルツ増殖の法則の考えられる原因: 腫瘍増殖の機構モデリング」 。 アクタ フィジカ ポロニカ A . 138 (6): 854–862 。 Bibcode : 2020AcPPA.138..854F 。 土井 : 10.12693/APhysPolA.138.854 。