ロジスティック方程式

ロジスティック方程式は、生物の個体数の変化の様子を表す数理モデルの一種である。ある単一種の生物が一定環境内で増殖するようなときに、その生物の個体数の変動を予測できる。人間の場合でいえば、人口の変動を表すモデルである。

1838年にベルギーの数学者ピエール＝フランソワ・フェルフルストによって、ロジスティック方程式は最初に発案された。フェルフルストは、1798年に発表されて大きな反響を呼んだトマス・ロバート・マルサスの『人口論』の不自然な点を解消するために、このモデルを考案した。マルサスは『人口論』で、人口は原理的に指数関数的に増加することを指摘したスチュアート。しかし、実際には環境や資源は限られているため、人口の増加にはいずれブレーキがかかると考えるのが自然である。人口が増えるに連れて人口増加率は低減し、人口はどこかで飽和すると考えられる。ロジスティック方程式はこの点を取り入れて、生物の個体数増殖をモデル化したものである。フェルフルスト以後には、アメリカの生物学者レイモンド・パールが式を普及させた。

具体的には、ロジスティック方程式は

という微分方程式で表される。//N// は個体数、//t// は時間、//dN/////dt// が個体数の増加率を意味する。//r// は内的自然増加率、//K// は環境収容力と呼ばれる定数である。個体数が増えて環境収容力に近づくほど、個体数増加率が減っていくというモデルになっている。

式の解はS字型の曲線を描き、個体数は最終的には環境収容力の値に収束する。この曲線や解の関数はロジスティック曲線やロジスティック関数として知られる。方程式の名称は、ロジスティック式やロジスティックモデル、ロジスティック微分方程式と表記される場合もある。発案者の名からVerhulst方程式、発案者と普及者の名からVerhulst-Pearl方程式とも呼ばれるティーメ。

ロジスティック方程式は、個体群生態学あるいは個体群動態論における数理モデルとしては入門的なものとして位置づけられ、より複雑な現象に対応する基礎を与える。数学分野としては、微分方程式論や力学系理論の初等的な話題としても取り上げられるHirsch et al.。

1. 生物の個体数のモデル
-------------

生物の個体数の変動については古くから興味を持たれ、研究が行われてきたマレー。フィボナッチ数の発見に繋がったレオナルド・フィボナッチのウサギの個体数の問題が、おそらく最も古い個体数の数理モデルといわれるスチュアート。

生物の個体数の増え方に関する研究は、個体群生態学の分野に属する日本数理生物学会。ここで、個体群とは簡単には、ある領域に生息している単一の種の個体の集まりのことを指す日本数理生物学会。

この個体群の大きさの指標としては、個体群内の総個体数が使用される大串。個体数の代わりに、領域の単位面積当たりの個体数である個体群密度や単位面積当たりの生物の総重量である生物量が、個体群サイズとして適切な指標となる場合もある寺本。人間でいえば、これらの指標は人口や人口密度に相当する瀬野。

1.1. マルサスモデル
------------

多くの生物では、親は多くの子孫を作るので、それがそのまま生き残ると仮定すれば、あっという間に莫大な個体数となる。ねずみ算など、数学的小話の種である巌佐。まずはこのような単純なものが、生物個体数の増加モデルとして考えられる。

ある個体群において、時刻 //t// に個体数が //N// 体が存在しているとする。実際の生物個体数は不連続な値をとるものであるが、数学的扱いを簡便にするために、個体数は連続な値をとるものとすることがしばしば行われる山口。実際の生物でいえば、個体数が多かったり各個体の世代が重なったりしていれば、このような近似も妥当性を帯びてくるブラウン。個体数を連続な値とすれば、個体数の増加率は //N// の時間微分 //dN/////dt// で表すことができるハーバーマン。

さらに話を単純化するために、個体は環境を出入りしないという状況を想定する渡辺。この場合、個体の出生と死亡という2つの要因のみによって個体数は増減する渡辺。個体群の出生率が死亡率を上回っていれば、個体数は増え続けるということになるマレー。さらに簡略化するために出生率と死亡率を常に一定であるとするマレー。個体数当たりの出生率を //b//、個体数当たりの死亡率を //d// とすれば、個体数の増加率は差し引きした //b// − //d// に個体数 //N// を掛け合わせた値となる寺本。よって個体数増加率 //dN/////dt// は

という微分方程式で表される巌佐。ここで //m// は比例定数であり、//m// = //b// − //d// である。

このような式で表される個体数増加は //t// の指数関数となり、人間でいえば、あっという間に人口爆発を引き起こすことになる。このような個体群成長のモデルは、生物個体の増加が幾何級数的であることを最初に指摘したトマス・ロバート・マルサスに因んでマルサスモデルと呼ばれる人口研究会。比例定数 //m// もマルサスの名からマルサス係数と呼ばれ、単位は一個体当たりの増加率となる。

しかし、このモデルは現実と違いすぎるマレー。現実の生物は、限られた環境下で生息しており、個体数が多くなると、各個体にとって必要な資源が得にくくなる巌佐。そこに生息できる個体数には上限があると見るのが自然である日本生態学会。つまり、個体数が多くなると、その増加にブレーキがかかるものと想像される日本生態学会。このような一種内での資源の取り合いは種内競争と呼ばれ、生物における競争関係の一種である巌佐。

2. ロジスティック方程式
-------------

上記のようにマルサスモデルは非現実的な面を持つ。個体数が多くなると増加率が抑えられることを表現するために、個体数 //N// が増加するにつれて増加率 //m// が減少するモデルが考えられる寺本。また、個体数がある上限を超えたら増加率は負となり、個体数は減少に向かうと考えられるStrogatz。これらの点を簡単に表せば、比例定数 //m// を

と置ける巌佐。すなわち、//m// の値は個体数がゼロに限りなく近いときに最大値で、その後は //N// の値の増加に比例して //m// の値は減少するというモデルである。これをマルサスモデルに代入して、次の微分方程式を得ることができる。

この微分方程式をロジスティック方程式と呼ぶ。個体群成長モデルの一種としてロジスティックモデルとも呼ばれる。この微分方程式は、数学的には //n// = 2 のベルヌーイの微分方程式にも該当する寺本。

ロジスティック方程式の //K// は環境収容力と呼ばれ、その環境が維持できる個体数を意味する巌佐。//r// の単位は上記のマルサス係数と同じく一個体当たりの増加率だが瀬野、特に内的自然増加率と呼ばれ、その生物が実現する可能性のある最大増加率を示している瀬野。通常のロジスティック方程式では、//K// と //r// は時間に関わらず一定とみなし、正の定数と考えるコーエン。

2.1. ロジスティック効果
--------------

マルサスモデルからロジスティック方程式へ拡張したときに行ったことは、個体群生態学における密度効果を取り入れたことに相当する渡辺。上記では //N// を個体数として説明したが、ロジスティック方程式では有限な環境を前提にしているので、//N// は単位面積当たりの個体数である個体群密度でもある瀬野。個体群密度がその個体群自身の変動に影響を与えることは、密度効果という名称で呼ばれる日本数理生物学会。特にロジスティック方程式では、個体群密度が高くなると増加率に負の効果を与える種類の密度効果となっており、これをロジスティック効果と呼ぶ人口研究会。

ロジスティック方程式では個体群密度増加に比例して増加率が一方的に低下することを想定したが、密度増加によって増加率が上昇する場合も考えられるティーメ。例えば、ある程度は密度が高くないと交尾の相手が見つけるのが困難となって、結果として増加率が低下する場合などである日本数理生物学会。よって、個体密度が低い内は個体群密度増加によって増加率が上昇する種類の密度効果も考えられ、このような種類の密度効果をアリー効果と呼ぶ大串。

2.2. 個体数と増加率の関係
---------------

ロジスティック方程式における個体数増加率 //dN/////dt// と個体数 //N// の関係に着目すれば、この関係は初等教育でも習う二次関数そのものとなっており、//dN/////dt// と //N// のグラフは放物線を描く寺本。方程式を解析的に解かなくとも、//N// と //dN/////dt// がどのような変化を起こすのかを、以下のようにグラフから読み解くこともできる。

まず、//N// = 0 と //N// = //K// のとき、//dN/////dt// = 0 となる。すなわち、いくら時間が経過しても個体数は増加も減少もしない状態となる。このような状態を定常状態や平衡状態と呼ぶマレー。//N// の値が 0 < //N// < //K// の範囲にあるとき、 //dN/////dt// の値は様々だが、値が正なのか負なのかで言えば、正の値であることがわかるHirsch et al.。//N// の値が //K// < //N// となると、//dN/////dt// は同じように負の値であるHirsch et al.。言い換えれば、個体数が環境収容力内では常に個体数は増加するが、環境収容力を超えると個体数は減少へ転ずる、ということであるHirsch et al.。

個体数増加率 //dN/////dt// の変化をさらに細かく見てみる。//N// が 0 から //K///2 まで増えると、その間 //dN/////dt// の値も増加し続ける。//N// = //K///2 は放物線の頂点であり、ここで //dN/////dt// は極大値を迎えるStrogatz。極大値は、//N// = //K///2 を式に代入して //dN/////dt// = //rK///4 である寺本。//N// = //K///2 を超えると //dN/////dt// は減少し始め、//N// = //K// で 0 となるStrogatz。このような数値の変化から読み取れることの一つは、個体数が環境収容力のちょうど半分になったときに個体増加率は最大となる点であるStrogatz。したがって、もし個体数の変化がロジスティック方程式に従うとしたら、増加率が最大になるときの個体数に注目することで環境収容力、すなわち最大個体数を予測できることになるスチュアート。

3. 式の解
------

3.1. ロジスティック曲線
--------------

ロジスティック方程式は非線形の微分方程式だが、標準的な微分方程式の解法である変数分離法を利用して解くことができる山口。時間 //t// = 0 における初期個体数を //N//0 とすると、//t// の関数として以下の解が得られるスチュアート一例として以下のように解くことができる。//N// の値の範囲を 0 < //N// < //K// に限定して解く方法とSfn。

ここで //e// はネイピア数である。分母・分子を //N//0//ert// で割り、次のような形でも示される。

この解の関数をロジスティック関数、この解によって描かれる曲線をロジスティック曲線と呼ぶ人口研究会。この曲線に従う個体群成長は、ロジスティック成長やロジスティック増殖とも呼ばれる寺本。関数は //t// → ∞ の極限で //N// → //K// となり、マルサスモデルと異なり発散しないことが確認できるHirsch et al.。

3.2. 曲線の形状
----------

横軸を //t//、縦軸を //N// とした平面上にロジスティック関数のグラフを描くと、曲線が描かれる。この曲線は前述のとおりにロジスティック曲線と呼ばれる。初期個体数が3つの範囲 //N//0 < 0, 0 < //N//0 < //K//, //K// < //N//0 のどれに該当するかによって、曲線の形状は大きく異なってくるHirsch et al.。ただし、//N//0 < 0 の範囲は負の個体数というものを意味するので、生物のモデルとしてはあまり意味がないHirsch et al.。時間 //t// = 0 から //t// → ∞ の極限までのロジスティック曲線の様相は、それぞれの //N//0 の値ごとに、以下のようになっている。

まず //N//0 が環境収容力の半分以下の場合、初期状態の点  から始まる曲線は、ゆっくりと右肩上がりに登っていく山口。//t// が増加するにつれて、曲線の傾きは増加していき、曲線は加速度的に立ち上がっていくStrogatz。しかし、ある時点で曲線は変曲点を迎え、傾きの増加は止むブラウン。その後は、傾きは減少しだし、曲線は横倒しになっていく山口。そして最終的には、傾きは0になり、曲線は水平な直線となるStrogatz。結局、曲線は、変曲点前では下に凸の曲線、変曲点後では上に凸の曲線となっており、全体としてアルファベットのSのような形を描くStrogatz。このため、S字型曲線やシグモイド曲線という名称でも呼ばれるスチュアート。間にある変曲点は個体数増加率が最大となる点で、前述の //dN/////dt// と //N// のグラフの頂点に相当するStrogatz。変曲点における個体数は前述のとおり //N// = //K///2 で、このときの時間は //t// = ln ///r// である人口研究会。ここで ln は自然対数である。最終的に //t// → ∞ で漸近する水平な直線は //N// = //K// の直線であり、時間が経過すると最終的には、個体数は環境収容力の値に収束するということであるStrogatz。

初期個体数が //N//0 = //K///2 の場合は、曲線は最初から変曲点から始まる。//K///2 < //N//0 < //K// のときは最初から変曲点を過ぎた曲線になるStrogatz。 初期個体数が環境収容力に一致している場合、//N//0 = //K// のときは、その値のまま一定となるマレー。//N//0 = 0 のときも同様に、//N// = 0 のままであるマレー。

次に、初期個体数が環境収容力を上回っているとき、すなわち //N//0 > //K// の場合は、この場合の曲線はS字型ではなく、全体として下に凸の曲線となるStrogatz。//N// は //N//0 から単調に減少しつづけ、この場合も、時間経過に従って //K// に収束していく寺本。

以上をまとめると、//N//0 > 0 であれば、どんな初期個体数であっても、個体数は最終的に常に環境収容力の値に収束していくということであるStrogatz。あるいは、//N//0 = 0 であれば、個体数は 0 のままということであるStrogatz。

最後に、生物個体数のモデルとしては無意味であるが //N//0 < 0 の場合も見てみると、この場合 //N// は時間発展に従って減少し続け、有限時間内で −∞ へ発散する曲線を描くHirsch et al.。

3.3. 平衡状態の安定性
-------------

上記で、//N// = 0 および //N// = //K// のときはいくら時間が経過しても個体数 //N// は増加も減少もしないことから、これらの状態を平衡状態や定常状態と呼ぶことを説明した。平衡状態では、//N// = 0 または //N// = //K// という一点に留まり続ける。数学の力学系分野では、このような点を不動点や平衡点と呼ぶStrogatz。平衡状態には安定な平衡状態と不安定な平衡状態がある巌佐。安定な平衡状態とは、 その平衡状態の点から少しずれたとしても、時間が経過すれば平衡状態へ戻り、収束することを意味している巌佐。また、不安定な平衡状態とは、平衡状態の点から少しずれたとき、時間経過すると平衡状態とのズレはどんどん大きくなっていき、平衡状態に戻らないことを意味しているマレー。ロジスティック方程式の場合は、//N// = //K// 時の平衡状態が安定、//N// = 0 時の平衡状態が不安定となっているStrogatz。すなわち、初期個体数 //N//0 が //K// または 0 であれば、時間経過によらず常に同じ値を取り続けることは同じだが、//N//0 が平衡状態から少しずれたときの挙動は正反対となる巌佐。

この安定・不安定の様子は、ロジスティック曲線の傾きをベクトル場として表すことで読み取ることができるHirsch et al.。時間経過に従って、全ての解は、これらのベクトルの矢印に沿って動いていく。初期個体数が //N//0 > 0 であれば、//t// → ∞ で //N// は //K// に収束し、//N//0 < 0 であれば、//t// → ∞ で //N// は −∞ に発散することが分かるHirsch et al.。

あるいは、上記で説明した個体数 //N// と増加率 //dN/////dt// の関係曲線からも、安定か不安定かの判別が可能である巌佐。//N// = //K// の点の右側に点があるとき、//dN/////dt// の値は負なので、//N// は減少していき、//K// に近づくことになる。//N// = //K// の点の左側に点があるときは、//dN/////dt// は正なので、//N// は増加していき、同じく //K// に近づくことになるHirsch et al.。//N// = 0 の点についても、左右にずれたときの //dN/////dt// の値の正負から、0 の点から離れていくことが理解できるHirsch et al.。

あるいは、安定性理論における線形安定性解析の考えにもとづいて、より一般的に安定性を判別することもできる。//dN/////dt// = //f// 、その //N// による微分を //d//)///dN// = //f//   ′、平衡状態の点を //Ne// と置くとする。このとき、//f//   ′ < 0 ならば //Ne// は安定な平衡点で、//f//   ′ > 0 ならば //Ne// は不安定な平衡点であると判別できるマレー。ロジスティック方程式の場合は、

なので、

となり、//f//   ′ = −//r// < 0, //f//   ′ = //r// > 0 となることが確認できるStrogatz。

4. 生物学的前提条件
-----------

実際の生物の個体数増殖においてロジスティック方程式が成り立ち、ロジスティック曲線がその増殖データに上手く当てはまるには、次のような生物学的条件が前提として挙げられる。

環境内には単一の種か、あるいは同等とみなせる種のみが存在する瀬野。
対象の生物の各世代は連続的に重なっているマレー。すなわち、連続的に子が生まれ、親と子が共存する期間が存在する山口。
個体は一定の大きさの環境内に常に存在する。すなわち、環境から移出したり、外部から移入が無いコーエン。
環境の大きさは変わらず、一定状態が保たれるコーエン。
個体群のために、食糧や資源が一定して供給されるStrogatz。ショウジョウバエや真正細菌といった、微生物や単純な生物を一定環境で増殖させた場合は、上記の条件に近く、ロジスティック方程式によって個体数変化の正確な予測ができる巌佐。しかし、例えば鹿や鳥類などのような、一定環境のもとで増殖する設定が成立しない個体群成長には、ロジスティック方程式を適用することはできない巌佐。

環境を整えた飼育実験によって、ロジスティック曲線に当てはまる個体数増殖のデータを得ることはできるが、上記の生物学的条件を実験上で整えることはいつも簡単というわけではない山口。増殖を抑える原因となる老廃物を定期的に取り除く、といった配慮も必要となる山口。

5. 実際のデータへの適用例
--------------

5.1. 実験生物
---------

いくつかの微生物や小型の昆虫の飼育実験で、ロジスティック曲線がよく当てはまる個体数増加や個体密度増加実験のデータが得られている。例として以下のようなものがある。特に、ゾウリムシや酵母菌は条件さえ整えればロジスティック曲線に沿った増加をほとんどの場合で示し、高校レベルの教科書にも載る定番でもある渡辺。

キイロショウジョウバエ
ゾウリムシ山口
大腸菌
タマミジンコ
出芽酵母、分裂酵母スチュアート File:Drosophila melanogaster - side .jpg|キイロショウジョウバエ File:Paramecium.jpg|ゾウリムシ File:EscherichiaColi NIAID.jpg|大腸菌 File:Moina sp.jpg|タマミジンコ File:S cerevisiae under DIC microscopy.jpg|出芽酵母

一方、ロジスティック曲線に当てはまるデータは得られなかったものとしては、次のような生物の実験がある。これらの実験では、時間経過後も個体数は一定に収束せず、周期的変動が繰り返されたり、大きなゆらぎが続く個体群変動となった山口。

ミバエ
コクヌストモドキ
マメゾウムシ File:Bactrocera dorsalis.jpg|ミバエ File:Tribolium castaneum.jpg|コクヌストモドキ File:Chrząszcz na ziarnach fasoli.jpg|マメゾウムシ

5.1.1. パールのキイロショウジョウバエ飼育実験
--------------------------

ロジスティック曲線を普及させたことで知られるレイモンド・パールは、ローウェル・リードと共にキイロショウジョウバエの飼育実験を行い、この曲線を実証した。ロジスティック曲線が上手く適合する実験の具体的様子の例として、内田俊郎の著作をもとにしてパールらの実験を簡単に説明する内田。

パールが用意した環境は小さな牛乳瓶で、供給する餌にはバナナを磨り潰して寒天で固めイーストを少し振りかけたものを使用した内田。牛乳瓶の中にハエと餌を入れ、温度などの環境条件を一定にし、一定時間間隔でハエの個体数を調べた内田。
実験としては3種類の実験が行われた。
1つ目では、餌を始めに入れた後に餌を補給しなかった内田。このため、個体数が増加して一定となった後、急激に減少してほぼ全滅状態となった内田。
2つ目では、一定時間間隔で餌の継ぎ足しを行い、一定状態が保たれる結果が得られた内田。
3つ目では、一定時間間隔で新しい餌の入った瓶へハエを移し替え、食糧条件だけでなく、その他の環境条件も一定に保った内田。この結果でも一定状態が保たれ、ロジスティック曲線が当てはまるデータが得られた内田。パールの元へ留学していた寺尾新もこのハエの個体群成長研究を行った。それによれば、ロジスティック成長の特徴である個体群密度上昇にともなう個体数増加率の低下は、死亡率の上昇よりも出生率の低下によって起こっていたKingsland。

5.2. 野外生物
---------

野外環境では、前提条件となるような環境が保持されることはほぼ無いため、ある個体群がロジスティック曲線が当てはまるような増加の仕方を示すことは少ない。自然界では環境条件は常に変化し、個体群変動のパターンも様々となる大串。

ロジスティック曲線によく当てはまる個体数増加が確認できた例として、パナマ熱帯雨林でのハキリアリの1つの巣における個体数増加結果がある。理由としては、天敵がいないこと、雨量・温度の気象条件が安定していることなどにより、ロジスティックモデルの前提条件に近い環境であったことによるものと考えられている。他の野外生物でロジスティック曲線に合致した例としては、アメリカ・アラスカ州のセントポール島におけるキタオットセイの個体数増加の結果があるレーヴンほか。植物の場合では、アイルランドのスルツェイ島で観測されたコケの成長の例がある。新規に露出した岩表面上へのコケの定着・広がり方が、ロジスティック曲線に当てはまる観測データを見せた。

5.3. 人口成長
---------

式を発案したフェルフルストは人口の成長の様子を表すためにロジスティック方程式を発案したKingsland。式を普及させたパールとリードも、ロジスティック方程式を使った最初の個体群成長研究は人口成長に対するものであったKingsland。彼らは共に、当時までの人口統計をもとにしてアメリカ合衆国の将来の人口を予測したが、どちらの予測も実際の人口成長を言い当てることはできなかったコーエン。パールとリードの結果では、1700年から1940年までの値は曲線によく合致していた。彼らが当てはめたロジスティック曲線では人口はその後飽和に向かうはずだったが、実際にはそれを裏切り、1940年以後もアメリカの人口は急増状態が続いた山口。

さらにパールは、当時の推定世界人口をもとに世界人口の上限値の推定を行った。1924年と1936年、パールはそれぞれ別の研究者とともに推定を行い、その値を発表した。それらの上限推定値は、前者では20億人、後者では26億人という値で、どちらも実際とはかけ離れたものとなったコーエン。

6. 生物学的・人口学的位置付け
----------------

ロジスティック方程式は、非常に簡単な生物学的意味からモデルを導くことができる山口。//r// と //K// の2つのパラメータに種の特性に関わる議論を集約して、とても簡明なモデルを構成している寺本。また、式の特徴である個体数密度の上昇が増加率を抑えるロジスティック効果は、個体群生態学における基本原理ともいわれる人口研究会。個体数が少ない内は指数関数的に増殖し、個体数が増えてくると増加が止むという現象自体は、正確に前提条件に当てはまらないような個体群成長であっても、広く認められる現象であり、この一般的傾向をロジスティック方程式は上手く表しているとも評される巌佐。

ただし、一見してロジスティック曲線のような個体群成長を示すデータであっても、そのデータに上手く曲線あてはめできる数理モデルは数多く存在する山口。ロジスティック方程式のみが唯一当てはまるということはまずない山口。この式が個体群成長の「普遍則」のように受け止められるのは誤解であると、数理生物学者のジェイムズ・D・マレーや応用数学者のスティーヴン・ストロガッツは指摘しているマレー。

人口予測に関しても、人口学者のジョエル・E・コーエンは「ロジスティック曲線は短期的な予測に関しては、他の連続でなめらかな曲線と比べて特に劣っていることもないが、長期的な予測に関しても格別に秀でているわけでもない」と評しているコーエン。式を普及させたレイモンド・パールは、ある期間の人口成長にロジステック曲線が適用できる条件として、人口成長に影響を与える新しい要素がその期間中に現れないことを挙げている。しかし、このような前提条件を人口という複雑な現象に課すのは困難である点を経済学者のアルバート・B・ウルフや人口学者のジョージ・ハンドリー・ニブスなどから批判されているKingsland。2010年代現在、将来人口推計にはコーホート要因法の使用が主流となっている人口研究会。ロジスティック曲線のような関数を過去の人口データに重ねて将来の人口を予測するという単純な方法は、現在ではほとんど行われていない。

以上のように、ロジスティック方程式が個体群成長の「普遍則」というわけではないが、個体群成長モデルにおける基礎的なアイデアを有しており、より複雑な現象に対応する様々なモデルへ拡張されたり、その考え方が取り入れられたりするマレー。個体群成長のモデルの中で「出発点」として位置づけされる渡辺。

7. 名称の由来
--------

フェルフルストは、1845年の論文で、"Nous donnerons le nom de //logistique// à la courbe"と述べ、ロジスティック方程式の解による曲線を //logistique// と名付けた人口研究会。これが、式が"ロジスティック"方程式、その解曲線が"ロジスティック"曲線と呼ばれる由来である人口研究会。しかし、フェルフルストは //logistique// という語を使った理由を説明しなかったので、それ以上の由来は分かっていない瀬野。

//logistique// と名付けられた理由のいくつかの推測は存在する。ベルギー王国陸軍士官学校の数学教授のHugo Pastijnは、実際の理由は不明と断った上で、

陸軍大学に勤めていたフェルフルストも馴染みが有ったであろう「兵站」の意味と関連付けて //logistique// と名付けたのではないか
フェルフルストのモデルでも扱われる人口のための限られた資源と関連させて、「住居」を意味するフランス語の //logis// から名付けたのではないかと、ありえそうな理由を2点ほど推測している。また、19世紀当時のフランスでは、//logistique// には「計算に巧みな」「計算の技巧」といった意味での用例があった点も指摘されている寺本。

8. モデルの拡張・応用
------------

既に述べたとおり、ロジスティック方程式を基本にすえて、様々なモデルが提案されてきた。以下では、そのようなモデルの拡張・応用の例を説明する。

8.1. 捕獲の影響
----------

人間が資源として利用するための捕獲や収穫は、その種を絶滅させる可能性もあるほどの大きな影響を持っている寺本。漁業分野では、水産資源を獲りつくさないように資源・漁業管理する必要性が認識されている北原。持続可能な漁業のためには、人間による漁獲量が漁獲対象の自然増加量を上回らないようにする必要がある。漁獲量と自然増加量が一致するとき資源は一定に保たれるので、このときの漁獲量を持続生産量と呼ぶ。さらに、可能な持続生産量の中でも最大のものを最大持続生産量と呼び、漁獲基準の一つの目安とされている北原。

この最大持続生産量の値をロジスティック方程式を利用して定量化するモデルを、ジェーファーのプロダクションモデルなどと呼ぶ。漁獲量を //Y// とすれば、次のように、ロジスティック方程式で表される個体数増加率から //Y// を差し引いた値が実際の増加率となる寺本。

//dN/////dt// が 0 のときが資源一定状態なので、このときの漁獲量が持続生産量を示している。//Y// を含まないときの //dN/////dt// の最大値は前述のとおり //rK///4 である。これに釣り合う漁獲量が最大持続生産量であるから、このモデルでは最大持続生産量を //rK///4 と得ることができる日本数理生物学会。

漁獲量 //Y// を単純な一定値とせずに、個体数に比例するようなモデルも考えられる。例えば、出漁する漁船の数が一定とすれば、捕獲の成果は生息している個体数に比例すると考える方が適当である寺本。//q// と //E// を定数として //q// と //E// と //N// を掛け合わせたもので漁獲量を表せば、個体数増加率は

となる。//q// は漁具効率、//E// は漁獲努力量と呼ばれる。このモデルの場合は、//qE// を内的自然増加率の半分 //r///2 となるようにすれば、漁獲量を最大持続生産量にすることができる。ただし、以上のモデルは現実をかなり単純化したモデルである。環境の変化や他の生物との相互作用など、現実には様々な要因が関係しているため、多数の相互作用がある実際の生態系では成り立たない日本数理生物学会。実際の最大持続生産量の決定にはより高度な手法も使用されている。

8.2. 2種存在する場合
-------------

ロジスティック方程式は環境内に1種のみが存在するときのモデルだが、実際の環境では複数以上の種が生息している巌佐。複数の種が存在するとき、それぞれの種の間には競争や相利共生、捕食-被食などの関係が存在して、それぞれの個体数が互いの個体数増加率に影響を与えるマレー。その中でも特に、環境内に競争関係にある2種が存在する場合にロジスティック方程式を拡張させたものとして、以下のロトカ・ヴォルテラの競争方程式が知られる大串。

係数の //N//1, //r//1, //K//1 は種1の個体数、内的自然増加率、環境収容力である。同様に、//N//2, //r//2, //K//2 は種2の個体数、内的自然増加率、環境収容力である。さらに、//a//12 が種2が種1に与える影響を、//a//21 が種1が種2に与える影響を表し、競争係数と呼ばれる日本生態学会。この式はアメリカの数学者アルフレッド・ロトカとイタリアの数学者ヴィト・ヴォルテラによって独立に考案された日本生態学会。

ロトカ・ヴォルテラの競争式では、それぞれの係数の値がある範囲内のときのみ2種が共存するが、それ以外の場合にはどちらかの種が絶滅する結果に至る巌佐。この結果はゲオルギー・ガウゼの競争排除則を裏付ける一例となっているハーバーマン。

8.3. 時間遅れの考慮
------------

ロジスティック方程式では、ある時刻の個体数 //N// が同時刻の個体数増加率 //dN/////dt// に瞬間的に影響を与えるというモデルになっている日本生態学会。しかし、妊娠期間や性成熟までの期間などが存在するため瞬間的に影響が出るというのは非現実的でもあるマレー。よって、モデルの中に影響の時間遅れを含ませることが考えられる。遅延時間を //T// とすると、ロジスティック方程式に時間遅れの効果を取り込んだモデルとして

がよく用いられるマレー。この式はジョージ・イヴリン・ハッチンソンが発案したためHutchinson方程式とも呼ばれる日本数理生物学会。このモデルでは、ロジスティック方程式におけるブレーキ効果の部分に、現時点での個体数 //N// ではなく、時間 //T// だけ前の時点での個体数 //N// が入力されている。

時間遅れを持つロジスティック方程式でも //N// = 0 または //N// = //K// が平衡状態であることに変わりはない日本数理生物学会。しかし個体数が環境収容力 //K// に達しても、//T// 時間前における個体数は //K// よりも小さいか大きいので、増加率は0とならない。そのため、個体数は環境収容力を通り過ぎてしまう。環境収容力を上回った個体数が継続すると、増加率は個体数を環境収容力に収束させる方向に働く。しかし、それによって個体数が環境収容力に戻っても、再度同じ現象が起き、また環境収容力を通り過ぎる。このように平衡状態を行き過ぎたり戻り過ぎたりしながら個体数が振動する現象が、このモデルでは起こり得るハーバーマン。より詳細にいえば、解の振る舞いは //rT// の値によって変化する。//rT// が //π///2 を超えるとホップ分岐を起こし、解は平衡状態を回るリミットサイクルとなる日本数理生物学会。周期変動を実際に起こすヒツジキンバエの実験データに対して、ロバート・メイがこの式の当てはめを行って良好な結果を得ているマレー。

8.4. 離散時間モデル
------------

ロジスティック方程式では、時間 //t// を連続な実数として個体数変動をモデル化した。しかし、世代の交代が同期的に起こり、世代の重なりがないようなときには、時間を飛び飛びの時間間隔でモデル化する方が妥当である寺本。ロジスティック方程式型の離散時間モデルにはいくつかの種類があるが、一例として次のような差分方程式がある巌佐。

ここで、//n// は世代で、//n// = 1世代, 2世代, 3世代,... といったような飛び飛びの時間間隔を意味している。//Nn// は、//n// 世代における個体数 //N// を意味している巌佐。上式と数学的には等価だが、ロジスティック写像と呼ばれる、次の形式での差分方程式もよく知られている。

これらの差分方程式はロジスティック方程式と一見似ているが、解の様相は全く異なり、個体数の変動はロジスティック方程式よりも遥かに複雑な振る舞いを見せるようになるHirsch et al.。//r//が小さい内は、これらの解はロジスティック方程式と同じように安定な平衡状態に収束する巌佐。//r// が大きくなってくると、個体数は多くなったり少なくなったりを交互に繰り返すようになる。さらに //r// が大きくなると、カオスと呼ばれる非周期的で極めて複雑な振る舞いを起こすようになる巌佐。

また、京都大学の森下正明が発案した次のような差分方程式がある。

ここで、//Δt// を差分時間間隔として、//a// と //b// は

である。通常、差分化を行うと元の方程式の解と誤差が生じる山口。しかしこの方程式では誤差を全く生じさせない。得られる解は離散的だが、その解はロジスティック方程式の解と一致し、解を //N//-//t// 平面上に描けば、ロジスティック曲線上に正確にプロットされる山口。

8.5. 生物個体数以外での例
---------------

本来の導入目的であった生物の個体数の変動以外にも、ロジスティックモデルがしばし使用される。興味対象の何かの変量が時間発展とともにS字型の曲線を描くようなときに、この式がよく当てはまる巌佐。水産資源管理の例では、生物の体の大きさの成長曲線にロジスティック曲線を当てはめることがある北原。また、人間の集団の中で無形なものが広まる様子を表すのにもロジスティックモデルが使われることがある。例えば、新技術の社会・産業全体への普及ブラウン、ある集団の中での噂の拡散がある。

また、時間変化ではないが、統計学においてはロジスティック関数と同形式の累積分布関数 //f// を持つ連続確率分布が用いられている。これをロジスティック分布と呼ぶ。人工ニューラルネットワークの研究で使われるシグモイド関数の一つとしてもロジステック関数が利用されている。

9. 他形式
------

上記ではロジスティック方程式を

と表したが、これ以外の表現もある。いずれも数学的には等価だが、その導出過程における生態学的意味づけは様々である瀬野。 //k// = //r// / //K// と置いて、ロジスティック方程式は

とも表される瀬野。//k// はVerhulst-Pearl係数や種内競争係数と呼ばれる寺本。個体群密度が増加率を減少させる影響の強さを //k// が表しているといえる瀬野。

他には、変数を //N// = //N/////K// と置きなおして、すなわち個体数ではなく環境収容力に対する個体数の割合を変数として

という形式もあるHirsch et al.。

非線形のロジスティック関数を扱いやすくするために線形の対数関数に変換する、フィッシャ・プライ変換と呼ばれる次のような変換もある。

ここで //FP// = //N// とすると、ロジスティック関数のパラメータとの関係は //K// = 1, //r// = //b//, //N//0 = //ea/// である。

10. 歴史
------

10.1. フェルフルストによる発表
------------------

ベルギー・ブリュッセルの陸軍大学の数学者であったピエール＝フランソワ・フェルフルストによって、ロジスティック方程式は発表されたコーエン。18世紀になると、トマス・ロバート・マルサスが出版した『人口論』に関心が高まっていた山口。マルサスモデルの説明で述べたように、マルサスは人口が指数関数的に成長していくモデルを発表し、その帰結として社会が飢饉の発生など破滅的状況を迎えることを予測した人口研究会。このセンセーショナルな予測は衝撃を与え、当時およびマルサス死後も長く続く論争を引き起こした人口研究会。「近代統計学の父」と呼ばれるアドルフ・ケトレーも、マルサスのモデルに関心を持ち、人口増減モデルについて論じた人口研究会。ケトレーは流体の抵抗をヒントにして、人口増加率への抵抗は人口増加率自体の二乗に比例すると考えた山口。

ケトレーから教えを受けたこともあり、友人でもあったフェルフルストは、ケトレー自身からケトレーのモデルに関する研究を勧められた瀬野。ケトレーの考えをもとにして、人口が人口自体によって増加する一方で、人口増加を抑制する何らかの機構が働く数学的なモデルを思案した人口研究会。1838年、フェルフルストは、"//Notice sur la loi que la population poursuit dans son accroissement//"という題で研究成果を発表し、この論文の中でロジスティック方程式が提案された山口。この論文の中でフェルフルストが実際に提案した式は、

という形であった瀬野。//p// は人口を意味する。フェルフルストは人口自体の二乗によって人口増加率の減少効果を表現し、上記の //φ// を導入した瀬野。当時はこの式の価値を認めるものはほとんどなく、彼の死亡時の告知にも、彼の業績として取り上げられなかった山口。

10.2. 式の再発見と論争と普及
-----------------

フェルフルスト発表の後、生物の個体群成長に関する実験などで、同じ式が独自にあちこちで使われ始めたが、フェルフルストの名が挙がることはなかった山口。1908年には生理学者のブレイルスフォード・ロバートソンが、動物、植物、人間といった生物の個体成長を同形式の曲線で記述したKingsland。ロバートソンはフェルフルストの発表を知らなかったが、同じ曲線を用い、さらに偶然にもケトレーのデータを使用している山口。この時点で同じ曲線が化学物質における自己触媒反応の過程を表すのに使われていたので、ロバートソンは曲線のことを自己触媒的曲線と呼んでいた寺本。

1920年、ジョンズ・ホプキンス大学のレイモンド・パールとローウェル・リードが、ロジスティック方程式と同形式のモデルを用いてアメリカ合衆国の人口増加について論じた。この研究も、フェルフルストにより先に発表されていたことを知らずに行われたコーエン。翌年の1921年には、これがすでに80年近く前にフェルフルストによって発見されたことをパールらも認めた山口。これによってパールらもロジスティック曲線という名称を使うようになり、やっとフェルフルストの名がこの式に結びつくことになる山口。これ以降、生物学ではロジスティック曲線という名称が定着した寺本。

パールはショウジョウバエの個体群成長の実験を行い、この式を実証した内田。1924年と1925年にも、アメリカ、スウェーデン、フランスなどの様々な国勢調査の人口統計にロジスティック曲線のあてはめを行い、よく一致することを示したマレー。このような積み重ねた証拠をもとに、パールは個体群全般がロジステック曲線に沿って成長することを強く確信し、ロジステック曲線が「法則」であると主張したKingsland。当時、パールとリードはこの式の価値を「控え目にいっても、それはケプラーの惑星の楕円運動法則に匹敵するものであるといってもよいように思われる」と自身らで評価しているコーエン。ロジスティック曲線は、経験的なものというよりも、個体群成長全般において普遍性を持つ法則であり、成長の長期的傾向の予測も可能にするとパールは考えていたKingsland。パールは、この式が個体群成長における普遍則であるという持論を広め、ロジスティック方程式の普及に大きく貢献することになるマレー。このため、ロジスティック曲線にはパールの名が題されることもあるコーエン。

一方で、パールの自説の展開には多くの批判も呼び、1940年にパールが死去するまで論争が続いたKingsland。経済学者のA. B. ウルフ、人口学者のジョージ・ハンドリー・ニブス、統計学者のエドウィン・ウィルソンなどがパールの主張に批判を加えているKingsland。動物学者のジェームズ・グレイ 、ランスロット・ホグベン、遺伝学者のシューアル・ライトからは、他のS字型曲線を使っても個体群成長のデータに当てはめができるので、ロジスティック曲線を使う必然性が欠けていることについて指摘を受けているKingsland。

ロジスティック曲線の有効性を支持し、その普及を担った人たちもいる。イギリスの統計学者ウドニー・ユールは、パールの理論を1924年のイギリスの学会で発表している。ユールはロジスティック曲線は長期予測には適用できないと考えており、その点を強調したが、基本的にはパールの研究を支持していたKingsland。アルフレッド・ロトカも、ロジスティック曲線の有効性を理解し、ロジスティック方程式について1925年の自書の中で一章を与えて説明したKingsland。ただしロトカは、ロジスティック方程式は実現象の近似の一種であるという考えを保っていたKingsland。ロシアのゲオルギー・ガウゼも、近似の一種と受け止めながらも、ロジスティック方程式が同種の集団の中での生存競争を定量的に表すことができると述べているKingsland。1934年、ガウゼは微生物の実験によってロジスティック方程式の検証を行い、この検証はロジスティックモデルを個体群動態論における古典的理論の一つとして確固たるものとした日本生態学会。モデルの限界には注意が払われながらもロジスティック方程式の受容は広まっていき、1940年代後半には個体群解析における一般的な道具として確立したKingsland。

10.3. ロジスティック方程式からの発展
---------------------

その後は、より現実的な個体群変動を表すことができるように、ロジスティック方程式を修正したモデルが提案されてきた。1948年には、ジョージ・イヴリン・ハッチンソンが時間遅れの影響をロジスティック方程式に導入した研究を行ったKingsland。ロジスティック方程式の前提条件を満たすような環境であっても、個体数が一定に収束せず、多くなったり少なくなったりをいつまでも繰り返すような生物実験の結果も得られた山口。京都大学の内田俊郎と藤井宏一がヨツモンマメゾウムシの培養実験でそのような結果を得たことを1953年に発表している巌佐。内田らは、ロジスティック方程式をもとにした差分方程式でこの結果を分析し、個体数の振動を再現した山口。

ロジスティック方程式における //r// は個体群密度がとても低いときの増加率で表しており、密度が低いときにどれだけ素早く繁殖できるかを意味している日本生態学会。また、//K// はその環境下で生存できる個体数あるいは個体群密度の上限を示す日本生態学会。1967年、ロバート・マッカーサーとエドワード・オズボーン・ウィルソンは、この //r// と //K// に着目して、島における生物個体群の定着と絶滅に関する理論を発案した木元。彼らの理論によれば、ある生物の島への定着が成功するには大きな //r// を持つことが重要であり、絶滅の回避には大きな //K// を持つことが重要であるとし、それぞれの方向へ淘汰されることを r淘汰、K淘汰と呼んだ木元。この説はr-K戦略説と呼ばれ、生物の生活史の進化に種内競争の観点から説明を与えた日本生態学会。

物理学から数理生態学へ転向してきたロバート・メイも個体群変動の問題に取り組んだ山口。メイはロジスティック方程式の離散化を行い、その式の解は、通常のロジスティック方程式の解とは全く異なる、現在ではカオスと呼ばれる非常に複雑な振る舞いを起こすことを示したスチュアート。この結果は1974年と1975年に発表され、大きな反響を得ると共に、その後のカオス理論の隆盛に大きく寄与することになるスチュアート。

11. 脚注
------

11.1. 注釈
--------

11.2. 出典
--------

12. 参考文献
--------

※文献内の複数個所に亘って参照したものを示す。

13. 関連項目
--------

フェルミ分布関数
増殖曲線

14. 外部リンク
---------

フェルフルストの原論文
last=Verhulst - Google ブックス
first= Pierre-François - Göttinger Digitalisierungszentrum