(57)【特許請求の範囲】
【請求項１】
野生動物の個体群動態について事前分布を用いて事後分布の導出を行うベイズ推定により前記野生動物の個体群動態を推定する野生動物の個体群動態推定装置であって、
前記野生動物の個体群動態の前記事前分布を作成する事前分布作成手段と、
作成された前記事前分布と、前記野生動物の個体群動態に関する既知のデータとを用いて、前記野生動物の個体群動態の前記事後分布を作成する事後分布作成手段とを備え、
前記事前分布作成手段は、
前記事前分布を規定する期待値、分散および尖度を独立して調整可能な式（Ａ）に示す第１の関数と、
前記事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、
前記期待値、前記分散および前記尖度が前記野生動物の個体群動態に対応するように調整され、前記事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する前記事前分布を作成するように構成されている、野生動物の個体群動態推定装置。
【数１】
0005992369.tif 000007
【請求項２】
前記第２の関数は、前記事前分布を非対称の分布形状に変換可能であり、
前記事前分布作成手段は、前記事前分布の上限値および下限値の前記一方の値において連続的な分布形状を有するように前記事前分布の上限値および下限値の前記一方が規定された状態の前記第２の関数を用いて、前記事前分布の分布形状を非対称形状にするように構成されている、請求項１に記載の野生動物の個体群動態推定装置。
【請求項３】
前記事後分布作成手段は、前記野生動物の個体群動態に関する既知のデータが記憶された記憶部から前記既知のデータを読み出し、読み出した前記既知のデータと、前記期待値、前記分散および前記尖度が調整されることにより作成された前記事前分布とを用いて、前記野生動物の個体群動態の前記事後分布を作成するように構成されている、請求項１または２に記載の野生動物の個体群動態推定装置。
【請求項４】
前記野生動物の個体群動態は、野生動物の自然増加率を含む、請求項１〜３のいずれか１項に記載の野生動物の個体群動態推定装置。
【請求項５】
前記第２の関数は、前記事前分布の上限値を規定可能な式（Ｂ）に示す関数である、請求項１〜４のいずれか１項に記載の野生動物の個体群動態推定装置。
【数２】
0005992369.tif 000008
ここで、Ｒは、事前分布の確率変数、ＲMAXは、上限値をそれぞれ表す。
【請求項６】
コンピュータを、
野生動物の個体群動態について事前分布を用いて事後分布を導出するベイズ推定において、前記野生動物の個体群動態の前記事前分布を作成する事前分布作成手段と、
作成された前記事前分布と、前記野生動物の個体群動態に関する既知のデータとを用いて、前記野生動物の個体群動態の前記事後分布を作成する事後分布作成手段として機能させ、
前記事前分布作成手段は、
前記事前分布を規定する期待値、分散および尖度を独立して調整可能な式（Ａ）に示す第１の関数と、
前記事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、
前記期待値、前記分散および前記尖度が前記野生動物の個体群動態に対応するように調整され、前記事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する前記事前分布を作成するように構成されている、野生動物の個体群動態推定プログラム。
【数３】
0005992369.tif 000009
【請求項７】
野生動物の個体群動態について事前分布を用いて事後分布の導出を行うベイズ推定により前記野生動物の個体群動態を推定する野生動物の個体群動態推定方法であって、
前記野生動物の個体群動態の前記事前分布を作成するステップと、
前記事前分布と、前記野生動物の個体群動態に関する既知のデータとを用いて、前記野生動物の個体群動態の前記事後分布を作成するステップとを備え、
前記事前分布を作成するステップは、
前記事前分布を規定する期待値、分散および尖度を独立して調整可能な式（Ａ）に示す第１の関数と、
前記事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、
前記期待値、前記分散および前記尖度が前記野生動物の個体群動態に対応するように調整され、前記事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する前記事前分布を作成するステップを含む、野生動物の個体群動態推定方法。
【数４】
0005992369.tif 000010
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、野生動物の個体群動態推定装置、野生動物の個体群動態推定プログラムおよび野生動物の個体群動態推定方法に関する。
【背景技術】
【０００２】
従来から、シカやイノシシなどの野生動物による農作物被害が深刻化しており、野生動物の保全管理を適切に行うことが望まれている。現状では、都道府県や市町村において、野生動物の自然増加率や個体数等の推定に基づいて、必要な捕獲数の決定等、野生動物の保全管理に関する意思決定が行われている。
【０００３】
また、従来、野生動物の自然増加率や個体数等（個体群動態）の推定を行う際に、野生動物の個体群動態についての事前分布を用いて事後分布の導出（推定）を行うベイズ推定により野生動物の個体群動態を推定する方法が提案されている（たとえば、非特許文献１参照）。
【０００４】
上記非特許文献１には、ベイズ推定により、ニホンジカの自然増加率や個体数等（個体群動態）についての事前分布と既知のデータとに基づいて事後分布を導出することによって、ニホンジカの個体群動態を推定する方法（野生動物の個体群動態推定方法）が開示されている。この野生動物の個体群動態推定方法では、自然増加率および個体数についての事前分布を正規分布として設定し、その事前分布を用いて事後分布を作成している。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】「「ニホンジカの個体群動態の推定と将来予測（兵庫県本州部２０１１年）」、兵庫ワイルドライフレポート１号（２０１２）Ｐ１〜１６、兵庫県森林動物研究センター」
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記非特許文献１の野生動物の個体群動態推定方法では、ニホンジカの自然増加率および個体数（個体群動態）について事前分布が正規分布として設定されるので、たとえば、野生動物の個体群動態の確率密度分布において、期待値近傍の確率密度が正規分布よりも広い範囲にわたって一様に高いと想定される場合には、正規分布に従った事前分布では、その想定に沿った分布形状に精度よく設定することができないという不都合がある。この場合には、事前分布の設定精度が低下するので、事前分布に基づく事後分布の精度が低下し、その結果、野生動物の個体群推定の精度が低下するという問題点がある。
【０００７】
この発明は、上記のような課題を解決するためになされたものであり、この発明の１つの目的は、事前分布の設定精度を高めることにより、事後分布の精度を向上させて野生動物の個体群推定の精度を向上させることが可能な野生動物の個体群動態推定装置、野生動物の個体群動態推定プログラムおよび野生動物の個体群動態推定方法を提供することである。
【課題を解決するための手段】
【０００８】
この発明の第１の局面による野生動物の個体群動態推定装置は、野生動物の個体群動態について事前分布を用いて事後分布の導出を行うベイズ推定により野生動物の個体群動態を推定する野生動物の個体群動態推定装置であって、野生動物の個体群動態の事前分布を作成する事前分布作成手段と、作成された事前分布と、野生動物の個体群動態に関する既知のデータとを用いて、野生動物の個体群動態の事後分布を作成する事後分布作成手段とを備え、事前分布作成手段は、事前分布を規定する期待値、分散および尖度を独立して調整可能な式（Ａ）に示す第１の関数と、事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、期待値、分散および尖度が野生動物の個体群動態に対応するように調整され、事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する事前分布を作成するように構成されている。
【数４】
0005992369.tif 000002
【０００９】
この発明の第１の局面による野生動物の個体群動態推定装置では、上記のように構成することによって、事前分布を正規分布として設定する場合とは異なり、野生動物の個体群動態の期待値近傍の確率密度が正規分布とは異なる分布形状を有すると想定される場合でも、野生動物の個体群動態に対応するように期待値、分散および尖度が調整された第１の関数を用いることにより、その想定に沿った分布形状に事前分布を精度よく設定することができる。すなわち、野生動物の個体群動態に関する知見を十分に活かして想定される分布形状に事前分布を精度よく設定することができる。その結果、野生動物の個体群動態の事前分布の設定精度を高めることにより、事後分布の精度を向上させることができるので、野生動物の個体群推定の精度を向上させることができる。また、たとえば、ニホンジカの場合、生態学的な特性および調査結果から得られる妊娠率や産仔数等の情報に基づいて、自然増加率（個体群動態）の上限値および下限値を想定することが可能であり、このような場合に、ニホンジカの自然増加率の上限値および下限値に対応する事前分布の上限値および下限値を規定可能な第２の関数を用いれば、事前分布をより精度よく設定することができる。
【００１２】
この場合、好ましくは、第２の関数は、事前分布を非対称の分布形状に変換可能であり、事前分布作成手段は、事前分布の上限値および下限値の一方の値において連続的な分布形状を有するように事前分布の上限値および下限値の一方が規定された状態の第２の関数を用いて、事前分布の分布形状を非対称形状にするように構成されている。このように構成すれば、野生動物の生態学的な特性や調査結果に基づいて、野生動物の個体群動態の上限値および下限値に加えて、野生動物の個体群動態の確率密度分布が非対称形状であると想定することが可能な場合に、野生動物の個体群動態の上限値（下限値）に対応するように事前分布の上限値（下限値）が規定され、かつ、事前分布の分布形状を非対称形状に変換可能な第２の関数を用いることにより、事前分布をさらに精度よく設定することができる。
【００１３】
上記第１の局面による野生動物の個体群動態推定装置において、好ましくは、事後分布作成手段は、野生動物の個体群動態に関する既知のデータが記憶された記憶部から既知のデータを読み出し、読み出した既知のデータと、期待値、分散および尖度が調整されることにより作成された事前分布とを用いて、野生動物の個体群動態の事後分布を作成するように構成されている。このように構成すれば、事後分布作成手段により、記憶部に記憶された野生動物の個体群動態に関する既知のデータが自動的に読み出されるので、容易に、読み出した既知のデータと期待値、分散および尖度が調整された事前分布とに基づいて事後分布を精度よく作成することができる。
【００１４】
上記第１の局面による野生動物の個体群動態推定装置において、好ましくは、野生動物の個体群動態は、野生動物の自然増加率を含む。この場合、事前分布の期待値、分散および尖度が野生動物の自然増加率に対応するように調整された状態の第１の関数を用いて、野生動物の自然増加率について事前分布の設定精度を高めて事後分布の精度を向上させることができる。これにより、野生動物の保全管理に関する意思決定において重要なパラメータである自然増加率の推定を精度よく行うことができる。
上記第１の局面による野生動物の個体群動態推定装置において、好ましくは、第２の関数は、事前分布の上限値を規定可能な式（Ｂ）に示す関数である。
【数５】
0005992369.tif 000003
ここで、Ｒは、事前分布の確率変数、ＲMAXは、上限値をそれぞれ表す。
【００１５】
この発明の第２の局面による野生動物の個体群動態推定プログラムは、コンピュータを、野生動物の個体群動態について事前分布を用いて事後分布を導出するベイズ推定において、野生動物の個体群動態の事前分布を作成する事前分布作成手段と、作成された事前分布と、野生動物の個体群動態に関する既知のデータとを用いて、野生動物の個体群動態の事後分布を作成する事後分布作成手段として機能させ、事前分布作成手段は、事前分布を規定する期待値、分散および尖度を独立して調整可能な上式（Ａ）に示す第１の関数と、事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、期待値、分散および尖度が野生動物の個体群動態に対応するように調整され、事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する事前分布を作成するように構成されている。
【００１６】
この発明の第２の局面による野生動物の個体群動態推定プログラムでは、上記のように構成することによって、事前分布を正規分布として設定する場合とは異なり、野生動物の個体群動態の期待値近傍の確率密度が正規分布とは異なる分布形状を有すると想定される場合でも、野生動物の個体群動態に対応するように期待値、分散および尖度が調整された関数を用いることにより、その想定に沿った分布形状に事前分布を精度よく設定することができる。すなわち、野生動物の個体群動態に関する知見を十分に活かして想定される分布形状に事前分布を精度よく設定することができる。その結果、野生動物の個体群動態の事前分布の設定精度を高めることにより、事後分布の精度を向上させることができるので、野生動物の個体群推定の精度を向上させることができる。また、たとえば、ニホンジカの場合、生態学的な特性および調査結果から得られる妊娠率や産仔数等の情報に基づいて、自然増加率（個体群動態）の上限値および下限値を想定することが可能であり、このような場合に、ニホンジカの自然増加率の上限値および下限値に対応する事前分布の上限値および下限値を規定可能な第２の関数を用いれば、事前分布をより精度よく設定することができる。
【００１７】
この発明の第３の局面による野生動物の個体群動態推定方法は、野生動物の個体群動態について事前分布を用いて事後分布の導出を行うベイズ推定により野生動物の個体群動態を推定する野生動物の個体群動態推定方法であって、野生動物の個体群動態の事前分布を作成するステップと、事前分布と、野生動物の個体群動態に関する既知のデータとを用いて、野生動物の個体群動態の事後分布を作成するステップとを備え、事前分布を作成するステップは、事前分布を規定する期待値、分散および尖度を独立して調整可能な上式（Ａ）に示す第１の関数と、事前分布の上限値又は下限値の少なくとも一方を含む第２の関数と、を用いて、期待値、分散および尖度が野生動物の個体群動態に対応するように調整され、事前分布の上限値又は下限値の少なくとも一方において連続的な分布形状を有する事前分布を作成するステップを含む。
【００１８】
この発明の第３の局面による野生動物の個体群動態推定方法では、上記のように構成することによって、事前分布を正規分布として設定する場合とは異なり、野生動物の個体群動態の期待値近傍の確率密度が正規分布とは異なる分布形状を有すると想定される場合でも、野生動物の個体群動態に対応するように期待値、分散および尖度が調整された第１の関数を用いることにより、その想定に沿った分布形状に事前分布を精度よく設定することができる。すなわち、野生動物の個体群動態に関する知見を十分に活かして想定される分布形状に事前分布を精度よく設定することができる。その結果、野生動物の個体群動態の事前分布の設定精度を高めることにより、事後分布の精度を向上させることができるので、野生動物の個体群推定の精度を向上させることができる。また、たとえば、ニホンジカの場合、生態学的な特性および調査結果から得られる妊娠率や産仔数等の情報に基づいて、自然増加率（個体群動態）の上限値および下限値を想定することが可能であり、このような場合に、ニホンジカの自然増加率の上限値および下限値に対応する事前分布の上限値および下限値を規定可能な第２の関数を用いれば、事前分布をより精度よく設定することができる。
【発明の効果】
【００２０】
本発明によれば、上記のように、事前分布の設定精度を高めることにより、事後分布の精度を向上させて野生動物の個体群推定の精度を向上させることができる。
【図面の簡単な説明】
【００２１】
【図１】本発明の一実施形態による個体群動態推定装置により行われるベイズ推定の概要を説明するための図である。
【図２】本発明の一実施形態による個体群動態推定装置の構成を示したブロック図である。
【図３】本発明の一実施形態による個体群動態推定装置により行われる第１の関数を用いた事前分布の作成について説明するための図である。
【図４】本発明の一実施形態による個体群動態推定装置により行われる第２の関数を用いた事前分布の作成について説明するための図である。
【図５】本発明の一実施形態による個体群動態推定装置により行われる個体群動態の推定処理について説明するためのフローチャートである。
【発明を実施するための形態】
【００２２】
以下、本発明を具体化した実施形態を図面に基づいて説明する。
【００２３】
まず、図１〜図４を参照して、本発明の一実施形態による個体群動態推定装置１００の構成について説明する。
【００２４】
本発明の一実施形態による個体群動態推定装置１００は、野生動物の自然増加率や個体数、捕獲率等の個体群動態について、事前分布を用いて事後分布の導出を行うベイズ推定により野生動物の個体群動態を推定する装置である。たとえば、野生動物の個体数の過程モデルを以下に示す式（１）のように仮定した場合、各年の個体数Ｎｔおよび自然増加率ｒは、未知の値（パラメータ）である。このようなパラメータになる個体群動態（個体数Ｎｔや自然増加率ｒ）について、個体群動態推定装置１００は、図１に示すように、個体群動態の事前分布と、捕獲数や糞塊密度、目撃効率（狩猟者から得られた目撃情報）等の野生動物の個体群動態に関する既知のデータとを用いて、マルコフ連鎖モンテカルロ法（ＭＣＭＣ）により、事後分布を導出して個体群動態の推定を行うように構成されている。
【００２５】
【数１】
0005992369.tif 000004
ここで、Ｎｔは、ｔ年における個体数、ｒは、自然増加率、ｃｔは、ｔ年における捕獲数をそれぞれ表す。
【００２６】
なお、本実施形態の個体群動態推定装置１００では、個体群動態の推定に際して、狩猟時の捕獲数（狩猟捕獲数）や有害駆除による捕獲数（有害捕獲数）、狩猟登録者数あたりの銃猟による捕獲数などの変数を用いて作成された確率モデルと、上記過程モデルとを用いて、マルコフ連鎖モンテカルロ法による階層ベイズ推定を行う。以下、個体群動態推定装置１００の具体的な構成について説明する。
【００２７】
個体群動態推定装置１００は、図２に示すように、ＣＰＵ１と、ハードディスク２と、読出装置３と、通信部４と、出力部５とから主として構成されている。ＣＰＵ１、ハードディスク２、読出装置３、通信部４、および出力部５は、バス６によって接続されている。なお、ＣＰＵ１は、本発明の「事前分布作成手段」および「事後分布作成手段」の一例であり、ハードディスク２は、本発明の「記憶部」の一例である。また、個体群動態推定装置１００は、本発明の「コンピュータ」の一例である。
【００２８】
ＣＰＵ１は、個体群動態推定装置１００の各部を制御可能に構成されている。また、ＣＰＵ１は、ハードディスク２にインストールされた後述の個体群動態推定プログラム１０ａを実行することにより、野生動物の個体群動態の推定処理を実施するように構成されている。
【００２９】
ハードディスク２には、図示しないオペレーティングシステムや個体群動態推定プログラム１０ａなど、ＣＰＵ１に実行させるための種々のコンピュータプログラムがインストールされている。また、ハードディスク２には、捕獲数や糞塊密度、目撃効率等の野生動物の個体群動態に関する既知のデータも格納されている。既知のデータは、図示しないキーボードを用いて打ち込むことによりハードディスク２に保存してもよいし、Ｅｔｈｅｒｎｅｔ（登録商標）インタフェースからなる通信部４を介して他の機器から取り込んでもよい。
【００３０】
読出装置３は、ディスクドライブにより構成されており、ＣＤ−ＲＯＭ等の可搬型記録媒体１０に記録された情報を読み出すことが可能である。個体群動態推定装置１００は、読出装置３により、可搬型記録媒体１０に記録された個体群動態推定プログラム１０ａを読み出してハードディスク２にインストールすることが可能である。また、個体群動態推定プログラム１０ａを、Ｅｔｈｅｒｎｅｔ（登録商標）インタフェースからなる通信部４を介して受信するようにしてもよい。
【００３１】
出力部５は、表示装置２０に接続されており、後述の個体群動態の推定処理において、事前分布および事後分布それぞれのグラフの情報や推定結果の情報などを表示装置２０に出力可能に構成されている。
【００３２】
ここで、本実施形態では、ＣＰＵ１は、個体群動態の事前分布を規定する尖度、期待値および分散が独立して調整された状態の関数を用いて、野生動物の個体群動態の事前分布を作成する機能を有している。具体的には、ＣＰＵ１は、図３に示すように、以下の式（２）に示す第１の関数において事前分布の尖度、期待値および分散がそれぞれ対象の野生動物の個体群動態に対応するようにユーザにより調整された状態で、第１の関数を用いて事前分布を作成する。なお、第１の関数は、尖度、期待値および分散を独立して調整可能な確率密度関数（誤差分布または指数べき分布）である。また、第１の関数は、本発明の「関数」の一例である。
【００３３】
【数２】
0005992369.tif 000005
【００３４】
すなわち、上記式（２）に示す第１の関数において、ユーザが対象の野生動物の個体群動態に対応するように尖度、期待値および分散を独立して調整することによって、ＣＰＵ１は、尖度、期待値および分散が調整された状態の上記第１の関数を用いて、事野生動物の個体群動態に対応する前分布を作成する。
【００３５】
また、ＣＰＵ１は、尖度、期待値および分散が調整された事前分布（図３参照）に対して、図４に示すように、事前分布の上限値において連続的な分布形状を有するように事前分布の上限値が規定された状態で、野生動物の個体群動態の事前分布を作成するように構成されている。具体的には、ＣＰＵ１は、以下の式（３）に示す第２の関数（変換式）において、野生動物の個体群動態に関する数値の上限値に対応するとともに事前分布の上限値において連続的な分布形状を有するように事前分布の上限値が規定された状態で、上記第１の関数および第２の関数を用いて、野生動物の個体群動態の事前分布を作成する。第２の関数は、野生動物の個体群動態に関する数値の上限値に対応する事前分布の上限値を規定可能な関数である。
【００３６】
【数３】
0005992369.tif 000006
ここで、Ｒは、推定する確率変数（たとえば、自然増加率）、ＲＭＡＸは、上限値、ｐは、上記式（２）の第１の関数に従う変数をそれぞれ表す。
【００３７】
たとえば、ニホンジカの自然増加率（個体群動態）は、生態学的な特性および調査結果から得られる妊娠率や産仔数等の情報に基づいて、上限値（たとえば、１．５）および下限値（たとえば、１．０）を想定することが可能である。この場合、第２の関数のＲＭＡＸに、ニホンジカの自然増加率の上限値（たとえば、１．５）を適用すれば、事前分布の上限値において連続的な分布形状を有した状態で、ニホンジカの自然増加率の上限値に対応する値になるように事前分布の上限値を規定することができる。なお、上記第２の関数では、事前分布の下限値は規定されない。
【００３８】
また、上記第２の関数は、事前分布の上限値を規定しながら、野生動物の個体群動態に対応するように事前分布を非対称の分布形状に変換可能な関数である。ＣＰＵ１は、上記のように事前分布の上限値が規定された状態の第２の関数を用いて、事前分布の分布形状を非対称形状にするように構成されている。
【００３９】
野生動物の個体群動態は、生態学的な特性や調査結果に基づいて、確率密度分布が期待値に対して非対称形状（たとえば、図４に示すように、上限値側に偏った分布形状）になると想定することができる場合がある。この場合には、野生動物の個体群動態に対応するように事前分布を非対称の分布形状に変換可能な上記第２の関数を用いることは、特に有効である。
【００４０】
また、ＣＰＵ１は、ハードディスク２から既知のデータ（捕獲数や糞塊密度、目撃効率等）を読み出すように構成されている。ＣＰＵ１は、読み出した既知のデータと、尖度、期待値および分散が調整されることにより作成された事前分布とを用いて、野生動物の個体群動態の事後分布を作成することにより個体群動態の推定を行うように構成されている。具体的には、ＣＰＵ１は、既知のデータと事前分布とを用いて、マルコフ連鎖モンテカルロ法により、野生動物の個体群動態の事後分布を作成するように構成されている。
【００４１】
次に、図１および図３〜図５を参照して、個体群動態推定装置１００のＣＰＵ１により実行される個体群動態の推定処理について説明する。
【００４２】
まず、ステップＳ１において、ＣＰＵ１は、ユーザにより野生動物の個体群動態に対応するように尖度、期待値および分散が調整された上記式（２）に示す第１の関数を用いて、図３に示すように、野生動物の個体群動態の事前分布を作成する。その後、ステップＳ２において、ＣＰＵ１は、第１の関数により調整された事前分布に対して、図４に示すように、事前分布の上限値が規定され、かつ、事前分布を非対称形状に変換可能な上記式（３）に示す第２の関数を用いて、上限値が規定された非対称形状の事前分布を作成する。すなわち、ＣＰＵ１は、上記ステップＳ１およびＳ２の動作により、第１の関数および第２の関数の両方を用いて、事前分布を作成する。
【００４３】
その後、ＣＰＵ１は、ステップＳ３において、ハードディスク２から既知のデータを読み出す。ＣＰＵ１は、ステップＳ４において、読み出した既知のデータと、作成された事前分布とを用いて、図１に示すように、マルコフ連鎖モンテカルロ法により、対象の個体群動態の事後分布を作成し、ステップＳ５において、個体群動態の推定結果（グラフデータや数値データ）を表示装置２０に出力する制御を行う。
【００４４】
本実施形態では、上記のように、事前分布を規定する期待値、分散および尖度のうち、少なくとも尖度を独立して調整可能な第１の関数において事前分布の尖度が野生動物の個体群動態に対応するように調整された状態で、第１の関数を用いて野生動物の個体群動態の事前分布を作成するＣＰＵ１を設けることによって、事前分布を正規分布として設定する場合とは異なり、野生動物の個体群動態の期待値近傍の確率密度が正規分布とは異なる分布形状を有すると想定される場合でも、野生動物の個体群動態に対応するように尖度が調整された第１の関数を用いることにより、その想定に沿った分布形状に事前分布を精度よく設定することができる。その結果、野生動物の個体群動態の事前分布の設定精度を高めることにより、事後分布の精度を向上させることができるので、野生動物の個体群推定の精度を向上させることができる。
【００４５】
また、本実施形態では、事前分布の尖度、期待値および分散が野生動物の個体群動態に対応するように調整された状態の第１の関数を用いて、野生動物の個体群動態の事前分布を作成するようにＣＰＵ１を構成する。これにより、野生動物の個体群動態の想定される確率密度分布に沿った分布形状に事前分布をより精度よく設定することができるので、事後分布の精度をより向上させることができる。
【００４６】
また、本実施形態では、事前分布の上限値を規定可能な第２の関数において、野生動物の個体群動態に関する数値の上限値に対応するとともに事前分布の上限値において連続的な分布形状を有するように事前分布の上限値が規定された状態で、第１の関数および第２の関数を用いて、野生動物の個体群動態の事前分布を作成するようにＣＰＵ１を構成する。これにより、たとえば、ニホンジカの場合、生態学的な特性および調査結果から得られる妊娠率や産仔数等の情報に基づいて、自然増加率（個体群動態）の上限値（たとえば、１．５）を想定することが可能であり、このような場合に、ニホンジカの自然増加率の上限値に対応する事前分布の上限値を規定可能な第２の関数を用いれば、事前分布をより精度よく設定することができる。
【００４７】
また、本実施形態では、事前分布の上限値において連続的な分布形状を有するように事前分布の上限値が規定された状態の第２の関数を用いて、事前分布の分布形状を非対称形状にするようにＣＰＵ１を構成する。これにより、野生動物の生態学的な特性や調査結果に基づいて、野生動物の個体群動態の上限値に加えて、野生動物の個体群動態の確率密度分布が非対称形状であると想定することが可能な場合に、野生動物の個体群動態の上限値に対応するように事前分布の上限値が規定され、かつ、事前分布の分布形状を非対称形状に変換可能な第２の関数を用いることにより、事前分布をさらに精度よく設定することができる。
【００４８】
また、本実施形態では、野生動物の個体群動態に関する既知のデータが記憶されたハードディスク２から既知のデータを読み出し、読み出した既知のデータと、尖度が調整されることにより作成された事前分布とを用いて、野生動物の個体群動態の事後分布を作成するようにＣＰＵ１を構成する。これにより、ＣＰＵ１により、ハードディスク２に記憶された野生動物の個体群動態に関する既知のデータが自動的に読み出されるので、容易に、読み出した既知のデータと尖度が調整された事前分布とに基づいて事後分布を精度よく作成することができる。
【００４９】
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれる。
【００５０】
たとえば、上記実施形態では、本発明の第１の関数の一例として、期待値、分散および尖度を調整可能な上記式（２）に示す第１の関数を示したが、本発明はこれに限られない。本発明の第１の関数は、期待値、分散および尖度のうち、少なくとも尖度を独立して調整可能であれば、上記式（２）に示す第１の関数以外の関数であってもよい。
【００５１】
また、上記実施形態では、本発明の第２の関数の一例として、上限値を規定する上記式（３）に示す第２の関数を示したが、本発明はこれに限られない。本発明の第２の関数は、下限値を規定する関数や上限値および下限値の両方を規定する関数など、上記式（３）に示す第２の関数以外の関数であってもよい。
【００５２】
また、上記実施形態では、事前分布と既知のデータとを用いて、マルコフ連鎖モンテカルロ法により事後分布を作成（導出）する例を示したが、本発明はこれに限られない。本発明では、事前分布と既知のデータとを用いて、マルコフ連鎖モンテカルロ法以外の方法により事後分布を作成（導出）してもよい。
【符号の説明】
【００５３】
１  ＣＰＵ（事前分布作成手段、事後分布作成手段）
２  ハードディスク（記憶部）
１０ａ  個体群動態推定プログラム
１００  個体群動態推定装置（コンピュータ）
【図１】
0005992369.tif 000011
【図２】
0005992369.tif 000012
【図３】
0005992369.tif 000013
【図４】
0005992369.tif 000014
【図５】
0005992369.tif 000015