Re: Interpreting the eigen value of a population matrix (2nd try)
[I've cc'ed this to the list because I think that it may be of value to
other useRs, at least for archiving purposes.]
Sorry, I did that work a long time ago in XLispStat (before I switched
to R). The delta method and the bootstrap method for standard errors for
eigenvalues are described in Caswell (2001) Matrix Population Models,
2nd Ed., Chapter 12.
On Wed, 2007-08-29 at 13:31 -0400, Anouk Simard wrote:
> Thanks a lot for your help!
> Following you advice I have a couple more questions. I you could help me
> with it I would really appreciate although I would understand if you don't
> have time for it.
> I do have SE estimate on the survival and reproductive rate at each age
> class but I don't know how to incorporate this in a matrix model in order
> to do the bootstrap. So I was just wondering if you have any example of
> codes for it. I am not very familiar with R coding for this and I would
> very appeciate any help!
> Thank for your time again,
> rom: Simon Blomberg <s.blomberg1 <at> uq.edu.au>
> Subject: Re: Interpreting the eigen value of a population matrix (2nd try)
> Newsgroups: gmane.comp.lang.r.general
> Date: 2007-08-29 13:10:59 GMT (3 hours and 31 minutes ago)
> To get a confidence interval on lambda, you need to have measures of
> variability in the elements of the
> transition matrix. If you have that, you can use a parametric bootstrap to
> get approximate confidence
> intervals. I have done this, and it seems to work. Alternatively, you
> could calculate a Bayesian
> posterior density for lambda using the Bayesian melding methods developed
> by Adrian Raftery et al., and
> calculate an HPD interval from that. I've done that too. It's slightly
> more difficult, however.
> Simon Blomberg, BSc (Hons), PhD, MAppStat.
> Lecturer and Consultant Statistician
> Faculty of Biological and Chemical Sciences
> The University of Queensland
> St. Lucia Queensland 4072
> T: +61 7 3365 2506
> email: S.Blomberg1_at_uq.edu.au
Simon Blomberg, BSc (Hons), PhD, MAppStat.
Lecturer and Consultant Statistician
Faculty of Biological and Chemical Sciences
The University of Queensland
St. Lucia Queensland 4072
Room 320 Goddard Building (8)
T: +61 7 3365 2506
1. I will NOT analyse your data for you.
2. Your deadline is your problem.
The combination of some data and an aching desire for
an answer does not ensure that a reasonable answer can
be extracted from a given body of data. - John Tukey.