Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- states <- c(1, 2)
- alpha <- c(1, 1)/2
- mat <- matrix(c(1/2, 1/2, 1/2, 1/2), nrow = 2, ncol = 2)
- nextX <- function(X, pMat)
- {
- #set.seed(1)
- probVec <- vector()
- if(X == states[1])
- {
- probVec <- pMat[1,]
- }
- if(X==states[2])
- {
- probVec <- pMat[2,]
- }
- return(sample(states, 1, replace=TRUE, prob=probVec))
- }
- steps <- function(alpha1, mat1, n1)
- {
- X0 <- sample(states, 1, replace=TRUE, prob=alpha1)
- if(n1 <=0)
- {
- return (X0)
- }
- else
- {
- vec <- vector(mode="numeric", length=n1)
- for (i in 1:n1)
- {
- X0 <- nextX(X0, mat1)
- vec[i] <- X0
- }
- return (vec)
- }
- }
- # steps(alpha1=alpha, mat1=mat, n1=5)
- simulate <- function(alpha1, mat1, n1)
- {
- for (i in 1:n1)
- {
- vec <- steps(alpha1, mat1, 5)
- print(vec)
- }
- }
- simulate(alpha, mat, 100)
- > simulate(alpha, mat, 100)
- [1] 1 2 2 2 2
- [1] 2 1 2 2 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 2 2 2 2
- [1] 2 1 1 2 2
- [1] 1 1 1 1 1
- [1] 2 2 1 2 2
- [1] 1 1 1 1 1
- [1] 2 2 1 1 2
- [1] 2 2 1 2 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 2 2 2
- [1] 1 1 1 1 1
- [1] 2 1 2 2 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 1 2 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 2 1 2 2
- [1] 2 2 2 1 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 1 2 2 1
- [1] 1 2 2 2 2
- [1] 1 1 2 2 2
- [1] 1 2 2 1 2
- [1] 1 1 1 1 1
- [1] 2 2 1 2 2
- [1] 2 2 2 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 2 2
- [1] 1 2 1 1 2
- [1] 2 2 1 1 1
- [1] 1 1 1 1 1
- [1] 2 2 2 2 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 2 2 2 2
- [1] 2 1 1 2 2
- [1] 1 1 1 1 1
- [1] 1 2 1 2 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 1 2 1 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 1 2 1 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 2 1 1 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 2 1 1 2
- [1] 1 1 1 1 1
- [1] 1 2 1 1 2
- [1] 1 1 1 1 1
- [1] 2 1 1 2 1
- [1] 1 1 1 1 1
- [1] 2 1 2 2 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 1 2 2 2 2
- [1] 1 1 1 1 1
- [1] 2 2 2 1 2
- [1] 2 2 2 1 1
- [1] 1 1 2 2 2
- [1] 1 1 1 1 1
- [1] 2 2 2 1 2
- [1] 1 1 1 1 1
- [1] 2 2 1 2 2
- [1] 2 2 1 2 1
- [1] 1 1 1 1 1
- [1] 2 2 2 2 2
- [1] 1 1 1 1 1
- [1] 1 2 1 2 2
- [1] 1 1 1 1 1
- [1] 2 1 1 2 1
- [1] 2 2 2 2 1
- [1] 2 2 2 2 2
- [1] 1 1 1 1 1
- [1] 2 2 2 1 1
- [1] 2 2 2 2 2
- [1] 1 1 1 1 1
- [1] 1 1 1 1 1
- [1] 2 1 2 2 1
- [1] 2 2 1 1 1
- [1] 1 1 1 1 1
- [1] 2 2 1 2 2
- [1] 2 1 2 2 2
- [1] 1 1 1 1 1
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement