Nothing
#function to extract matrices
extractMatrices <- function(mcObj) {
require(matlab)
mcObj <- canonicForm(object = mcObj)
#get the indices of transient and absorbing
transIdx <- which(states(mcObj) %in% transientStates(mcObj))
absIdx <- which(states(mcObj) %in% absorbingStates(mcObj))
#get Q, R and I
Q <- as.matrix(mcObj@transitionMatrix[transIdx,transIdx])
R <- as.matrix(mcObj@transitionMatrix[transIdx,absIdx])
I <- as.matrix(mcObj@transitionMatrix[absIdx, absIdx])
#get fundamental matrix
N <- solve(eye(size(Q)) - Q)
#final absorbion probabilities
NR <- N %*% R
#return
out <- list(
canonicalForm = mcObj,
Q = Q,
R = R,
I = I,
N=N,
NR=NR
)
return(out)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.