Ошибка в «indepTest» в алгоритме P C для проверки условной независимости

Question

Я использую функцию алгоритма P C, в которой условная независимость является одним из атрибутов. Обнаружена ошибка в следующем коде. Обратите внимание, что здесь «данные» - это данные, которые я использовал, а 1,6,2 в gaussCItest - это позиции узлов в моей матрице смежности 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)