利用bnlearn函数进行贝叶斯推断时利用cpquery封装函数出错

liufayue

2360

收藏 2014-09-14

悬赏 20 个论坛币未解决

在贝叶斯网络中，利用bnlearn进行推断时出错。
r代码（完整代码，在安装bnlearn包后可直接复制到r中运行）：####################################
library(bnlearn)
data(learning.test)
fitted = bn.fit(hc(learning.test), learning.test)

bn.function <- function(network, evidence_data) {
  a <- NULL
  b <- nrow(evidence_data)
  for (i in 1:b) {
evi <- paste("(", names(evidence_data), "=='",
            sapply(evidence_data[i,], as.character), "')",
            sep = "", collapse = " & ")
a <- cpquery(network, (C=='c'), eval(parse(text=evi))) # 应该是这条指令有误
  }
  return(a)
}

test <- bn.function(fitted, learning.test)
[size=12.8000001907349px]########################################

报错信息为
Error in parse(text = evi) : object 'evi' not found

但是如果把cpquery函数单独拿出来执行，则没有问题
[size=12.8000001907349px]这是网上的相关解释，但是还是不能解决问题
[size=12.8000001907349px] http://r.789695.n4.nabble.com/Us ... loop-td4677927.html