我有以下状态的数据框：愤怒、平静和疲倦。下面的数据框提供了从一种状态转换到另一种状态的个别情况。
pre<-cbind(c(rep("tired",100),rep("angry",100),rep("calm",100)))
post<-cbind(c(rep("tired",50),rep("angry",70),rep("calm",100),rep("tired",80)))
df<-cbind(pre,post)
df<-as.data.frame(df)
colnames(df)<-c("pre","post")
我想要实现的是为三个状态构建马尔可夫链图，也称为“游乐场”，如下所示：
我将如何在 R 中做到这一点？
提前致谢！
【问题讨论】：
标签： r markov-chains
我们可以使用markovchain 包，使用diagram 中的plot 方法。
让我们首先从您的df 计算转移概率矩阵
states<-c("tired","angry","calm")
probsCase<-function(i,j){
sum(as.character(df$pre)==states[i] & as.character(df$post)==states[j])/sum(as.character(df$pre)==states[i])
}
transitionMatrix<-outer(1:3,1:3,Vectorize(probsCase))
colnames(transitionMatrix)<-states
rownames(transitionMatrix)<-states
现在使用markovchain 来初始化和绘制矩阵
library(markovchain)
markovChain <- new("markovchain", states=states, transitionMatrix=transitionMatrix)
plot(markovChain,package="diagram")
编辑：
如果你安装markovchain包有问题，我们其实可以不用，直接使用diagram包，只需要transitionMatrix
library(diagram)
plotmat(transitionMatrix,relsize=0.7)
然后您可以使用plotmat 中的选项根据您的喜好调整外观
【讨论】：
markovchain 时遇到问题，特别是它的依赖项RCurl。我有 I R 版本 3.2.3，在 Crunchbang Linux Waldorf 上运行。你是怎么安装的？你至少可以添加你的答案来发布情节的样子吗？再次感谢。
plotmat 中的relsize 参数。对于第二个，您需要重命名转换矩阵的行和列。编辑将发挥作用
box.col=c("red","green","blue") 会起作用。再次感谢。