R中iGraph的random_walk函数未在吸收态终止的问题

嘿，这个问题我之前也碰到过！核心原因就是你的游走逻辑里没加吸收态的终止判断——不管是手动写循环还是用现成工具包，默认都不会自动在吸收态停下，得自己明确设置终止条件才行。下面分两种常见场景给你解决方案：

一、如果是手动实现的随机游走 #

假设你写了类似这样的基础循环（没处理吸收态）：

# 示例转移矩阵（状态1-5，2和5为吸收态）
trans_mat <- matrix(
  c(0.5, 0.5, 0, 0, 0,
    0, 1, 0, 0, 0,
    0.3, 0, 0.4, 0.3, 0,
    0, 0, 0.5, 0, 0.5,
    0, 0, 0, 0, 1),
  nrow = 5, byrow = TRUE,
  dimnames = list(c(1,2,3,4,5), c(1,2,3,4,5))
)

current_state <- 1
visited <- c(current_state)

# 错误写法：无终止条件，碰到吸收态会无限循环
while(TRUE) {
  current_state <- sample(colnames(trans_mat), 1, prob = trans_mat[current_state, ])
  visited <- c(visited, current_state)
}

这种代码会一直跑下去，因为吸收态（比如2）的转移概率永远指向自己，但循环没判断要不要停。

修正方法：在循环里加入吸收态的终止判断：

current_state <- 1
visited <- c(current_state)

# 手动指定吸收态，或者自动识别（见下方）
absorbing_states <- c("2", "5")

# 当当前状态不是吸收态时继续游走
while(!(current_state %in% absorbing_states)) {
  current_state <- sample(colnames(trans_mat), 1, prob = trans_mat[current_state, ])
  visited <- c(visited, current_state)
}

# 输出所有访问过的顶点
print(visited)

如果不想手动指定吸收态，还可以从转移矩阵自动识别（更灵活）：

# 自动判断：行中仅对角线元素为1的状态就是吸收态
absorbing_states <- rownames(trans_mat)[
  apply(trans_mat, 1, function(row) {
    state_idx <- which(rownames(trans_mat) == names(row))
    all(row == 0) && row[state_idx] == 1
  })
]

二、如果用的是markovchain包 #

这个包自带马尔可夫链工具，但默认生成序列是固定长度，不会自动在吸收态终止。比如错误用法：

library(markovchain)
mc <- new("markovchain", transitionMatrix = trans_mat)
# 生成固定10步的序列，碰到吸收态也会重复输出该状态
seq <- markovchainSequence(n = 10, markovchain = mc, t0 = 1)

正确做法：封装一个函数，生成序列直到碰到吸收态：

library(markovchain)
mc <- new("markovchain", transitionMatrix = trans_mat)

# 用包自带函数识别吸收态
absorbing_states <- absorbingStates(mc)

generate_absorbing_walk <- function(mc, start_state) {
  current <- start_state
  visited <- c(current)
  # 未到达吸收态就继续游走
  while(!(current %in% absorbing_states)) {
    current <- rmarkovchain(n = 1, object = mc, t0 = current)
    visited <- c(visited, current)
  }
  return(visited)
}

# 生成完整游走路径
walk_sequence <- generate_absorbing_walk(mc, start_state = "1")
print(walk_sequence)

额外提醒 #

一定要确认你的转移矩阵定义正确：吸收态的转移行必须满足“从该状态出发只能留在自己”（比如状态2的转移行是c(0,1,0,0,0)），否则就算加了终止判断，逻辑也会出错。

内容的提问来源于stack exchange，提问作者J. Doe