<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;">
<div dir="ltr"><font style="font-size:12pt" face="Calibri,sans-serif" size="3" color="#000000">
<div>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB">Good morning IPMusers,</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB"><br>
thank you Dr. Jongejans for answered me so quickly, and sorry for my delay in answered you.
</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB"><br>
</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB">I use popbio with matrices, so I think it is not the best solution for IPMs (but remain A solution). However, I’m not actually able to create my own function, I will study more.</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB"><br>
</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB">About Dr. Kendall’s solution, I’m searching other methods to help me with my resampling of each year of data, because it is my first time with this type of analyses. Maybe the package “boot”
can do the bootstrapping? Is that package correct to follow also the first suggestion of Joachim Topper?</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB"> </span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB">Hence, at the moment, I’ll study the resampling/bootstrapping methods and I'll wait for other solutions, if exists, so I thanks to all for your suggestion, I agree with Dr. Jongejans and
I hope that others IPMusers share with me any possible idea for resolve my problem.</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB"><br>
</span></p>
<p class="MsoNormal"><span style="mso-ansi-language:EN-GB" lang="EN-GB">Elena</span><font style="font-size:12pt" face="Calibri,sans-serif" size="3" color="#000000"></font><br>
</p>
</div>
</font><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"> (Italia)<br>
</font><i>E-mail:</i> elenasulis@hotmail.it <br>
<i>Skype:</i> elena.sulis<br>
<br>
<br>
<div>
<hr id="stopSpelling">
Subject: Re: [Ipmpack-users] help with IPMpack<br>
To: elenasulis@hotmail.it; ipmpack-users@lists.r-forge.r-project.org<br>
From: e.jongejans@science.ru.nl<br>
Date: Wed, 27 Jul 2016 11:23:35 +0200<br>
<br>
Dear Elena,<br>
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>
<br>
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>
<br>
I hope this helps,<br>
Eelke<br>
<br>
<br>
<br>
<br>
<div class="ecxmoz-cite-prefix">On 20-Jul-16 1:48 PM, Elena S. wrote:<br>
</div>
<blockquote>
<div dir="ltr"><font style="font-size:12pt" face="Calibri,sans-serif" size="3" color="#000000"><br id="ecxFontBreak">
</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="ecxmoz-txt-link-abbreviated" href="mailto:elenasulis@hotmail.it">
elenasulis@hotmail.it</a> <br>
<i>Skype:</i> elena.sulis<br>
</div>
<br>
<fieldset class="ecxmimeAttachmentHeader"></fieldset> <br>
<pre>_______________________________________________
IPMpack-users mailing list
<a class="ecxmoz-txt-link-abbreviated" href="mailto:IPMpack-users@lists.r-forge.r-project.org">IPMpack-users@lists.r-forge.r-project.org</a>
<a class="ecxmoz-txt-link-freetext" href="http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/ipmpack-users" target="_blank">http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/ipmpack-users</a>
</pre>
</blockquote>
<br>
</div>
</div>
</div>
</body>
</html>