[ADMB Users] making a pin file
dave fournier
otter at otter-rsch.com
Wed Aug 5 09:00:49 PDT 2009
the easiest way to make a pin file is to run the model with the command
line option
-maxfn 0
and copy the output par file (*.p01) in this case to *.pin and edit it.
On another note your random effects will cause trouble
because cvec can be <0
// compute vector of c values
▒
cvec = c + sigma_c*(Z*u);
▒
// power-Ricker
I think you want something like exp(cvec)
or the admb function mfexp(cvec) which is bounded above for very large
values to help avoid overflow.
▒
prob =
1/(1/(elem_prod(mfexp(cvec),elem_prod(pow(size/d,g),exp(-size/d))))+h*init▒
▒
--
David A. Fournier
P.O. Box 2040,
Sidney, B.C. V8l 3S3
Canada
Phone/FAX 250-655-3364
http://otter-rsch.com
More information about the Users
mailing list