<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Dear Elena,</p>
<p>the function stochGrowthRateSampleList only returns the mean of
the per-time-step population growth rates in a stochastic
simulation. If you want to know the distribution of these
per-time-step population growth rates you could write your own
version of this function and make sure that that distribution is
stored as well. On the other hand if you want a CI for the
resulting mean stochastic population growth rate, you code write
code in which stochGrowthRateSampleList is called multiple times,
and analyse the resulting distribution of mean stochastic
population growth rates. If other people have other solutions,
please share them! <br>
</p>
<p>Please also note the description of the stochGrowthRateSampleList
function, that reads "Note that the function stoch.growth.rate in
the package popbio does this more efficiently and with more useful
output; but may fail for large bin numbers.".<br>
</p>
<p>I hope this helps,</p>
<p>Eelke<br>
</p>
<p><br>
</p>
<br>
<div class="moz-cite-prefix">On 20-Jul-16 1:48 PM, Elena S. wrote:<br>
</div>
<blockquote cite="mid:DUB127-W831364287C23FD15264FC3A0080@phx.gbl"
type="cite">
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style>
<div dir="ltr"><font style="font-size:12pt;"
face="Calibri,sans-serif" color="#000000" size="3"><br
id="FontBreak">
</font>Hi,<br>
<br>
I try to calculate CI (confidence intervals) for my stochastic
population growth rate. I started with 4 IPM (four years of
study), then I run the following function:<br>
<br>
exp(stochGrowthRateSampleList(list(IPM9,IPM10,IPM11,IPM12),nRunIn = 100,
tMax = 5000))<br>
<br>
or the following ones, based on the IPMpack manual (Metcalf et
al. 2014):<br>
<br>
exp(stochGrowthRateSampleList(nRunIn=100,tMax=5000,listIPMmatrix=list(IPM9,IPM10,IPM11,IPM12),<br>
listPmatrix=list(Pmatrix9,Pmatrix10,Pmatrix11,Pmatrix12), <br>
listFmatrix=list(Fmatrix9,Fmatrix10,Fmatrix11,Fmatrix12),<br>
seedList=NULL,<br>
densDep=FALSE))<br>
<br>
and I obtained a stochastic lambda of 1.25, but, at the moment,
I don't know how to calculate the confidence intervals.<br>
<br>
Anyone can help me?<br>
<br>
Thanks in advance,<br>
<br>
best regards,<br>
<br>
Elena<br>
<br>
<br>
<br>
_________________________________________________________<br>
Dott.ssa Elena Sulis, PhD<br>
<font style="font-size:10pt;" size="2">Centro Conservazione
Biodiversitą (CCB)</font><br>
<font style="font-size:10pt;" size="2">Sezione Botanica -
Dipartimento di Scienze della Vita e dell'Ambiente (DISVA)</font><font
style="font-size:10pt;" size="2"><br>
</font><font style="font-size:10pt;" size="2">Universitą degli
Studi di Cagliari</font><font style="font-size:10pt;" size="2"><br>
</font><font style="font-size:10pt;" size="2">V.le Sant'Ignazio
da Laconi, 11-13</font><font style="font-size:10pt;" size="2"><br>
</font><font style="font-size:10pt;" size="2">09123 Cagliari
(ITALIA)</font><font style="font-size:10pt;" size="2"><br>
</font><font style="font-size:10pt;" size="2">Tel:
+39-0706753509</font><br>
<i>E-mail:</i> <a class="moz-txt-link-abbreviated" href="mailto:elenasulis@hotmail.it">elenasulis@hotmail.it</a> <br>
<i>Skype:</i> elena.sulis<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
IPMpack-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:IPMpack-users@lists.r-forge.r-project.org">IPMpack-users@lists.r-forge.r-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/ipmpack-users">http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/ipmpack-users</a>
</pre>
</blockquote>
<br>
</body>
</html>