Simulating a sample path using continuous time markov chain
I am trying to simulate a sample path using continuous time markov chain. This is what I've done:
set.seed(183427)
require(ECctmc)
# rates
r1 <- 1 # 1->2
r2 <- 0.75 # 2->3
r3 <- 0.5 # 3->1
r4 <- 0.5 # 3-> 2
Q <- matrix(c(-r1, r1, 0, 0, -r2, r2, r3, r4, -(r3+r4)),nrow = 3, byrow = TRUE)
# sample path
path <- sample_path(a=1, b = 3,t0=0, t1=5, Q=Q)
plot(stepfun(x=path[1:(nrow(path)-1),"time"], y = path[,"state"]), xlim = c(0,5), xlab = "Time", ylab = "State", main = "Sample path")
This works fine and I get a graph as an output where I have State on y axis and time on x-axis. However, in this case I have to specify my ending state b. I want to modify this code such that I don't have to specify the ending state and just give a rate matrix, initial distribution and the length of time interval T as an input. But I don't know how to do that. Help would be appreciated.
r markov-chains
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
add a comment |
I am trying to simulate a sample path using continuous time markov chain. This is what I've done:
set.seed(183427)
require(ECctmc)
# rates
r1 <- 1 # 1->2
r2 <- 0.75 # 2->3
r3 <- 0.5 # 3->1
r4 <- 0.5 # 3-> 2
Q <- matrix(c(-r1, r1, 0, 0, -r2, r2, r3, r4, -(r3+r4)),nrow = 3, byrow = TRUE)
# sample path
path <- sample_path(a=1, b = 3,t0=0, t1=5, Q=Q)
plot(stepfun(x=path[1:(nrow(path)-1),"time"], y = path[,"state"]), xlim = c(0,5), xlab = "Time", ylab = "State", main = "Sample path")
This works fine and I get a graph as an output where I have State on y axis and time on x-axis. However, in this case I have to specify my ending state b. I want to modify this code such that I don't have to specify the ending state and just give a rate matrix, initial distribution and the length of time interval T as an input. But I don't know how to do that. Help would be appreciated.
r markov-chains
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
add a comment |
I am trying to simulate a sample path using continuous time markov chain. This is what I've done:
set.seed(183427)
require(ECctmc)
# rates
r1 <- 1 # 1->2
r2 <- 0.75 # 2->3
r3 <- 0.5 # 3->1
r4 <- 0.5 # 3-> 2
Q <- matrix(c(-r1, r1, 0, 0, -r2, r2, r3, r4, -(r3+r4)),nrow = 3, byrow = TRUE)
# sample path
path <- sample_path(a=1, b = 3,t0=0, t1=5, Q=Q)
plot(stepfun(x=path[1:(nrow(path)-1),"time"], y = path[,"state"]), xlim = c(0,5), xlab = "Time", ylab = "State", main = "Sample path")
This works fine and I get a graph as an output where I have State on y axis and time on x-axis. However, in this case I have to specify my ending state b. I want to modify this code such that I don't have to specify the ending state and just give a rate matrix, initial distribution and the length of time interval T as an input. But I don't know how to do that. Help would be appreciated.
r markov-chains
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
I am trying to simulate a sample path using continuous time markov chain. This is what I've done:
set.seed(183427)
require(ECctmc)
# rates
r1 <- 1 # 1->2
r2 <- 0.75 # 2->3
r3 <- 0.5 # 3->1
r4 <- 0.5 # 3-> 2
Q <- matrix(c(-r1, r1, 0, 0, -r2, r2, r3, r4, -(r3+r4)),nrow = 3, byrow = TRUE)
# sample path
path <- sample_path(a=1, b = 3,t0=0, t1=5, Q=Q)
plot(stepfun(x=path[1:(nrow(path)-1),"time"], y = path[,"state"]), xlim = c(0,5), xlab = "Time", ylab = "State", main = "Sample path")
This works fine and I get a graph as an output where I have State on y axis and time on x-axis. However, in this case I have to specify my ending state b. I want to modify this code such that I don't have to specify the ending state and just give a rate matrix, initial distribution and the length of time interval T as an input. But I don't know how to do that. Help would be appreciated.
r markov-chains
r markov-chains
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
edited Dec 20 '18 at 8:11
Rohan Nadagouda
257212
257212
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
asked Nov 20 '18 at 19:25
Shafa HaiderShafa Haider
468
468
add a comment |
add a comment |
0
active
oldest
votes
Your Answer
StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
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
},
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%2fstackoverflow.com%2fquestions%2f53400157%2fsimulating-a-sample-path-using-continuous-time-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 Stack Overflow!
- 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.
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%2fstackoverflow.com%2fquestions%2f53400157%2fsimulating-a-sample-path-using-continuous-time-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
