Glosario de métodos

Academia: los métodos de descubrimiento causal detrás de nuestro trabajo.

Una referencia técnica breve para responsables de bioestadística y equipos clínicos: qué hace cada algoritmo, qué supuestos exige y cuándo tiene sentido usarlo.

PC (Peter-Clark)

Algoritmo de descubrimiento causal basado en tests de independencia condicional que reconstruye la estructura de un grafo a partir de datos observacionales.

El algoritmo PC parte de un grafo completo no dirigido y elimina progresivamente las aristas entre pares de variables que resultan condicionalmente independientes dado algún subconjunto de las demás variables. Cada arista eliminada representa la ausencia de una relación causal directa detectable con los datos disponibles. El resultado es un grafo esqueleto sobre el que se orientan después las aristas restantes aplicando reglas locales (por ejemplo, la identificación de estructuras en V o colisionadores) hasta donde la evidencia lo permite.

El resultado no es un DAG único sino un CPDAG (grafo parcialmente dirigido de una clase de equivalencia de Markov): un conjunto de aristas dirigidas donde la dirección causal es identificable sin ambigüedad, y aristas no dirigidas donde varias orientaciones son igualmente compatibles con los datos observados. PC no fuerza una dirección donde los datos no la determinan.

X ⊥⊥ Y | Z → eliminar arista X—Y del esqueleto

Los supuestos clave son fidelidad (faithfulness) —las independencias observadas en los datos reflejan la estructura del grafo y no cancelaciones casuales de parámetros—, suficiencia causal —no hay confusores no medidos entre las variables incluidas— y tests de independencia condicional fiables, que en la práctica requieren tamaños muestrales razonables cuando el número de variables crece.

Lo usamos como primer paso exploratorio cuando el equipo tiene múltiples variables candidatas y quiere una hipótesis estructural de partida, antes de refinarla con conocimiento de dominio.

GES (Greedy Equivalence Search)

Búsqueda voraz sobre el espacio de clases de equivalencia de Markov que optimiza directamente un criterio de bondad de ajuste penalizado.

GES explora el espacio de CPDAGs en dos fases. En la fase de avance añade aristas una a una, aceptando en cada paso el cambio que más mejora una puntuación penalizada (habitualmente BIC), hasta que ninguna adición mejora la puntuación. En la fase de retroceso elimina aristas siguiendo la misma lógica voraz, refinando el resultado de la primera fase.

A diferencia de PC, que decide arista a arista mediante tests de independencia, GES evalúa la calidad global de cada estructura candidata con una puntuación estadística. Esto tiende a ser más robusto frente a errores acumulados de tests individuales cuando el modelo subyacente se ajusta razonablemente bien al criterio de puntuación elegido.

Bajo los supuestos habituales (fidelidad, suficiencia causal y una familia de modelos donde la puntuación es consistente), GES converge asintóticamente a la clase de equivalencia de Markov correcta. El resultado, igual que en PC, es un CPDAG: la dirección de algunas aristas puede quedar sin determinar.

Lo usamos cuando interesa una puntuación global interpretable (BIC) y el conjunto de variables es de tamaño moderado a grande, donde encadenar muchos tests de independencia condicional —como hace PC— sería más frágil.

LiNGAM (Linear Non-Gaussian Acyclic Model)

Modelo lineal acíclico no gaussiano que identifica la dirección causal explotando la asimetría estadística de los residuos, no solo la estructura de independencias.

LiNGAM asume que cada variable es una combinación lineal de sus causas directas más un término de ruido, y que ese ruido no es gaussiano. Esa no-gaussianidad es la pieza clave: en un modelo lineal gaussiano, causa y efecto son estadísticamente indistinguibles por simetría, pero con residuos no gaussianos la dirección causal correcta deja una firma estadística asimétrica detectable, típicamente mediante análisis de componentes independientes (ICA).

A diferencia de PC y GES, que se detienen en una clase de equivalencia, LiNGAM puede identificar un DAG completamente dirigido a partir de datos puramente observacionales, sin necesitar restricciones adicionales del investigador para resolver la orientación de las aristas.

El precio de esa identificabilidad más fuerte es un supuesto más restrictivo: linealidad de las relaciones y no gaussianidad de todos los términos de ruido salvo, como mucho, uno. Cuando estos supuestos no se cumplen razonablemente, la dirección estimada puede ser incorrecta con confianza aparente.

Es útil cuando hay evidencia de que las relaciones son aproximadamente lineales y las distribuciones de las variables (o de los residuos tras un ajuste preliminar) se alejan claramente de la normalidad, algo frecuente en datos biológicos con colas pesadas o asimetría.

FCI (Fast Causal Inference)

En preparación

Extensión de PC que admite confusores latentes no medidos, sin exigir el supuesto de suficiencia causal.

Intervención y do-calculus

En preparación

Cálculo formal para estimar el efecto de intervenciones hipotéticas a partir de un grafo causal y datos puramente observacionales.

Mediación causal

En preparación

Descomposición de un efecto total en vías directas e indirectas a través de variables mediadoras.

Comparativa de métodos disponibles

Supuestos, complejidad y robustez de los métodos de descubrimiento causal disponibles.
MétodoSupuestosComplejidadRobustez
PC (Peter-Clark)
  • Suficiencia causal (sin confusores no medidos)
  • Fidelidad (faithfulness)
  • Aciclicidad
MediaSensible a errores de test en grafos densos o alta dimensión
GES (Greedy Equivalence Search)
  • Suficiencia causal
  • Fidelidad
  • Score descomponible (p.ej. BIC)
AltaMás estable que PC en muestras finitas, pero computacionalmente costoso
LiNGAM (Linear Non-Gaussian Acyclic Model)
  • Linealidad
  • Ruido no gaussiano
  • Suficiencia causal
  • Aciclicidad
BajaIdentificación fuerte si se cumple no-gaussianidad; falla si el ruido es gaussiano