[ADMB Users] "xxglobal.tmp" error with tpl2rem
Steve Martell
s.martell at fisheries.ubc.ca
Tue Aug 4 08:35:16 PDT 2009
Ben,
Try using the scripts that Arni posted at:
http://admb-project.org/community/editing-tools/admb-ide/scripts-linux.zip/view?searchterm=Scripts
Steve
On 3-Aug-09, at 10:09 PM, Bolker,Benjamin Michael wrote:
>
> I'm just beginning to use AD Model Builder.
> I have successfully (after some struggles) created a few different
> likelihood models, with results matching previous model fits in R.
> Now I'm trying to add a random effect. Versions of the TPL file
> below compiled OK with tpl2cpp before I added the random effect;
> now, with tpl2rem, I'm getting
>
> cat: xxalloc4.tmp: No such file or directory
> cat: xxalloc5.tmp: No such file or directory
> Error executing command cat xxglobal.tmp xxhtop.tmp header.tmp
> xxalloc1.tmp xxalloc2.tmp xxalloc3.tmp xxalloc4.tmp
> xxalloc5.tmp xxtopm.tmp xxalloc6.tmp > mccoypred4.cpp
>
> There was an issue with similar keywords listed on the ADMB site,
> but it
> was listed as "resolved".
>
> I'm running on Ubuntu Linux 9.04, with a recent (?? can't find
> version number ??)
> version of ADMB. (Maybe 9.0.202 from the google code page?)
>
> I'd be happy to provide any other diagnostics/information if that
> will
> help solve / diagnose ...
>
> sincerely
> Ben Bolker
>
>
>
> -------------------------------------
> DATA_SECTION
>
> init_int nobs // # of observations
> init_int nblock // # of blocks
> init_vector killed(1,nobs) // # killed per trial
> init_vector size(1,nobs) // size of individuals
> init_vector initial(1,nobs) // starting density (# individuals)
> init_matrix Z(1,nobs,1,nblock) // random-effects model matrix
>
> PARAMETER_SECTION
>
> init_number c // *mean* c value
> init_number d
> init_bounded_number h(0.0,1.0,1)
> init_number g
> init_bounded_number sigma_c(0.00001,1.0,1)
> random_effects_vector u(1,nblock)
> vector prob(1,nobs)
> vector cvec(1,nobs)
> objective_function_value f
>
> PROCEDURE_SECTION
> // compute vector of c values
> cvec = c + sigma_c*(Z*u)
> // power-Ricker
> prob = 1/(1/(elem_prod(cvec,elem_prod(pow(size/d,g),exp(-size/d))))
> +h*initial);
> // binomial negative log-likelihood
> f -= sum( log_comb(initial,killed)+
> elem_prod(killed,log(prob))+
> elem_prod(initial-killed,log(1-prob)));
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
Steve Martell
s.martell at fisheries.ubc.ca
More information about the Users
mailing list