我想绘制一个转换矩阵,但每个状态都需要 2 列。
我的矩阵是:
> R
0 30 60 90 <NA>
0 0.75 0.37 0.17 0.07 0.97
30 0.15 0.40 0.32 0.02 0.02
60 0.00 0.20 0.19 0.05 0.01
90 0.00 0.00 0.03 0.52 0.00
NA 0.10 0.03 0.29 0.35 0.00
因此,从状态 0 开始,剩下 75%,15% 变为 30,依此类推。
问题是我不想要以下情节:
library(diagram)
plotmat(R)
相反,我想要每个状态 2 列......根据这个答案 https://stackoverflow.com/questions/16112162/graph-flow-chart-of-transition-from-states我需要创建一个 10x10 的表......有没有其他方法可以做到同样的事情而不必创建这样的表?
我的想法是无需修改原始表即可获得此图:
为此,我使用以下代码转换了原始矩阵:
L<-matrix(nrow = 10, ncol = 10, byrow = TRUE, data = 0)
for (i in 1:(nrow(R))){
for (j in 1: ncol(R))
{L[i*2,j*2-1]<-R[i,j]
}}
rownames(L)<-c('0','0', '30','30','60','60','90','90','NA','NA')
plotmat(L[1:6,1:6])
Thanks