Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bug in example speed parameterization #25

Open
camrinbraun opened this issue Jul 7, 2020 · 1 comment
Open

Bug in example speed parameterization #25

camrinbraun opened this issue Jul 7, 2020 · 1 comment
Labels

Comments

@camrinbraun
Copy link
Owner

camrinbraun commented Jul 7, 2020

From @fmasbervejillo here

In @paulgatti example code:

Line 320:
# pars sigma1 (resident), sigma2 (migrant), p11 switch from behavior 1 to 1 (remain resident), p22

Line 321 you define these initial parameters:
pars.init=c(2,.2,.6,.8)

My guess is that the first parameter is, following the logic of the former method, the migration speed (not the resident speed), and the second parameter is the fraction of that speed that is attributed to the resident speed (hence the .2?). Is this correct? otherwise, the resident speed would be 10 fold the migratory speed.

@camrinbraun camrinbraun added the bug label Jul 7, 2020
@camrinbraun
Copy link
Owner Author

from @paulgatti:
The idea here is to estimate the two movement parameters. I can’t recall if HMMoce forces sigma1=10*sigma2 or if it is a starting guess. In case, the intention is indeed to force that, only one (sigma1) should be estimated as the two are not independant. In fact, the code could be easily changed to estimate only 3 parameter (sigma1 and the transition probability) and thus likely save some computation time. Anyway, I am still wondering if it is better to estimate or fix swimming speed. In the case where it is possible to have an empirical estimate of swimming speed, I would likely use it (as @camrinbraun did with sharks equipped with spot tags) and avoid the tricky issue of parameter estimation which can be time consuming and I would rather test the sensitivity by altering a bit this value. Also, for my own case, I am now using a very common trick in parameter optimization, which consist in transforming the parameters within plus and minus infinity. This usually help to avoid being stuck in a local optimum (it is not a guarantee but it can help). I can provide you with that (if you are interested).

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant