我正在尝试在网络上执行 Pearson's X2 测试的两个版本。
这是我正在运行的代码：
library(bnlearn)
library(Rgraphviz)
library(gRain)
library(graph)
library(grid)
library(snow)
dag <- empty.graph(names(alarm))
modelstring(dag) <- paste("[HIST|LVF][CVP|LVV][PCWP|LVV][HYP][LVV|HYP:LVF]",
"[LVF][STKV|HYP:LVF][ERLO][HRBP|ERLO:HR][HREK|ERCA:HR][ERCA][HRSA|ERCA:HR]",
"[ANES][APL][TPR|APL][ECO2|ACO2:VLNG][KINK][MINV|INT:VLNG][FIO2]",
"[PVS|FIO2:VALV][SAO2|PVS:SHNT][PAP|PMB][PMB][SHNT|INT:PMB][INT]",
"[PRSS|INT:KINK:VTUB][DISC][MVS][VMCH|MVS][VTUB|DISC:VMCH]",
"[VLNG|INT:KINK:VTUB][VALV|INT:VLNG][ACO2|VALV][CCHL|ACO2:ANES:SAO2:TPR]",
"[HR|CCHL][CO|HR:STKV][BP|CO:TPR]", sep = "")
alarm.gs <- gs(alarm, test = "x2")
alarm.mc <- gs(alarm, test = "mc-x2", B=10000)
x2 测试运行，但 mc-x2 测试停止。
我希望有人能帮助解决它。