さあ!諸君!捕鯨問題だ!

Yahoo! Japan 掲示板トピックビューアー

[ << 最初のページ | < 前のページ | メッセージリスト | 掲示板表示 | [ メッセージ # ] | 次のページ > | 最後のページ >> ]

「水産資源学を語る」解説(5)

投稿者: legal_guardian01 投稿日時: 2012/05/13 12:26 投稿番号: [58139 / 62227]
P(t+1)=P(t)−C(t)+1.4184μP(t){1−(P(t)/P(0))(P(t)/P(0))}

では、計算してみようという話に移る。今、我々は三つの数値を持っているP(t+1)、P(t)、C(t)だ。μとP(0)は定数だが未知数だ。そこで、μとP(0)を適当な値に仮定して、P(t)、C(t)を入れ計算すると、計算されたP(t+1)がでてくる。ここで注意が必要なのは、この数値が予め分かっているP(t+1)とは異なることだ。

原典の「水産資源学を語る」では、この周辺の表記が省略されすぎていて不親切だといわざるを得ない。もう一度いうと、持っている数値P(t)、C(t)と適当に仮定した数値μとP(0)から、P(t+1)を出し、それを現在資源量P(T)としている。(Msg.32124)
P(T)のTが大文字であることに注目してもらいたい。これはP(t)ではないのだ。ここが誤解を生むポイントとなっている。

整理しよう。適当に仮定した数値μとP(0)から、P(T)がでた。これとは別にP(t+1)という数値がある。前者は計算値、後者は観測値だ。このP(T)を使うと、

L(T)=3μ{(P(T)/P(0))−0.54}P(T)

からT年の捕獲限度量L(T)がでる。分かるだろうか、Tとはt+1を示している。

ここでは、重要な説明が省略されている。ここまでの計算で、適当に仮定した数値μとP(0)から、L(T)が出た訳だが、それとともにP(t+1)とP(T)から尤度を計算している。このことを「この値を観測された資源量と比較し、差の小さい場合には仮定したμとP(0)が正しい値である可能性が高いとして重みをを重くし、逆に大きく違っている時には軽くする。」と表現しているわけだが、これだと尤度函数を使って推定しているとは読めない。もともと啓蒙書なので、そこまで要求する読者は想定していないだろうから、仕方ないのかもしれない。だが、統計学を分かっている人ほど、この部分の記述が実際には何をやっているか分からないのだと思う。

尤度(ゆうど)というのは、出てきた数値がどれくらい尤(もっと)もらしいのかを数値化したもの。尤度函数を使って算出する。この本にはそこまで記載していないのでどの尤度函数を使うのか分からないが、報告書などを読むと種々の尤度函数を使って検討しているらしい。

例えれば、コインを3回投げて2回表が出たとき、表が出る確率を2/3としたら、それがどのくらい真値にあっているのか(尤もらしいのか)を示すということだ。勘どころがある人なら、「3回投げて2回表が出たとき」と「30回投げて20回表が出たとき」の尤度が異なることがわかる。

もう一度整理しよう。
適当に仮定した数値μとP(0)から、L(T)と尤度がでた。

次に行なうのは「μとP(0)の組を指定された範囲内に満遍なく配置して各組のL(T)を求め、」だ。最初に仮定したμとP(0)を少しづつ変え、それに対するL(T)と尤度を算出する。

「これらの未知数に対しては、あらかじめ常識的に考えられる範囲を設定しておく。」という意味は、常識的に妥当でない範囲を算出しても良いんだけど、尤度の低いL(T)をいくら出してもその後の計算にその数値は採用されないから無駄だということだ。最近の計算機はパーソナルタイプでも相当な計算が出来るが、P(t+1)、P(t)、C(t)は複数年度分のものがあるわけだし、高速計算できる計算機だって余計な演算は排除したい所だ。

ここで、適当に仮定した数値μとP(0)の複数の組み合わせから、その組み合わせに対応するL(T)と尤度がでたわけだ。次にこの尤度を正規化する。簡単にいえば全ての尤度を足すと1になるようにするということだ。

次に、(μ、P(0)、L(T)、尤度)の組をL(T)の大きさに対してソートする。L(T)の小さい順に尤度を足していく。この尤度の和が40%になったときのL(T)を限度量とする。

この手順の記載が「μとP(0)の組を指定された範囲内に満遍なく配置して各組のL(T)を求め、これを大きさの順に配列し、小さいほうからそれぞれのL(T)に与えられた重みの累積和を計算し、これが重み全体の和の40%くらいになるL(T)を実際の限度量として採用する。」ということになる。
[ << 最初のページ | < 前のページ | メッセージリスト | 掲示板表示 | [ メッセージ # ] | 次のページ > | 最後のページ >> ]

Yahoo! Japan 掲示板 アーカイヴ

[検索ページ] (中東) (東亜) (捕鯨 / 捕鯨詳細)