[R] Markov chain simulation
Kaptue Tchuente, Armel
armel.kaptue at sdstate.edu
Tue Jan 28 19:23:44 CET 2014
Thank you for your reply.
Indeed, there are many ways to simulate a markov chain sequence.
For instance, if I assume that n=20, rnd_occ1<-c(0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1) and rnd_occ2<-c(0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0) are two markov chain sequences having the transition matrix TransitionMatrix<- matrix(c(0.7, 0.3, 0.4, 0.6),byrow=TRUE, nrow=2)
For the first sequence the state "1" has 9 occurences while for the second sequence, the state "1" has 6 occurences.
What I'm looking for it is an algorithm or a library to simulate efficiently such a markov chain sequence with for instance 12 ooccurences of the state '1'
From: Kehl Dániel [kehld at ktk.pte.hu]
Sent: Tuesday, January 28, 2014 8:00 AM
To: Kaptue Tchuente, Armel; R-help at r-project.org
Subject: RE: [R] Markov chain simulation
Hi, see inline
Feladó: r-help-bounces at r-project.org [r-help-bounces at r-project.org] ; meghatalmazó: Armel KAPTUE [armel.kaptue at sdstate.edu]
Küldve: 2014. január 28. 8:16
To: R-help at r-project.org
Tárgy: [R] Markov chain simulation
I'm wonder if in R there is a way to simulate a discrete Markov chains with
a specific number of occurence of state knowing the transition matrixway.
Indeed there are many ways to do that in R!
For example, how to simualte a markov chain of length n with p occurences
(p<n) of the sate '0' for a transition matrix defined by:
I am not sure what you mean here?! If you give the length of the chain, depending on the transition matrix you are going to get some occurences of state '0', I do not think you can specify n, p and the transition matrix all at once!
TransitionMatrix<- matrix(c(0.7, 0.3, 0.4, 0.6),byrow=TRUE, nrow=2)
colnames(TransitionMatrix) <- c('0','1')
row.names(TransitionMatrix) <- c('0','1')
Please try googling first, there are a lot of good examples on your problem, among the first ones:
[[alternative HTML version deleted]]
R-help at r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
More information about the R-help