[ADMB Users] multinomial weights and -gh
Jeff Laake
jefflaake at gmail.com
Fri Mar 1 15:29:37 PST 2013
So I'm confused. In the binomial example (binomial_w2.tpl) it shows
multiplication of weight and likelihood value but in the distance sampling
example I did with Hans the weighting was all behind the scenes to use -gh.
Can someone hep me with the example cjsre.tpl at:
https://docs.google.com/folder/d/0B77g1ScdUwVeWkNDaFA5OVp2SjA/edit?usp=sharing
which I modeled after what Hans helped me with. You'll see code
multiplying weights that is commented out.
The weights here are frequencies of the occurrence of a particular capture
history.
If you rename the file yes_cjsre.dat to cjsre.dat and run it then it will
work giving a likelihood about 333. But if you rename the example
no_cjsre.dat to cjsre.dat and run it, it will give a value of about 110.
Both should be in the 300s. The difference is that in the no_cjsre.dat
there are two random effects per call whereas in yes_cjsre.dat it is on
only one of the parameters. Clearly the weights aren't being applied
properly when there are random effects on both parameters.
I found the following comment in the code df1bg2h
// !!! This only works for one random efect in each separable call
// at present.
Am I out of luck at present?
--jeff
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20130301/b6b407bd/attachment.html>
More information about the Users
mailing list