<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Dries,<div class=""><br class=""></div><div class=""> Thanks for sending this very detailed explanation. I’ll try to answer to each of your points below:</div><div class=""><br class=""><div apple-content-edited="true" class="">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><blockquote type="cite" style="orphans: auto; widows: auto;" class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;">(i) <b class="">Recapture probability</b> is around 10 %. We however did not check for presence/absence of the bird every day after fledging. Usually on the 5th/6th day (range 3-9). This thus underestimates recapture probability dramatically. We noted on the BaSTA website that future versions would include covariates on recapture and recovery probabilities. We could for example add "search effort" as covariate, stating the specific days when we actually looked for that individual. Is there a possibility of trying this out already ?</div></div></blockquote><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><br class=""></div>Fernando: Unfortunately we have not yet been able to implement this.</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><br class=""></div><div style="text-align: start; text-indent: 0px; orphans: 2; widows: 2;" class=""><blockquote type="cite" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: auto; widows: auto;" class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;">(ii) Our <b class="">confidence intervals</b> on the mortality rates just after fledging are big. To that extent that it obscures the real shape of the mortality function when added to the graphs. <a href="https://www.dropbox.com/sh/kk1wztsftfljbw5/AAC4OHKT-pNXepioaR4qf78oa?dl=0" class="">See PlotFancyBaSTA with CI and regular plot without CI</a>. Because we only searched for the juveniles some days after fledging, the first resightings are only done on day 3 and later. Is it possible that the particular sampling strategy is on the origin of this big CI's ? </div></div></blockquote><div class=""><div dir="ltr" class=""><div class="" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12.8px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div><div class="">Fernando: By looking at your plots, particularly the PlotFancy, I suspect that not all your parameters had converged… Did you have a chance at checking the traces?.</div><div class=""><br class=""></div><div class=""><blockquote type="cite" style="orphans: auto; widows: auto;" class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;">(iii) Based on the lowest DIC from running all basic models (no covariates), we selected the WEIBULL model with a <b class="">Makeham</b> shape. Adding the Makeham constant avoids creating immortal individuals when having a declining mortality shape (i.e. when 0 < b0 < 1) by setting a "base mortality". How do we interpret this value when having other shapes (b0 <= 1) ? Why is this in these cases not (partially) embedded within the scale parameter (b1) ?</div></div></blockquote><div class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;"><br class=""></div></div></div></div><div class="">Fernando: Actually embedding the third parameter (c) into b1 does not solve the problem of starting mortality at 0, inherent to the Weibull model, while the intended effect of the Makeham constant is to shift up or down the mortality. If it’s embebed in b1, then you are only changing mortality proportionally.</div><div class=""><br class=""></div><div class=""><blockquote type="cite" style="orphans: auto; widows: auto;" class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;">(iv) The <b class="">serial autocorrelation</b> of most coefficient estimates appears to be rather high (> 0.4). Would you recommend readjusting the thinning argument or other model settings ? Or does this not really effect our estimates and outcome ? </div></div></blockquote></div><div class=""><br class=""></div><div class="">Fernando: Yes, I would recommend running the model for more iterations with a large thinning argument. If you run it in parallel then you can run more chains for less iterations but with higher thinning, that’ll be even better and it would save you quite a bit of time.</div><div class="" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12.8px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div><div class=""><blockquote type="cite" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: auto; widows: auto;" class=""><div dir="ltr" class=""><div class="" style="font-size: 12.8px;">NB : we noticed a discrepancy between the default values in the package versus the ones listed in the BaSTA vignette (niter = 11000, burnin = 1001 and thinning = 20 vs. niter = 50000, burnin = 5001 and thinning = 50). Which settings to you recommend as default values ?</div></div></blockquote><div class=""><div dir="ltr" class=""><div class="" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12.8px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div><div class="" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12.8px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div><div class="">Fernando: Thanks for spotting this! We’ll check the settings.</div><div class="" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12.8px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class=""></div></div></div></div></div></div>Best,</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Fernando<br class="Apple-interchange-newline"><br class=""></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class=""><b class="">Fernando Colchero</b></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Assistant Professor</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Department of Mathematics and Computer Science</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Max-Planck Odense Center on the Biodemography of Aging</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class=""><b class=""><i class=""><br class=""></i></b></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Tlf. +45 65 50 23 24</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Email <a href="mailto:colchero@imada.sdu.dk" class="">colchero@imada.sdu.dk</a></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Web <a href="http://www.sdu.dk/staff/colchero" class="">www.sdu.dk/staff/colchero</a></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Pers. web <a href="http://www.sdu.dk/staff/colchero" class="">www.colchero.com</a></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class="">Adr. Campusvej 55, 5230, Odense, Dk</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-size: 12px; orphans: 2; widows: 2;" class=""><b class="">University of Southern Denmark</b></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On 16 Sep 2015, at 15:13, Dries Van de Loock <<a href="mailto:dries.vandeloock@ugent.be" class="">dries.vandeloock@ugent.be</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class=""><span style="font-size:12.8px" class="">Dear BaSTA community, </span><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">I'm involved in the same project as Diederik Strubbe on post-fledging survival of a tropical passerine (see <a href="http://lists.r-forge.r-project.org/pipermail/basta-users/2015-August/thread.html" target="_blank" class="">archive</a>). </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Nests were monitored until fledging and juveniles were relocated every 5 to 6 days after fledging. </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Our input matrix is formatted accordingly : the birthyear is considered the point of fledging and is set to "1" for all individuals. The resighting matrix starts at day 1 and is "0" for all individuals on that day. Day 2 in the matrix is thus the first day after fledging. <a href="https://www.dropbox.com/sh/kk1wztsftfljbw5/AAC4OHKT-pNXepioaR4qf78oa?dl=0" class="">See screenshot</a> from RStudio. </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">When running the models, we bumped into following questions :</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">(i) <b class="">Recapture probability</b> is around 10 %. We however did not check for presence/absence of the bird every day after fledging. Usually on the 5th/6th day (range 3-9). This thus underestimates recapture probability dramatically. We noted on the BaSTA website that future versions would include covariates on recapture and recovery probabilities. We could for example add "search effort" as covariate, stating the specific days when we actually looked for that individual. Is there a possibility of trying this out already ?</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">(ii) Our <b class="">confidence intervals</b> on the mortality rates just after fledging are big. To that extent that it obscures the real shape of the mortality function when added to the graphs. <a href="https://www.dropbox.com/sh/kk1wztsftfljbw5/AAC4OHKT-pNXepioaR4qf78oa?dl=0" class="">See PlotFancyBaSTA with CI and regular plot without CI</a>. Because we only searched for the juveniles some days after fledging, the first resightings are only done on day 3 and later. Is it possible that the particular sampling strategy is on the origin of this big CI's ? </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">(iii) Based on the lowest DIC from running all basic models (no covariates), we selected the WEIBULL model with a <b class="">Makeham</b> shape. Adding the Makeham constant avoids creating immortal individuals when having a declining mortality shape (i.e. when 0 < b0 < 1) by setting a "base mortality". How do we interpret this value when having other shapes (b0 <= 1) ? Why is this in these cases not (partially) embedded within the scale parameter (b1) ?</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">(iv) The <b class="">serial autocorrelation</b> of most coefficient estimates appears to be rather high (> 0.4). Would you recommend readjusting the thinning argument or other model settings ? Or does this not really effect our estimates and outcome ? </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">NB : we noticed a discrepancy between the default values in the package versus the ones listed in the BaSTA vignette (niter = 11000, burnin = 1001 and thinning = 20 vs. niter = 50000, burnin = 5001 and thinning = 50). Which settings to you recommend as default values ?</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">We performed all analysis using BaSTA 1.9.4</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Your feedback and suggestions are greatly appreciated, </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Best wishes, </div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Dries</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">-----</div><div style="font-size:12.8px" class=""><div class="">Dries Van de Loock</div><div class="">PhD student</div><div class="">Terrestrial Ecology Unit</div><div class="">Department of Biology</div><div class="">Ghent University</div><div class="">KL Ledeganckstraat 35</div><div class="">B-9000 Ghent, Belgium</div><div class="">Phone: <a href="tel:%2B32%20%280%299%20265%2050%2039" value="+3292655039" target="_blank" class="">+32 (0)9 265 50 39</a></div><div class=""><a href="http://www.ecology.ugent.be/terec/" target="_blank" class="">http://www.ecology.ugent.be/terec/</a></div></div></div>
_______________________________________________<br class="">Basta-users mailing list<br class=""><a href="mailto:Basta-users@lists.r-forge.r-project.org" class="">Basta-users@lists.r-forge.r-project.org</a><br class="">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/basta-users<br class=""></div></blockquote></div><br class=""></div></body></html>