Computing the Scaling Limit of a Nonnegative Markov ChainMarkov Chain PerturbationUnderstanding a Markov ChainThe expected time until reaching a specified set in a Markov chainMarkov chain with dynamic higher ordersMarkov chain of transition probabilitiesRequirements for approximating general stochastic processes as the limit of a sequence of Markov Chains?Showing transience condition for this continuous markov chainMarkov Chain: Balance equations for a failing systemWhat's the transition semigroup of the Markov chain generated by the Metropolis-Hastings algorithm?Reversal of an Autoregressive Cauchy Markov Chain

A link redirect to http instead of https: how critical is it?

Is it normal that my co-workers at a fitness company criticize my food choices?

Why does Bach not break the rules here?

Why doesn't the EU now just force the UK to choose between referendum and no-deal?

In a future war, an old lady is trying to raise a boy but one of the weapons has made everyone deaf

How to explain that I do not want to visit a country due to personal safety concern?

How could a scammer know the apps on my phone / iTunes account?

How to terminate ping <dest> &

Do the common programs (for example: "ls", "cat") in Linux and BSD come from the same source code?

Do I need to be arrogant to get ahead?

A Cautionary Suggestion

PTIJ: Who should I vote for? (21st Knesset Edition)

How to write cleanly even if my character uses expletive language?

Define, (actually define) the "stability" and "energy" of a compound

Why is the President allowed to veto a cancellation of emergency powers?

Is a party consisting of only a bard, a cleric, and a warlock functional long-term?

Welcoming 2019 Pi day: How to draw the letter π?

Python if-else code style for reduced code for rounding floats

Awsome yet unlucky path traversal

Co-worker team leader wants to inject his friend's awful software into our development. What should I say to our common boss?

Dice rolling probability game

What do Xenomorphs eat in the Alien series?

Look at your watch and tell me what time is it. vs Look at your watch and tell me what time it is

What has been your most complicated TikZ drawing?



Computing the Scaling Limit of a Nonnegative Markov Chain


Markov Chain PerturbationUnderstanding a Markov ChainThe expected time until reaching a specified set in a Markov chainMarkov chain with dynamic higher ordersMarkov chain of transition probabilitiesRequirements for approximating general stochastic processes as the limit of a sequence of Markov Chains?Showing transience condition for this continuous markov chainMarkov Chain: Balance equations for a failing systemWhat's the transition semigroup of the Markov chain generated by the Metropolis-Hastings algorithm?Reversal of an Autoregressive Cauchy Markov Chain













1












$begingroup$


Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':



  1. From $x^t$, move to $x^t+1/2 = x^t + y^t$, where $y^t sim textGamma(alpha h, 1)$ .

    • To fix notation, by this I mean $p(y) propto y^alpha h-1 exp(-y)$.


  2. From $x^t+1/2$, move to $x^t+1 = x^t cdot z^t$, where $z^t sim textBeta(alpha, alpha h)$.

    • Again, for clarity, by this I mean $p(z) propto z^alpha-1 (1-z)^alpha h-1$.


I use $K_h$ to denote the composite Markov kernel which takes me from $x^t$ to $x^t+1$.



I know that this chain has $textGamma(alpha, 1)$ as an invariant measure, and moreover, that the chain is reversible with respect to this measure.



I want to compute the behaviour as $h to 0^+$ of this chain. Initially, I thought that the chain would admit a diffusion limit (something like the CIR process), but after carrying out some simulations, it appears more likely that it's something like a jump-diffusion, or even a pure jump process.



Coarse Resolution SimulationModerate Resolution SimulationFine Resolution Simulation



Anyhow, I'm not sure how I should go about i) identifying a limiting process, and ii) proving rigorously that it is the true limit of these Markov kernels. I can identify that if I define



$$B(x;h) triangleq mathbfE_K_h [x^t+1 - x^t | x^t = x]$$



$$V(x;h) triangleq textbfVar_K_h [x^t+1 - x^t | x^t = x], $$



then



beginalign
B(x;h) &= h (alpha - x ) + o(h) \
V(x;h) &= h left alpha^2 + fracalphaalpha + 1 x^2 right + o(h).
endalign



I initially though that this would mean the limiting process would be



$$dX_t = ( alpha - X_t) dt + sqrtalpha^2 + fracalphaalpha + 1 x^2 dW_t$$



but this both i) doesn't have $textGamma(alpha, 1)$ as a stationary measure, and ii) doesn't account for the jump behaviour I observe in the simulations.



Any advice would be well-received.










share|cite|improve this question











$endgroup$
















    1












    $begingroup$


    Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':



    1. From $x^t$, move to $x^t+1/2 = x^t + y^t$, where $y^t sim textGamma(alpha h, 1)$ .

      • To fix notation, by this I mean $p(y) propto y^alpha h-1 exp(-y)$.


    2. From $x^t+1/2$, move to $x^t+1 = x^t cdot z^t$, where $z^t sim textBeta(alpha, alpha h)$.

      • Again, for clarity, by this I mean $p(z) propto z^alpha-1 (1-z)^alpha h-1$.


    I use $K_h$ to denote the composite Markov kernel which takes me from $x^t$ to $x^t+1$.



    I know that this chain has $textGamma(alpha, 1)$ as an invariant measure, and moreover, that the chain is reversible with respect to this measure.



    I want to compute the behaviour as $h to 0^+$ of this chain. Initially, I thought that the chain would admit a diffusion limit (something like the CIR process), but after carrying out some simulations, it appears more likely that it's something like a jump-diffusion, or even a pure jump process.



    Coarse Resolution SimulationModerate Resolution SimulationFine Resolution Simulation



    Anyhow, I'm not sure how I should go about i) identifying a limiting process, and ii) proving rigorously that it is the true limit of these Markov kernels. I can identify that if I define



    $$B(x;h) triangleq mathbfE_K_h [x^t+1 - x^t | x^t = x]$$



    $$V(x;h) triangleq textbfVar_K_h [x^t+1 - x^t | x^t = x], $$



    then



    beginalign
    B(x;h) &= h (alpha - x ) + o(h) \
    V(x;h) &= h left alpha^2 + fracalphaalpha + 1 x^2 right + o(h).
    endalign



    I initially though that this would mean the limiting process would be



    $$dX_t = ( alpha - X_t) dt + sqrtalpha^2 + fracalphaalpha + 1 x^2 dW_t$$



    but this both i) doesn't have $textGamma(alpha, 1)$ as a stationary measure, and ii) doesn't account for the jump behaviour I observe in the simulations.



    Any advice would be well-received.










    share|cite|improve this question











    $endgroup$














      1












      1








      1





      $begingroup$


      Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':



      1. From $x^t$, move to $x^t+1/2 = x^t + y^t$, where $y^t sim textGamma(alpha h, 1)$ .

        • To fix notation, by this I mean $p(y) propto y^alpha h-1 exp(-y)$.


      2. From $x^t+1/2$, move to $x^t+1 = x^t cdot z^t$, where $z^t sim textBeta(alpha, alpha h)$.

        • Again, for clarity, by this I mean $p(z) propto z^alpha-1 (1-z)^alpha h-1$.


      I use $K_h$ to denote the composite Markov kernel which takes me from $x^t$ to $x^t+1$.



      I know that this chain has $textGamma(alpha, 1)$ as an invariant measure, and moreover, that the chain is reversible with respect to this measure.



      I want to compute the behaviour as $h to 0^+$ of this chain. Initially, I thought that the chain would admit a diffusion limit (something like the CIR process), but after carrying out some simulations, it appears more likely that it's something like a jump-diffusion, or even a pure jump process.



      Coarse Resolution SimulationModerate Resolution SimulationFine Resolution Simulation



      Anyhow, I'm not sure how I should go about i) identifying a limiting process, and ii) proving rigorously that it is the true limit of these Markov kernels. I can identify that if I define



      $$B(x;h) triangleq mathbfE_K_h [x^t+1 - x^t | x^t = x]$$



      $$V(x;h) triangleq textbfVar_K_h [x^t+1 - x^t | x^t = x], $$



      then



      beginalign
      B(x;h) &= h (alpha - x ) + o(h) \
      V(x;h) &= h left alpha^2 + fracalphaalpha + 1 x^2 right + o(h).
      endalign



      I initially though that this would mean the limiting process would be



      $$dX_t = ( alpha - X_t) dt + sqrtalpha^2 + fracalphaalpha + 1 x^2 dW_t$$



      but this both i) doesn't have $textGamma(alpha, 1)$ as a stationary measure, and ii) doesn't account for the jump behaviour I observe in the simulations.



      Any advice would be well-received.










      share|cite|improve this question











      $endgroup$




      Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':



      1. From $x^t$, move to $x^t+1/2 = x^t + y^t$, where $y^t sim textGamma(alpha h, 1)$ .

        • To fix notation, by this I mean $p(y) propto y^alpha h-1 exp(-y)$.


      2. From $x^t+1/2$, move to $x^t+1 = x^t cdot z^t$, where $z^t sim textBeta(alpha, alpha h)$.

        • Again, for clarity, by this I mean $p(z) propto z^alpha-1 (1-z)^alpha h-1$.


      I use $K_h$ to denote the composite Markov kernel which takes me from $x^t$ to $x^t+1$.



      I know that this chain has $textGamma(alpha, 1)$ as an invariant measure, and moreover, that the chain is reversible with respect to this measure.



      I want to compute the behaviour as $h to 0^+$ of this chain. Initially, I thought that the chain would admit a diffusion limit (something like the CIR process), but after carrying out some simulations, it appears more likely that it's something like a jump-diffusion, or even a pure jump process.



      Coarse Resolution SimulationModerate Resolution SimulationFine Resolution Simulation



      Anyhow, I'm not sure how I should go about i) identifying a limiting process, and ii) proving rigorously that it is the true limit of these Markov kernels. I can identify that if I define



      $$B(x;h) triangleq mathbfE_K_h [x^t+1 - x^t | x^t = x]$$



      $$V(x;h) triangleq textbfVar_K_h [x^t+1 - x^t | x^t = x], $$



      then



      beginalign
      B(x;h) &= h (alpha - x ) + o(h) \
      V(x;h) &= h left alpha^2 + fracalphaalpha + 1 x^2 right + o(h).
      endalign



      I initially though that this would mean the limiting process would be



      $$dX_t = ( alpha - X_t) dt + sqrtalpha^2 + fracalphaalpha + 1 x^2 dW_t$$



      but this both i) doesn't have $textGamma(alpha, 1)$ as a stationary measure, and ii) doesn't account for the jump behaviour I observe in the simulations.



      Any advice would be well-received.







      probability-theory stochastic-processes markov-chains markov-process weak-convergence






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Mar 11 at 20:28







      πr8

















      asked Feb 14 at 21:53









      πr8πr8

      9,69331025




      9,69331025




















          0






          active

          oldest

          votes











          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',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader:
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          ,
          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%2f3113240%2fcomputing-the-scaling-limit-of-a-nonnegative-markov-chain%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes















          draft saved

          draft discarded
















































          Thanks for contributing an answer to Mathematics Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid


          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.

          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3113240%2fcomputing-the-scaling-limit-of-a-nonnegative-markov-chain%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Moe incest case Sentencing See also References Navigation menu"'Australian Josef Fritzl' fathered four children by daughter""Small town recoils in horror at 'Australian Fritzl' incest case""Victorian rape allegations echo Fritzl case - Just In (Australian Broadcasting Corporation)""Incest father jailed for 22 years""'Australian Fritzl' sentenced to 22 years in prison for abusing daughter for three decades""RSJ v The Queen"

          John Burke, 9th Earl of Clanricarde References Navigation menuA General and heraldic dictionary of the peerage and baronetage of the British EmpireLeigh Rayment's Peerage Pages

          Football at the 1986 Brunei Merdeka Games Contents Teams Group stage Knockout stage References Navigation menu"Brunei Merdeka Games 1986".