R 创建三状态马尔可夫链图
我有以下三种状态的数据框:愤怒、平静和疲惫。下面的数据框提供了一种状态转换为另一种状态的个别情况R 创建三状态马尔可夫链图,r,markov-chains,R,Markov Chains,我有以下三种状态的数据框:愤怒、平静和疲惫。下面的数据框提供了一种状态转换为另一种状态的个别情况 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) co
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")
pre我们可以使用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
包时遇到问题,我们实际上不能使用它,而直接使用diagram
包,它只需要transitionMatrix
library(diagram)
plotmat(transitionMatrix,relsize=0.7)
然后,您可以使用plotmat
中的选项根据自己的喜好调整外观。感谢您的回复,我成功地访问了transitionMatrix,但在安装程序包markovchain
时遇到了问题,特别是它的依赖项RCurl
。我有IR3.2.3版,运行在Crunchbang Linux Waldorf上。你是怎么安装的?你能至少把你的答案加起来,贴出情节的样子吗?再次感谢,谢谢!它可以工作,但由于某些原因,绘图的右边缘和上边缘超出了R图形设备。有没有办法将数字1、2和3替换为绘图中的“疲惫”、“愤怒”、“平静”状态?对于第一种状态,请更改plotmat
中的relsize
参数。第二步,需要重命名转换矩阵的行和列。编辑将播放一些技巧,有没有一种方法可以为绘图中的树状态圆添加颜色?@我想出来了,对于其他人来说,添加box.col=c(“红色”、“绿色”、“蓝色”)
会起作用。再次感谢。
library(diagram)
plotmat(transitionMatrix,relsize=0.7)