コンテンツにスキップ

英文维基 | 中文维基 | 日文维基 | 草榴社区

「ロジスティック方程式」の版間の差分

出典: フリー百科事典『ウィキペディア(Wikipedia)』
削除された内容 追加された内容
人口動態学について
Cewbot (会話 | 投稿記録)
m 解消済み仮リンクアルフレッド・ロトカを内部リンクに置き換えます (今回のBot作業のうち66.3%が完了しました)
 
(21人の利用者による、間の99版が非表示)
1行目: 1行目:
[[File:Logistic curve, r=1, K=100, N0=1.png|thumb|270px|ロジスティック方程式の解曲線(ロジスティック曲線)の一例。S字の形を描き、環境収容力に収束する。]]
{{参照方法|date=2013年6月22日 (土) 00:12 (UTC)}}
[[File:Drosophila melanogaster laboratory culture-vial.jpg|thumb|230px|培養容器内の[[キイロショウジョウバエ]]。ロジスティック曲線に当てはまる個体数増加が確認された例である。]]
'''ロジスティック方程式'''は、[[個体群生態学]]において、[[個体群]]成長のモデルとして考案された[[微分方程式]]である。その後、[[カオス理論]]の出発点の一つともなり、現在では、[[生態学]]のみならず、多くの分野で応用が行われている。
'''ロジスティック方程式'''(ロジスティックほうていしき、英語:logistic equation{{Sfn|Strogatz|2015|p=25}})は、[[生物]]の個体数の変化の様子を表す[[数理モデル]]の一種である。ある単一種の生物が一定環境内で増殖するようなときに、その生物の個体数([[個体群]]サイズ)の変動を予測できる。人間の場合でいえば、[[人口]]の変動を表すモデルである。


1838年にベルギーの数学者[[ピエール=フランソワ・フェルフルスト]](Pierre-François Verhulst)によって、ロジスティック方程式は最初に発案された。フェルフルストは、1798年に発表されて大きな反響を呼んだ[[トマス・ロバート・マルサス]]の『[[人口論]]』の不自然な点を解消するために、このモデルを考案した<ref name="大澤2007">{{cite book| 和書 |author=大澤光 |title=社会システム工学の考え方 |year=2007 |edition=初版 |publisher=オーム社 |isbn =978-4-274-06675-7 |pages=193-194}}</ref>。マルサスは『人口論』で、人口は原理的に[[指数関数]]的に増加することを指摘した{{Sfn|スチュアート|2012|p=78}}。しかし、実際には環境や資源は限られているため、人口の増加にはいずれブレーキがかかると考えるのが自然である。人口が増えるに連れて人口増加率は低減し、人口はどこかで飽和すると考えられる。ロジスティック方程式はこの点を取り入れて、生物の個体数増殖をモデル化したものである。フェルフルスト以後には、アメリカの生物学者[[レイモンド・パール]](Raymond Pearl)が式を普及させた。
== 個体群増加のモデル ==
[[生物]]の個体数の増え方に関する研究は、[[個体群生態学]]の分野に属する。[[人口]]推計や、[[害虫]]発生の予想などの応用的側面もあり、古くから研究が行われた。多くの生物では、実際に生存するより遙かに多くの子孫を作り、それがそのまま生き残れば、あっという間に莫大な個体数となる。[[ねずみ算]]など、数学的小話の種である。しかし、これでは現実とは違いすぎる。そのため、実際の[[個体群成長]]を扱うためには、より現実的な[[数理モデル|数学モデル]]が必要となる。


具体的には、ロジスティック方程式は
ただし、一般に生物個体数は整数の値をとるものであり、多くの場合、繁殖は特定の時期に行われるので、個体数増加は段階的な形を取る。しかし、数学的扱いを簡便にするために、その増加も個体数も連続した値をとるものと見なして扱うことが多い。
:<math>\frac{dN}{dt}\ = r N \left(1- \frac{N}{K} \right)</math>
という[[微分方程式]]で表される。''N'' は個体数、''t'' は[[時間]]、''dN''/''dt'' が個体数の増加率を意味する。''r'' は[[内的自然増加率]]、''K'' は[[環境収容力]]と呼ばれる定数である。個体数が増えて環境収容力に近づくほど、個体数増加率が減っていくというモデルになっている。


式の解(個体数と時間の関係)はS字型の曲線を描き、個体数は最終的には環境収容力の値に収束する。この曲線や解の関数は'''ロジスティック曲線'''や'''ロジスティック関数'''として知られる。方程式の名称は、ロジスティック式やロジスティックモデル、ロジスティック微分方程式と表記される場合もある<ref name="巌佐2008">{{cite book |和書 |author=巌佐庸 |title=生命の数理 |publisher=共立出版 |date=2008-02-25 |edition=初版 |isbn=978-4-320-05662-6 |pages=2-3 }}</ref><ref name="西村欣也"/><ref name="アリグッド_92">{{cite book |和書 |author=K.T.アリグッド; T.D.サウアー; J.A.ヨーク |translator =星野高志・阿部巨仁・黒田拓・松本和宏 |others=津田一郎(監訳) |editor=シュプリンガー・ジャパン |title=カオス 第2巻 力学系入門 |publisher=丸善出版 |year=2012 |isbn=978-4-621-06279-1 |page=92}}</ref>。発案者の名からVerhulst方程式、発案者と普及者の名からVerhulst-Pearl方程式とも呼ばれる{{Sfn|ティーメ|2006|p=38}}。
通常、親が作る子孫の数は、ほぼ一定であるから、増加率を''r'' とすれば、個体数''N'' の個体群における時間に対する(絶対)増加率は
:<math>\frac{dN}{dt}\ = rN</math>
で表される。これは指数曲線になって、あっという間に[[人口爆発]]を引き起こす。この様な個体群成長の型を、生物個体(人口)の増加が[[幾何級数]]的であることを最初に指摘した[[トマス・ロバート・マルサス]]にちなんで'''マルサス的成長'''と呼ぶこともある。


ロジスティック方程式は、[[個体群生態学]]あるいは[[個体群動態論]]における数理モデルとしては入門的なものとして位置づけられ、より複雑な現象に対応する基礎を与える<ref name="西村欣也"/>。数学分野としては、[[微分方程式|微分方程式論]]や[[力学系|力学系理論]]の初等的な話題としても取り上げられる<ref>{{cite book |和書 |author=稲岡毅 |title=基礎からの微分方程式―実例でよくわかる |publisher=森北出版 |year=2012 |isbn=978-4-627-07671-6 |pages=22-23}}</ref>{{Sfn|Hirsch et al.|2007|pp=4&ndash;7}}。
しかし、現実の生物は、ある特定の[[環境]]下で生活しており、そこに生活できる個体数には上限があると見るのが自然である。つまり、個体数が多くなると、その増加に[[ブレーキ]]がかかるものと想像される。そこで、そのような、現実の個体数変化を説明するためには、次のような性質の式が必要になる。
*個体数''N'' = 0 では、増加率は''r'' = 0 になる。
*個体数''N'' が増加するにつれ、増加率''r'' は減少する。
*環境の収容可能個体数に限度があるから、その数を''K'' とすれば、''N'' = ''K'' のとき、増加率は''r'' = 0 になる。


== 生物の個体数のモデル ==
== ロジスティック方程式の表現 ==
[[File:FibonacciRabbit.svg|thumb|310px|フィボナッチによるウサギのつがいの増殖問題]]
ロジスティック方程式は、1838年に[[ピエール=フランソワ・フェルフルスト]](ベルハルストとも)が、人口増加を説明するモデルとして考案した<ref>彼が[[兵站学]]([[ロジスティクス]])教官であったためロジスティックと命名したといわれる</ref>。その後、独自に同様の式を提示した個体群生態学者などもおり、次第に、個体群モデルの基礎となった。ロジスティック方程式は、上記の条件をすべて備えている。
{{see also|個体群動態論}}
生物の個体数の変動については古くから興味を持たれ、研究が行われてきた{{Sfn|マレー|2014|p=1}}。[[フィボナッチ数]]の発見に繋がった[[レオナルド・フィボナッチ]]の[[ウサギ]]の個体数の問題が、おそらく最も古い個体数の[[数理モデル]]といわれる{{Sfn|スチュアート|2012|p=333}}。


[[生物]]の個体数の増え方に関する研究は、[[個体群生態学]]の分野に属する{{Sfn|日本数理生物学会(編)|2008|p=61}}。ここで、[[個体群]]とは簡単には、ある領域に生息している単一の[[種 (分類学)|種]]の個体の集まりのことを指す{{Sfnm|日本数理生物学会(編)|2008|1p=181|寺本|1997|2p=2}}。
ロジスティック方程式は、次の式である。
:<math>\frac{dN}{dt}\ = r \left( \frac{K - N}{K} \right) N</math>


この個体群の大きさ(個体群サイズ)の指標としては、個体群内の総個体数が使用される{{Sfn|大串|2014|p=49}}。個体数の代わりに、領域の単位面積当たりの個体数である個体群密度や単位面積当たりの生物の総重量である[[生物量]]が、個体群サイズとして適切な指標となる場合もある{{Sfn|寺本|1997|pp=2&ndash;3}}。人間でいえば、これらの指標は[[人口]]や[[人口密度]]に相当する{{Sfn|瀬野|2007|p=1}}。
ここで、''K'' は[[環境収容力]]、つまり、その環境における個体数の定員である。''r'' は(相対)内的増加率で、その生物が実現する可能性のある、最大の増加率である。実際の増加率''r'' (''K'' - ''N'' )/''K'' は個体数''N'' が環境収容力''K'' に近づくにつれて減少し、''N'' = ''K'' ならば増加率は 0 である。''N'' > ''K'' だと、増加率は負となり、個体数が''K'' になるまで減少する。


=== マルサスモデル ===
また、ここで''k'' = ''r'' /''K'' と置けば、
[[File:Malthusian growth curves.png|thumb|280px|マルサスモデルによる個体数増加曲線の様子。赤色が ''m'' = 4、紫色が ''m'' = 2、藍色が ''m'' = 1。いずれも最初は ''N'' =''1'' だが、その後の急激な成長が見て取れる。]]
:<math>\frac{dN}{dt}\ = N(r - kN)</math>
{{see also|マルサスモデル}}
と書ける。この場合、''k'' は、一個体の増加によって増加率が減少する率を現す。つまり、個体群密度の増加が増加率にブレーキをかけるので、これを[[密度効果]]という。
多くの生物では、親は多くの子孫を作るので、それがそのまま生き残ると仮定すれば、あっという間に莫大な個体数となる。[[ねずみ算]]など、数学的小話の種である{{Sfn|巌佐|1990|p=2}}。まずはこのような単純なものが、生物個体数の増加モデルとして考えられる。


ある個体群において、時刻 ''t'' に個体数が ''N'' 体が存在しているとする。実際の生物個体数は不連続な値([[整数]])をとるものであるが、数学的扱いを簡便にするために、個体数は[[連続 (数学)|連続]]な値([[実数]])をとるものとする(1.5体といったような値も含める)ことがしばしば行われる{{Sfn|山口|1992|p=59}}。実際の生物でいえば、個体数が多かったり各個体の世代が重なったりしていれば、このような近似も妥当性を帯びてくる{{Sfn|ブラウン|2012|p=30}}<ref name="Silvertown"/>。個体数を連続な値とすれば、個体数の増加率は ''N'' の[[時間微分]] ''dN''/''dt'' で表すことができる{{Sfn|ハーバーマン|1992|p=13}}。
== ロジスティック方程式の解 ==
[[Image:SigmoidFunction.png|thumb|シグモイド関数(ロジスティック関数)の例。]]
初期値''N'' (0) が 0 < ''N'' (0) < ''K'' を満たす場合、ロジスティック方程式で表される微分方程式の解は
:<math>N(t) = K \, \varsigma_1 (r K (t_0 - t)) = \frac{K}{1 + \exp(r K (t_0 - t))}</math>
となる(<math>\varsigma_1</math>は標準[[シグモイド関数]]、''t''<sub>0</sub> は初期値で決まる任意性)。<!--この関数は、<math>\lim_{t \rightarrow -\infty} N \rightarrow 0</math> で始まり、初めはゆっくりと、つぎに急激に、最後にまたゆっくりと増加し、<math>\lim_{t \rightarrow \infty} N \rightarrow K</math> で終わる。-->この解は'''ロジスティック関数'''と呼ばれる(描く曲線は'''ロジスティック曲線'''と呼ばれる)。
<!--
<math>r = 1</math>、<math>K = 1</math>、<math>t_0 = 0</math> のとき、解は[[シグモイド関数]]になる。逆に、一般のロジスティック関数の解をシグモイド関数 <math>\mathrm{P}(x)</math> を使って


さらに話を単純化するために、個体は環境を出入りしないという状況を想定する{{Sfn|渡辺|2012|pp=50&ndash;51}}。この場合、個体の出生と死亡という2つの要因のみによって個体数は増減する{{Sfn|渡辺|2012|pp=50&ndash;51}}。個体群の[[出生率]]が[[死亡率]]を上回っていれば、個体数は増え続けるということになる{{Sfn|マレー|2014|p=1}}。さらに簡略化するために出生率と死亡率を常に一定であるとする{{Sfn|マレー|2014|p=1}}。個体数当たりの出生率を ''b''、個体数当たりの死亡率を ''d'' とすれば、個体数の増加率は差し引きした ''b'' &minus; ''d'' に個体数 ''N'' を掛け合わせた値となる{{Sfn|寺本|1997|pp=4&ndash;5}}。よって個体数増加率 ''dN''/''dt'' は
:<math>N = K \, \mathrm{P}(r K (t_0 - t)) \,</math>
:<math>\frac{dN}{dt}\ = mN</math>
という[[微分方程式]]で表される{{Sfn|巌佐|2015|p=17}}。ここで ''m'' は比例定数であり、''m'' = ''b'' &minus; ''d'' である<ref name="Silvertown"/>。


このような式で表される個体数増加は ''t'' の[[指数関数]]となり、人間でいえば、あっという間に[[人口爆発]]を引き起こすことになる<ref>{{cite book |和書 |author=山本将史 |title=Excelで学ぶ微分積分 |publisher=オーム社 |year=2012 |isbn=978-4-274-06888-1 |page=264}}</ref>。このような個体群成長のモデルは、生物個体(人口)の増加が[[幾何級数]]的であることを最初に指摘した[[トマス・ロバート・マルサス]]に因んで[[マルサスモデル]]と呼ばれる{{Sfnm|人口研究会(編)|2010|1pp=281&ndash;282|マレー|2014|2p=38}}。比例定数 ''m'' もマルサスの名からマルサス係数と呼ばれ、[[物理単位|単位]]は一個体当たりの増加率となる<ref name="巌佐2008"/>。
と表すことができる。また、[[双曲線関数|双曲線正接関数]] <math>\tanh</math> を使って表すこともできる。
-->


しかし、このモデルは現実と違いすぎる{{Sfn|マレー|2014|pp=1&ndash;2}}。現実の生物は、限られた[[環境]]下で生息しており、個体数が多くなると、各個体にとって必要な資源が得にくくなる{{Sfn|巌佐|1990|p=4}}。そこに生息できる個体数には上限があると見るのが自然である{{Sfn|日本生態学会(編)|2004|p=125}}。つまり、個体数が多くなると、その増加にブレーキがかかるものと想像される{{Sfn|日本生態学会(編)|2004|p=123}}。このような一種内での資源の取り合いは'''種内競争'''と呼ばれ、生物における[[競争 (生物)|競争関係]]の一種である{{Sfn|巌佐|2015|p=21}}。
ロジスティック関数は非線形だが、次の変換によって線形の扱いやすい関数にすることができる。これはフィッシャ・プライ変換(Fisher-Pry transform)と呼ばれる<ref name=watanabe>{{cite|和書 |author=渡辺千仭 |title=技術経済システム |publisher=創成社 |year=2007 |isbn=978-4-7944-3089-2 |pages=87}}</ref>。
:<math>\begin{align}
& \ln FP(t)=rK(t_0-t), \\
& FP(t):=\frac{N(t)/K}{1-N(t)/K}
\end{align}</math>


==ロジスティック方程式==
== 生物学的解釈 ==
[[File:Logistic equation m vs N.svg|thumb|個体数 ''N'' と一個体当たり個体群増加率 ''m'' の関係]]
ロジスティック方程式そのものは、生物学的には、かなりありえない仮定に基づいている。
上記のようにマルサスモデルは非現実的な面を持つ。個体数が多くなると増加率が抑えられることを表現するために、個体数 ''N'' が増加するにつれて増加率 ''m'' が減少するモデルが考えられる{{Sfn|寺本|1997|p=8}}。また、個体数がある上限を超えたら増加率は負となり、個体数は減少に向かうと考えられる{{Sfn|Strogatz|2015|p=25}}。これらの点を簡単に表せば、比例定数 ''m'' を
*まず、個体数の増加が連続的に生じること。多くの生物では、特定の時期にのみ増加が起こる。[[昆虫]]など、世代が重ならないものでは、個体数増加は世代を追って段階的に生じる。
:<math>m=r \left(1- \frac{N}{K} \right)</math>
*個体数増加は増加率を抑制するが、親個体も子の個体も、同じだけの率で抑制に関わる。多くの生物では、親子では大きさが異なるので、このようなことはありえない。昆虫では、親と子では生活の場が異なるものも多い。
と置ける{{Sfn|巌佐|1990|p=4}}。すなわち、''m'' の値は個体数がゼロに限りなく近いときに最大値で、その後は ''N'' の値の増加に比例して ''m'' の値は減少するというモデルである<ref name="渡辺守">{{cite book |和書 |author=渡辺守 |title=昆虫の保全生態学 |date=2007-12-20 |edition=初版 |publisher=東京大学出版会 |isbn =978-4-13-062215-8 |pages =39-40}}</ref>。これをマルサスモデルに代入して、次の微分方程式を得ることができる。
*個体数の増加は、その瞬間に増加率に影響を与える。もちろん現実には、瞬間ということはあり得ないにせよ、親子で大きさが異なったり、昆虫など、親と子では生活の場そのものが異なる場合もあり、個体数の増加が増加率に影響するまでに、かなりの時間が必要と思われる例が少なくない。
:<math>\frac{dN}{dt}\ = r N \left(1- \frac{N}{K} \right)</math>
この微分方程式を'''ロジスティック方程式'''と呼ぶ<ref name="Silvertown"/>。個体群成長モデルの一種として'''ロジスティックモデル'''とも呼ばれる<ref name="西村欣也">{{cite book |和書 |author=西村欣也 |title=生態学のための数理的方法―考えながら学ぶ個体群生態学 |publisher=文一総合出版 |year=2012 |isbn=978-4-8299-6520-7 |pages=168-169}}</ref>。この微分方程式は、数学的には ''n'' = 2 の[[ベルヌーイの微分方程式]]にも該当する{{Sfn|寺本|1997|pp=10&ndash;11}}。


ロジスティック方程式の ''K'' は[[環境収容力]]と呼ばれ、その環境が維持できる個体数を意味する{{Sfn|巌佐|1990|p=4}}。''r'' の単位は上記のマルサス係数と同じく一個体当たりの増加率だが{{Sfnm|瀬野|2007|1pp=11, 13&ndash;14|コーエン|1998|2p=112}}、特に[[内的自然増加率]]と呼ばれ、その生物が実現する可能性のある最大増加率を示している{{Sfn|瀬野|2007|p=14}}。通常のロジスティック方程式では、''K'' と ''r'' は時間に関わらず一定とみなし、正の[[定数]]と考える{{Sfnm|コーエン|1998|1p=112|マレー|2014|2p=2}}。
したがって、ロジスティック方程式を単純に適用できるのは、ほとんど大きさに差のない形で増殖し、始終増えている[[細菌]]や、世代が完全に重なって、繁殖期がはっきりしない[[ヒト]]のようなものに限られるともいわれる。しかし、実際には様々な生物の個体群研究において、ロジスティック方程式は個体数変化の基本的モデルとして利用され、多くの成果が得られている。

=== ロジスティック効果 ===
{{see also|密度効果}}
マルサスモデルからロジスティック方程式へ拡張したときに行ったことは、個体群生態学における[[密度効果]]を取り入れたことに相当する{{Sfn|渡辺|2012|p=51}}。上記では ''N'' を個体数として説明したが、ロジスティック方程式では有限な環境を前提にしているので、''N'' は単位面積当たりの個体数である個体群密度でもある{{Sfn|瀬野|2007|p=20}}。個体群密度がその個体群自身の変動に影響を与えることは、密度効果という名称で呼ばれる{{Sfn|日本数理生物学会(編)|2008|p=4}}。特にロジスティック方程式では、個体群密度が高くなると増加率に負の効果を与える種類の密度効果となっており、これを'''ロジスティック効果'''と呼ぶ<ref>{{cite book |和書 |author=ミンモ・イアネリ、稲葉寿、國谷紀良 |title=人口と感染症の数理―年齢構造ダイナミクス入門 |year=2014 |edition=初版 |publisher=東京大学出版会 |isbn =978-4-13-061309-5 |page =52}}</ref>{{Sfn|人口研究会(編)|2010|p=307}}。

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

===個体数と増加率の関係===
[[File:Logistic equation dNdt vs N.svg|thumb|280px|縦軸が ''dN''/''dt''、横軸が ''N'' のグラフ。ロジスティック方程式における、''dN''/''dt'' と ''N'' の関係が示されている。]]
ロジスティック方程式における個体数増加率 ''dN''/''dt'' と個体数 ''N'' の関係に着目すれば、この関係は初等教育でも習う[[二次関数]]そのものとなっており、''dN''/''dt'' と ''N'' の[[グラフ (関数)|グラフ]]は[[放物線]]を描く{{Sfn|寺本|1997|p=22}}。方程式を解析的に解かなくとも、''N'' と ''dN''/''dt'' がどのような変化を起こすのかを、以下のようにグラフから読み解くこともできる。

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

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

== 式の解 ==
=== ロジスティック曲線 ===
[[File:Malthusian growth vs logistic growth.png|thumb|280px|マルサスモデルによる指数関数的増加曲線(赤)とロジスティック曲線(青)の比較]]
ロジスティック方程式は[[非線形]]の微分方程式だが、標準的な微分方程式の解法である[[変数分離法]]を利用して解くことができる{{Sfn|山口|1992|pp=62&ndash;65}}。時間 ''t'' = 0 における初期個体数を ''N''<sub>0</sub> とすると、''t'' の関数として以下の解が得られる{{Sfn|スチュアート|2012|p=417}}{{efn|一例として以下のように解くことができる。''N'' の値の範囲を 0 < ''N'' < ''K'' に限定して解く方法と{{Sfn|山口|1992|pp=62&ndash;65}}、特に限定せずに解く方法がある<ref>{{cite book |和書 |author=小寺忠・長谷川健二 |title=工学系学生のための常微分方程式 |publisher=森北出版 |year=2006 |edition=第2版 |isbn=978-4-627-07452-1 |pages=44-46}}</ref>。ここでは範囲を限定しない解き方を示す。まずロジスティック方程式を[[変数分離]]変形して
:<math>\frac{K}{N(K-N)}dN = r dt</math>
を得る。さらに左辺を[[部分分数分解]]すれば
:<math>\left( \frac{1}{N}+\frac{1}{K-N} \right)dN = r dt</math>
となる。両辺を積分して
:<math>\ln \vert N \vert-\ln \vert K-N \vert =rt +C</math>
:<math>\ln \left | \frac{N}{K-N} \right \vert =rt +C</math>
となり、ここで ''C'' は[[積分定数]]である。両辺の指数をとり、絶対値を外せば
:<math>\frac{N}{K-N} =\pm e^{rt +C}</math>
となる。''t'' &#x3D; 0 のときの ''N'' の値を ''N''<sub>0</sub> で表せば
:<math>\frac{N_0}{K-N_0} =\pm e^{C}</math>
なので、これを上式に代入して
:<math>\frac{N}{K-N} =\frac{N_0}{K-N_0} e^{rt}</math>
となる。式を整理して
:<math>N=\frac{N_0 K e^{rt}}{K-N_0+N_0 e^{rt}} </math>
となる。
}}。
:<math>N=\frac{N_0 K e^{rt}}{K-N_0+N_0 e^{rt}} </math>
ここで ''e'' は[[ネイピア数]]である。分母・分子を ''N''<sub>0</sub>''e<sup>rt</sup>'' で割り、次のような形でも示される<ref name="巌佐2008"/>。
:<math>N=\frac{K}{1+(K/N_0-1)e^{-rt}} </math>
この解の関数を'''ロジスティック関数'''(英語:logistic function)<ref>{{cite book |和書 | author=和田光平 | year = 2015 | title = 人口統計学の理論と推計への応用 | publisher = オーム社 | edition = 初版 | isbn = 978-4274217166 |pages =31-34}}</ref>、この解によって描かれる曲線を'''ロジスティック曲線'''(英語:logistic curve)と呼ぶ{{Sfn|人口研究会(編)|2010|p=307}}。この曲線に従う個体群成長は、ロジスティック成長やロジスティック増殖とも呼ばれる<ref>{{cite book |和書 |editor= 日本生物物理学会 |others=巌佐庸(担当編集委員) |title=数理生態学 |series=シリーズ・ニューバイオフィジックス (10) |year=1997 |edition=初版 |publisher=共立出版 |isbn =4-320-05473-3 |page =33}}</ref>{{Sfn|寺本|1997|p=10}}。関数は ''t'' &rarr; &infin; の[[極限]]で ''N'' &rarr; ''K'' となり、マルサスモデルと異なり発散しないことが確認できる{{Sfn|Hirsch et al.|2007|p=5}}。

=== 曲線の形状 ===
[[File:Logistic curve in some initial conditions.png|thumb|270px|right|いくつかの ''N''<sub>0</sub> から始まるロジスティック曲線。''N'' > 0 の範囲では、時間の経過に従って ''N'' は ''K'' に収束する。]]
[[File:Logistic curve overall view.png|thumb|270px|right|時間と個体数が負の場合も含めたロジスティック曲線の全体図。縦軸を ''N''/''K''、横軸を ''rt'' として無次元化している。]]
横軸を ''t''、縦軸を ''N'' とした平面上にロジスティック関数の[[グラフ (関数)|グラフ]]を描くと、曲線が描かれる。この曲線は前述のとおりにロジスティック曲線と呼ばれる。初期個体数が3つの範囲 ''N''<sub>0</sub> < 0, 0 < ''N''<sub>0</sub> < ''K'', ''K'' < ''N''<sub>0</sub> のどれに該当するかによって、曲線の形状は大きく異なってくる{{Sfn|Hirsch et al.|2007|pp=5&ndash;6}}。ただし、''N''<sub>0</sub> < 0 の範囲は負の個体数というものを意味するので、生物のモデルとしてはあまり意味がない{{Sfn|Hirsch et al.|2007|p=5}}。時間 ''t'' = 0 から ''t'' &rarr; &infin; の[[極限]]までのロジスティック曲線の様相は、それぞれの ''N''<sub>0</sub> の値ごとに、以下のようになっている。

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

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

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

以上をまとめると、''N''<sub>0</sub> > 0 であれば(個体が存在してさえいれば)、どんな初期個体数であっても、個体数は最終的に常に環境収容力の値に収束していくということである{{Sfn|Strogatz|2015|p=26}}。あるいは、''N''<sub>0</sub> = 0 であれば(個体が存在してなければ)、個体数は 0 のままということである{{Sfn|Strogatz|2015|p=26}}。

最後に、生物個体数のモデルとしては無意味であるが ''N''<sub>0</sub> < 0 の場合も見てみると、この場合 ''N'' は時間発展に従って減少し続け、有限時間内で &minus;&infin; へ発散する曲線を描く{{Sfn|Hirsch et al.|2007|p=5}}。

===平衡状態の安定性===
上記で、''N'' = 0 および ''N'' = ''K'' のときはいくら時間が経過しても個体数 ''N'' は増加も減少もしないことから、これらの状態を平衡状態や定常状態と呼ぶことを説明した。平衡状態では、''N'' = 0 または ''N'' = ''K'' という一点に留まり続ける。数学の[[力学系]]分野では、このような点を[[不動点]]や平衡点と呼ぶ{{Sfnm|Strogatz|2015|1p=19|Hirsch et al.|2007|2p=2}}。平衡状態には'''安定な平衡状態'''と'''不安定な平衡状態'''がある{{Sfn|巌佐|2015|p=23}}。安定な平衡状態とは、 その平衡状態の点から少しずれたとしても、時間が経過すれば平衡状態へ戻り、収束することを意味している{{Sfn|巌佐|2015|p=23}}。また、不安定な平衡状態とは、平衡状態の点から少しずれたとき、時間経過すると平衡状態とのズレはどんどん大きくなっていき、平衡状態に戻らないことを意味している{{Sfn|マレー|2014|p=5}}。ロジスティック方程式の場合は、''N'' = ''K'' 時の平衡状態が安定、''N'' = 0 時の平衡状態が不安定となっている{{Sfn|Strogatz|2015|pp=25&ndash;26}}。すなわち、初期個体数 ''N''<sub>0</sub> が ''K'' または 0 であれば、時間経過によらず常に同じ値を取り続けることは同じだが、''N''<sub>0</sub> が平衡状態から少しずれたときの挙動は正反対となる{{Sfn|巌佐|2015|pp=22&ndash;23}}。

[[File:Logistic curve vector field.png|thumb|270px|right|ロジスティック曲線とその傾きのベクトル場の様子]]
この安定・不安定の様子は、ロジスティック曲線の[[傾き (数学)|傾き]]を[[ベクトル場]]として表すことで読み取ることができる{{Sfn|Hirsch et al.|2007|pp=5&ndash;6}}。時間経過に従って、全ての解は、これらのベクトルの矢印に沿って動いていく<ref name="アリグッド_94">{{cite book |和書 |author=K.T.アリグッド; T.D.サウアー; J.A.ヨーク |translator =星野高志・阿部巨仁・黒田拓・松本和宏 |others=津田一郎(監訳) |editor=シュプリンガー・ジャパン |title=カオス 第2巻 力学系入門 |publisher=丸善出版 |year=2012 |isbn=978-4-621-06279-1 |page=94}}</ref>。初期個体数が ''N''<sub>0</sub> > 0 であれば、''t'' &rarr; &infin; で ''N'' は ''K'' に収束し、''N''<sub>0</sub> < 0 であれば、''t'' &rarr; &infin; で ''N'' は &minus;&infin; に発散することが分かる{{Sfn|Hirsch et al.|2007|pp=5&ndash;6}}。

[[File:Stability of logistic equation.png|thumb|270px|right|''N'' と ''dN''/''dt'' の関係曲線。''N'' 軸と曲線の交点が平衡状態の点で、右が安定な点、左が不安定な点である。]]
あるいは、上記で説明した個体数 ''N'' と増加率 ''dN''/''dt'' の関係曲線からも、安定か不安定かの判別が可能である{{Sfn|巌佐|2015|p=24}}。''N'' = ''K'' の点の右側に点があるとき、''dN''/''dt'' の値は負なので、''N'' は減少していき、''K'' に近づくことになる。''N'' = ''K'' の点の左側に点があるときは、''dN''/''dt'' は正なので、''N'' は増加していき、同じく ''K'' に近づくことになる{{Sfn|Hirsch et al.|2007|p=6}}。''N'' = 0 の点についても、左右にずれたときの ''dN''/''dt'' の値の正負から、0 の点から離れていくことが理解できる{{Sfn|Hirsch et al.|2007|p=6}}。

あるいは、[[安定性理論]]における線形安定性解析の考えにもとづいて、より一般的に安定性を判別することもできる。''dN''/''dt'' = ''f''(''N'') 、その ''N'' による微分を ''d''(''f''(''N''))/''dN'' = ''f''&thinsp;&prime;(''N'')、平衡状態の点を ''N<sub>e</sub>'' と置くとする。このとき、''f''&thinsp;&prime;(''N<sub>e</sub>'') < 0 ならば ''N<sub>e</sub>'' は安定な平衡点で、''f''&thinsp;&prime;(''N<sub>e</sub>'') > 0 ならば ''N<sub>e</sub>'' は不安定な平衡点であると判別できる{{Sfnm|マレー|2014|1p=5|Strogatz|2015|2p=28}}。ロジスティック方程式の場合は、
:<math>f(N)=\frac{dN}{dt}\ = r N \left( 1-\frac{N}{K} \right) </math>
なので、
:<math>f'(N)=r \left( 1-\frac{2N}{K} \right)</math>
となり、''f''&thinsp;&prime;(''K'') = &minus;''r'' < 0, ''f''&thinsp;&prime;(0) = ''r'' > 0 となることが確認できる{{Sfn|Strogatz|2015|p=29}}。

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

*環境内には単一の種か、あるいは同等とみなせる種のみが存在する{{Sfn|瀬野|2007|pp=5, 21&ndash;22}}。
*対象の生物の各世代(親子)は連続的に重なっている{{Sfn|マレー|2014|p=37}}。すなわち、連続的に子が生まれ、親と子が共存する期間が存在する{{Sfn|山口|1992|p=71}}。
*個体は一定の大きさの環境内に常に存在する。すなわち、環境から移出したり、外部から移入が無い{{Sfn|コーエン|1998|p=112}}。(用語としては'''閉じた個体群'''とも呼ばれる{{Sfn|ティーメ|2006|p=7}})
*環境の大きさは変わらず、一定状態が保たれる{{Sfn|コーエン|1998|p=112}}。
*個体群のために、食糧や資源が一定して供給される{{Sfnm|Strogatz|2015|1p=27|コーエン|1998|2p=112}}。

[[ショウジョウバエ]]や[[真正細菌]]といった、[[微生物]]や単純な生物を一定環境で増殖させた場合は、上記の条件に近く、ロジスティック方程式によって個体数変化の正確な予測ができる{{Sfnm|巌佐|2015|1p=24|コーエン|1998|2p=116|Strogatz|2015|3p=27}}。しかし、例えば[[鹿]]や[[鳥類]]などのような、一定環境のもとで増殖する設定が成立しない個体群成長には、ロジスティック方程式を適用することはできない{{Sfn|巌佐|2015|p=24}}。

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

==実際のデータへの適用例==
===実験生物===
[[File:Gause's experiment and fitted logistic curves (single species).svg|thumb|center|550px|ソ連・ロシアの生物学者[[ゲオルギー・ガウゼ]]による2種の[[酵母]](''Saccharomyces cerevisiae'', ''Schizosaccharomyces kefir'')の個体群サイズ成長実験データ<ref name="Gause1932">{{cite journal |last=Gause |first=G. F. |year=1932 |title=Experimental Studies on the Struggle for Existence |url=http://jeb.biologists.org/content/9/4/389 |journal=Journal of Experimental Biology |publisher=The Company of Biologists Ltd |volume=9 |issue=4 |pages=389-402 |issn=1477-9145}}</ref>と、それらのデータに対して[[曲線あてはめ|フィッティング]]させたロジスティック曲線を示した図。縦軸は菌全体の体積、横軸は時間を示している<ref name="Gause1932"/>。このガウゼの実験はロジスティック曲線がよく当てはまった個体群成長実験としてよく知られる{{Sfn|スチュアート|2012|pp=335&ndash;336}}。]]

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

*[[キイロショウジョウバエ]]{{Sfn|内田|1972|p=29}}
*[[ゾウリムシ]]{{Sfn|山口|1992|pp=67&ndash;68}}<ref>G. F. Gause. {{Google books|v01OToAhJboC|The Struggle for Existence|page=106}}</ref>
*[[大腸菌]]{{Sfn|山口|1992|pp=67&ndash;69}}
*[[タマミジンコ科|タマミジンコ]]{{Sfn|巌佐|1990|p=6}}
*[[出芽酵母]]、[[分裂酵母]]<small>(''Saccharomyces cerevisiae'', ''Schizosaccharomyces kefir'')</small>{{Sfn|スチュアート|2012|pp=335&ndash;336}}<ref>{{Cite book |last=Gause |first=G. F. |url=https://books.google.co.jp/books?id=v01OToAhJboC&lpg=PP1&hl=ja&pg=PA77#v=onepage&q&f=false |edition=Dover Phoenix Editions |year=2003 (original 1934) |title= The Struggle for Existence |publisher=Dover Pubns |page=77 |isbn=0-486-49520-5 }}</ref>
<gallery mode="packed-hover" height="140">
File:Drosophila melanogaster - side (aka).jpg|[[キイロショウジョウバエ]]
File:Paramecium.jpg|[[ゾウリムシ]]
File:E. coli Bacteria (7316101966).jpg|[[大腸菌]]
File:Moina sp.jpg|[[タマミジンコ科|タマミジンコ]]
File:S cerevisiae under DIC microscopy.jpg|[[出芽酵母]]
</gallery>

一方、ロジスティック曲線に当てはまるデータは得られなかったものとしては、次のような生物の実験がある。これらの実験では、時間経過後も個体数は一定に収束せず、周期的変動が繰り返されたり、大きなゆらぎが続く個体群変動となった{{Sfnm|山口|1992|1pp=69&ndash;71|Strogatz|2015|2p=27}}。
*[[ミバエ]]{{Sfn|Strogatz|2015|p=27}}
*[[コクヌストモドキ]]{{Sfn|Strogatz|2015|p=27}}
*[[マメゾウムシ]]{{Sfn|山口|1992|pp=69&ndash;71}}
<gallery mode="packed-hover" height="140">
File:Bactrocera dorsalis.jpg|[[ミバエ]]([[ミカンコミバエ]])
File:Tribolium castaneum.jpg|[[コクヌストモドキ]]
File:Chrząszcz na ziarnach fasoli.jpg|[[マメゾウムシ]](インゲンマメゾウムシ)
</gallery>

====パールのキイロショウジョウバエ飼育実験====
ロジスティック曲線を普及させたことで知られる[[レイモンド・パール]]は、{{仮リンク|ローウェル・リード|en|Lowell Reed}}と共に[[キイロショウジョウバエ]]の飼育実験を行い、この曲線を実証した。ロジスティック曲線が上手く適合する実験の具体的様子の例として、[[内田俊郎]]の著作をもとにしてパールらの実験を簡単に説明する{{Sfn|内田|1972|pp=27&ndash;29}}。
*パールが用意した環境は小さな牛乳瓶で、供給する餌にはバナナを磨り潰して寒天で固め[[出芽酵母|イースト]]を少し振りかけたものを使用した{{Sfn|内田|1972|pp=27&ndash;28}}。牛乳瓶の中にハエと餌を入れ、温度などの環境条件を一定にし、一定時間間隔でハエの個体数を調べた{{Sfn|内田|1972|p=27}}。
*実験としては3種類の実験が行われた。
*1つ目では、餌を始めに入れた後に餌を補給しなかった{{Sfn|内田|1972|p=28}}。このため、個体数が増加して一定となった後、急激に減少してほぼ全滅状態となった{{Sfn|内田|1972|p=28}}。
*2つ目では、一定時間間隔で餌の継ぎ足しを行い、一定状態が保たれる結果が得られた{{Sfn|内田|1972|p=29}}。
*3つ目では、一定時間間隔で新しい餌の入った瓶へハエを移し替え、食糧条件だけでなく、その他の環境条件も一定に保った{{Sfn|内田|1972|p=29}}。この結果でも一定状態が保たれ、ロジスティック曲線が当てはまるデータが得られた{{Sfn|内田|1972|p=29}}。

パールの元へ留学していた[[寺尾新 (動物学者)|寺尾新]]もこのハエの個体群成長研究を行った。それによれば、ロジスティック成長の特徴である個体群密度上昇にともなう個体数増加率の低下は、[[死亡率]]の上昇よりも[[出生率]]の低下によって起こっていた{{Sfn|Kingsland|1982|p=44}}。

===野外生物===
[[File:Leafcutter ant.jpg|thumb|180px|[[ハキリアリ]]]]
[[File:Northern fur seal callorhinus ursinus.jpg|thumb|180px|[[キタオットセイ]]]]
野外環境では、前提条件となるような環境が保持されることはほぼ無いため、ある個体群がロジスティック曲線が当てはまるような増加の仕方を示すことは少ない<ref name="伊藤1994">{{cite book |和書 |author=伊藤嘉昭 |title=生態学と社会―経済・社会系学生のための生態学入門 |year=1994 |edition=初版 |publisher=東海大学出版会 |isbn =4-486-01272-0 |pages =48&ndash;51}}</ref>。自然界では環境条件は常に変化し、個体群変動のパターンも様々となる{{Sfn|大串|2014|pp=53&ndash;54}}。

ロジスティック曲線によく当てはまる個体数増加が確認できた例として、[[パナマ]]熱帯雨林での[[ハキリアリ]]の1つの巣における個体数増加結果がある<ref name="伊藤1994"/>。理由としては、天敵がいないこと、雨量・温度の気象条件が安定していることなどにより、ロジスティックモデルの前提条件に近い環境であったことによるものと考えられている<ref name="伊藤1994"/>。他の野外生物でロジスティック曲線に合致した例としては、アメリカ・[[アラスカ州]]のセントポール島における[[キタオットセイ]]({{Snamei|Callorhinus ursinus}})の個体数増加の結果がある{{Sfn|レーヴンほか|2007|p=1151}}。植物の場合では、アイルランドの[[スルツェイ島]]で観測された[[コケ植物|コケ]]の成長の例がある。新規に露出した岩表面上へのコケの定着・広がり方が、ロジスティック曲線に当てはまる観測データを見せた<ref name="Silvertown">{{cite book |和書 |translator= 河野昭一・高田壮典・大原雅|author=Jonathan W. Silvertown |title=植物の個体群生体学 第2版 |year=1992 |edition=初版 |publisher=東海大学出版会 |isbn =4-486-01157-0 |pages =49-50}}</ref>。

===人口成長===
[[File:World-Population-1800-2100.png|thumb|180px|世界人口のグラフ]]
式を発案したフェルフルストは[[人口]]の成長の様子を表すためにロジスティック方程式を発案した{{Sfn|Kingsland|1982|p=30}}。式を普及させたパールとリードも、ロジスティック方程式を使った最初の個体群成長研究は人口成長に対するものであった{{Sfn|Kingsland|1982|p=31}}。彼らは共に、当時までの人口統計をもとにして[[アメリカ合衆国]]の将来の人口を予測したが、どちらの予測も実際の人口成長を言い当てることはできなかった{{Sfn|コーエン|1998|p=113&ndash;115}}。パールとリードの結果では、1700年から1940年までの値は曲線によく合致していた。彼らが当てはめたロジスティック曲線では人口はその後飽和に向かうはずだったが、実際にはそれを裏切り、1940年以後もアメリカの人口は急増状態が続いた{{Sfn|山口|1992|pp=58&ndash;59}}。

さらにパールは、当時の推定[[世界人口]]をもとに世界人口の上限値(環境収容力 ''K'')の推定を行った。1924年と1936年、パールはそれぞれ別の研究者とともに推定を行い、その値を発表した。それらの上限推定値は、前者では20億人、後者では26億人という値で、どちらも実際とはかけ離れたものとなった{{Sfn|コーエン|1998|p=115}}。

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

ただし、一見してロジスティック曲線のような個体群成長を示すデータであっても、そのデータに上手く[[曲線あてはめ]]できる数理モデルは数多く存在する{{Sfnm|山口|1992|1p=66|スチュアート|2012|2p=336}}。ロジスティック方程式のみが唯一当てはまるということはまずない{{Sfn|山口|1992|p=66}}。この式が個体群成長の「普遍則」のように受け止められるのは誤解であると、数理生物学者の{{仮リンク|ジェイムズ・D・マレー|en|James D. Murray}}や応用数学者の{{仮リンク|スティーヴン・ストロガッツ|en|Steven Strogatz}}は指摘している{{Sfnm|マレー|2014|1p=3|Strogatz|2015|2p=27}}。

人口予測に関しても、人口学者の{{仮リンク|ジョエル・E・コーエン|en|Joel E. Cohen}}は「ロジスティック曲線は短期的な予測に関しては、他の連続でなめらかな曲線と比べて特に劣っていることもないが、長期的な予測に関しても格別に秀でているわけでもない」と評している{{Sfn|コーエン|1998|p=116より引用}}。式を普及させたレイモンド・パールは、ある期間の人口成長にロジステック曲線が適用できる条件として、人口成長に影響を与える新しい要素がその期間中に現れないことを挙げている。しかし、このような前提条件を人口という複雑な現象に課すのは困難である点を経済学者の{{仮リンク|アルバート・B・ウルフ|en|Albert B. Wolfe}}や人口学者の{{仮リンク|ジョージ・ハンドリー・ニブス|en|George Handley Knibbs}}などから批判されている{{Sfn|Kingsland|1982|p=36}}。2010年代現在、将来人口推計には[[コーホート要因法]]の使用が主流となっている{{Sfnm|人口研究会(編)|2010|1p=74|コーエン|1998|2p=144}}。ロジスティック曲線のような関数を過去の人口データに重ねて将来の人口を予測するという単純な方法は、現在ではほとんど行われていない<ref>{{cite book |和書 |author=山重慎二・加藤久和・小黒一正 |title=人口動態と政策―経済学的アプローチへの招待 |date=2013-09-10 |edition=第1版 |publisher=日本評論社 |isbn =978-4-535-55750-5 |page =28}}</ref>。

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

==名称の由来==
フェルフルストは、1845年の論文で、"Nous donnerons le nom de ''logistique'' à la courbe"(参考訳:私たちはその曲線にロジスティック "''logistique''" という名前を与える)と述べ<ref>{{cite journal|first= Pierre-François |last=Verhulst |year= 1845| title = Recherches mathématiques sur la loi d'accroissement de la population | journal = Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Bruxelles |volume = 18| pages = 1–42 | url = http://gdz.sub.uni-goettingen.de/dms/load/img/?PPN=PPN129323640_0018&DMDID=dmdlog7 }}のp. 8より引用</ref>、ロジスティック方程式の解による曲線を ''logistique'' と名付けた{{Sfn|人口研究会(編)|2010|p=307}}。これが、式が"ロジスティック"方程式、その解曲線が"ロジスティック"曲線と呼ばれる由来である{{Sfnm|人口研究会(編)|2010|1p=307|瀬野|2007|2p=20}}。しかし、フェルフルストは ''logistique'' という語を使った理由を説明しなかったので、それ以上の由来は分かっていない<ref name="MacTutor"/>{{Sfn|瀬野|2007|p=20}}<ref name="大澤2007"/>。

''logistique'' と名付けられた理由のいくつかの推測は存在する。ベルギー王国陸軍士官学校の数学教授のHugo Pastijnは、実際の理由は不明と断った上で、
*陸軍大学に勤めていたフェルフルストも馴染みが有ったであろう「[[兵站]]」の意味と関連付けて ''logistique'' と名付けたのではないか
*フェルフルストのモデルでも扱われる人口のための限られた資源と関連させて、「[[住居]]」を意味するフランス語の ''logis'' から名付けたのではないか
と、ありえそうな理由を2点ほど推測している<ref name="MacTutor">{{Cite web |author=O'Connor, John J; Robertson, Edmund F |url=http://www-history.mcs.st-andrews.ac.uk/Biographies/Verhulst.html |title=Pierre François Verhulst |year=2014 |month=1 |work=MacTutor History of Mathematics archive |publisher=University of St Andrews |accessdate=2015-12-17}} 出典での引用元は、{{Cite book |last=Pastijn |first=Hugo |chapter=Chaotic Growth with the Logistic Model of P.-F. Verhulst |chapterurl=http://link.springer.com/chapter/10.1007/3-540-32023-7_1 |year=2006 |pages=3-11 |title=The logistic map and the route to chaos: From the beginnings to modern applications, Understanding Complex Systems |publisher=Springer Berlin Heidelberg |doi=10.1007/3-540-32023-7_1 |isbn=978-3-540-28366-9}}</ref>。また、19世紀当時のフランスでは、''logistique'' には「計算に巧みな」「計算の技巧」といった意味での用例があった点も指摘されている{{Sfnm|寺本|1997|1p=10|人口研究会(編)|2010|2p=307}}。

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

===捕獲の影響===
人間が資源として利用するための捕獲や収穫は、その種を絶滅させる可能性もあるほどの大きな影響を持っている{{Sfn|寺本|1997|p=22}}。漁業分野では、水産資源を獲りつくさないように資源・漁業管理する必要性が認識されている{{Sfn|北原(編)|2003|pp=1&ndash;2}}。持続可能な漁業のためには、人間による漁獲量が漁獲対象の自然増加量を上回らないようにする必要がある<ref name="海洋保全生態学">{{cite book |和書 |editor=白山義久・桜井泰憲・古谷研・中原裕幸・松田裕之・加々美康彦 |title=海洋保全生態学 |year=2012 |publisher=講談社 |isbn =978-4-06-155231-9 |pages=172-174}}</ref>。漁獲量と自然増加量が一致するとき資源は一定に保たれるので、このときの漁獲量を持続生産量と呼ぶ<ref name="海洋保全生態学"/>。さらに、可能な持続生産量の中でも最大のものを{{仮リンク|最大持続生産量|en|Maximum sustainable yield}}(MSY)と呼び、漁獲基準の一つの目安とされている{{Sfn|北原(編)|2003|p=244}}。

この最大持続生産量の値をロジスティック方程式を利用して定量化するモデルを、ジェーファーのプロダクションモデルなどと呼ぶ<ref name="海洋保全生態学"/>。漁獲量(漁獲速度)を ''Y'' とすれば、次のように、ロジスティック方程式で表される個体数増加率(自然増加率)から ''Y'' を差し引いた値が実際の増加率となる{{Sfn|寺本|1997|p=22}}。
:<math>\frac{dN}{dt}\ = rN \left( 1- \frac{N}{K} \right) -Y</math>
''dN''/''dt'' が 0 のときが資源一定状態なので、このときの漁獲量が持続生産量を示している<ref name="海洋保全生態学"/>。''Y'' を含まないときの ''dN''/''dt'' の最大値は前述のとおり ''rK''/4 である。これに釣り合う漁獲量が最大持続生産量であるから、このモデルでは最大持続生産量を ''rK''/4 と得ることができる{{Sfn|日本数理生物学会(編)|2008|p=187}}。

漁獲量 ''Y'' を単純な一定値とせずに、個体数に比例するようなモデルも考えられる。例えば、出漁する漁船の数が一定とすれば、捕獲の成果は生息している個体数に比例すると考える方が適当である{{Sfn|寺本|1997|p=23}}。''q'' と ''E'' を定数として ''q'' と ''E'' と ''N'' を掛け合わせたもので漁獲量を表せば、個体数増加率は
:<math>\frac{dN}{dt}\ = r N \left( 1- \frac{N}{K} \right)-qEN</math>
となる<ref name="水産ハンドブック"/>。''q'' は漁具効率、''E'' は漁獲努力量と呼ばれる<ref name="水産ハンドブック">{{cite book |和書 |editor= 島一雄・關文威・前田昌調・木村伸吾・佐伯宏樹・桜本和美・末永芳美・長野章・森永勤・八木信行・山中英明 |title=最新 水産ハンドブック |year=2012 |publisher=講談社サイエンティフィク |isbn =978-4-06-153736-1 |page=114}}</ref>。このモデルの場合は、''qE'' を内的自然増加率の半分 ''r''/2 となるようにすれば、漁獲量を最大持続生産量にすることができる<ref name="海洋保全生態学"/>。ただし、以上のモデルは現実をかなり単純化したモデルである。環境の変化や他の生物との相互作用など、現実には様々な要因が関係しているため、多数の相互作用がある実際の生態系では成り立たない{{Sfn|日本数理生物学会(編)|2008|p=188}}。実際の最大持続生産量の決定にはより高度な手法も使用されている<ref name="海洋保全生態学"/>。

===2種存在する場合===
[[File:Comportements competition lv 2 sp.svg|thumb|250px|ロトカ・ヴォルテラの競争式では係数の値がある範囲内のときのみ2種が共存し(図の3)、それ以外ではどちらかが絶滅する(図の1, 2, 4){{Sfn|ハーバーマン|1992|pp=130&ndash;136}}。]]
{{see also|競争 (生物)|ロトカ・ヴォルテラの競争方程式}}
ロジスティック方程式は環境内に1種のみが存在するときの(あるいは1種とみなせるときの)モデルだが、実際の環境では複数以上の種が生息している{{Sfn|巌佐|1990|p=13}}。複数の種が存在するとき、それぞれの種の間には[[競争 (生物)|競争]]や[[相利共生]]、[[捕食-被食関係|捕食-被食]]などの関係が存在して、それぞれの個体数が互いの個体数増加率に影響を与える{{Sfn|マレー|2014|p=65}}。その中でも特に、環境内に競争関係にある2種が存在する場合にロジスティック方程式を拡張させたものとして、以下の[[ロトカ・ヴォルテラの競争方程式]]が知られる{{Sfn|大串|2014|p=66}}。
:<math>\frac{dN_1}{dt} = r_1 N_1 \left(1 - \frac{N_1 + a_{12} N_2}{K_1} \right)</math>
:<math>\frac{dN_2}{dt} = r_2 N_2 \left(1 - \frac{N_2 + a_{21} N_1}{K_2} \right)</math>

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

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

===時間遅れの考慮===
[[File:Lucilia cuprina (15848887484).jpg|thumb|[[ヒツジキンバエ]]。[[ロバート・メイ]]が時間遅れを持つロジスティック方程式でこのハエの個体数変動の解析を行った。]]
ロジスティック方程式では、ある時刻の個体数 ''N''(''t'') が同時刻の個体数増加率 ''dN''(''t'')/''dt'' に瞬間的に影響を与えるというモデルになっている{{Sfn|日本生態学会(編)|2004|p=125}}。しかし、妊娠期間や性成熟までの期間などが存在するため瞬間的に影響が出るというのは非現実的でもある{{Sfn|マレー|2014|p=12}}。よって、モデルの中に影響の時間遅れを含ませることが考えられる。遅延時間を ''T'' とすると、ロジスティック方程式に時間遅れの効果を取り込んだモデルとして
:<math>\frac{dN(t)}{dt} = r N(t) \left(1 - \frac{N(t-T)}{K} \right)</math>
がよく用いられる{{Sfn|マレー|2014|p=12}}。この式は[[ジョージ・イヴリン・ハッチンソン]]が発案したためHutchinson方程式とも呼ばれる<ref>{{cite book |editor=Arino, O.; Hbid, M.L.; Dads, E. Ait |author=Ruan, S. |title=Delay Differential Equations and Applications: Proceedings of the NATO Advanced Study Institute held in Marrakech, Morocco, 9-21 September 2002 |publisher=Springer Netherlands |chapter=DELAY DIFFERENTIAL EQUATIONS IN SINGLE SPECIES DYNAMICS |year=2006 |edition=1 |doi=10.1007/1-4020-3647-7 |pages=479}}</ref>{{Sfn|日本数理生物学会(編)|2008|p=67}}。このモデルでは、ロジスティック方程式におけるブレーキ効果の部分に、現時点での個体数 ''N''(''t'') ではなく、時間 ''T'' だけ前の時点での個体数 ''N''(''t'' &minus; ''T'') が入力されている。

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

===離散時間モデル===
[[File:Discrete logistic equation-time evolutions.svg|thumb|280px|離散時間型モデルの場合の個体数の変化の様子。いずれも ''K'' と ''N''<sub>0</sub> の値は同じだが、青が ''r'' = 0.6、赤が ''r'' = 2.1、緑が ''r'' = 3 のときを示している。]]
{{see also|ロジスティック写像}}
ロジスティック方程式では、時間 ''t'' を連続な実数として個体数変動をモデル化した。しかし、世代の交代が同期的に起こり、世代の重なりがないようなときには、時間を飛び飛びの時間間隔(離散時間)でモデル化する方が妥当である{{Sfn|寺本|1997|p=11}}。ロジスティック方程式型の離散時間モデルにはいくつかの種類があるが、一例として次のような[[差分方程式]]がある{{Sfn|巌佐|1990|p=49}}。
:<math> N_{n+1}=N_{n}+rN_n \left( 1-\frac{N_n}{K} \right) </math>
ここで、''n'' は世代で、''n'' = 1世代, 2世代, 3世代,... といったような飛び飛びの時間間隔を意味している。''N<sub>n</sub>'' は、''n'' 世代における個体数 ''N'' を意味している{{Sfn|巌佐|1990|p=49}}。上式と数学的には等価だが、[[ロジスティック写像]]と呼ばれる、次の形式での差分方程式もよく知られている<ref name="ウィロックス"/>。
:<math> x_{n+1}=ax_n(1-x_n) </math>
これらの差分方程式はロジスティック方程式と一見似ているが、解の様相は全く異なり、個体数の変動はロジスティック方程式よりも遥かに複雑な振る舞いを見せるようになる{{Sfn|Hirsch et al.|2007|p=347}}。''r''(または ''a'')が小さい内は、これらの解はロジスティック方程式と同じように安定な平衡状態に収束する{{Sfnm|巌佐|1990|1p=51|マレー|2014|2p=44}}。''r'' が大きくなってくると、個体数は多くなったり少なくなったりを交互に繰り返すようになる。さらに ''r'' が大きくなると、[[カオス理論|カオス]]と呼ばれる非周期的で極めて複雑な振る舞いを起こすようになる{{Sfnm|巌佐|1990|1p=53|マレー|2014|2pp=44&ndash;47}}。

また、京都大学の[[森下正明]]が発案した次のような差分方程式がある<ref>{{cite journal |author=Masaaki Morisita |title=The fitting of the logistic equation to the rate of increase of population density |publisher=Springer-Verlag |journal=Researches on Population Ecology |volume= 7 |issue 1 |year=1965 |month=6 |doi=10.1007/BF02518815 |pages=52&ndash;55}}</ref>。
:<math> N_{n+1}=\frac{(1+a)N_n}{1+bN_n} </math>
ここで、''&Delta;t'' を差分時間間隔として、''a'' と ''b'' は
:<math> a=e^{r \Delta t}-1 </math>
:<math> b=\frac{a}{K} </math>
である。通常、差分化を行うと元の方程式の解と誤差が生じる{{Sfn|山口|1992|p=72}}。しかしこの方程式では誤差を全く生じさせない<ref name="ウィロックス">{{cite journal |和書 |author=ウィロックス, ラルフ |title=自然現象の離散・超離散系によるモデル化 |publisher=九州大学応用力学研究所 |journal=応用力学研究所研究集会報告 No.22AO-S8 |year=2011 |month=3 |doi=10.15017/23387 |page=3}}</ref>。得られる解は離散的だが、その解はロジスティック方程式の解と一致し、解を ''N''-''t'' 平面上に描けば、ロジスティック曲線上に正確にプロットされる{{Sfn|山口|1992|p=73}}。

===生物個体数以外での例===
本来の導入目的であった生物の個体数の変動以外にも、ロジスティックモデルがしばし使用される。興味対象の何かの変量が時間発展とともにS字型の曲線を描くようなときに、この式がよく当てはまる{{Sfn|巌佐|1990|p=5}}。水産資源管理の例では、生物の体の大きさの成長曲線にロジスティック曲線を当てはめることがある{{Sfn|北原(編)|2003|pp=220&ndash;221}}。また、人間の集団の中で無形なものが広まる様子を表すのにもロジスティックモデルが使われることがある。例えば、新技術の社会・産業全体への普及{{Sfn|ブラウン|2012|pp=45&ndash;46}}<ref name="watanabe"/>、ある集団の中での[[噂]]の拡散<ref>{{cite book |和書 |author=橋本洋志・石井千春 |title=Scilab/Scicosで学ぶシミュレーションの基礎―自然・社会現象から、経済・金融、システム制御まで |publisher=オーム社 |year=2008 |isbn=978-4-274-20487-6 |pages=103&ndash;105}}</ref>がある。

また、時間変化ではないが、統計学においてはロジスティック関数と同形式の[[累積分布関数]] ''f''(''x'') を持つ[[連続確率分布]]が用いられている。これを[[ロジスティック分布]]と呼ぶ<ref>{{cite book| 和書 |author=白石高章 |title=統計科学の基礎―データと確率の結びつきがよくわかる数理 |year=2012 |edition=第1版 |publisher=日本評論社 |isbn =978-4-535-78700-1 |page=159}}</ref>。[[人工ニューラルネットワーク]]の研究で使われる[[シグモイド関数]]の一つとしてもロジステック関数が利用されている<ref>{{cite book |和書 |author=岡谷貴之 |title=深層学習 |publisher=講談社サイエンティフィク |series=機械学習プロフェッショナルシリーズ |date=2015-04-07 |isbn=978-4-06-152902-1 |page=10}}</ref>。

==他形式==
上記ではロジスティック方程式を
:<math>\frac{dN}{dt}\ = rN \left( 1-\frac{N}{K} \right) </math>
と表したが、これ以外の表現もある。いずれも数学的には等価だが、その導出過程における生態学的意味づけは様々である{{Sfn|瀬野|2007|pp=i, 25&ndash;26}}。 ''k'' = ''r'' / ''K'' と置いて、ロジスティック方程式は
:<math>\frac{dN}{dt}\ = N(r - kN)</math>
とも表される{{Sfnm|瀬野|2007|1p=22|ブラウン|2012|2p=32|ハーバーマン|1992|3p=41}}。''k'' は'''Verhulst-Pearl係数'''や'''種内競争係数'''と呼ばれる{{Sfn|寺本|1997|p=9}}。個体群密度が増加率を減少させる影響の強さを ''k'' が表しているといえる{{Sfn|瀬野|2007|p=21}}。

他には、変数を ''N'' = ''N''/''K'' と置きなおして、すなわち個体数ではなく環境収容力に対する個体数の割合を変数として
:<math>\frac{dN}{dt}\ = r N(1-N)</math>
という形式もある{{Sfn|Hirsch et al.|2007|p=4}}。

[[非線形]]のロジスティック関数を扱いやすくするために[[線型性|線形]]の対数関数に変換する、フィッシャ・プライ変換(英語:Fisher-Pry transform)と呼ばれる次のような変換もある<ref name="watanabe">{{cite book|和書 |author=渡辺千仭 |title=技術経済システム |publisher=創成社 |year=2007 |edition=初版 |isbn=978-4-7944-3089-2 |pages=84&ndash;87}}</ref>。
:<math>FP=\frac{1}{1+e^{-bt-a}} </math>
:<math>bt + a = \ln \frac{FP}{1-FP}</math>
ここで ''FP'' = ''N'' とすると、ロジスティック関数のパラメータとの関係は ''K'' = 1, ''r'' = ''b'', ''N''<sub>0</sub> = ''e<sup>a</sup>''/(1 + ''e<sup>a</sup>'') である。


== 歴史 ==
== 歴史 ==
===フェルフルストによる発表===
この式を最初に発表したのは[[ピエール=フランソワ・フェルフルスト]]である。彼は1838年に始まる数本の論文で人口増加について論じた中で、それを表す式としてこれを提案した。当時はこの価値を認めるものはほとんどなく、彼の死亡時の告知にも、彼の業績として取り上げられなかった。
[[File:Pierre Francois Verhulst.jpg|thumb|upright|[[ピエール=フランソワ・フェルフルスト]] (''Pierre-François Verhulst'')]]
[[ベルギー]]・[[ブリュッセル]]の陸軍大学の数学者であった[[ピエール=フランソワ・フェルフルスト]]によって、ロジスティック方程式は発表された{{Sfn|コーエン|1998|pp=112&ndash;113}}。18世紀になると、[[トマス・ロバート・マルサス]]が出版した『[[人口論]]』に関心が高まっていた{{Sfn|山口|1992|p=54}}。[[マルサスモデル]]の説明で述べたように、マルサスは人口が指数関数的に成長していくモデルを発表し、その帰結として社会が飢饉の発生など破滅的状況を迎えることを予測した{{Sfn|人口研究会(編)|2010|pp=280&ndash;282}}。このセンセーショナルな予測は衝撃を与え、当時およびマルサス死後も長く続く論争を引き起こした{{Sfn|人口研究会(編)|2010|p=282}}。「近代統計学の父」と呼ばれる[[アドルフ・ケトレー]]も、マルサスのモデルに関心を持ち、人口増減モデルについて論じた{{Sfn|人口研究会(編)|2010|p=315}}。ケトレーは[[抗力|流体の抵抗]]をヒントにして、人口増加率への抵抗は人口増加率自体の二乗に比例すると考えた{{Sfnm|山口|1992|1p=55|寺本|1997|2p=10}}。


ケトレーから教えを受けたこともあり、友人でもあったフェルフルストは、ケトレー自身からケトレーのモデルに関する研究を勧められた{{Sfn|瀬野|2007|p=20}}。ケトレーの考えをもとにして、人口が人口自体によって増加する一方で、人口増加を抑制する何らかの機構が働く数学的なモデルを思案した{{Sfn|人口研究会(編)|2010|p=307}}。1838年、フェルフルストは、"''Notice sur la loi que la population poursuit dans son accroissement''"という題で研究成果を発表し、この論文の中でロジスティック方程式が提案された{{Sfn|山口|1992|p=56}}。この論文の中でフェルフルストが実際に提案した式は、
その後1910年代から生物の個体群成長に関する実験などが行われる中で、この式は独自にあちこちで使われ始めたが、フェルフルストの名が挙がることはなかった。1920年、[[パール]]らがアメリカ合衆国の人口増加について論じ、[[ショウジョウバエ]]の実験個体群の成長を研究したとき、やはりこの式を使い、翌年にこれがすでに90年近く前にフェルフルストによって発見されたことを認めた。これによって、やっと彼の名がこの式に結びついた。
:<math>\frac{dp}{dt}\ =mp-\phi(p) </math>
:<math>\phi(p)=np^2 </math>
という形であった<ref name="Verhulst1838">{{cite journal |last=Verhulst |first=Pierre-François |year=1838 |title=Notice sur la loi que la population suit dans son accroissement |url=https://books.google.co.jp/books?id=8GsEAAAAYAAJ&pg=PA113&hl=ja&source=gbs_toc_r&cad=4#v=onepage&q&f=false |journal=Correspondance mathématique et physique |volume=10 |pages=113-121 }}のp.115</ref>{{Sfn|瀬野|2007|p=18}}。''p'' は人口を意味する<ref name="Verhulst1838"/>。フェルフルストは人口自体の二乗によって人口増加率の減少効果を表現し、上記の ''φ''(''p'') を導入した{{Sfn|瀬野|2007|p=18}}。当時はこの式の価値を認めるものはほとんどなく、彼の死亡時の告知にも、彼の業績として取り上げられなかった{{Sfn|山口|1992|p=57}}。


===式の再発見と論争と普及===
上記のように、一般の生物に当てはめるには難しい点もあるが、これ以降、実験室や野外での生物の個体数変動を扱う基礎モデルとして、この式は広く認められるようになった。
[[File:Raymond Pearl (cropped).jpg|thumb|right|upright|レイモンド・パール (Raymond Pearl)]]
フェルフルスト発表の後、生物の個体群成長に関する実験などで、同じ式が独自にあちこちで使われ始めたが、フェルフルストの名が挙がることはなかった{{Sfn|山口|1992|p=58}}。1908年には生理学者の{{仮リンク|ブレイルスフォード・ロバートソン|en|Thorburn Brailsford Robertson}}が、動物、植物、人間といった生物の個体成長を同形式の曲線で記述した{{Sfn|Kingsland|1982|p=31}}。ロバートソンはフェルフルストの発表を知らなかったが、同じ曲線を用い、さらに偶然にもケトレーのデータを使用している{{Sfn|山口|1992|p=58}}。この時点で同じ曲線が化学物質における[[自己触媒反応]]の過程を表すのに使われていたので、ロバートソンは曲線のことを自己触媒的曲線(autocatalytic curve)と呼んでいた{{Sfn|寺本|1997|p=10}}。


1920年、[[ジョンズ・ホプキンス大学]]の[[レイモンド・パール]]と{{仮リンク|ローウェル・リード|en|Lowell Reed}}が、ロジスティック方程式と同形式のモデルを用いてアメリカ合衆国の人口増加について論じた<ref name="大澤2007"/>。この研究も、フェルフルストにより先に発表されていたことを知らずに行われた{{Sfn|コーエン|1998|p=113}}。翌年の1921年には、これがすでに80年近く前にフェルフルストによって発見されたことをパールらも認めた{{Sfn|山口|1992|p=58}}。これによってパールらもロジスティック曲線という名称を使うようになり、やっとフェルフルストの名がこの式に結びつくことになる{{Sfn|山口|1992|p=58}}。これ以降、生物学ではロジスティック曲線という名称が定着した{{Sfn|寺本|1997|p=10}}。
また、ロジスティック方程式における''r'' はその種が実現できる最大の相対増加率であり、これが大きい方が素早く増殖できる可能性がある。また、''K'' はその環境下で生存できる個体数上限を示す。[[島嶼生物学]]の分野で、マッカーサーとウィルソンは[[島]]における生物個体群の定着と[[絶滅]]を論じ、定着の成功には大きな''r'' を持つことが重要であり、絶滅の回避には大きな''K'' を持つことが重要であるとし、それぞれを'''r淘汰'''、'''K淘汰'''と呼んだ。これが[[r-K戦略説]]、ひいては[[生活史戦略]]論の始まりとなった。


パールは[[ショウジョウバエ]]の個体群成長の実験を行い、この式を実証した{{Sfn|内田|1972|pp=28&ndash;30}}。1924年と1925年にも、アメリカ、スウェーデン、フランスなどの様々な国勢調査の人口統計にロジスティック曲線のあてはめを行い、よく一致することを示した{{Sfnm|マレー|2014|1pp=4&ndash;5|コーエン|1998|2p=113}}。このような積み重ねた証拠をもとに、パールは個体群全般がロジステック曲線に沿って成長することを強く確信し、ロジステック曲線が「法則」であると主張した{{Sfn|Kingsland|1982|pp=33&ndash;34}}。当時、パールとリードはこの式の価値を「控え目にいっても、それは[[ケプラーの法則|ケプラーの惑星の楕円運動法則]]に匹敵するものであるといってもよいように思われる」と自身らで評価している{{Sfn|コーエン|1998|pp=113&ndash;114より引用。出典での引用元は、Pearl, Raymond, and Lowell J. Reed. 1924 The growth of human population. In ''Studies in human biology'', ed. Raymond Pearl. Baltimore: Williams and Wilkins, pp. 584-637.のp.585より}}。ロジスティック曲線は、経験的なものというよりも、個体群成長全般において普遍性を持つ法則であり、成長の長期的傾向の予測も可能にするとパールは考えていた{{Sfn|Kingsland|1982|pp=32, 35, 42}}。パールは、この式が個体群成長における普遍則であるという持論を広め、ロジスティック方程式の普及に大きく貢献することになる{{Sfnm|マレー|2014|1pp=3&ndash;4|人口研究会(編)|2010|2p=323}}。このため、ロジスティック曲線にはパールの名が題されることもある{{Sfn|コーエン|1998|p=115}}。
また、ロジスティック方程式を[[差分方程式]]にすると、''K'' の値の取り方次第で、個体数は''N'' &rarr;''K'' に安定する場合もあるが、''K'' の上下2つの値の間を行き来したり、あるいは4つの値の間を行き来する場合もある。[[内田俊郎]]らによる実験室内での昆虫個体群の研究によると、この現象には実例があり、その原因は個体数の増加が増加率に影響する時間差である。なお、[[ロバート・メイ]]はこの式をさらに追求して、非周期的にあらゆる値をとる場合にまでいたるさまざまな形が出現することをコンピュータ・[[シミュレーション]]によって示し、これに対して'''カオス的''' (chaotic) という言葉を当てたのが[[カオス理論]]の始まりの一つである。


一方で、パールの自説の展開には多くの批判も呼び、1940年にパールが死去するまで論争が続いた{{Sfn|Kingsland|1982|p=32}}。経済学者のA. B. ウルフ、人口学者の{{仮リンク|ジョージ・ハンドリー・ニブス|en|George Handley Knibbs}}、統計学者の{{仮リンク|エドウィン・ウィルソン|en|Edwin Bidwell Wilson}}などがパールの主張に批判を加えている{{Sfn|Kingsland|1982|pp=35&ndash;36, 38&ndash;39}}。動物学者の{{仮リンク|ジェームズ・グレイ (動物学者)|label=ジェームズ・グレイ|en|James Gray (zoologist)}}、[[ランスロット・ホグベン]]、遺伝学者の[[シューアル・ライト]]からは、他のS字型曲線を使っても個体群成長のデータに当てはめができるので、ロジスティック曲線を使う必然性が欠けていることについて指摘を受けている{{Sfn|Kingsland|1982|p=40}}。
==二重ロジスティック関数==

定義は以下の通り。[[ガウス関数]]を変形した物である。
ロジスティック曲線の有効性を支持し、その普及を担った人たちもいる。イギリスの統計学者{{仮リンク|ウドニー・ユール|en|Udny Yule}}は、パールの理論を1924年のイギリスの学会で発表している。ユールはロジスティック曲線は長期予測には適用できないと考えており、その点を強調したが、基本的にはパールの研究を支持していた{{Sfn|Kingsland|1982|p=34}}。[[アルフレッド・ロトカ]]も、ロジスティック曲線の有効性を理解し、ロジスティック方程式について1925年の自書の中で一章を与えて説明した{{Sfn|Kingsland|1982|pp=41&ndash;42}}。ただしロトカは、ロジスティック方程式は実現象の近似の一種であるという考えを保っていた{{Sfn|Kingsland|1982|p=43}}。ロシアの{{仮リンク|ゲオルギー・ガウゼ|en|Georgy Gause}}も、近似の一種と受け止めながらも、ロジスティック方程式が同種の集団の中での生存競争を定量的に表すことができると述べている{{Sfn|Kingsland|1982|pp=45&ndash;46}}。1934年、ガウゼは微生物の実験によってロジスティック方程式の検証を行い、この検証はロジスティックモデルを[[個体群動態論]]における古典的理論の一つとして確固たるものとした{{Sfn|日本生態学会(編)|2004|p=241}}。モデルの限界には注意が払われながらもロジスティック方程式の受容は広まっていき、1940年代後半には個体群解析における一般的な道具として確立した{{Sfn|Kingsland|1982|p=49}}。
:<math> y = \mbox{sgn}(x-d) \, \left[1-\exp\left\{-\left(\frac{x-d}{s}\right)^2\right\}\right]</math>

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

ロジスティック方程式における ''r'' は個体群密度がとても低いときの増加率で表しており、密度が低いときにどれだけ素早く繁殖できるかを意味している{{Sfn|日本生態学会(編)|2004|p=62}}。また、''K'' はその環境下で生存できる個体数あるいは個体群密度の上限を示す{{Sfnm|日本生態学会(編)|2004|1p=62|レーヴンほか|2007|2p=1150}}。1967年、[[ロバート・マッカーサー]]と[[エドワード・オズボーン・ウィルソン]]は、この ''r'' と ''K'' に着目して、[[島]]における生物個体群の定着と[[絶滅]]に関する理論を発案した{{Sfn|木元|1979|pp=108&ndash;109}}。彼らの理論によれば、ある生物の島への定着が成功するには大きな ''r'' を持つことが重要であり、絶滅の回避には大きな ''K'' を持つことが重要であるとし、それぞれの方向へ[[選択 (進化)|淘汰]]されることを '''r淘汰'''、'''K淘汰'''と呼んだ{{Sfn|木元|1979|pp=116&ndash;117}}。この説は[[r-K戦略説]]と呼ばれ、生物の[[生活史 (生物)|生活史]]の進化に種内競争の観点から説明を与えた{{Sfn|日本生態学会(編)|2004|pp=61, 64}}。

物理学から数理生態学へ転向してきた[[ロバート・メイ]]も個体群変動の問題に取り組んだ{{Sfn|山口|1992|p=76}}。メイはロジスティック方程式の離散化を行い、その式の解は、通常のロジスティック方程式の解とは全く異なる、現在では[[カオス理論|カオス]]と呼ばれる非常に複雑な振る舞いを起こすことを示した{{Sfn|スチュアート|2012|pp=342&ndash;343}}。この結果は1974年と1975年に発表され、大きな反響を得ると共に、その後の[[カオス理論]]の隆盛に大きく寄与することになる{{Sfn|スチュアート|2012|p=342}}<ref>{{cite book |和書 |title=カオスはこうして発見された |author1=ティェンイェン・リー |author2=ジェームス・A・ヨーク |others =ラルフ・エイブラハム、ヨシスケ・ウエダ (編) 稲垣耕作、赤松則男(訳) |publisher=共立出版 |year=2002 |edition=初版 |isbn= 4-320-03418-X |chapter=第10章 区間上のカオスを探索する |pages= 169-170}}</ref>。


== 脚注 ==
== 脚注 ==
=== 注釈 ===
{{reflist}}
{{notelist}}
=== 出典 ===
{{Reflist|2}}


== 参考文献 ==
== 参考文献 ==
※文献内の複数個所に亘って参照したものを示す。
*[[山口昌哉]]、『カオスとフラクタル』(1986)、[[講談社]]([[ブルーバックス]])、ISBN 4-06-132652-X
*{{cite book ja-jp
*[[内田俊郎]]、『動物の人口論』(1972)、[[NHKブックス]]([[NHK出版|日本放送出版協会]])、ISBN 9784140011645
|author= 山口 昌哉
|title=カオスとフラクタル―非線形の不思議
|url = https://bookclub.kodansha.co.jp/product?item=0000138469
|series=ブルーバックス
|publisher= 講談社
|year=1992
|edition=第17刷
|isbn=4-06-132652-X
|ref={{Sfnref|山口|1992}}
}}
*{{cite book ja-jp
|author= 内田 俊郎
|title=動物の人口論―過密・過疎の生態をみる
|series=NHKブックス164
|publisher=日本放送出版協会
|year=1972
|id=1345-001164-6023
|isbn=9784140011645
|ref={{Sfnref|内田|1972}}
}}
*{{cite book ja-jp
|author= 木元 新作
|title= 南の島の生きものたち―島の生物地理学
|url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320006959
|year=1979
|edition=初版
|publisher=共立出版
|id=1345-472380-1371
|isbn =978-4320006959
|series=科学ブックス38
|ref={{Sfnref|木元|1979}}
}}
*{{Cite book ja-jp
|author= 瀬野 裕美
|title=数理生物学―個体群動態の数理モデリング入門
|url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320056565
|year=2007
|edition=初版
|publisher=共立出版
|isbn =978-4-320-05656-5
|ref ={{Sfnref|瀬野|2007}}
}}
*{{cite book ja-jp
|author=巌佐 庸
|title=数理生物学入門―生物社会のダイナミックスを探る
|year=1990
|edition=初版
|publisher=HBJ出版局
|isbn =4-8337-6011-8
|ref ={{Sfnref|巌佐|1990}}
}}
*{{cite book ja-jp
|editor= 日本生態学会
|others=巌佐 庸・舘田 英典(担当編集委員)
|author=巌佐 庸
|chapter = 第2章 人口増殖と環境収容力
|pages=17&ndash;27
|title= 集団生物学
|url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320057449
|series=シリーズ 現代の生態学 1
|year=2015
|edition= 初版
|publisher=共立出版
|isbn =978-4-320-05744-9
|ref ={{Sfnref|巌佐|2015}}
}}
*{{cite book ja-jp
|author=寺本 英
|editor =川崎 廣吉・重定 南奈子・中島 久男・東 正彦・山村 則男
|title= 数理生態学
|year=1997
|edition=初版
|publisher=朝倉書店
|isbn =4-254-17100-5
|ref ={{Sfnref|寺本|1997}}
}}
*{{cite book ja-jp
|author= 大串 隆之
|chapter = 3章 昆虫の個体群と群集
|pages=49&ndash;98
|title= 昆虫生態学
|url = https://www.asakura.co.jp/detail.php?book_code=42039
|year=2014
|edition=初版
|publisher=朝倉書店
|isbn =978-4-254-42039-5
|ref ={{Sfnref|大串|2014}}
}}
*{{cite book ja-jp
|author=渡辺 守
|title= 生態学のレッスン ―身近な言葉から学ぶ
|url = http://www.utp.or.jp/book/b306454.html
|publisher= 東京大学出版会
|year=2012
|isbn=978-4-13-063334-5
|ref ={{Sfnref|渡辺|2012}}
}}
*{{cite book ja-jp
|author=ジョエル・E・コーエン
|title=新「人口論」―生態学的アプローチ
|url = https://shop.ruralnet.or.jp/b_no=01_4540970569/
|translator = 重定 南奈子・瀬野 裕美・高須 夫悟
|publisher= 農山漁村文化協会
|year=1998
|edition=初版
|isbn=4-540-97056-9
|ref={{Sfnref|コーエン|1998}}
}}
*{{cite book ja-jp
|others =三村 昌泰(総監修)、瀬野 裕美ほか(監修)
|author=ジェームス・D・マレー
|translator =勝瀬 一登・吉田 雄紀・青木 修一郎・宮嶋 望・半田 剛久・山下 博司
|title=マレー数理生物学入門
|year=2014
|edition=初版
|publisher=丸善出版
|isbn =978-4-621-08674-2
|ref ={{Sfnref|マレー|2014}}
}}
*{{cite book ja-jp
|author=イアン・スチュアート
|translator =水谷 淳
|title=数学で生命の謎を解く
|url = https://www.sbcr.jp/product/4797369694/
|year=2012
|edition=初版
|publisher=ソフトバンククリエイティブ
|isbn =978-4-7973-6969-4
|ref ={{Sfnref|スチュアート|2012}}
}}
*{{cite book ja-jp
|author=ホルスト R. ティーメ
|others =斉藤 保久(監訳)
|title= 生物集団の数学(上)―人口学,生態学,疫学へのアプローチ
|url = https://www.nippyo.co.jp/shop/book/2932.html
|year=2006
|edition=第1版
|publisher=日本評論社
|isbn =4-535-78418-3
|ref ={{Sfnref|ティーメ|2006}}
}}
*{{cite book ja-jp
|author=R. ハーバーマン
|translator =稲垣 宣生
|title= 生態系の微分方程式
|url = https://www.gensu.jp/product/%e7%94%9f%e6%85%8b%e7%b3%bb%e3%81%ae%e5%be%ae%e5%88%86%e6%96%b9%e7%a8%8b%e5%bc%8f%e2%80%95%e5%80%8b%e4%bd%93%e7%be%a4%e6%88%90%e9%95%b7%e3%81%ae%e6%95%b0%e5%ad%a6%e3%83%a2%e3%83%87%e3%83%ab/
|year=1992
|edition=初版
|publisher= 現代数学社
|isbn =4-7687-0307-0
|ref ={{Sfnref|ハーバーマン|1992}}
}}
*{{cite book ja-jp
|author=M. ブラウン
|title= 微分方程式 上―その数学と応用
|translator = 一樂 重雄・河原 正治・河原 雅子・一樂 祥子
|url = https://www.maruzen-publishing.co.jp/item/b294278.html
|editor=シュプリンガー・ジャパン
|publisher=丸善出版
|year=2012
|isbn=978-4-621-06196-1
|ref={{Sfnref|ブラウン|2012}}
}}
*{{cite book ja-jp
|author=P. レーヴン; G. ジョンソン; J. ロソス; S. シンガー
|others =R/J Biology 翻訳委員会(監訳)
|title= レーヴン・ジョンソン 生物学 下(原書第7版)
|year=2007
|edition=第7版
|publisher= 培風館
|isbn =978-4-563-07797-6
|ref={{Sfnref|レーヴンほか|2007}}
}}
*{{cite book ja-jp
|author=Steven H. Strogatz
|translator=田中 久陽・中尾 裕也・千葉 逸人
|title =ストロガッツ 非線形ダイナミクスとカオス―数学的基礎から物理・生物・化学・工学への応用まで
|url = https://www.maruzen-publishing.co.jp/item/b294857.html
|publisher= 丸善出版
|year=2015
|isbn=978-4-621-08580-6
|ref={{Sfnref|Strogatz|2015}}
}}
*{{cite book ja-jp
|author=Morris W. Hirsch; Stephen Smale; Robert L. Devaney
|translator=桐木 紳・三波 篤朗・谷川 清隆・辻井 正人
|title=力学系入門 原著第3版―微分方程式からカオスまで
|url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320018471
|publisher= 共立出版
|edition=3版
|year=2007
|isbn=9784320111363
|ref={{Sfnref|Hirsch et al.|2007}}
}}
*{{cite book ja-jp
|others=瀬野 裕美(責任編集)
|editor=日本数理生物学会
|title=「数」の数理生物学
|url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320056756
|publisher=共立出版
|series=シリーズ 数理生物学要論 巻1
|year=2008
|edition=初版
|isbn=978-4-320-05675-6
|ref={{Sfnref|日本数理生物学会(編)|2008}}
}}
*{{cite book ja-jp
| editor= 人口研究会
| year = 2010
| title = 現代人口辞典
| url = http://www.harashobo.co.jp/book/b369510.html
| location = 東京都
| publisher = 原書房
| edition = 初版
| isbn = 978-4-562-09140-9
|ref ={{Sfnref|人口研究会(編)|2010}}
}}
*{{cite book ja-jp
|editor= 日本生態学会
|title= 生態学入門
|year=2004
|edition=初版
|publisher= 東京化学同人
|isbn =4-8079-0598-8
|ref ={{Sfnref|日本生態学会(編)|2004}}
}}
*{{cite book ja-jp
|editor=北原 武
|title= 水産資源管理学
|publisher=成山堂書店
|year=2003
|isbn=4-425-82991-3
|ref ={{Sfnref|北原(編)|2003}}
}}
*{{cite journal
|author= Sharon Kingsland
|title=The Refractory Model: The Logistic Curve and the History of Population Ecology
|year=1982
|month=March
|journal=The Quarterly Review of Biology
|volume=57
|issue=1
|pages=29&ndash;52
|publisher=The University of Chicago Press
|doi =10.1086/412574
|jstor =2825134
|ref ={{Sfnref|Kingsland|1982}}
}}


== 関連項目 ==
== 関連項目 ==
* [[フェルミ分布関数]]
*[[個体群動態学]]
* [[増殖曲線]]
*[[競争 (生物)|競争]]

*[[人口動態学]]
==外部リンク==
*[[ロトカ=ヴォルテラの方程式]]
{{Commonscat|Logistic functions|ロジスティック関数}}
*フェルフルストの原論文
**{{cite journal |last=Verhulst |first=Pierre-François |year=1838 |title=Notice sur la loi que la population suit dans son accroissement |url=https://books.google.co.jp/books?id=8GsEAAAAYAAJ&pg=PA113&hl=ja&source=gbs_toc_r&cad=4#v=onepage&q&f=false |journal=Correspondance mathématique et physique |volume=10 |pages=113-121 }} - [[Google ブックス]]
**{{cite journal|first= Pierre-François |last=Verhulst |year= 1845| title = Recherches mathématiques sur la loi d'accroissement de la population | journal = Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Bruxelles |volume = 18| pages = 1–42 | url = http://gdz.sub.uni-goettingen.de/dms/load/img/?PPN=PPN129323640_0018&DMDID=dmdlog7}} - Göttinger Digitalisierungszentrum
*{{Kotobank|ロジスティック方程式|2=}}
*{{Kotobank|ロジスティック曲線|2=}}
*{{MathWorld|title=Logistic Equation|urlname=LogisticEquation}}


{{Good article}}
{{デフォルトソート:ろしすていつくしき}}
{{デフォルトソート:ろしすていつくほうていしき}}
[[Category:生態学]]
[[Category:複雑系]]
[[Category:個体群生態学]]
[[Category:力学系]]
[[Category:力学系]]
[[Category:統計学]]
[[Category:関数]]
[[Category:カオス理論]]
[[Category:カオス理論]]
[[Category:数理生物学]]
[[Category:数理生物学]]
[[Category:人口統計学]]
{{Math-stub}}
[[Category:微分方程式]]

2024年7月19日 (金) 21:09時点における最新版

ロジスティック方程式の解曲線(ロジスティック曲線)の一例。S字の形を描き、環境収容力に収束する。
培養容器内のキイロショウジョウバエ。ロジスティック曲線に当てはまる個体数増加が確認された例である。

ロジスティック方程式(ロジスティックほうていしき、英語:logistic equation[1])は、生物の個体数の変化の様子を表す数理モデルの一種である。ある単一種の生物が一定環境内で増殖するようなときに、その生物の個体数(個体群サイズ)の変動を予測できる。人間の場合でいえば、人口の変動を表すモデルである。

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

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

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

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

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

生物の個体数のモデル

[編集]
フィボナッチによるウサギのつがいの増殖問題

生物の個体数の変動については古くから興味を持たれ、研究が行われてきた[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 とすれば、個体数の増加率は差し引きした bd に個体数 N を掛け合わせた値となる[23]。よって個体数増加率 dN/dt

という微分方程式で表される[24]。ここで m は比例定数であり、m = bd である[20]

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

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

ロジスティック方程式

[編集]
個体数 N と一個体当たり個体群増加率 m の関係

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

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

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

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

ロジスティック効果

[編集]

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

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

個体数と増加率の関係

[編集]
縦軸が dN/dt、横軸が N のグラフ。ロジスティック方程式における、dN/dtN の関係が示されている。

ロジスティック方程式における個体数増加率 dN/dt と個体数 N の関係に着目すれば、この関係は初等教育でも習う二次関数そのものとなっており、dN/dtNグラフ放物線を描く[45]。方程式を解析的に解かなくとも、NdN/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]

ここで eネイピア数である。分母・分子を N0ert で割り、次のような形でも示される[4]

この解の関数をロジスティック関数(英語:logistic function)[53]、この解によって描かれる曲線をロジスティック曲線(英語:logistic curve)と呼ぶ[42]。この曲線に従う個体群成長は、ロジスティック成長やロジスティック増殖とも呼ばれる[54][55]。関数は t → ∞ の極限NK となり、マルサスモデルと異なり発散しないことが確認できる[56]

曲線の形状

[編集]
いくつかの N0 から始まるロジスティック曲線。N > 0 の範囲では、時間の経過に従って NK に収束する。
時間と個体数が負の場合も含めたロジスティック曲線の全体図。縦軸を 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/dtN のグラフの頂点に相当する[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]NN0 から単調に減少しつづけ、この場合も、時間経過に従って 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]。すなわち、初期個体数 N0K または 0 であれば、時間経過によらず常に同じ値を取り続けることは同じだが、N0 が平衡状態から少しずれたときの挙動は正反対となる[65]

ロジスティック曲線とその傾きのベクトル場の様子

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

NdN/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]。ロジスティック方程式の場合は、

なので、

となり、f ′(K) = −r < 0, f ′(0) = r > 0 となることが確認できる[69]

生物学的前提条件

[編集]

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

  • 環境内には単一の種か、あるいは同等とみなせる種のみが存在する[70]
  • 対象の生物の各世代(親子)は連続的に重なっている[71]。すなわち、連続的に子が生まれ、親と子が共存する期間が存在する[72]
  • 個体は一定の大きさの環境内に常に存在する。すなわち、環境から移出したり、外部から移入が無い[73]。(用語としては閉じた個体群とも呼ばれる[74])
  • 環境の大きさは変わらず、一定状態が保たれる[73]
  • 個体群のために、食糧や資源が一定して供給される[75]

ショウジョウバエ真正細菌といった、微生物や単純な生物を一定環境で増殖させた場合は、上記の条件に近く、ロジスティック方程式によって個体数変化の正確な予測ができる[76]。しかし、例えば鹿鳥類などのような、一定環境のもとで増殖する設定が成立しない個体群成長には、ロジスティック方程式を適用することはできない[67]

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

実際のデータへの適用例

[編集]

実験生物

[編集]
ソ連・ロシアの生物学者ゲオルギー・ガウゼによる2種の酵母Saccharomyces cerevisiae, Schizosaccharomyces kefir)の個体群サイズ成長実験データ[78]と、それらのデータに対してフィッティングさせたロジスティック曲線を示した図。縦軸は菌全体の体積、横軸は時間を示している[78]。このガウゼの実験はロジスティック曲線がよく当てはまった個体群成長実験としてよく知られる[79]

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

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

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

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

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

2種存在する場合

[編集]
ロトカ・ヴォルテラの競争式では係数の値がある範囲内のときのみ2種が共存し(図の3)、それ以外ではどちらかが絶滅する(図の1, 2, 4)[124]

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

係数の N1, r1, K1 は種1の個体数、内的自然増加率、環境収容力である。同様に、N2, r2, K2 は種2の個体数、内的自然増加率、環境収容力である。さらに、a12 が種2が種1に与える影響を、a21 が種1が種2に与える影響を表し、競争係数と呼ばれる[128]。この式はアメリカの数学者アルフレッド・ロトカとイタリアの数学者ヴィト・ヴォルテラによって独立に考案された[128]

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

時間遅れの考慮

[編集]
ヒツジキンバエロバート・メイが時間遅れを持つロジスティック方程式でこのハエの個体数変動の解析を行った。

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

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

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

離散時間モデル

[編集]
離散時間型モデルの場合の個体数の変化の様子。いずれも KN0 の値は同じだが、青が r = 0.6、赤が r = 2.1、緑が r = 3 のときを示している。

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

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

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

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

ここで、Δt を差分時間間隔として、ab

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

生物個体数以外での例

[編集]

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

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

他形式

[編集]

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

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

とも表される[155]kVerhulst-Pearl係数種内競争係数と呼ばれる[156]。個体群密度が増加率を減少させる影響の強さを k が表しているといえる[157]

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

という形式もある[158]

非線形のロジスティック関数を扱いやすくするために線形の対数関数に変換する、フィッシャ・プライ変換(英語:Fisher-Pry transform)と呼ばれる次のような変換もある[150]

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

歴史

[編集]

フェルフルストによる発表

[編集]
ピエール=フランソワ・フェルフルスト (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]。この論文の中でフェルフルストが実際に提案した式は、

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

式の再発見と論争と普及

[編集]
レイモンド・パール (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年、ロバート・マッカーサーエドワード・オズボーン・ウィルソンは、この rK に着目して、における生物個体群の定着と絶滅に関する理論を発案した[190]。彼らの理論によれば、ある生物の島への定着が成功するには大きな r を持つことが重要であり、絶滅の回避には大きな K を持つことが重要であるとし、それぞれの方向へ淘汰されることを r淘汰K淘汰と呼んだ[191]。この説はr-K戦略説と呼ばれ、生物の生活史の進化に種内競争の観点から説明を与えた[192]

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

脚注

[編集]

注釈

[編集]
  1. ^ 一例として以下のように解くことができる。N の値の範囲を 0 < N < K に限定して解く方法と[50]、特に限定せずに解く方法がある[52]。ここでは範囲を限定しない解き方を示す。まずロジスティック方程式を変数分離変形して
    を得る。さらに左辺を部分分数分解すれば
    となる。両辺を積分して
    となり、ここで C積分定数である。両辺の指数をとり、絶対値を外せば
    となる。t = 0 のときの N の値を N0 で表せば
    なので、これを上式に代入して
    となる。式を整理して
    となる。

出典

[編集]
  1. ^ a b Strogatz 2015, p. 25.
  2. ^ a b c 大澤光『社会システム工学の考え方』(初版)オーム社、2007年、193-194頁。ISBN 978-4-274-06675-7 
  3. ^ スチュアート 2012, p. 78.
  4. ^ a b c 巌佐庸『生命の数理』(初版)共立出版、2008年2月25日、2-3頁。ISBN 978-4-320-05662-6 
  5. ^ a b c d 西村欣也『生態学のための数理的方法―考えながら学ぶ個体群生態学』文一総合出版、2012年、168-169頁。ISBN 978-4-8299-6520-7 
  6. ^ K.T.アリグッド; T.D.サウアー; J.A.ヨーク 著、星野高志・阿部巨仁・黒田拓・松本和宏 訳、シュプリンガー・ジャパン 編『カオス 第2巻 力学系入門』津田一郎(監訳)、丸善出版、2012年、92頁。ISBN 978-4-621-06279-1 
  7. ^ ティーメ 2006, p. 38.
  8. ^ 稲岡毅『基礎からの微分方程式―実例でよくわかる』森北出版、2012年、22-23頁。ISBN 978-4-627-07671-6 
  9. ^ Hirsch et al. 2007, pp. 4–7.
  10. ^ a b c マレー 2014, p. 1.
  11. ^ スチュアート 2012, p. 333.
  12. ^ 日本数理生物学会(編) 2008, p. 61.
  13. ^ 日本数理生物学会(編) 2008, p. 181; 寺本 1997, p. 2.
  14. ^ 大串 2014, p. 49.
  15. ^ 寺本 1997, pp. 2–3.
  16. ^ 瀬野 2007, p. 1.
  17. ^ 巌佐 1990, p. 2.
  18. ^ 山口 1992, p. 59.
  19. ^ ブラウン 2012, p. 30.
  20. ^ a b c d Jonathan W. Silvertown 著、河野昭一・高田壮典・大原雅 訳『植物の個体群生体学 第2版』(初版)東海大学出版会、1992年、49-50頁。ISBN 4-486-01157-0 
  21. ^ ハーバーマン 1992, p. 13.
  22. ^ a b 渡辺 2012, pp. 50–51.
  23. ^ 寺本 1997, pp. 4–5.
  24. ^ 巌佐 2015, p. 17.
  25. ^ 山本将史『Excelで学ぶ微分積分』オーム社、2012年、264頁。ISBN 978-4-274-06888-1 
  26. ^ 人口研究会(編) 2010, pp. 281–282; マレー 2014, p. 38.
  27. ^ マレー 2014, pp. 1–2.
  28. ^ a b c 巌佐 1990, p. 4.
  29. ^ a b 日本生態学会(編) 2004, p. 125.
  30. ^ 日本生態学会(編) 2004, p. 123.
  31. ^ 巌佐 2015, p. 21.
  32. ^ 寺本 1997, p. 8.
  33. ^ 渡辺守『昆虫の保全生態学』(初版)東京大学出版会、2007年12月20日、39-40頁。ISBN 978-4-13-062215-8 
  34. ^ 寺本 1997, pp. 10–11.
  35. ^ 瀬野 2007, pp. 11, 13–14; コーエン 1998, p. 112.
  36. ^ 瀬野 2007, p. 14.
  37. ^ コーエン 1998, p. 112; マレー 2014, p. 2.
  38. ^ 渡辺 2012, p. 51.
  39. ^ a b c 瀬野 2007, p. 20.
  40. ^ a b 日本数理生物学会(編) 2008, p. 4.
  41. ^ ミンモ・イアネリ、稲葉寿、國谷紀良『人口と感染症の数理―年齢構造ダイナミクス入門』(初版)東京大学出版会、2014年、52頁。ISBN 978-4-13-061309-5 
  42. ^ a b c d e f 人口研究会(編) 2010, p. 307.
  43. ^ ティーメ 2006, p. 66.
  44. ^ 大串 2014, p. 57.
  45. ^ a b c d 寺本 1997, p. 22.
  46. ^ a b c マレー 2014, p. 3.
  47. ^ a b c d e Hirsch et al. 2007, p. 6.
  48. ^ a b c d e f g h i j k l Strogatz 2015, p. 26.
  49. ^ a b スチュアート 2012, p. 335.
  50. ^ a b 山口 1992, pp. 62–65.
  51. ^ スチュアート 2012, p. 417.
  52. ^ 小寺忠・長谷川健二『工学系学生のための常微分方程式』(第2版)森北出版、2006年、44-46頁。ISBN 978-4-627-07452-1 
  53. ^ 和田光平『人口統計学の理論と推計への応用』(初版)オーム社、2015年、31-34頁。ISBN 978-4274217166 
  54. ^ 日本生物物理学会 編『数理生態学』巌佐庸(担当編集委員)(初版)、共立出版〈シリーズ・ニューバイオフィジックス (10)〉、1997年、33頁。ISBN 4-320-05473-3 
  55. ^ a b c d 寺本 1997, p. 10.
  56. ^ a b c Hirsch et al. 2007, p. 5.
  57. ^ a b c Hirsch et al. 2007, pp. 5–6.
  58. ^ 山口 1992, pp. 65–66.
  59. ^ ブラウン 2012, p. 34.
  60. ^ a b 山口 1992, p. 65.
  61. ^ Strogatz 2015, p. 19; Hirsch et al. 2007, p. 2.
  62. ^ a b 巌佐 2015, p. 23.
  63. ^ マレー 2014, p. 5.
  64. ^ Strogatz 2015, pp. 25–26.
  65. ^ 巌佐 2015, pp. 22–23.
  66. ^ K.T.アリグッド; T.D.サウアー; J.A.ヨーク 著、星野高志・阿部巨仁・黒田拓・松本和宏 訳、シュプリンガー・ジャパン 編『カオス 第2巻 力学系入門』津田一郎(監訳)、丸善出版、2012年、94頁。ISBN 978-4-621-06279-1 
  67. ^ a b 巌佐 2015, p. 24.
  68. ^ マレー 2014, p. 5; Strogatz 2015, p. 28.
  69. ^ Strogatz 2015, p. 29.
  70. ^ 瀬野 2007, pp. 5, 21–22.
  71. ^ マレー 2014, p. 37.
  72. ^ 山口 1992, p. 71.
  73. ^ a b コーエン 1998, p. 112.
  74. ^ ティーメ 2006, p. 7.
  75. ^ Strogatz 2015, p. 27; コーエン 1998, p. 112.
  76. ^ 巌佐 2015, p. 24; コーエン 1998, p. 116; Strogatz 2015, p. 27.
  77. ^ a b c 山口 1992, p. 66.
  78. ^ a b Gause, G. F. (1932). “Experimental Studies on the Struggle for Existence”. Journal of Experimental Biology (The Company of Biologists Ltd) 9 (4): 389-402. ISSN 1477-9145. http://jeb.biologists.org/content/9/4/389. 
  79. ^ a b スチュアート 2012, pp. 335–336.
  80. ^ 渡辺 2012, p. 53.
  81. ^ a b c d 内田 1972, p. 29.
  82. ^ 山口 1992, pp. 67–68.
  83. ^ G. F. Gause. The Struggle for Existence, p. 106, - Google ブックス
  84. ^ 山口 1992, pp. 67–69.
  85. ^ 巌佐 1990, p. 6.
  86. ^ Gause, G. F. (2003 (original 1934)). The Struggle for Existence (Dover Phoenix Editions ed.). Dover Pubns. p. 77. ISBN 0-486-49520-5. https://books.google.co.jp/books?id=v01OToAhJboC&lpg=PP1&hl=ja&pg=PA77#v=onepage&q&f=false 
  87. ^ 山口 1992, pp. 69–71; Strogatz 2015, p. 27.
  88. ^ a b Strogatz 2015, p. 27.
  89. ^ a b 山口 1992, pp. 69–71.
  90. ^ 内田 1972, pp. 27–29.
  91. ^ 内田 1972, pp. 27–28.
  92. ^ 内田 1972, p. 27.
  93. ^ a b 内田 1972, p. 28.
  94. ^ Kingsland 1982, p. 44.
  95. ^ a b c 伊藤嘉昭『生態学と社会―経済・社会系学生のための生態学入門』(初版)東海大学出版会、1994年、48–51頁。ISBN 4-486-01272-0 
  96. ^ 大串 2014, pp. 53–54.
  97. ^ レーヴンほか 2007, p. 1151.
  98. ^ Kingsland 1982, p. 30.
  99. ^ a b Kingsland 1982, p. 31.
  100. ^ コーエン 1998, p. 113–115.
  101. ^ 山口 1992, pp. 58–59.
  102. ^ a b コーエン 1998, p. 115.
  103. ^ 寺本 1997, p. iii.
  104. ^ 巌佐 2015, p. 25.
  105. ^ 山口 1992, p. 66; スチュアート 2012, p. 336.
  106. ^ マレー 2014, p. 3; Strogatz 2015, p. 27.
  107. ^ コーエン 1998, p. 116より引用.
  108. ^ Kingsland 1982, p. 36.
  109. ^ 人口研究会(編) 2010, p. 74; コーエン 1998, p. 144.
  110. ^ 山重慎二・加藤久和・小黒一正『人口動態と政策―経済学的アプローチへの招待』(第1版)日本評論社、2013年9月10日、28頁。ISBN 978-4-535-55750-5 
  111. ^ マレー 2014, p. 4.
  112. ^ 渡辺 2012, p. 56.
  113. ^ 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. http://gdz.sub.uni-goettingen.de/dms/load/img/?PPN=PPN129323640_0018&DMDID=dmdlog7. のp. 8より引用
  114. ^ 人口研究会(編) 2010, p. 307; 瀬野 2007, p. 20.
  115. ^ a b O'Connor, John J; Robertson, Edmund F (2014年1月). “Pierre François Verhulst”. MacTutor History of Mathematics archive. University of St Andrews. 2015年12月17日閲覧。 出典での引用元は、Pastijn, Hugo (2006). “Chaotic Growth with the Logistic Model of P.-F. Verhulst”. The logistic map and the route to chaos: From the beginnings to modern applications, Understanding Complex Systems. Springer Berlin Heidelberg. pp. 3-11. doi:10.1007/3-540-32023-7_1. ISBN 978-3-540-28366-9. http://link.springer.com/chapter/10.1007/3-540-32023-7_1 
  116. ^ 寺本 1997, p. 10; 人口研究会(編) 2010, p. 307.
  117. ^ 北原(編) 2003, pp. 1–2.
  118. ^ a b c d e f 白山義久・桜井泰憲・古谷研・中原裕幸・松田裕之・加々美康彦 編『海洋保全生態学』講談社、2012年、172-174頁。ISBN 978-4-06-155231-9 
  119. ^ 北原(編) 2003, p. 244.
  120. ^ 日本数理生物学会(編) 2008, p. 187.
  121. ^ 寺本 1997, p. 23.
  122. ^ a b 島一雄・關文威・前田昌調・木村伸吾・佐伯宏樹・桜本和美・末永芳美・長野章・森永勤・八木信行・山中英明 編『最新 水産ハンドブック』講談社サイエンティフィク、2012年、114頁。ISBN 978-4-06-153736-1 
  123. ^ 日本数理生物学会(編) 2008, p. 188.
  124. ^ ハーバーマン 1992, pp. 130–136.
  125. ^ 巌佐 1990, p. 13.
  126. ^ マレー 2014, p. 65.
  127. ^ 大串 2014, p. 66.
  128. ^ a b 日本生態学会(編) 2004, p. 133.
  129. ^ 巌佐 1990, p. 15.
  130. ^ ハーバーマン 1992, p. 134.
  131. ^ a b マレー 2014, p. 12.
  132. ^ Ruan, S. (2006). “DELAY DIFFERENTIAL EQUATIONS IN SINGLE SPECIES DYNAMICS”. In Arino, O.; Hbid, M.L.; Dads, E. Ait. Delay Differential Equations and Applications: Proceedings of the NATO Advanced Study Institute held in Marrakech, Morocco, 9-21 September 2002 (1 ed.). Springer Netherlands. pp. 479. doi:10.1007/1-4020-3647-7 
  133. ^ 日本数理生物学会(編) 2008, p. 67.
  134. ^ 日本数理生物学会(編) 2008, p. 68.
  135. ^ ハーバーマン 1992, pp. 46–47.
  136. ^ 日本数理生物学会(編) 2008, pp. 70–71.
  137. ^ マレー 2014, pp. 13–14.
  138. ^ 寺本 1997, p. 11.
  139. ^ a b 巌佐 1990, p. 49.
  140. ^ a b ウィロックス, ラルフ「自然現象の離散・超離散系によるモデル化」『応用力学研究所研究集会報告 No.22AO-S8』、九州大学応用力学研究所、2011年3月、3頁、doi:10.15017/23387 
  141. ^ Hirsch et al. 2007, p. 347.
  142. ^ 巌佐 1990, p. 51; マレー 2014, p. 44.
  143. ^ 巌佐 1990, p. 53; マレー 2014, pp. 44–47.
  144. ^ Masaaki Morisita (6 1965). “The fitting of the logistic equation to the rate of increase of population density”. Researches on Population Ecology (Springer-Verlag) 7: 52–55. doi:10.1007/BF02518815. 
  145. ^ 山口 1992, p. 72.
  146. ^ 山口 1992, p. 73.
  147. ^ 巌佐 1990, p. 5.
  148. ^ 北原(編) 2003, pp. 220–221.
  149. ^ ブラウン 2012, pp. 45–46.
  150. ^ a b 渡辺千仭『技術経済システム』(初版)創成社、2007年、84–87頁。ISBN 978-4-7944-3089-2 
  151. ^ 橋本洋志・石井千春『Scilab/Scicosで学ぶシミュレーションの基礎―自然・社会現象から、経済・金融、システム制御まで』オーム社、2008年、103–105頁。ISBN 978-4-274-20487-6 
  152. ^ 白石高章『統計科学の基礎―データと確率の結びつきがよくわかる数理』(第1版)日本評論社、2012年、159頁。ISBN 978-4-535-78700-1 
  153. ^ 岡谷貴之『深層学習』講談社サイエンティフィク〈機械学習プロフェッショナルシリーズ〉、2015年4月7日、10頁。ISBN 978-4-06-152902-1 
  154. ^ 瀬野 2007, pp. i, 25–26.
  155. ^ 瀬野 2007, p. 22; ブラウン 2012, p. 32; ハーバーマン 1992, p. 41.
  156. ^ 寺本 1997, p. 9.
  157. ^ 瀬野 2007, p. 21.
  158. ^ Hirsch et al. 2007, p. 4.
  159. ^ コーエン 1998, pp. 112–113.
  160. ^ 山口 1992, p. 54.
  161. ^ 人口研究会(編) 2010, pp. 280–282.
  162. ^ 人口研究会(編) 2010, p. 282.
  163. ^ 人口研究会(編) 2010, p. 315.
  164. ^ 山口 1992, p. 55; 寺本 1997, p. 10.
  165. ^ 山口 1992, p. 56.
  166. ^ a b Verhulst, Pierre-François (1838). “Notice sur la loi que la population suit dans son accroissement”. Correspondance mathématique et physique 10: 113-121. https://books.google.co.jp/books?id=8GsEAAAAYAAJ&pg=PA113&hl=ja&source=gbs_toc_r&cad=4#v=onepage&q&f=false. のp.115
  167. ^ a b 瀬野 2007, p. 18.
  168. ^ 山口 1992, p. 57.
  169. ^ a b c d 山口 1992, p. 58.
  170. ^ コーエン 1998, p. 113.
  171. ^ 内田 1972, pp. 28–30.
  172. ^ マレー 2014, pp. 4–5; コーエン 1998, p. 113.
  173. ^ Kingsland 1982, pp. 33–34.
  174. ^ コーエン 1998, pp. 113–114より引用。出典での引用元は、Pearl, Raymond, and Lowell J. Reed. 1924 The growth of human population. In Studies in human biology, ed. Raymond Pearl. Baltimore: Williams and Wilkins, pp. 584-637.のp.585より.
  175. ^ Kingsland 1982, pp. 32, 35, 42.
  176. ^ マレー 2014, pp. 3–4; 人口研究会(編) 2010, p. 323.
  177. ^ Kingsland 1982, p. 32.
  178. ^ Kingsland 1982, pp. 35–36, 38–39.
  179. ^ Kingsland 1982, p. 40.
  180. ^ Kingsland 1982, p. 34.
  181. ^ Kingsland 1982, pp. 41–42.
  182. ^ Kingsland 1982, p. 43.
  183. ^ Kingsland 1982, pp. 45–46.
  184. ^ 日本生態学会(編) 2004, p. 241.
  185. ^ a b Kingsland 1982, p. 49.
  186. ^ 巌佐 1990, p. 50.
  187. ^ 山口 1992, pp. 73–75.
  188. ^ 日本生態学会(編) 2004, p. 62.
  189. ^ 日本生態学会(編) 2004, p. 62; レーヴンほか 2007, p. 1150.
  190. ^ 木元 1979, pp. 108–109.
  191. ^ 木元 1979, pp. 116–117.
  192. ^ 日本生態学会(編) 2004, pp. 61, 64.
  193. ^ 山口 1992, p. 76.
  194. ^ スチュアート 2012, pp. 342–343.
  195. ^ スチュアート 2012, p. 342.
  196. ^ ティェンイェン・リー、ジェームス・A・ヨーク「第10章 区間上のカオスを探索する」『カオスはこうして発見された』ラルフ・エイブラハム、ヨシスケ・ウエダ (編) 稲垣耕作、赤松則男(訳)(初版)、共立出版、2002年、169-170頁。ISBN 4-320-03418-X 

参考文献

[編集]

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

関連項目

[編集]

外部リンク

[編集]