[Basta-users] out$lifetable errors

Owen Jones jonesor at gmail.com
Mon Nov 7 11:18:20 CET 2011


Hi Cliff,

The immediate problem appears to be a problem in the MakeLifeTable function, which BaSTA calls to produce a life table. The error message is saying that all of the DeathAges that are predicted by basta, and sent to the MakeLifeTable function, are "NA". Therefore, it appears that basta is somehow failing to predict any ages at death.

What does the rest of "out" look like? (e.g. summary(out))
Does it allow you to produce plots of traces/trajectories?

Cheers,
Owen






On 7 Nov 2011, at 11:03, Cliff Cunningham wrote:

> Dear BaSTa,
> 	After several runs I have typed  out$lifetable and got the following
> 
> > out$lifetable
> NULL
> 
> 	I should report I got several warning messages (see bottom) when I ran the following: 
> Survival parameters converged appropriately. 
> DIC was calculated
> .Warning messages:
> 1: In max(DeathAges, na.rm = TRUE) :
>   no non-missing arguments to max; returning -Inf
> 2: In max(wa) : no non-missing arguments to max; returning -Inf
> 3: In MakeLifeTable(x, ax = 0.5, n = 1) : NAs introduced by coercion
> 
> 
> My code was as follows, 
> 
> out <- basta(object = Males.dat$newData, studyStart = 1979, studyEnd = 1982, burnin=1000 ,niter=10000, nsim=4, parallel=TRUE, ncpus=4)
> 
> 	In this one I had all the birthdays at 1978, and first captures in 1979.
> 
> 
> 
> 
> Here are the errors from the 
> Multiple simulations finished.
> Total MCMC computing time: 3.91 minutes
> 
> Thanks so much!!
> 
>    	Cliff Cunningham
> 
> 
>         Professor
> 	Biology Department
>     	Duke University
>         130 Science Drive
> 	Durham NC 27708
> 	Phone 919-660-7356
> 
> http://www.biology.duke.edu/cunningham
> 
> 
> 
> 
> _______________________________________________
> Basta-users mailing list
> Basta-users at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/basta-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/basta-users/attachments/20111107/16a9b937/attachment.htm>


More information about the Basta-users mailing list