我有以下带有状态的数据框:生气、平静和疲倦。下面的数据框提供了从一种状态转换为另一种状态的个别情况。
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 中做到这一点?
提前致谢!
最佳答案
我们可以使用 markovchain
包,使用 plot
中的 diagram
方法。
让我们首先从您的 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
中的选项根据您的喜好调整外观关于r - 创建三态马尔可夫链图,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/36574814/