Stochastic Process Simulation

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
0
down vote

favorite












I am trying to understand some basic stochastic process simulation for population growth. Let us say, a population grows exponentially at a constant rate $R$, given by the equation,
$$
P(t+1) = R P(t)
$$
In the book "A Biologist's Guide to Mathematical Modeling in Ecology and Evolution" the authors generate various realizations of this model by taking a value of the rate parameter (like $R = 1.2$) and sampling from a Poisson distribution with mean $RP(t)$ to get $P(t+1)$. In reality, in some cases, we may not know the rate parameter $R$. Is there a way to estimate $R$ if we know the mean value of $P(t)$ and $P(t+1)$ at some point $t$, by such simulations? Thanks.







share|cite|improve this question























    up vote
    0
    down vote

    favorite












    I am trying to understand some basic stochastic process simulation for population growth. Let us say, a population grows exponentially at a constant rate $R$, given by the equation,
    $$
    P(t+1) = R P(t)
    $$
    In the book "A Biologist's Guide to Mathematical Modeling in Ecology and Evolution" the authors generate various realizations of this model by taking a value of the rate parameter (like $R = 1.2$) and sampling from a Poisson distribution with mean $RP(t)$ to get $P(t+1)$. In reality, in some cases, we may not know the rate parameter $R$. Is there a way to estimate $R$ if we know the mean value of $P(t)$ and $P(t+1)$ at some point $t$, by such simulations? Thanks.







    share|cite|improve this question





















      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      I am trying to understand some basic stochastic process simulation for population growth. Let us say, a population grows exponentially at a constant rate $R$, given by the equation,
      $$
      P(t+1) = R P(t)
      $$
      In the book "A Biologist's Guide to Mathematical Modeling in Ecology and Evolution" the authors generate various realizations of this model by taking a value of the rate parameter (like $R = 1.2$) and sampling from a Poisson distribution with mean $RP(t)$ to get $P(t+1)$. In reality, in some cases, we may not know the rate parameter $R$. Is there a way to estimate $R$ if we know the mean value of $P(t)$ and $P(t+1)$ at some point $t$, by such simulations? Thanks.







      share|cite|improve this question











      I am trying to understand some basic stochastic process simulation for population growth. Let us say, a population grows exponentially at a constant rate $R$, given by the equation,
      $$
      P(t+1) = R P(t)
      $$
      In the book "A Biologist's Guide to Mathematical Modeling in Ecology and Evolution" the authors generate various realizations of this model by taking a value of the rate parameter (like $R = 1.2$) and sampling from a Poisson distribution with mean $RP(t)$ to get $P(t+1)$. In reality, in some cases, we may not know the rate parameter $R$. Is there a way to estimate $R$ if we know the mean value of $P(t)$ and $P(t+1)$ at some point $t$, by such simulations? Thanks.









      share|cite|improve this question










      share|cite|improve this question




      share|cite|improve this question









      asked Aug 2 at 3:48









      user2167741

      1019




      1019




















          1 Answer
          1






          active

          oldest

          votes

















          up vote
          1
          down vote



          accepted










          Let $p(t) = logP(t), r = log R$. Then the equation becomes $p(t+1)=r+p(t)$, and you can turn that into something like $p(t+t_0)-p(t_0)=rt$. You can run a standard linear regression on the latter.






          share|cite|improve this answer





















          • Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
            – user2167741
            Aug 2 at 11:35











          • Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
            – user2167741
            Aug 2 at 14:26










          Your Answer




          StackExchange.ifUsing("editor", function ()
          return StackExchange.using("mathjaxEditing", function ()
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          );
          );
          , "mathjax-editing");

          StackExchange.ready(function()
          var channelOptions =
          tags: "".split(" "),
          id: "69"
          ;
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function()
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled)
          StackExchange.using("snippets", function()
          createEditor();
          );

          else
          createEditor();

          );

          function createEditor()
          StackExchange.prepareEditor(
          heartbeatType: 'answer',
          convertImagesToLinks: true,
          noModals: false,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          );



          );








           

          draft saved


          draft discarded


















          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2869699%2fstochastic-process-simulation%23new-answer', 'question_page');

          );

          Post as a guest






























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          up vote
          1
          down vote



          accepted










          Let $p(t) = logP(t), r = log R$. Then the equation becomes $p(t+1)=r+p(t)$, and you can turn that into something like $p(t+t_0)-p(t_0)=rt$. You can run a standard linear regression on the latter.






          share|cite|improve this answer





















          • Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
            – user2167741
            Aug 2 at 11:35











          • Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
            – user2167741
            Aug 2 at 14:26














          up vote
          1
          down vote



          accepted










          Let $p(t) = logP(t), r = log R$. Then the equation becomes $p(t+1)=r+p(t)$, and you can turn that into something like $p(t+t_0)-p(t_0)=rt$. You can run a standard linear regression on the latter.






          share|cite|improve this answer





















          • Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
            – user2167741
            Aug 2 at 11:35











          • Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
            – user2167741
            Aug 2 at 14:26












          up vote
          1
          down vote



          accepted







          up vote
          1
          down vote



          accepted






          Let $p(t) = logP(t), r = log R$. Then the equation becomes $p(t+1)=r+p(t)$, and you can turn that into something like $p(t+t_0)-p(t_0)=rt$. You can run a standard linear regression on the latter.






          share|cite|improve this answer













          Let $p(t) = logP(t), r = log R$. Then the equation becomes $p(t+1)=r+p(t)$, and you can turn that into something like $p(t+t_0)-p(t_0)=rt$. You can run a standard linear regression on the latter.







          share|cite|improve this answer













          share|cite|improve this answer



          share|cite|improve this answer











          answered Aug 2 at 6:56









          ertl

          427110




          427110











          • Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
            – user2167741
            Aug 2 at 11:35











          • Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
            – user2167741
            Aug 2 at 14:26
















          • Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
            – user2167741
            Aug 2 at 11:35











          • Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
            – user2167741
            Aug 2 at 14:26















          Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
          – user2167741
          Aug 2 at 11:35





          Brilliant! So you mean to say we can simulate from a Poisson for $p(t+t_0)$ and $p(t_0)$ and do a regression. But the $log$ of a Poisson will not be a Poisson right? I suppose we can still sample from Poisson and take a log.
          – user2167741
          Aug 2 at 11:35













          Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
          – user2167741
          Aug 2 at 14:26




          Also, the above equation should be $p(t+t_0)-p(t) = rt_0$.
          – user2167741
          Aug 2 at 14:26












           

          draft saved


          draft discarded


























           


          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2869699%2fstochastic-process-simulation%23new-answer', 'question_page');

          );

          Post as a guest













































































          Comments

          Popular posts from this blog

          What is the equation of a 3D cone with generalised tilt?

          Relationship between determinant of matrix and determinant of adjoint?

          Color the edges and diagonals of a regular polygon