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
$begingroup$
Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':
- 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)$.
- 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.
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
$endgroup$
add a comment |
$begingroup$
Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':
- 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)$.
- 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.
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
$endgroup$
add a comment |
$begingroup$
Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':
- 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)$.
- 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.
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
$endgroup$
Fix $alpha >0$, and for $h > 0$, consider the Markov kernel $K_h$ derived by composing the following two `moves':
- 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)$.
- 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.
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
probability-theory stochastic-processes markov-chains markov-process weak-convergence
edited Mar 11 at 20:28
πr8
asked Feb 14 at 21:53
πr8πr8
9,69331025
9,69331025
add a comment |
add a comment |
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
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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