在概率图模型中,cptable中表示条件概率表,使用之前需要调用一些包,代码如下:
library("gRbase")
library("gRain")
library("Rgraphviz")
machine_val <-c("working", "broken")
light_bulb_val <-c("good", "bad")
machine_prob <-c(99, 1)
light_bulb_prob <-c(99, 1, 60, 40)
M <-cptable(~machine, values = machine_prob, levels = machine_val)
L <-cptable(~light_bulb |machine, values = light_bulb_prob, levels = light_bulb_val)
plist <-compileCPT(list(M, L))
plist
plist$machine
plist$light_bulb输出结果:
> plist
P( machine )
P( light_bulb | machine )
> plist$machine
machine
working broken
0.99 0.01
> plist$light_bulb
machine
light_bulb working broken
good 0.99 0.6
bad 0.01 0.4
问题成功解决!
更多内容请移步主页查看!
点赞+关注,下次不迷路!