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

          Lowndes Grove History Architecture References Navigation menu32°48′6″N 79°57′58″W / 32.80167°N 79.96611°W / 32.80167; -79.9661132°48′6″N 79°57′58″W / 32.80167°N 79.96611°W / 32.80167; -79.9661178002500"National Register Information System"Historic houses of South Carolina"Lowndes Grove""+32° 48' 6.00", −79° 57' 58.00""Lowndes Grove, Charleston County (260 St. Margaret St., Charleston)""Lowndes Grove"The Charleston ExpositionIt Happened in South Carolina"Lowndes Grove (House), Saint Margaret Street & Sixth Avenue, Charleston, Charleston County, SC(Photographs)"Plantations of the Carolina Low Countrye

          random experiment with two different functions on unit interval Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 00:00UTC (8:00pm US/Eastern)Random variable and probability space notionsRandom Walk with EdgesFinding functions where the increase over a random interval is Poisson distributedNumber of days until dayCan an observed event in fact be of zero probability?Unit random processmodels of coins and uniform distributionHow to get the number of successes given $n$ trials , probability $P$ and a random variable $X$Absorbing Markov chain in a computer. Is “almost every” turned into always convergence in computer executions?Stopped random walk is not uniformly integrable

          How should I support this large drywall patch? Planned maintenance scheduled April 23, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?How do I cover large gaps in drywall?How do I keep drywall around a patch from crumbling?Can I glue a second layer of drywall?How to patch long strip on drywall?Large drywall patch: how to avoid bulging seams?Drywall Mesh Patch vs. Bulge? To remove or not to remove?How to fix this drywall job?Prep drywall before backsplashWhat's the best way to fix this horrible drywall patch job?Drywall patching using 3M Patch Plus Primer