<div dir="ltr"><div>Fernando,</div><div>Thanks very much for your help! Running the simulations definitely helps to learn about the models. And yes, your explanation of the Gompertz bathtub clarifies what the parameters represent.</div><div>Thanks,</div><div>Ken</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 2, 2018 at 3:55 AM, Fernando Colchero <span dir="ltr"><<a href="mailto:colchero@imada.sdu.dk" target="_blank">colchero@imada.sdu.dk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="-ms-word-wrap: break-word;">Hi Ken,<div><br></div><div> Sorry (again) for the late reply. I attached a code to simulate BaSTA like data for a Gompertz model, let me know if this helps. </div><div><br></div><div> About the interpretation of the parameters, here is a summary:</div><div><br></div><div> exp(a0) + c + exp(b0) provides the baseline mortality at age x, commonly it is dominated by exp(a0) + c. The parameter a1 is always preceded by a negative sign and it controls the speed at which infant mortality declines. The b1 parameter controls the speed at which adult mortality increases, so it can give you a measure of the ageing rate. Let me know if this clarifies a bit.<br><div>
<div style="color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><br class="m_-5567661300835477783Apple-interchange-newline">Best,</div><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><br></div><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"> Fernando</div><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"></div></div></div></div></div><br><div style="-ms-word-wrap: break-word;"><div><div><div style="color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"></div><div style="font:18px/normal Helvetica;color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="font-size:12px"><b>Fernando Colchero</b></div><div style="font-size:12px">Associate Professor</div><div style="font-size:12px">Department of Mathematics and Computer Science</div><div style="font-size:12px">Max-Planck Odense Center on the Biodemography of Aging</div><div style="font-size:12px"><b><i><br></i></b></div><div style="font-size:12px">Tlf. <a href="tel:+45%2065%2050%2023%2024" target="_blank" value="+4565502324">+45 65 50 23 24</a></div><div style="font-size:12px">Email <a href="mailto:colchero@imada.sdu.dk" target="_blank">colchero@imada.sdu.dk</a></div><div style="font-size:12px">Web <a href="http://www.sdu.dk/staff/colchero" target="_blank">www.sdu.dk/staff/colchero</a></div><div style="font-size:12px">Pers. web <a href="http://www.sdu.dk/staff/colchero" target="_blank">www.colchero.com</a></div><div style="font-size:12px">Adr. Campusvej 55, 5230, Odense, Dk</div><div style="font-size:12px"><br></div><div style="font-size:12px"><b>University of Southern Denmark</b></div></div></div>
</div>
<br><div><blockquote type="cite"><div>On 27 Feb 2018, at 18:49, Ken Honeycutt <<a href="mailto:rkenhoneycutt@gmail.com" target="_blank">rkenhoneycutt@gmail.com</a>> wrote:</div><br class="m_7112313900607288793Apple-interchange-newline"><div><div dir="ltr"><div>Fernando,</div><div><br></div><div>Thanks for your response and for sending along the latest BaSTA version, especially the bit about using "recaptTrans." I tried that and it worked well. </div><div><br></div><div><span>I do have a couple more questions.</span></div><div><span> </span></div><div><span>First, is your R-code to simulate populations (as in your 2012 J of Animal Ecology paper) publicly available somewhere? I'm thinking simulating populations would help with confidence in my results.</span></div><div><span><br></span></div><div><span>Second, I'm having a little trouble interpreting the meanings of parameters from the Gompertz bathtub model (I believe there are alpha0, alpha1, beta0, beta1, and c). I see these are described mathematically in the BaSTA MME supporting information, but I'm having trouble visualizing the meanings. Specifically, I'm wondering at what values would these parameters represent a bathtub shape for mortality.</span></div><div><span><br></span></div><div><span>Thanks,</span></div><div><span>Ken</span></div><span><p><br></p><p><br></p><p><br></p><p><br></p><span><p><br></p><span><p><br></p><span><p><br></p><div><br></div></span></span></span></span><div><span><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Feb 24, 2018 at 12:58 AM, Fernando Colchero <span dir="ltr"><<a href="mailto:colchero@imada.sdu.dk" target="_blank">colchero@imada.sdu.dk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div>Hi Ken,<div><br></div><div> Sorry for the late reply. You can use months, just bearing in mind that the parameters correspond then to the distribution of ages at death in months, not in years. There is an argument that you can use for the recapture probability changes during the months you didn’t sample, which is “recaptTrans”. That argument takes a vector of times during which recapture probabilities have changed. </div><div><br></div><div> About the lack of information on birth and death, in principle BaSTA can handle it, but of course, the less information, the higher the uncertainty in the estimation. We have tested the package on datasets without birth and death information and it works, but your setup is quite particular… Please let us know if you need more assistance. Also, I would recommend that you install the latest version of BaSTA. It is not in cran, but you can donwnload it from:</div><div><br></div><div> <a href="https://github.com/fercol/BaSTA/blob/master/latestBuild/BaSTA_1.9.5.tar.gz" target="_blank">https://github.com/fercol/Ba<wbr>STA/blob/master/latestBuild/Ba<wbr>STA_1.9.5.tar.gz</a></div><div><div>
<div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="font:16px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><br class="m_7112313900607288793m_2390374109760895622Apple-interchange-newline">Best,</div><div style="font:16px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><br></div><div style="font:16px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Fernando<br class="m_7112313900607288793m_2390374109760895622Apple-interchange-newline">______________________________<wbr>___</div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><b><br></b></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><b>Fernando Colchero</b></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Associate Professor</div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Department of Mathematics and Computer Science</div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Biodemography Unit</div><div style="text-transform:none;line-height:normal;text-indent:0px;letter-spacing:normal;font-family:Helvetica;font-size:12px;font-variant:normal;word-spacing:0px;white-space:normal"><br></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Tlf. <a href="tel:+45%2065%2050%2023%2024" target="_blank" value="+4565502324">+45 65 50 23 24</a></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Email <a href="mailto:colchero@imada.sdu.dk" target="_blank">colchero@imada.sdu.dk</a></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Web <a href="http://www.sdu.dk/staff/colchero" target="_blank">www.sdu.dk/staff/colchero</a></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Pers. web <span class="m_7112313900607288793m_2390374109760895622Apple-converted-space"> </span><a href="http://www.sdu.dk/staff/colchero" target="_blank">www.colchero.com</a></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal">Adr. <a href="https://maps.google.com/?q=Campusvej+55,+5230,+Odense,+Dk&entry=gmail&source=g" target="_blank">Campusvej 55, 5230, Odense, Dk</a></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><br></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><b>University of Southern Denmark</b></div><div style="font:12px/normal Helvetica;text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:normal"><div style="font-size:16px">______________________________<wbr>___</div><div><br></div></div></div></div></div></div></div></div><br class="m_7112313900607288793m_2390374109760895622Apple-interchange-newline"><br class="m_7112313900607288793m_2390374109760895622Apple-interchange-newline">
</div>
<br><div><blockquote type="cite"><div><div class="m_7112313900607288793h5"><div>On 23 Feb 2018, at 22:36, Ken Honeycutt <<a href="mailto:rkenhoneycutt@gmail.com" target="_blank">rkenhoneycutt@gmail.com</a>> wrote:</div><br class="m_7112313900607288793m_2390374109760895622Apple-interchange-newline"></div></div><div><div><div class="m_7112313900607288793h5"><div dir="ltr">Hi BaSTA users,<div><br></div><div>I just noticed the partial capture history I posted for example did not appear as I expected; and the text continues after many many rows of 0s and 1s. Here's the capture history again in csv. Hopefully this is better.</div><div><br></div><div>Thanks,</div><div>Ken</div><div><br></div><div><div>ID,Birth,Death,6,7,8,9,...,17,<wbr>18,19,20</div><div>1,0,0,1,0,1,0,0,0,1,1,1</div><div>2,0,0,0,1,0,0,0,0,0,1,0</div><div>3,0,0,1,0,1,0,0,0,0,0,1</div><div>4,0,0,0,1,0,0,0,0,1,0,1</div><div>…,0,0,1,0,1,0,0,0,0,1,1</div><div>1544,0,0,0,1,0,0,0,0,0,1,0</div></div></div></div></div>
______________________________<wbr>_________________<br>Basta-users mailing list<br><a href="mailto:Basta-users@lists.r-forge.r-project.org" target="_blank">Basta-users@lists.r-forge.r-pr<wbr>oject.org</a><br><a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/basta-users" target="_blank">https://lists.r-forge.r-projec<wbr>t.org/cgi-bin/mailman/listinfo<wbr>/basta-users</a><br></div></blockquote></div><br></div></div></blockquote></div><br></div>
</div></blockquote></div><br></div></div><br></blockquote></div><br></div>