我在R中处理一个简单的问题。这是我的代码:
library(bnlearn)dag <- model2network("[Location][Quality][Cost|Location:Quality][NoPeople|Location:Cost]")plot(dag)quality.values <- factor(c("Good", "Normal", "Bad"))location.values <- factor(c("Good", "Bad"))cost.values <- factor(c("High", "Low"))nopeople.values <- factor(c("High", "Low"))quality.prob <- array(c(0.3, 0.5, 0.2), dim = 3, dimnames = list(quality = quality.values))location.prob <- array(c(0.6, 0.4), dim = 2, dimnames = list(location = location.values))cost.prob <- array(c(0.8, 0.2, 0.6, 0.4, 0.1, 0.9, 0.6, 0.4, 0.6, 0.4, 0.05, 0.95), dim = c(2, 3, 2), dimnames = list(cost = cost.values, quality = quality.values, location = location.values))nopeople.prob <- array(c(0.6, 0.4, 0.8, 0.2, 0.1, 0.9, 0.6, 0.4), dim = c(2, 2, 2), dimnames = list(NoPeople = nopeople.values, cost = cost.values, location = location.values))condProbTable <- list(Location = location.prob, Quality = quality.prob, Cost = cost.prob, NoPeople = nopeople.prob)bn <- custom.fit(dag, condProbTable)
但当我尝试运行这段代码时,得到的是一个错误Error in check.dnode.vs.spec(dist[[cpd]], old = fitted[[cpd]]$parents, : wrong dimensions for node Cost
。我不知道哪里出了问题。我刚开始学习R,所以任何帮助都会非常有用。
供参考,我试图构建的是这个:
谢谢!
回答:
原来问题出在大小写上。做如下修改
dag <- model2network("[location][quality][cost|location:quality][nopeople|location:cost]")
就可以解决问题了。