我正在使用R中的gRain包创建一个贝叶斯网络。当试图编译条件概率表时,我得到错误'表维度不匹配'。有问题的表采用A | B形式,其中A可以取三个可能的值,B可以取两个。我基于六种可能的组合在cptable定义中输入了总共12个值。任何帮助将非常感激,因为我检查和研究无济于事,无法看出出了什么问题。这是我的代码。
bucket_by_sequence_length
错误如下:
# define levels
lh <- c("low", "high")
lmh <- c("low", "medium", "high")
# specify the conditional probability tables
eh <- cptable(~eh, values=c(0.2, 0.8), levels=lh)
inf.oil.eh <- cptable(~inf | oil:eh, values = c(0.9, 0.1, 0.1, 0.9, 0.1, 0.9, 0.01, 0.99), levels=lh)
bp.oil <- cptable(~bp | oil, values=c(0.9, 0.1, 0.1, 0.9, 0, 1, 0.1, 0.9, 0.4, 0.6, 0.5, 0.5), levels=lmh)
oil.eh <- cptable(~oil | eh, values=c(0.9, 0.1, 0.05, 0.95), levels=lh)
rt.inf.eh <-cptable(~rt | inf:eh, values=c(0.9, 0.1, 0.1, 0.9, 0.1, 0.9, 0.01, 0.99), levels=lh)
# compile the tables
plist <- compileCPT(list(eh, oil.eh, inf.oil.eh, bp.oil, rt.inf.eh))
答案 0 :(得分:1)
您为bp|oil
的条件概率表指定了太多值。由于oil
有2个等级且bp
3,因此您需要6个条件概率,但行中有12个
bp.oil <- cptable(~bp | oil, values=c(0.9, 0.1, 0.1, 0.9, 0, 1, 0.1, 0.9, 0.4, 0.6, 0.5, 0.5), levels=lmh)