[ADMB Users] Reading in random effects start values from .pin file
Ben Stevenson
bcs5 at st-andrews.ac.uk
Wed Nov 6 04:31:32 PST 2013
Hi ADMBers,
A while ago I noticed that all elements of a random_effects_vector seem to
start out at 0, regardless of what appears in the .pin file. I have
recently discovered that using the -noinit option results in behaviour that
one would expect, i.e., random effects start values are read in from the
.pin file.
However, when you have a SEPARABLE_FUNCTION, this persists even with the
-noinit option.
Hans suggested a nice workaround, whereby you read in a vector of offsets
(which can be set to your desired start values) from the .dat file, and use
v as your random effect, where
v = u + offset;
and u is the random_effects_vector which has elements with start value 0.
Can anyone shed any light on why start values in the .pin file for a
random_effects_vector do seem to do what one would expect? The ADMB-RE
manual suggests the .pin file initialisation should work.
I have attached a minimum working example of the behaviour I am seeing. The
compiled executable won't optimise anything, but start values for the
random effects are printed out. I get a bunch of 0s, even though the .pin
file provides the start values {1, 2, 3, 4, 5}.
Cheers,
Ben
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20131106/6edf7667/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.dat
Type: application/octet-stream
Size: 14 bytes
Desc: not available
URL: <http://lists.admb-project.org/pipermail/users/attachments/20131106/6edf7667/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.pin
Type: application/octet-stream
Size: 21 bytes
Desc: not available
URL: <http://lists.admb-project.org/pipermail/users/attachments/20131106/6edf7667/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.tpl
Type: application/vnd.groove-tool-template
Size: 316 bytes
Desc: not available
URL: <http://lists.admb-project.org/pipermail/users/attachments/20131106/6edf7667/attachment.bin>
More information about the Users
mailing list