r - 用于条件独立测试的 PC 算法中的“indepTest”错误

Question

我正在使用 PC 算法函数，其中条件独立是属性之一。面临以下代码中的错误。请注意，这里的“数据”是我一直在使用的数据，1,6,2ingaussCItest是数据的邻接矩阵 x 和 y 中的节点位置。

代码：

library(pcalg)   
suffstat <- list(C = cor(data), n = nrow(data))    
pc.data <- pc(suffstat,
              indepTest=gaussCItest(1,6,2,suffstat),  
              p=ncol(data),alpha=0.01)

错误：

indepTest(x, y, nbrs[S], suffStat) 中的错误：
找不到函数“indepTest”

Answer 1

下面是运行的代码。去掉了 gaussCItest 的参数作为它的一个函数，可以直接使用。

library(pcalg)   
suffstat <- list(C = cor(data), n = nrow(data))    
pc.data <- pc(suffstat,indepTest=gaussCItest, p=ncol(data),alpha=0.01)