sequence<-c("a", "b", "a", "a", "a", "a", "b", "a", "b", "a", "b", "a", "a", "b", "b", "b", "a")
mcFit<-markovchainFit(data=sequence,byrow=FALSE)
verifyMarkovProperty(sequence)
assessOrder(sequence)
assessStationarity(sequence,1)
divergenceTest(sequence, mcFit$estimate@transitionMatrix)Run the code above in your browser using DataLab