ロジスティック‐ほうていしき〔‐ハウテイシキ〕【ロジスティック方程式】
ロジスティック方程式
出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2024/05/01 01:58 UTC 版)
ロジスティック方程式（ロジスティックほうていしき、英語：logistic equation[1]）は、生物の個体数の変化の様子を表す数理モデルの一種である。ある単一種の生物が一定環境内で増殖するようなときに、その生物の個体数（個体群サイズ）の変動を予測できる。人間の場合でいえば、人口の変動を表すモデルである。
1838年にベルギーの数学者ピエール＝フランソワ・フェルフルスト（Pierre-François Verhulst）によって、ロジスティック方程式は最初に発案された。フェルフルストは、1798年に発表されて大きな反響を呼んだトマス・ロバート・マルサスの『人口論』の不自然な点を解消するために、このモデルを考案した[2]。マルサスは『人口論』で、人口は原理的に指数関数的に増加することを指摘した[3]。しかし、実際には環境や資源は限られているため、人口の増加にはいずれブレーキがかかると考えるのが自然である。人口が増えるに連れて人口増加率は低減し、人口はどこかで飽和すると考えられる。ロジスティック方程式はこの点を取り入れて、生物の個体数増殖をモデル化したものである。フェルフルスト以後には、アメリカの生物学者レイモンド・パール（Raymond Pearl）が式を普及させた。
具体的には、ロジスティック方程式は
フィボナッチによるウサギのつがいの増殖問題 「個体群動態論」も参照生物の個体数の変動については古くから興味を持たれ、研究が行われてきた[10]。フィボナッチ数の発見に繋がったレオナルド・フィボナッチのウサギの個体数の問題が、おそらく最も古い個体数の数理モデルといわれる[11]。
生物の個体数の増え方に関する研究は、個体群生態学の分野に属する[12]。ここで、個体群とは簡単には、ある領域に生息している単一の種の個体の集まりのことを指す[13]。
この個体群の大きさ（個体群サイズ）の指標としては、個体群内の総個体数が使用される[14]。個体数の代わりに、領域の単位面積当たりの個体数である個体群密度や単位面積当たりの生物の総重量である生物量が、個体群サイズとして適切な指標となる場合もある[15]。人間でいえば、これらの指標は人口や人口密度に相当する[16]。
マルサスモデル
マルサスモデルによる個体数増加曲線の様子。赤色が m = 4、紫色が m = 2、藍色が m = 1。いずれも最初は N =1 だが、その後の急激な成長が見て取れる。 「マルサスモデル」も参照多くの生物では、親は多くの子孫を作るので、それがそのまま生き残ると仮定すれば、あっという間に莫大な個体数となる。ねずみ算など、数学的小話の種である[17]。まずはこのような単純なものが、生物個体数の増加モデルとして考えられる。
ある個体群において、時刻 t に個体数が N 体が存在しているとする。実際の生物個体数は不連続な値（整数）をとるものであるが、数学的扱いを簡便にするために、個体数は連続な値（実数）をとるものとする（1.5体といったような値も含める）ことがしばしば行われる[18]。実際の生物でいえば、個体数が多かったり各個体の世代が重なったりしていれば、このような近似も妥当性を帯びてくる[19][20]。個体数を連続な値とすれば、個体数の増加率は N の時間微分 dN/dt で表すことができる[21]。
さらに話を単純化するために、個体は環境を出入りしないという状況を想定する[22]。この場合、個体の出生と死亡という2つの要因のみによって個体数は増減する[22]。個体群の出生率が死亡率を上回っていれば、個体数は増え続けるということになる[10]。さらに簡略化するために出生率と死亡率を常に一定であるとする[10]。個体数当たりの出生率を b、個体数当たりの死亡率を d とすれば、個体数の増加率は差し引きした b − d に個体数 N を掛け合わせた値となる[23]。よって個体数増加率 dN/dt は
個体数 N と一個体当たり個体群増加率 m の関係 上記のようにマルサスモデルは非現実的な面を持つ。個体数が多くなると増加率が抑えられることを表現するために、個体数 N が増加するにつれて増加率 m が減少するモデルが考えられる[32]。また、個体数がある上限を超えたら増加率は負となり、個体数は減少に向かうと考えられる[1]。これらの点を簡単に表せば、比例定数 m を
縦軸が dN/dt、横軸が N のグラフ。ロジスティック方程式における、dN/dt と N の関係が示されている。 ロジスティック方程式における個体数増加率 dN/dt と個体数 N の関係に着目すれば、この関係は初等教育でも習う二次関数そのものとなっており、dN/dt と N のグラフは放物線を描く[45]。方程式を解析的に解かなくとも、N と dN/dt がどのような変化を起こすのかを、以下のようにグラフから読み解くこともできる。
まず、N = 0 と N = K のとき、dN/dt = 0 となる。すなわち、いくら時間が経過しても個体数は増加も減少もしない状態となる。このような状態を定常状態や平衡状態と呼ぶ[46]。N の値が 0 < N < K の範囲にあるとき、 dN/dt の値は様々だが、値が正なのか負なのかで言えば、正の値であることがわかる[47]。N の値が K < N となると、dN/dt は同じように負の値である[47]。言い換えれば、個体数が環境収容力内では常に個体数は増加するが、環境収容力を超えると個体数は減少へ転ずる、ということである[47]。
個体数増加率 dN/dt の変化をさらに細かく見てみる。N が 0 から K/2 まで増えると、その間 dN/dt の値も増加し続ける。N = K/2 は放物線の頂点であり、ここで dN/dt は極大値を迎える[48]。極大値は、N = K/2 を式に代入して dN/dt = rK/4 である[45]。N = K/2 を超えると dN/dt は減少し始め、N = K で 0 となる[48]。このような数値の変化から読み取れることの一つは、個体数が環境収容力のちょうど半分になったときに個体増加率は最大となる点である[48]。したがって、もし個体数の変化がロジスティック方程式に従うとしたら、増加率が最大になるときの個体数に注目することで環境収容力、すなわち最大個体数を予測できることになる[49]。
式の解
ロジスティック曲線
マルサスモデルによる指数関数的増加曲線（赤）とロジスティック曲線（青）の比較 ロジスティック方程式は非線形の微分方程式だが、標準的な微分方程式の解法である変数分離法を利用して解くことができる[50]。時間 t = 0 における初期個体数を N0 とすると、t の関数として以下の解が得られる[51][注釈 1]。
いくつかの N0 から始まるロジスティック曲線。N > 0 の範囲では、時間の経過に従って N は K に収束する。 
時間と個体数が負の場合も含めたロジスティック曲線の全体図。縦軸を N/K、横軸を rt として無次元化している。 横軸を t、縦軸を N とした平面上にロジスティック関数のグラフを描くと、曲線が描かれる。この曲線は前述のとおりにロジスティック曲線と呼ばれる。初期個体数が3つの範囲 N0 < 0, 0 < N0 < K, K < N0 のどれに該当するかによって、曲線の形状は大きく異なってくる[57]。ただし、N0 < 0 の範囲は負の個体数というものを意味するので、生物のモデルとしてはあまり意味がない[56]。時間 t = 0 から t → ∞ の極限までのロジスティック曲線の様相は、それぞれの N0 の値ごとに、以下のようになっている。
まず N0 が環境収容力の半分以下（0 < N0 < K/2 ）の場合、初期状態の点 (t = 0, N = N0) から始まる曲線は、ゆっくりと右肩上がりに登っていく[58]。t が増加するにつれて、曲線の傾き（個体数増加率）は増加していき、曲線は加速度的に立ち上がっていく[48]。しかし、ある時点で曲線は変曲点を迎え、傾きの増加は止む[59]。その後は、傾きは減少しだし、曲線は横倒しになっていく[60]。そして最終的には、傾きは0になり、曲線は水平な直線となる[48]。結局、曲線は、変曲点前では下に凸の曲線、変曲点後では上に凸の曲線となっており、全体としてアルファベットのSのような形を描く[48]。このため、S字型曲線やシグモイド曲線という名称でも呼ばれる[49]。間にある変曲点は個体数増加率が最大となる点で、前述の dN/dt と N のグラフの頂点に相当する[48]。変曲点における個体数は前述のとおり N = K/2 で、このときの時間は t = ln (K/N0 - 1)/r である[42]。ここで ln は自然対数である。最終的に t → ∞ で漸近する水平な直線は N = K の直線であり、時間が経過すると最終的には、個体数は環境収容力の値に収束するということである[48]。
初期個体数が N0 = K/2 の場合は、曲線は最初から変曲点から始まる。K/2 < N0 < K のときは最初から変曲点を過ぎた曲線になる[48]。 初期個体数が環境収容力に一致している場合、N0 = K のときは、その値のまま一定となる[46]。N0 = 0 のときも同様に、N = 0 のままである[46]。
次に、初期個体数が環境収容力を上回っているとき、すなわち N0 > K の場合は、この場合の曲線はS字型ではなく、全体として下に凸の曲線となる[48]。N は N0 から単調に減少しつづけ、この場合も、時間経過に従って K に収束していく[55]。
以上をまとめると、N0 > 0 であれば（個体が存在してさえいれば）、どんな初期個体数であっても、個体数は最終的に常に環境収容力の値に収束していくということである[48]。あるいは、N0 = 0 であれば（個体が存在してなければ）、個体数は 0 のままということである[48]。
最後に、生物個体数のモデルとしては無意味であるが N0 < 0 の場合も見てみると、この場合 N は時間発展に従って減少し続け、有限時間内で −∞ へ発散する曲線を描く[56]。
平衡状態の安定性
上記で、N = 0 および N = K のときはいくら時間が経過しても個体数 N は増加も減少もしないことから、これらの状態を平衡状態や定常状態と呼ぶことを説明した。平衡状態では、N = 0 または N = K という一点に留まり続ける。数学の力学系分野では、このような点を不動点や平衡点と呼ぶ[61]。平衡状態には安定な平衡状態と不安定な平衡状態がある[62]。安定な平衡状態とは、 その平衡状態の点から少しずれたとしても、時間が経過すれば平衡状態へ戻り、収束することを意味している[62]。また、不安定な平衡状態とは、平衡状態の点から少しずれたとき、時間経過すると平衡状態とのズレはどんどん大きくなっていき、平衡状態に戻らないことを意味している[63]。ロジスティック方程式の場合は、N = K 時の平衡状態が安定、N = 0 時の平衡状態が不安定となっている[64]。すなわち、初期個体数 N0 が K または 0 であれば、時間経過によらず常に同じ値を取り続けることは同じだが、N0 が平衡状態から少しずれたときの挙動は正反対となる[65]。
ロジスティック曲線とその傾きのベクトル場の様子 この安定・不安定の様子は、ロジスティック曲線の傾きをベクトル場として表すことで読み取ることができる[57]。時間経過に従って、全ての解は、これらのベクトルの矢印に沿って動いていく[66]。初期個体数が N0 > 0 であれば、t → ∞ で N は K に収束し、N0 < 0 であれば、t → ∞ で N は −∞ に発散することが分かる[57]。
N と dN/dt の関係曲線。N 軸と曲線の交点が平衡状態の点で、右が安定な点、左が不安定な点である。 あるいは、上記で説明した個体数 N と増加率 dN/dt の関係曲線からも、安定か不安定かの判別が可能である[67]。N = K の点の右側に点があるとき、dN/dt の値は負なので、N は減少していき、K に近づくことになる。N = K の点の左側に点があるときは、dN/dt は正なので、N は増加していき、同じく K に近づくことになる[47]。N = 0 の点についても、左右にずれたときの dN/dt の値の正負から、0 の点から離れていくことが理解できる[47]。
あるいは、安定性理論における線形安定性解析の考えにもとづいて、より一般的に安定性を判別することもできる。dN/dt = f(N) 、その N による微分を d(f(N))/dN = f ′(N)、平衡状態の点を Ne と置くとする。このとき、f ′(Ne) < 0 ならば Ne は安定な平衡点で、f ′(Ne) > 0 ならば Ne は不安定な平衡点であると判別できる[68]。ロジスティック方程式の場合は、
ソ連・ロシアの生物学者ゲオルギー・ガウゼによる2種の酵母（Saccharomyces cerevisiae, Schizosaccharomyces kefir）の個体群サイズ成長実験データ[78]と、それらのデータに対してフィッティングさせたロジスティック曲線を示した図。縦軸は菌全体の体積、横軸は時間を示している[78]。このガウゼの実験はロジスティック曲線がよく当てはまった個体群成長実験としてよく知られる[79]。 いくつかの微生物や小型の昆虫の飼育実験で、ロジスティック曲線がよく当てはまる個体数増加や個体密度増加実験のデータが得られている。例として以下のようなものがある。特に、ゾウリムシや酵母菌は条件さえ整えればロジスティック曲線に沿った増加をほとんどの場合で示し、高校レベルの教科書にも載る定番でもある[80]。
- キイロショウジョウバエ[81]
- ゾウリムシ[82][83]
- 大腸菌[84]
- タマミジンコ[85]
- 出芽酵母、分裂酵母（Saccharomyces cerevisiae, Schizosaccharomyces kefir）[79][86]
一方、ロジスティック曲線に当てはまるデータは得られなかったものとしては、次のような生物の実験がある。これらの実験では、時間経過後も個体数は一定に収束せず、周期的変動が繰り返されたり、大きなゆらぎが続く個体群変動となった[87]。
-
-
-
マメゾウムシ（インゲンマメゾウムシ）
パールのキイロショウジョウバエ飼育実験
ロジスティック曲線を普及させたことで知られるレイモンド・パールは、ローウェル・リードと共にキイロショウジョウバエの飼育実験を行い、この曲線を実証した。ロジスティック曲線が上手く適合する実験の具体的様子の例として、内田俊郎の著作をもとにしてパールらの実験を簡単に説明する[90]。
- パールが用意した環境は小さな牛乳瓶で、供給する餌にはバナナを磨り潰して寒天で固めイーストを少し振りかけたものを使用した[91]。牛乳瓶の中にハエと餌を入れ、温度などの環境条件を一定にし、一定時間間隔でハエの個体数を調べた[92]。
- 実験としては3種類の実験が行われた。
- 1つ目では、餌を始めに入れた後に餌を補給しなかった[93]。このため、個体数が増加して一定となった後、急激に減少してほぼ全滅状態となった[93]。
- 2つ目では、一定時間間隔で餌の継ぎ足しを行い、一定状態が保たれる結果が得られた[81]。
- 3つ目では、一定時間間隔で新しい餌の入った瓶へハエを移し替え、食糧条件だけでなく、その他の環境条件も一定に保った[81]。この結果でも一定状態が保たれ、ロジスティック曲線が当てはまるデータが得られた[81]。
パールの元へ留学していた寺尾新もこのハエの個体群成長研究を行った。それによれば、ロジスティック成長の特徴である個体群密度上昇にともなう個体数増加率の低下は、死亡率の上昇よりも出生率の低下によって起こっていた[94]。
野外生物
ハキリアリ 
キタオットセイ 野外環境では、前提条件となるような環境が保持されることはほぼ無いため、ある個体群がロジスティック曲線が当てはまるような増加の仕方を示すことは少ない[95]。自然界では環境条件は常に変化し、個体群変動のパターンも様々となる[96]。
ロジスティック曲線によく当てはまる個体数増加が確認できた例として、パナマ熱帯雨林でのハキリアリの1つの巣における個体数増加結果がある[95]。理由としては、天敵がいないこと、雨量・温度の気象条件が安定していることなどにより、ロジスティックモデルの前提条件に近い環境であったことによるものと考えられている[95]。他の野外生物でロジスティック曲線に合致した例としては、アメリカ・アラスカ州のセントポール島におけるキタオットセイ（Callorhinus ursinus）の個体数増加の結果がある[97]。植物の場合では、アイルランドのスルツェイ島で観測されたコケの成長の例がある。新規に露出した岩表面上へのコケの定着・広がり方が、ロジスティック曲線に当てはまる観測データを見せた[20]。
人口成長
世界人口のグラフ 式を発案したフェルフルストは人口の成長の様子を表すためにロジスティック方程式を発案した[98]。式を普及させたパールとリードも、ロジスティック方程式を使った最初の個体群成長研究は人口成長に対するものであった[99]。彼らは共に、当時までの人口統計をもとにしてアメリカ合衆国の将来の人口を予測したが、どちらの予測も実際の人口成長を言い当てることはできなかった[100]。パールとリードの結果では、1700年から1940年までの値は曲線によく合致していた。彼らが当てはめたロジスティック曲線では人口はその後飽和に向かうはずだったが、実際にはそれを裏切り、1940年以後もアメリカの人口は急増状態が続いた[101]。
さらにパールは、当時の推定世界人口をもとに世界人口の上限値（環境収容力 K）の推定を行った。1924年と1936年、パールはそれぞれ別の研究者とともに推定を行い、その値を発表した。それらの上限推定値は、前者では20億人、後者では26億人という値で、どちらも実際とはかけ離れたものとなった[102]。
生物学的・人口学的位置付け
ロジスティック方程式は、非常に簡単な生物学的意味からモデルを導くことができる[77]。r と K の2つのパラメータに種の特性に関わる議論を集約して、とても簡明なモデルを構成している[103]。また、式の特徴である個体数密度の上昇が増加率を抑えるロジスティック効果は、個体群生態学における基本原理ともいわれる[42]。個体数が少ない内は指数関数的に増殖し、個体数が増えてくると増加が止むという現象自体は、正確に前提条件に当てはまらないような個体群成長であっても、広く認められる現象であり、この一般的傾向をロジスティック方程式は上手く表しているとも評される[104]。
ただし、一見してロジスティック曲線のような個体群成長を示すデータであっても、そのデータに上手く曲線あてはめできる数理モデルは数多く存在する[105]。ロジスティック方程式のみが唯一当てはまるということはまずない[77]。この式が個体群成長の「普遍則」のように受け止められるのは誤解であると、数理生物学者のジェイムズ・D・マレーや応用数学者のスティーヴン・ストロガッツは指摘している[106]。
人口予測に関しても、人口学者のジョエル・E・コーエンは「ロジスティック曲線は短期的な予測に関しては、他の連続でなめらかな曲線と比べて特に劣っていることもないが、長期的な予測に関しても格別に秀でているわけでもない」と評している[107]。式を普及させたレイモンド・パールは、ある期間の人口成長にロジステック曲線が適用できる条件として、人口成長に影響を与える新しい要素がその期間中に現れないことを挙げている。しかし、このような前提条件を人口という複雑な現象に課すのは困難である点を経済学者のアルバート・B・ウルフや人口学者のジョージ・ハンドリー・ニブスなどから批判されている[108]。2010年代現在、将来人口推計にはコーホート要因法の使用が主流となっている[109]。ロジスティック曲線のような関数を過去の人口データに重ねて将来の人口を予測するという単純な方法は、現在ではほとんど行われていない[110]。
以上のように、ロジスティック方程式が個体群成長の「普遍則」というわけではないが、個体群成長モデルにおける基礎的なアイデアを有しており、より複雑な現象に対応する様々なモデルへ拡張されたり、その考え方が取り入れられたりする[5][111]。個体群成長のモデルの中で「出発点」として位置づけされる[112]。
名称の由来
フェルフルストは、1845年の論文で、"Nous donnerons le nom de logistique à la courbe"（参考訳：私たちはその曲線にロジスティック "logistique" という名前を与える）と述べ[113]、ロジスティック方程式の解による曲線を logistique と名付けた[42]。これが、式が"ロジスティック"方程式、その解曲線が"ロジスティック"曲線と呼ばれる由来である[114]。しかし、フェルフルストは logistique という語を使った理由を説明しなかったので、それ以上の由来は分かっていない[115][39][2]。
logistique と名付けられた理由のいくつかの推測は存在する。ベルギー王国陸軍士官学校の数学教授のHugo Pastijnは、実際の理由は不明と断った上で、
- 陸軍大学に勤めていたフェルフルストも馴染みが有ったであろう「兵站」の意味と関連付けて logistique と名付けたのではないか
- フェルフルストのモデルでも扱われる人口のための限られた資源と関連させて、「住居」を意味するフランス語の logis から名付けたのではないか
と、ありえそうな理由を2点ほど推測している[115]。また、19世紀当時のフランスでは、logistique には「計算に巧みな」「計算の技巧」といった意味での用例があった点も指摘されている[116]。
モデルの拡張・応用
既に述べたとおり、ロジスティック方程式を基本にすえて、様々なモデルが提案されてきた。以下では、そのようなモデルの拡張・応用の例を説明する。
捕獲の影響
人間が資源として利用するための捕獲や収穫は、その種を絶滅させる可能性もあるほどの大きな影響を持っている[45]。漁業分野では、水産資源を獲りつくさないように資源・漁業管理する必要性が認識されている[117]。持続可能な漁業のためには、人間による漁獲量が漁獲対象の自然増加量を上回らないようにする必要がある[118]。漁獲量と自然増加量が一致するとき資源は一定に保たれるので、このときの漁獲量を持続生産量と呼ぶ[118]。さらに、可能な持続生産量の中でも最大のものを最大持続生産量（MSY）と呼び、漁獲基準の一つの目安とされている[119]。
この最大持続生産量の値をロジスティック方程式を利用して定量化するモデルを、ジェーファーのプロダクションモデルなどと呼ぶ[118]。漁獲量（漁獲速度）を Y とすれば、次のように、ロジスティック方程式で表される個体数増加率（自然増加率）から Y を差し引いた値が実際の増加率となる[45]。
ロトカ・ヴォルテラの競争式では係数の値がある範囲内のときのみ2種が共存し（図の3）、それ以外ではどちらかが絶滅する（図の1, 2, 4）[124]。 「競争 (生物)」および「ロトカ・ヴォルテラの競争方程式」も参照ロジスティック方程式は環境内に1種のみが存在するときの（あるいは1種とみなせるときの）モデルだが、実際の環境では複数以上の種が生息している[125]。複数の種が存在するとき、それぞれの種の間には競争や相利共生、捕食-被食などの関係が存在して、それぞれの個体数が互いの個体数増加率に影響を与える[126]。その中でも特に、環境内に競争関係にある2種が存在する場合にロジスティック方程式を拡張させたものとして、以下のロトカ・ヴォルテラの競争方程式が知られる[127]。
ヒツジキンバエ。ロバート・メイが時間遅れを持つロジスティック方程式でこのハエの個体数変動の解析を行った。 ロジスティック方程式では、ある時刻の個体数 N(t) が同時刻の個体数増加率 dN(t)/dt に瞬間的に影響を与えるというモデルになっている[29]。しかし、妊娠期間や性成熟までの期間などが存在するため瞬間的に影響が出るというのは非現実的でもある[131]。よって、モデルの中に影響の時間遅れを含ませることが考えられる。遅延時間を T とすると、ロジスティック方程式に時間遅れの効果を取り込んだモデルとして
離散時間型モデルの場合の個体数の変化の様子。いずれも K と N0 の値は同じだが、青が r = 0.6、赤が r = 2.1、緑が r = 3 のときを示している。 「ロジスティック写像」も参照ロジスティック方程式では、時間 t を連続な実数として個体数変動をモデル化した。しかし、世代の交代が同期的に起こり、世代の重なりがないようなときには、時間を飛び飛びの時間間隔（離散時間）でモデル化する方が妥当である[138]。ロジスティック方程式型の離散時間モデルにはいくつかの種類があるが、一例として次のような差分方程式がある[139]。
ピエール＝フランソワ・フェルフルスト (Pierre-François Verhulst) ベルギー・ブリュッセルの陸軍大学の数学者であったピエール＝フランソワ・フェルフルストによって、ロジスティック方程式は発表された[159]。18世紀になると、トマス・ロバート・マルサスが出版した『人口論』に関心が高まっていた[160]。マルサスモデルの説明で述べたように、マルサスは人口が指数関数的に成長していくモデルを発表し、その帰結として社会が飢饉の発生など破滅的状況を迎えることを予測した[161]。このセンセーショナルな予測は衝撃を与え、当時およびマルサス死後も長く続く論争を引き起こした[162]。「近代統計学の父」と呼ばれるアドルフ・ケトレーも、マルサスのモデルに関心を持ち、人口増減モデルについて論じた[163]。ケトレーは流体の抵抗をヒントにして、人口増加率への抵抗は人口増加率自体の二乗に比例すると考えた[164]。
ケトレーから教えを受けたこともあり、友人でもあったフェルフルストは、ケトレー自身からケトレーのモデルに関する研究を勧められた[39]。ケトレーの考えをもとにして、人口が人口自体によって増加する一方で、人口増加を抑制する何らかの機構が働く数学的なモデルを思案した[42]。1838年、フェルフルストは、"Notice sur la loi que la population poursuit dans son accroissement"という題で研究成果を発表し、この論文の中でロジスティック方程式が提案された[165]。この論文の中でフェルフルストが実際に提案した式は、
レイモンド・パール (Raymond Pearl) フェルフルスト発表の後、生物の個体群成長に関する実験などで、同じ式が独自にあちこちで使われ始めたが、フェルフルストの名が挙がることはなかった[169]。1908年には生理学者のブレイルスフォード・ロバートソンが、動物、植物、人間といった生物の個体成長を同形式の曲線で記述した[99]。ロバートソンはフェルフルストの発表を知らなかったが、同じ曲線を用い、さらに偶然にもケトレーのデータを使用している[169]。この時点で同じ曲線が化学物質における自己触媒反応の過程を表すのに使われていたので、ロバートソンは曲線のことを自己触媒的曲線（autocatalytic curve）と呼んでいた[55]。
1920年、ジョンズ・ホプキンス大学のレイモンド・パールとローウェル・リードが、ロジスティック方程式と同形式のモデルを用いてアメリカ合衆国の人口増加について論じた[2]。この研究も、フェルフルストにより先に発表されていたことを知らずに行われた[170]。翌年の1921年には、これがすでに80年近く前にフェルフルストによって発見されたことをパールらも認めた[169]。これによってパールらもロジスティック曲線という名称を使うようになり、やっとフェルフルストの名がこの式に結びつくことになる[169]。これ以降、生物学ではロジスティック曲線という名称が定着した[55]。
パールはショウジョウバエの個体群成長の実験を行い、この式を実証した[171]。1924年と1925年にも、アメリカ、スウェーデン、フランスなどの様々な国勢調査の人口統計にロジスティック曲線のあてはめを行い、よく一致することを示した[172]。このような積み重ねた証拠をもとに、パールは個体群全般がロジステック曲線に沿って成長することを強く確信し、ロジステック曲線が「法則」であると主張した[173]。当時、パールとリードはこの式の価値を「控え目にいっても、それはケプラーの惑星の楕円運動法則に匹敵するものであるといってもよいように思われる」と自身らで評価している[174]。ロジスティック曲線は、経験的なものというよりも、個体群成長全般において普遍性を持つ法則であり、成長の長期的傾向の予測も可能にするとパールは考えていた[175]。パールは、この式が個体群成長における普遍則であるという持論を広め、ロジスティック方程式の普及に大きく貢献することになる[176]。このため、ロジスティック曲線にはパールの名が題されることもある[102]。
一方で、パールの自説の展開には多くの批判も呼び、1940年にパールが死去するまで論争が続いた[177]。経済学者のA. B. ウルフ、人口学者のジョージ・ハンドリー・ニブス、統計学者のエドウィン・ウィルソンなどがパールの主張に批判を加えている[178]。動物学者のジェームズ・グレイ、ランスロット・ホグベン、遺伝学者のシューアル・ライトからは、他のS字型曲線を使っても個体群成長のデータに当てはめができるので、ロジスティック曲線を使う必然性が欠けていることについて指摘を受けている[179]。
ロジスティック曲線の有効性を支持し、その普及を担った人たちもいる。イギリスの統計学者ウドニー・ユールは、パールの理論を1924年のイギリスの学会で発表している。ユールはロジスティック曲線は長期予測には適用できないと考えており、その点を強調したが、基本的にはパールの研究を支持していた[180]。アルフレッド・ロトカも、ロジスティック曲線の有効性を理解し、ロジスティック方程式について1925年の自書の中で一章を与えて説明した[181]。ただしロトカは、ロジスティック方程式は実現象の近似の一種であるという考えを保っていた[182]。ロシアのゲオルギー・ガウゼも、近似の一種と受け止めながらも、ロジスティック方程式が同種の集団の中での生存競争を定量的に表すことができると述べている[183]。1934年、ガウゼは微生物の実験によってロジスティック方程式の検証を行い、この検証はロジスティックモデルを個体群動態論における古典的理論の一つとして確固たるものとした[184]。モデルの限界には注意が払われながらもロジスティック方程式の受容は広まっていき、1940年代後半には個体群解析における一般的な道具として確立した[185]。
ロジスティック方程式からの発展
その後は、より現実的な個体群変動を表すことができるように、ロジスティック方程式を修正したモデルが提案されてきた。1948年には、ジョージ・イヴリン・ハッチンソンが時間遅れの影響をロジスティック方程式に導入した研究を行った[185]。ロジスティック方程式の前提条件を満たすような環境であっても、個体数が一定に収束せず、多くなったり少なくなったりをいつまでも繰り返すような生物実験の結果も得られた[89]。京都大学の内田俊郎と藤井宏一がヨツモンマメゾウムシの培養実験でそのような結果を得たことを1953年に発表している[186]。内田らは、ロジスティック方程式をもとにした差分方程式でこの結果を分析し、個体数の振動を再現した[187]。
ロジスティック方程式における r は個体群密度がとても低いときの増加率で表しており、密度が低いときにどれだけ素早く繁殖できるかを意味している[188]。また、K はその環境下で生存できる個体数あるいは個体群密度の上限を示す[189]。1967年、ロバート・マッカーサーとエドワード・オズボーン・ウィルソンは、この r と K に着目して、島における生物個体群の定着と絶滅に関する理論を発案した[190]。彼らの理論によれば、ある生物の島への定着が成功するには大きな r を持つことが重要であり、絶滅の回避には大きな K を持つことが重要であるとし、それぞれの方向へ淘汰されることを r淘汰、K淘汰と呼んだ[191]。この説はr-K戦略説と呼ばれ、生物の生活史の進化に種内競争の観点から説明を与えた[192]。
物理学から数理生態学へ転向してきたロバート・メイも個体群変動の問題に取り組んだ[193]。メイはロジスティック方程式の離散化を行い、その式の解は、通常のロジスティック方程式の解とは全く異なる、現在ではカオスと呼ばれる非常に複雑な振る舞いを起こすことを示した[194]。この結果は1974年と1975年に発表され、大きな反響を得ると共に、その後のカオス理論の隆盛に大きく寄与することになる[195][196]。
脚注
注釈
ウィキメディア・コモンズには、ロジスティック関数に関連するカテゴリがあります。
- フェルフルストの原論文
- Verhulst, Pierre-François (1838). “Notice sur la loi que la population suit dans son accroissement”. Correspondance mathématique et physique 10: 113-121. - Google ブックス
- Verhulst, Pierre-François (1845). “Recherches mathématiques sur la loi d'accroissement de la population”. Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Bruxelles 18: 1–42. - Göttinger Digitalisierungszentrum
- 『ロジスティック方程式』 - コトバンク
- 『ロジスティック曲線』 - コトバンク
- Weisstein, Eric W. "Logistic Equation". mathworld.wolfram.com (英語).
