specify.Mesnil.gamma: 再生産関係推定でSR=HSにおいて同一尤度を持つ複数のパラメータセットがえられたとき、SR=Mesnilとしてcheck.SRfitを呼び出し、gammaを変えながら尤度が一意に定まる値を探す関数

View source: R/stock_recruit.r

specify.Mesnil.gammaR Documentation

再生産関係推定でSR=HSにおいて同一尤度を持つ複数のパラメータセットがえられたとき、SR=Mesnilとしてcheck.SRfitを呼び出し、gammaを変えながら尤度が一意に定まる値を探す関数

Description

再生産関係推定でSR=HSにおいて同一尤度を持つ複数のパラメータセットがえられたとき、SR=Mesnilとしてcheck.SRfitを呼び出し、gammaを変えながら尤度が一意に定まる値を探す関数

Usage

specify.Mesnil.gamma(resSR, n = 100, seed = 1, sigma = 5, gamma_ini = 10)

Arguments

resSR

fit.SRfit.SRregimeのオブジェクトで、かつSR=Mesnilのもの

n

初期値を変えてパラメータ推定する回数 (check.SRfitの引数)

seed

set.seedで使用するseed(check.SRfitの引数)

sigma

初期値を変えるときの生起乱数の標準偏差(check.SRfitの引数)

gamma_ini

SR=Mesnilの曲がり具合を決めるパラメータ デフォルトは10で、9,8,…,2,1,0.9,0.8,…と複数のパラメータが得られるまで探していく。


ichimomo/frasyr documentation built on May 3, 2024, 1:30 a.m.