马尔可夫链的R可视化|手工改变转移矩阵中的值

install.packages("igraph")
install.packages("ggplot2")
library(igraph)
library(ggplot2)

P <- matrix(c(0.7, 0.2, 0.1,
              0.3, 0.4, 0.3,
              0.2, 0.3, 0.5), nrow = 3, byrow = TRUE)

# 创建状态节点
states <- c("State 1", "State 2", "State 3")
g <- graph.empty(directed = TRUE)
g <- add.vertices(g, length(states), label = states)

# 添加转移边和权重
for (i in 1:nrow(P)) {
  for (j in 1:ncol(P)) {
    if (P[i, j] > 0) {
      g <- add.edges(g, c(i-1, j-1))
      E(g)$weight <- P[i, j]
    }
  }
}

# 绘制图形
plot(g, edge.label = E(g)$weight, vertex.label.cex = 1.5, vertex.size = 30)

# 修改转移矩阵
P[1, 2] <- 0.4

# 清除之前的图形
plot.new()

# 使用新的转移矩阵重新绘制图形
for (i in 1:nrow(P)) {
  for (j in 1:ncol(P)) {
    if (P[i, j] > 0) {
      g <- add.edges(g, c(i-1, j-1))
      E(g)$weight <- P[i, j]
    }
  }
}

plot(g, edge.label = E(g)$weight, vertex.label.cex = 1.5, vertex.size = 30)