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

          Solar Wings Breeze Design and development Specifications (Breeze) References Navigation menu1368-485X"Hang glider: Breeze (Solar Wings)"e

          Kathakali Contents Etymology and nomenclature History Repertoire Songs and musical instruments Traditional plays Styles: Sampradayam Training centers and awards Relationship to other dance forms See also Notes References External links Navigation menueThe Illustrated Encyclopedia of Hinduism: A-MSouth Asian Folklore: An EncyclopediaRoutledge International Encyclopedia of Women: Global Women's Issues and KnowledgeKathakali Dance-drama: Where Gods and Demons Come to PlayKathakali Dance-drama: Where Gods and Demons Come to PlayKathakali Dance-drama: Where Gods and Demons Come to Play10.1353/atj.2005.0004The Illustrated Encyclopedia of Hinduism: A-MEncyclopedia of HinduismKathakali Dance-drama: Where Gods and Demons Come to PlaySonic Liturgy: Ritual and Music in Hindu Tradition"The Mirror of Gesture"Kathakali Dance-drama: Where Gods and Demons Come to Play"Kathakali"Indian Theatre: Traditions of PerformanceIndian Theatre: Traditions of PerformanceIndian Theatre: Traditions of PerformanceIndian Theatre: Traditions of PerformanceMedieval Indian Literature: An AnthologyThe Oxford Companion to Indian TheatreSouth Asian Folklore: An Encyclopedia : Afghanistan, Bangladesh, India, Nepal, Pakistan, Sri LankaThe Rise of Performance Studies: Rethinking Richard Schechner's Broad SpectrumIndian Theatre: Traditions of PerformanceModern Asian Theatre and Performance 1900-2000Critical Theory and PerformanceBetween Theater and AnthropologyKathakali603847011Indian Theatre: Traditions of PerformanceIndian Theatre: Traditions of PerformanceIndian Theatre: Traditions of PerformanceBetween Theater and AnthropologyBetween Theater and AnthropologyNambeesan Smaraka AwardsArchivedThe Cambridge Guide to TheatreRoutledge International Encyclopedia of Women: Global Women's Issues and KnowledgeThe Garland Encyclopedia of World Music: South Asia : the Indian subcontinentThe Ethos of Noh: Actors and Their Art10.2307/1145740By Means of Performance: Intercultural Studies of Theatre and Ritual10.1017/s204912550000100xReconceiving the Renaissance: A Critical ReaderPerformance TheoryListening to Theatre: The Aural Dimension of Beijing Opera10.2307/1146013Kathakali: The Art of the Non-WorldlyOn KathakaliKathakali, the dance theatreThe Kathakali Complex: Performance & StructureKathakali Dance-Drama: Where Gods and Demons Come to Play10.1093/obo/9780195399318-0071Drama and Ritual of Early Hinduism"In the Shadow of Hollywood Orientalism: Authentic East Indian Dancing"10.1080/08949460490274013Sanskrit Play Production in Ancient IndiaIndian Music: History and StructureBharata, the Nāṭyaśāstra233639306Table of Contents2238067286469807Dance In Indian Painting10.2307/32047833204783Kathakali Dance-Theatre: A Visual Narrative of Sacred Indian MimeIndian Classical Dance: The Renaissance and BeyondKathakali: an indigenous art-form of Keralaeee

          Method to test if a number is a perfect power? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 00:00UTC (8:00pm US/Eastern)Detecting perfect squares faster than by extracting square rooteffective way to get the integer sequence A181392 from oeisA rarely mentioned fact about perfect powersHow many numbers such $n$ are there that $n<100,lfloorsqrtn rfloor mid n$Check perfect squareness by modulo division against multiple basesFor what pair of integers $(a,b)$ is $3^a + 7^b$ a perfect square.Do there exist any positive integers $n$ such that $lfloore^nrfloor$ is a perfect power? What is the probability that one exists?finding perfect power factors of an integerProve that the sequence contains a perfect square for any natural number $m $ in the domain of $f$ .Counting Perfect Powers