Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/69.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 如何将转换矩阵读入用于马尔可夫链建模的Cytoscape?_R_Format_Transition_Markov Chains_Cytoscape - Fatal编程技术网

R 如何将转换矩阵读入用于马尔可夫链建模的Cytoscape?

R 如何将转换矩阵读入用于马尔可夫链建模的Cytoscape?,r,format,transition,markov-chains,cytoscape,R,Format,Transition,Markov Chains,Cytoscape,我想读一个转换矩阵到Cytoscape中,用于马尔可夫链建模。举个例子,我从这篇文章中借用了以下转换矩阵(包括R代码): Cytoscape如何读取这样的过渡矩阵 多谢各位 转换矩阵的格式 这就是transitionMatrix的外观: tired angry calm tired 0.5 0.5 0.0 angry 0.0 0.2 0.8 calm 0.8 0.0 0.2 R代码 pre<-cbind(c(rep("tired",

我想读一个转换矩阵到Cytoscape中,用于马尔可夫链建模。举个例子,我从这篇文章中借用了以下转换矩阵(包括R代码):

Cytoscape如何读取这样的过渡矩阵

多谢各位

转换矩阵的格式

这就是transitionMatrix的外观:

       tired angry calm
 tired   0.5   0.5  0.0
 angry   0.0   0.2  0.8
 calm    0.8   0.0  0.2
R代码

 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")
 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

pre有一个名为aMatReader的应用程序,可以让你阅读邻接矩阵

--滑板车