La PizarraMatemática.

Máster y doctorado

Análisis topológico de datos

Una nube de puntos no viene con agujeros dibujados. Debemos decidir qué puntos conectamos y qué caras rellenamos. Construiremos objetos pequeños, calcularemos sus componentes y ciclos, y seguiremos cómo cambian al aumentar una distancia umbral. La escala, las caras y el cuerpo de coeficientes quedarán explícitos.

1. Conectar puntos y rellenar caras son decisiones diferentes

Un grafo contiene vértices y aristas. Un complejo simplicial puede contener además triángulos rellenos, tetraedros y símplices superiores, siempre con todas sus caras. Un contorno cerrado puede ser ciclo del grafo sin representar un agujero del complejo.

Rellenar la cara cambia el agujero aunque conserve el contorno
Un contorno, dos complejosSolo aristasCon carab₁=1b₁=0Un ciclo puede ser una frontera.Rips añade la cara si admitelas tres parejas del triángulo.Las caras forman parte del objeto.
Dos complejos con tres vértices y tres aristas formando un triángulo. A la izquierda solo se incluyen aristas y el ciclo no es frontera, con Betti uno igual uno. A la derecha se incluye también la cara triangular, coloreada, y ese mismo ciclo es su borde, por lo que Betti uno es cero. Ambos objetos tienen una componente conexa. El grafo de la izquierda no representa el Rips completo cuando las tres parejas de esos vértices ya están admitidas: en ese caso su definición exige incluir también la cara. El relleno indica un símplex real del modelo, no una elección ornamental.

La métrica, el ruido y la escala elegida influyen en la construcción. La homología describe el objeto construido bajo esas decisiones, no una figura única certificada por los puntos sin más información.

El mismo contorno, dos objetos

Tres vértices A,B,C y las aristas AB,BC,AC.

  1. Solo aristas

    El contorno triangular no está rellenado.

  2. Añadir cara ABC

    Ahora el contorno es el borde de una superficie incluida.

  3. Resultado

    El primero tiene un agujero de dimensión uno y el segundo no.

Resultado

No se determina el agujero mirando únicamente los segmentos exteriores.

Lo comprobamos

Una cara es parte de los datos del complejo, no solo color decorativo.

Componente no significa agujero

Compara tres puntos aislados con una cadena A−B−C.

  1. Aislados

    Hay tres componentes y ningún ciclo de aristas.

  2. Cadena

    Hay una componente y tampoco ciclo.

  3. Cerrar

    Añadir AC crea un ciclo del grafo; añadir además ABC lo rellena.

Resultado

Conectividad y homología de dimensión uno responden a preguntas distintas.

Lo comprobamos

Los índices cero y uno de Betti no se intercambian.

2. Rips añade un símplex cuando todas sus parejas están cerca

Fijamos r≥0 como distancia máxima entre pares. Vietoris–Rips incluye un conjunto de vértices como símplex si todas sus distancias son ≤r. Dos dan arista; tres dan cara triangular si cumplen las tres parejas; cuatro dan un tetraedro abstracto si cumplen las seis.

Algunos textos usan 2r o radios de bolas: cambiar convención cambia los números de nacimiento. Aquí r es distancia máxima y las igualdades se incluyen. Un complejo abstracto puede tener dimensión mayor que el espacio donde se midieron los puntos.

Un triángulo rectángulo de datos

A=(0,0), B=(1,0), C=(1,1).

  1. Distancias

    AB=BC=1 y AC=√2.
  2. r=1

    Aparece la cadena de dos aristas, sin ciclo.

  3. r=√2

    Llega AC y también la cara ABC: no nace un agujero de duración positiva.

Resultado

El grafo completo y el complejo completo no son el mismo objeto.

Lo comprobamos

En Rips no se añade la última arista dejando fuera su cara triangular.

Cuatro vértices de un cuadrado unitario

Lados uno y diagonales √2.

  1. r<1

    Cuatro puntos aislados.

  2. 1≤r<√2

    Los cuatro lados forman un ciclo sin diagonales ni caras.

  3. r≥√2

    Las seis aristas, las cuatro caras triangulares y el tetraedro abstracto están presentes.

Resultado

Hay una clase de dimensión uno solo en el intervalo intermedio.

Lo comprobamos

El tetraedro abstracto no afirma que los cuatro datos planos sean vértices no coplanarios en R³.

3. Aumentar r genera inclusiones, aunque puedan morir clases

Si todas las distancias de un símplex son ≤r, también son ≤s para s≥r. Por eso \operatorname{Rips}(r)\subseteq\operatorname{Rips}(s). Se añaden símplices; no se borran. Una clase puede morir porque una cara nueva la convierte en frontera.

Una nube finita tiene finitos umbrales de cambio, pero puede producir muchísimos símplices. Truncar dimensión ahorra cálculo y también puede impedir detectar fronteras necesarias. Para calcular H_{k} hacen falta los símplices de dimensión k+1 que correspondan.

Morir sin borrar aristas

El cuadrado al pasar de r=1 a r=√2.

  1. Conservar

    Los cuatro lados siguen estando.

  2. Añadir

    Entran diagonales y caras.

  3. Rellenar

    El ciclo de lados pasa a ser una frontera; desaparece su clase en H₁.

Resultado

Muere una clase de homología, no el contorno como cadena.

Lo comprobamos

Las inclusiones de complejos se mantienen.

Un truncamiento puede inventar agujeros

Cuatro vértices con todas las distancias admitidas.

  1. Solo grafo

    Seis aristas y cuatro vértices conexos dan tres ciclos independientes.

  2. Añadir caras

    H₁ se anula, pero la superficie de las cuatro caras tiene H₂ de dimensión uno.

  3. Añadir tetraedro

    Esa superficie se vuelve frontera y H₂ también se anula.

Resultado

La dimensión máxima de cálculo afecta qué resultados son correctos.

Lo comprobamos

Para H₁ bastan las caras de dimensión dos; para H₂ aquí hace falta la de dimensión tres.

4. El operador borde convierte dibujos en álgebra

Trabajaremos sobre F₂: coeficientes cero o uno y 1+1=0.

Una cadena suma símplices. El borde de AB es A+B y el de ABC es AB+AC+BC. Al aplicar borde dos veces, cada contribución se repite y se cancela: ∂²=0.

Sobre otros cuerpos se introducen orientaciones y signos coherentes. La identidad sigue diciendo que toda frontera es un ciclo. Esa inclusión es la que permite formar después un espacio cociente.

Calcular el borde del borde

Cara ABC en F₂.

  1. Primer borde

    AB+AC+BC.

  2. Segundo

    (A+B)+(A+C)+(B+C).

  3. Cancelar

    Cada vértice aparece dos veces: cero.

Resultado

El borde de una cara es una cadena cerrada.

Lo comprobamos

Cero es la cadena nula en este cuerpo, no una nueva posición geométrica.

Cerrar una cadena

Empieza con AB+BC.

  1. Borde

    A+B+B+C=A+C: no es ciclo.
  2. Añadir AC

    El borde adicional es A+C.

  3. Total

    A+C+A+C=0: AB+BC+AC sí es ciclo.
Resultado

Ser ciclo se comprueba con ∂, sin decidir todavía si está rellenado.

Lo comprobamos

Si ABC está incluido, ese ciclo es además una frontera.

5. Homología distingue ciclos de fronteras

Definimos H_k=\ker\partial_k/\operatorname{im}\partial_{k+1}.

Los elementos del núcleo son ciclos; la imagen son fronteras de símplices y cadenas superiores. El cociente identifica ciclos que difieren por una frontera. El número de Betti b_{k} es su dimensión.

Con n_{k} símplices de dimensión k, b_k=n_k-\operatorname{rank}\partial_k-\operatorname{rank}\partial_{k+1}. Usamos homología no reducida: b₀ cuenta componentes. No basta contar todos los ciclos del grafo ignorando las caras.

Triángulo sin y con cara

Tres vértices y tres aristas.

  1. Borde de aristas

    Su matriz tiene rango dos; el espacio de ciclos tiene dimensión 3−2=1.

  2. Sin cara

    La imagen de ∂₂ es cero: b₁=1.
  3. Con cara

    El borde ABC genera una imagen de dimensión uno: b₁=1−1=0.

Resultado

La cuenta explica el cambio del dibujo inicial.

Lo comprobamos

En ambos casos b₀=1 porque el objeto está conectado.

Un cuadrado sin caras

Cuatro vértices y cuatro aristas en ciclo.

  1. Rango

    ∂₁ tiene rango tres.

  2. Componentes

    b₀=4−3=1.
  3. Ciclos

    b₁=4−3−0=1.
Resultado

El perímetro aporta una clase independiente.

Lo comprobamos

La afirmación exige que no haya caras incluidas que lo rellenen.

6. Las barras siguen clases mediante las inclusiones

Las inclusiones de complejos inducen aplicaciones entre sus homologías. Sobre un cuerpo y para esta filtración finita, la persistencia se describe mediante intervalos. Una barra [b,d) nace en b y deja de estar viva en d; su longitud es d−b.

Las barras del cuadrado distinguen componentes y ciclo
Cuadrado de lado 11 ≤ r < √2: un cicloDiagonales de longitud √2H₀012rH₁1√2rH₁: [1,√2); longitud √2−1
Cuatro vértices de un cuadrado de lado uno, con lados unidos en el intervalo de escalas uno incluido a raíz de dos excluido. Abajo se dibuja homología no reducida sobre un cuerpo para Rips completo. En H cero hay cuatro barras que nacen en cero; tres mueren en uno cuando se conectan los lados y una sigue indefinidamente, indicada con flecha. En H uno hay una única barra de uno a raíz de dos, cerrada al nacer y abierta al morir. En raíz de dos entran diagonales y todas las caras, y el ciclo se rellena. La convención de r es distancia máxima, no radio de bolas. Los ejes de barras usan la misma escala.

Una clase que sigue al final puede tener muerte infinita o quedar censurada como no observada. No se inventa su muerte al terminar un archivo. Empates pueden cambiar representantes elegidos, sin cambiar las multiplicidades de intervalos pertinentes.

La barra del cuadrado

Cuadrado unitario con Rips completo.

  1. Nacimiento

    El ciclo aparece en r=1.
  2. Muerte

    Se rellena en r=√2.
  3. Longitud

    La barra H₁ es [1,√2) y su persistencia √2−1.

Resultado

En r=√2 esa clase ya no está viva.

Lo comprobamos

No se ha usado radio de bolas ni distancia dividida por dos.

Leer una barra y un final abierto

Compara [1/2,2) con [0,∞).

  1. Primera

    Longitud 2−1/2=3/2.
  2. En su muerte

    La primera no sigue viva en dos.

  3. Segunda

    No tiene una longitud finita que calculemos restando el último r explorado.

Resultado

Una clase esencial y una observación interrumpida deben distinguirse según el resultado disponible.

Lo comprobamos

La escala conserva sus unidades.

7. Las componentes se fusionan y el grafo no muestra todas las muertes

En H₀ cada punto aislado inicia una componente. Cuando una arista une dos componentes distintas, queda una menos; una clase sobrevive a la fusión y otra muere. Cuando sus nacimientos empatan se elige un representante coherente.

En H₁ hay que revisar también las caras. En un Rips completo, tres aristas de un triángulo implican su cara en ese mismo umbral. El cálculo solo del grafo puede informar ciclos que no existen como clases del complejo completo.

Tres puntos en una recta

Posiciones cero, dos y cinco.

  1. r<2

    Hay tres componentes.

  2. r=2

    Se unen cero y dos; en r=3 se une cinco mediante la arista dos a cinco.

  3. Barras

    H₀ tiene [0,2), [0,3) y [0,∞).

Resultado

La arista de longitud cinco posterior no vuelve a reducir componentes.

Lo comprobamos

La identidad del representante superviviente es una elección, no un punto geométrico privilegiado.

Un triángulo equilátero

Tres distancias iguales a uno.

  1. Antes de uno

    Tres componentes y ningún ciclo.

  2. En uno

    Aparecen las tres aristas y la cara a la vez.

  3. Homología

    Queda una componente; no hay una barra H₁ de duración positiva.

Resultado

Un contorno visible no garantiza persistencia de un agujero.

Lo comprobamos

Dibujar solo líneas omitiría parte de la definición de Rips.

8. Estabilidad de la construcción y significado estadístico son preguntas distintas

Existen teoremas de estabilidad con métricas e hipótesis precisas. Por ejemplo, sobre los mismos vértices, si todas las distancias cambian como máximo η, el diámetro de cada símplex cambia como máximo η. Hay inclusiones de una filtración en la otra desplazando la escala por η.

Eso controla perturbaciones de la construcción; no certifica que cada barra larga sea una señal científica ni cada corta ruido. Se estudian muestreo, unidades, truncamiento, estabilidad pertinente y relación con la pregunta. Los números de Betti tampoco identifican una forma única.

Comprobar el desplazamiento de escala

Las métricas d y d′ satisfacen |d−d′|≤η en cada pareja.

  1. Un símplex

    Si todos sus pares tienen d≤r, entonces d′≤r+η.
  2. Incluir

    Rips_d(r) está dentro de Rips_d′(r+η).

  3. A la inversa

    El mismo argumento intercambiando d y d′ da la otra inclusión.

Resultado

Este control concreto exige comparar los mismos vértices y la cota de distancias indicada.

Lo comprobamos

No se extiende sin condiciones a cualquier cambio de muestreo.

Una barra larga depende también de unidades

Multiplica todas las distancias del cuadrado por cien.

  1. Nacimiento

    Pasa de uno a cien.

  2. Muerte

    De √2 a 100√2.

  3. Persistencia

    También se multiplica por cien, sin cambiar la forma relativa.

Resultado

Longitud absoluta mayor no significa por sí sola mayor significado estadístico.

Lo comprobamos

La comparación necesita una escala o normalización declarada.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Dos puntos

Dos puntos están a distancia 3. En Rips con umbral r, ¿cuándo se fusionan sus componentes?

  1. Datos

    Dos vértices distintos a distancia tres.

  2. Convención

    r es distancia máxima permitida para una arista.

  3. Antes

    Si r<3, no hay arista y quedan dos componentes.
  4. Umbral

    En r=3 se cumple distancia≤r.
  5. Fusión

    La arista une las componentes y queda una.

  6. Comprobar

    No usamos r=3/2 porque no es un radio de bolas en esta convención.

Resultado

Se fusionan en r=3.

Lo comprobamos

Se ha fijado umbral de distancia, no radio de bolas.

EJEMPLO RESUELTO 2

Leer una barra

Una clase nace en 2 y muere en 7. Calcula persistencia.

  1. Nacimiento

    La escala inicial de la barra es dos.

  2. Muerte

    Su extremo de muerte es siete.

  3. Definir

    Persistencia es muerte menos nacimiento.

  4. Restar

    7−2=5.
  5. Unidad

    Cinco unidades de la escala elegida.

  6. Comprobar

    En la convención [2,7), la clase ya no está viva en siete.

Resultado

La persistencia es 5.

Lo comprobamos

La longitud usa las unidades de la escala.

Los símbolos que hemos usado

Símplex
Vértice, arista, triángulo relleno o su análogo de dimensión mayor.
Complejo simplicial
Colección de símplices que contiene todas las caras de cada uno.
Filtración
Familia de complejos anidados al crecer un parámetro.
Ciclo / frontera
Cadena con borde cero / cadena que es borde de una dimensión superior.
Hₖ
Espacio de homología: ciclos de dimensión k identificados si difieren en una frontera.
bₖ
Número de Betti: dimensión de Hₖ sobre un cuerpo de coeficientes.
Nacimiento / muerte
Escalas de aparición y desaparición de una clase en una filtración.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar