[Basta-users] simulation problem

[Colchero, Fernando]

Hi Felix,

   I think the best bet is that you get the BaSTA version on R forge. You can type on the R console:

   install.packages("BaSTA", repos = "http://R-Forge.R-project.org", type = "source")

   From the error message you're getting, I think what's happening is that the starting parameters were producing 0 for some individuals in the pdf of ages at death. With this new version that shouldn't happen. Please, let us know if this solves it, if not, we can try to reproduce the error message and find the bug.

   Best,

   Fernando

______________________________________________
 
Fernando Colchero

Max Planck Institute for Demographic Research
Konrad Zuse Str. 1, 18057, Rostock, Germany
Tel:  +49 (0) 381 2081 113
Fax: +49 (0) 381 2081 529
E-mail: Colchero at demogr.mpg.de
______________________________________________





On Feb 9, 2012, at 7:33 PM, Felix Zajitschek wrote:

> Dear Basta team,
> 
> I tried to use some of my own data, but couldn't get the simulation to run. I don't have a reproducible example right now, but I guess the error I get might stem from some basic error in how I set up my data for basta. DataCheck gives me correct values, but there seems to be a problem because one or more ages are estimated as being zero (see below). All individuals have at least two captures (i.e. at least a first capture and one recapture), some of these are on two consecutive days. Could that be a problem?
> 
> Best regards,
> Felix
> 
> > check <- DataCheck(btf.use5, studyStart = 36014, studyEnd = 39636, autofix = rep(1,7), silent = F)
> No problems were detected with the data.
> 
> *DataSummary*
> - Number of individuals         =     479
> - Number with known birth year  =       0
> - Number with known death year  =       0
> - Number with known birth
> AND death years                =       0
> 
> - Total number of detections
> in recapture matrix            =    1,163
> 
> - Earliest detection time       =   36014
> - Latest detection time         =   39636
> 
> 
> out <- basta(object = btf.use5, studyStart = 36014, studyEnd = 39636,
> + model = "GO", shape = "simple")
> No problems were detected with the data.
> 
> Simulation started...
> 
> MCMC stopped at step 3
> Pdf of ages at death equal to 0 for some individuals.
> Reduce jumps to avoid pdf of ages at death equal to 0.
> 
> -- 
> ____________________________________________
> Dr Felix Zajitschek
> Department of Animal Ecology
> Evolutionary Biology Centre
> Uppsala University
> Norbyvägen 18D
> Uppsala, 752 36
> Sweden
> Tel:   +46 18 471 2613
> Fax:   +46 18 471 6484
> 
> Homepage:
> felix.zajitschek.net
> 
> Ageing Research Group page:
> http://www.ebc.uu.se/Research/groups/ageing/
> ____________________________________________
> 
> _______________________________________________
> 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

----------
This mail has been sent through the MPI for Demographic Research.  Should you receive a mail that is apparently from a MPI user without this text displayed, then the address has most likely been faked. If you are uncertain about the validity of this message, please check the mail header or ask your system administrator for assistance.