[ADMB Users] Problems with RE model, using FUNCTIONs for integration

H. Skaug hskaug at gmail.com
Fri Nov 30 09:43:31 PST 2012


I recommend staying away from calling adromb() and any FUNCTION
at the current stage because it is under development. This
is especially so when you are new to ADMB. I have included
your trapz code directly into the PROCEDURE SECTION. That
is the way to go for now.

Note my print out
cout << "trapz - exact=" << Ht-ht*t(i,j)/r << endl;
which compares which stays zero for a long time, but then strange things
happen, resulting from some instability in you integration code.

My advice remains: develop the joint model with the analytical Weibull function.
As it stands now, you have dug into some numerical integration detail,
without checking that the model as a whole will work.

Your main argument against the analytical Weibull function is time varying
covariates in exp(). I am sure they are available at discrete time points,
in which case exp() can be linearly interpolated, and multiplied
with the weibull baseline there still is an analytical solution.


On Fri, Nov 30, 2012 at 5:33 PM, Mark Wolters <kramsretlow at gmail.com> wrote:
> Hi Jeff (and Dave),
> Yes, I followed that discussion with great interest as it happened,
> and I tested out the .tpl file that was posted (even commenting out
> what I thought were the redundant initialization lines).  I got the
> same results that you did--the best I could get it to do was to
> compile with warnings. In any event I didn't feel that any solution
> requiring post-facto modification of the .cpp file was ready for prime
> time (read: for use by a novice like me). So I was hoping to keep it
> simple and find a workaround.
> At the moment I'm teetering on the edge of dropping this work entirely
> and moving on to something else.  I'm very attracted by the potential
> of ADMB, but not sure if I have the time to take the deep dive into it
> that seems required.
> Thanks for your input,
> Mark
-------------- next part --------------
A non-text attachment was scrubbed...
Name: kidneyB.tpl
Type: application/octet-stream
Size: 2790 bytes
Desc: not available
URL: <http://lists.admb-project.org/pipermail/users/attachments/20121130/10ae9a0f/attachment.obj>

More information about the Users mailing list