<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Dear Basta,</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>I get the following error when I run this script:</div><div><br></div><div> Fem1978only <- basta(object = Males.dat$newData, studyStart = 1978, studyEnd = 1982, model="GO" ,burnin=1000 ,niter=10000, nsim=4,parallel=TRUE, ncpus=4)</div><div><br></div><div><div>Survival parameters converged appropriately. </div><div>DIC was calculated</div><div>.Warning messages:</div><div>1: In max(DeathAges, na.rm = TRUE) :</div><div> no non-missing arguments to max; returning -Inf</div><div>2: In max(wa) : no non-missing arguments to max; returning -Inf</div><div>3: In MakeLifeTable(x, ax = 0.5, n = 1) : NAs introduced by coercion</div></div><div>Full script and errors below</div><div><br></div><div>This prevented me from using the out$lifeTable function, but was able to use the</div><div>out$Xq function, and count the 50% deaths.</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span> I have 943 individuals of unknown age studied with the Gompertz, all marked in 1978, and followed till 1982</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>Of these 773 were only seen the first year of the study (1978), yet basta estimates the first deaths at 3 years of age </div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>The qx for those first deaths at 3 years old is is extremely high, but the qx for year 2 is zero.</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>I interpret this absence of deaths in year 2 as being a way to fit the Gompertz mortality curve to the data, starting at 0 mortlity, and going up. If the deaths are reported in year 3, then year 2 has 0 mortality, and goes up rapidly (the first qx is 0.95)</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>This seems incorrect to me: 773 were only seen in year 1, some must have died in year 2. As you can see the pi is 0.94, so you wouldn't expect such a long lag time to report the first deaths. </div><div><br></div><div>
<table border="0" cellpadding="0" cellspacing="0" width="520" style="border-collapse: collapse; width: 520pt; position: static; z-index: auto; ">
<!--StartFragment-->
<colgroup><col width="65" span="8" style="width:65pt">
</colgroup><tbody><tr height="15" style="height:15.0pt">
<td height="15" colspan="7" width="455" style="height:15.0pt;mso-ignore:colspan;
width:455pt"><span style="mso-spacerun:yes"> </span>Estimate<span style="mso-spacerun:yes"> </span>StdErr Lower95%CI Upper95%CI SerAutocor
UpdateRate PotScaleReduc</td>
<td width="65" style="width:65pt"></td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">b0</td>
<td align="right">-2.3926</td>
<td align="right">0.183197</td>
<td align="right">-2.7062</td>
<td align="right">-2.0079</td>
<td class="xl63">0.62626</td>
<td class="xl63">0.3755</td>
<td class="xl63">1.001</td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">b1</td>
<td align="right">0.6992</td>
<td align="right">0.052093</td>
<td align="right">0.5881</td>
<td align="right">0.7925</td>
<td class="xl64">0.60219</td>
<td class="xl65">0.3755</td>
<td class="xl66">1</td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">pi</td>
<td align="right">0.9438</td>
<td align="right">0.009978</td>
<td align="right">0.9223</td>
<td align="right">0.9617</td>
<td class="xl65">0.04917</td>
<td class="xl65">1</td>
<td class="xl66">1.001</td>
</tr>
<!--EndFragment-->
</tbody></table></div><div><br></div><div><br></div><div><br></div><div>When I run the individuals first marked in 1979 in the same population and followed through to 1982, I get a similar pattern of first deaths at age 3, even though 660 of the total 768 individuls were only seen in 1979 and never again, also with extremely high "pi" values of 94%</div><div><br></div><div>
<table border="0" cellpadding="0" cellspacing="0" width="520" style="border-collapse: collapse; width: 520pt; position: static; z-index: auto; ">
<!--StartFragment-->
<colgroup><col width="65" span="8" style="width:65pt">
</colgroup><tbody><tr height="15" style="height:15.0pt">
<td height="15" colspan="7" width="455" style="height:15.0pt;mso-ignore:colspan;
width:455pt"><span style="mso-spacerun:yes"> </span>Estimate<span style="mso-spacerun:yes"> </span>StdErr Lower95%CI Upper95%CI SerAutocor
UpdateRate PotScaleReduc</td>
<td width="65" style="width:65pt"></td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">b0</td>
<td align="right">-2.3926</td>
<td align="right">0.183197</td>
<td align="right">-2.7062</td>
<td align="right">-2.0079</td>
<td class="xl63">0.62626</td>
<td class="xl63">0.3755</td>
<td class="xl63">1.001</td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">b1</td>
<td align="right">0.6992</td>
<td align="right">0.052093</td>
<td align="right">0.5881</td>
<td align="right">0.7925</td>
<td class="xl64">0.60219</td>
<td class="xl65">0.3755</td>
<td class="xl66">1</td>
</tr>
<tr height="15" style="height:15.0pt">
<td height="15" style="height:15.0pt">pi</td>
<td align="right">0.9438</td>
<td align="right">0.00991978</td>
<td align="right">0.9223</td>
<td align="right">0.9617</td>
<td class="xl65">0.04917</td>
<td class="xl65">1</td>
<td class="xl66">1.001</td>
</tr>
<!--EndFragment-->
</tbody></table><br></div><div><br></div><div><br></div><div>Finally, when I combine the 1978 and 1979 marked individuals, all hell breaks loose. No errors are reported, but the </div><div><br></div><div>Fem78and79$lifeTable command only reveals the life table for the 768 individuals first observed in 1979, even though basta reports that all 1711 individuals were read, and I can see all the ages at death when I use the Fem78and79$Xq command</div><div><br></div><div>This gives individuals marked in 1978 their first deaths at age 3, and individuals marked in 1979 their first deaths at age 2.</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>I know this is complicated, but this is a dataset you may want to take a look at, as I consider the output difficult to interpret.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div style="font-size: 12px; "><span class="Apple-tab-span" style="white-space: pre; "> </span>All the best!</div><div style="font-size: 12px; "><br></div><div style="font-size: 12px; "> <span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre; "> </span>Cliff Cunningham</div><div style="font-size: 12px; "><br></div><div style="font-size: 12px; "><br class="Apple-interchange-newline"> Professor</div><div style="font-size: 12px; "><span class="Apple-tab-span" style="white-space: pre; "> </span>Biology Department</div><div style="font-size: 12px; "> <span class="Apple-tab-span" style="white-space: pre; "> </span>Duke University</div><div style="font-size: 12px; "> 130 Science Drive</div><div style="font-size: 12px; "><span class="Apple-tab-span" style="white-space: pre; "> </span>Durham NC 27708<br class="webkit-block-placeholder"></div><div style="font-size: 12px; "><span class="Apple-tab-span" style="white-space: pre; "> </span>Phone 919-660-7356<br></div><div style="font-size: 12px; "><br></div><div style="font-size: 12px; "><a href="http://www.biology.duke.edu/cunningham">http://www.biology.duke.edu/cunningham</a></div><div style="font-size: 12px; "><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br></body></html>