我有一系列事件:
library(markovchain)
sequence<-c("LHR - BA","BOS - BA","BOS - ZE","IAD - ZE","BOS - BA","LHR - BA",
"LGW - BA","TPA - BA","TPA - BA","LGW - BA","LHR - BA","BOS - BA",
"BOS - ZE","BOS - ZE","BOS - BA","LHR - BA","LHR - BA","BOS - BA",
"BOS - ZE","BOS - ZE","BOS - DL","ATL - DL","LHR - BA","BRU - BA")
使用这个序列,我得到一个使用以下函数的马尔可夫链:
sequenceMatr <- createSequenceMatrix(sequence, sanitize=FALSE)
mcFit <- markovchainFit(data=sequence, method="mle")
考虑下一个状态是"LHR - BA"
; 那么我如何以以下格式识别跨州的概率分布:
"LHR - BA", "BOS - ZE", "IAD - ZE", "BOS - BA", "LHR - BA", "LGW - BA", "TPA - BA"
"LHR - BA" 0.1 , 0.2 , 0.2 , 0.3 , 0.1 , 0.1 , 0.1