La idea más común de red aleatoria es que son grafos cuyas aristas están conectadas aleatoriamente con una determinada probabilidad.
En el terreno analítico son muy útiles para experimentar los escenarios de una red desconectada a una red totalmente conectada,
siendo importantes las fases intermedias que puede dar lugar a las famosas redes de mundo pequeño o redes complejas también llamadas
redes libres de escala (Watts y Strogatz, 1998). Esta sección tiene como objetivo describir las características centrales de las redes }
aleatorias y simular el modelo Erdos-Renyi en Netlogo.
El objetivo de esta sección es rescatar las características centrales de las redes aleatorias y simular la red de Erdos-Renyi en Netlogo.
Los primeros estudios importantes sobre la modelación de las redes aleatorias fueron iniciados por Rapoport (1949), una década más tarde,
Paul Erdos (1959) y Alfred Renyi (1960) estudiaron exhaustivamente estos grafos. Algunas características generales de la red Erdos-Renyi (ER)
es que un nuevo nodo se enlaza con igual probabilidad con el resto de la red, además, se dice que cuando los nodos y aristas adoptan un valor
finito corresponde a una distribución binomial, y cuando el número de nodos y aristas tienden a infinito, adoptan una distribución Poisson.
Keywords: distribución de grado, red aleatoria, red Erdos-Renyi, distribución Poisson, distribución de grado Poisson.
Es un grafo no dirigido cuyas aristas están conectadas aleatoriamente con una probabilidad \(p\), por tanto, siguen una distribución aleatoria.
Los grafos aleatorios son muy útiles para modelar sistemas muy grandes.
Ejemplos de redes aleatorias con N=100 y p = 0.03.
La red aleatoria Erdös-Rényi (1960)
La red aleatoria más representativa es la Erdós-Rényi. En esta red se tiene que un nuevo nodo que se enlaza con igual probabilidad con el
resto de la red, es decir, posee independencia estadística con el resto de nodos de la red (Newman y Watts Strogatz, 2002).
Si consideramos \(N\) nodos de una red sin conectar y distribuidos de forma aleatoria, podemos imaginar que en un instante inicial enlazamos
dos cualesquiera, de esta forma en pasos sucesivos vamos enlazando aleatoriamente de dos a dos nodos. Los nodos que se encuentren enlazados
se descartan. Si repetimos el proceso \(M\) veces eligiendo un par de nodos en cada turno al final habremos establecido como máximo de \(M\) enlaces
entre parejas de nodos. Si \(M\) es un valor pequeño con respecto al valor total de nodos muchos de los nodos estarán desconectados, mientras
que otros nodos estarán formando pequeñas islas.
Por el contrario, si \(M\) es grande en comparación con \(N\) el número total de nodos, es muy posible que casi todos los nodos estén enlazados entre sí.
Cuando se enlazan los nodos de esta forma aparecen propiedades de la distribución de grado p ya que posee propiedades de distribución Poisson.
Propiedades de la red ER
Para calcular la probabilidad distribución de grado \(p\) de que un nodo tenga \(k\) conexiones en la red aleatoria Erdös–Rényi, primero se intenta
calcular la probabilidad de que una pareja elegida al azar esté enlazada entre sí. Para ello se calcula el número total de posibles parejas en
una red de \(N\) nodos, a ese número total lo denominamos y su expresión es,
\[N_p=\displaystyle\binom{N}{2}=\frac{N(N-1)}{2}\]
como el número de parejas enlazadas por el modelo es M, se tiene por lo tanto la expresión analítica de la probabilidad
\[N_p=\frac{2M}{N(N-1)}\]
Si tomamos en la red generada un nodo particular al azar y le denominamos , el número de nodos enlazados a pares que contuvieran a sería \(N-1\),
ya que se puede enlazar con exactamente \(N-1\) nodos restantes de la red. Sin embargo, en los \(M\) enlaces generados, puede que no estuviera. Suponemos
entonces que estuviera en \(k\) de ellas. La probabilidad en este caso de que estuviera contenido en \(k\) parejas de las \(N-1\) posibles es,
\[N_p=\displaystyle\binom{N-1}{2}(p_c)^k=(1-p_c)^{N-1-k}\]
Esta fórmula corresponde a una distribución binomial para \(M\) y \(N\) de valor finito. Si se tiene en consideración ahora que la red empieza a crecer
hasta llegar a valores grandes del número de nodos (\(N\)) y de enlaces (\(M\)) hasta llegar al punto en que De esta forma se tiene que la cantidad,
\[ z=\frac{2M}{N}\]
permanece en valores completamente finitos y la distribución de grado \(P(k)\) se convierte en una distribución de Poisson de la forma,
\[P(k)= \frac{e^{-z}z^{k}}{k!}\]
que como se ha mencionado es una distribución de Poisson de media z. En los escritos posteriores del año 1960 Erdös y Rényi empezaron a estudiar
la dinámica de las redes en crecimiento profundizando en las transiciones de fase en las redes.
De hecho, siguiendo a Newman (2008), en la década de los años sesenta Erdos y Renyi demostraron que muchas propiedades del gráfico aleatorio
se pueden resolver exactamente tomando el límite de n grande. Demostraron que la red aleatoria posee una distribución de grado Poisson, en
donde la presencia o ausencia de aristas es independiente. De modo que la probabilidad de que un vértice tenga grado \(k\) es,
\[P(k)=\displaystyle\binom{n}{k}p^{k}(1-p)^{(n-k)}\cong\frac{z^{k}e^{-z}}{k!}\]
con la última igualdad aproximada volviéndose exacta en el límite de \(n\) grande y \(k\) fijo. Esta es la razón del nombre "Gráfico aleatorio de
Poisson".
Una vez instalada la versión más actual de NetLogo (por ejemplo la versión 6.0.4) en su computadora:
Abrir NetLogo
Ir a”Bibloteca de modelos”
Dar click en “IAMB”
Seleccionar “Chapter 5”
Seleccionar “Random Network”
Ahí encontrará la siguiente pantalla (Figura 1):
Figura 1: Variables de control del modelo ER
NetLogo contiene tres variables de control: el número de nodos, la cantidad de aristas y la probabilidad. Usted puede simular los escenarios
manipulando dichos parámetros.
El modelo ER presenta cuatro formas diferentes de crear redes aleatorias, dos de las cuales (las dos últimas) crean las variantes de las redes
clásicas de Erdős-Rényi. NetLogo contiene las cuatro formas de crear redes aleatorias, a saber:
WIRE1: Cada tortuga le pide a otra tortuga que se vincule a ella.
WIRE2: Elija una tortuga al azar y pídale que intente vincularse con otra tortuga al azar. Esto se hace NUM-NODES veces.
WIRE3: La clásica red aleatoria de Erdős – Rényi. Crea exactamente NUM-LINKS enlaces.
WIRE4: una variante de la clásica red aleatoria de Erdős-Rényi. Para cada par de tortugas, crea un vínculo entre ellas con probabilidad WIRING-PROB.
En la Figura 1 se pueden situar las cuatro maneras de generar redes aleatorias.
Para simular cualesquiera de los modelos mencionados: primero, tiene que diseñar el experimento, eligiendo el número de nodos, el número de aristas
y la probabilidad de enlace o conexión; segundo, debe cargar su experimento dando click en el botón “setup”; y tercero, tiene que elegir el método
de generación de la red, al dar click al método elegido aparecerá la red generada. Para efectos de este escrito se adopta el método Wire 3.
En la parte inferior izquierda del Interface también puede visualizar el grado máximo y el grado mínimo de la red recién creada. A partir de
estos elementos ya podemos realizar algunos experimentos de construcción de una red aleatoria.
Ejemplos de creación de una red en NetLogo
Utilizando el método Wire 3, con \(N=10\), \(k=6\) y \(p=0.02\) muestre que no es posible conformar una red, puede haber muchas instanciaciones,
diferentes en el número de enlaces y en los enlaces concretos. El grado máximo varía entre 2, 3 y 4; grado mínimo es 0.
Del Experimento 1, ¿Qué pasa si aumenta la probabilidad y la cantidad N de nodos? Verifique que tampoco se conforma una red.
Por ejemplo, \(N=77\) y \(p=0.02\), el grado máximo varía entre 3 y 8; y el grado mínimo es 0.
¿Qué pasa si aumenta el número de enlaces, por ejemplo, \(num-nodes=50\), \(num-links=100\) con el método wire3?, verifique que ahora sí
es posible conformar una red conexa. El grado máximo está entre 9 y 10, y el grado mínimo varía entre 1 y 0.
¿Qué pasa con la red si la cantidad de aristas es muy grande, por ejemplo, \( num-nodes=100\) y \(p=0.20\)? La red se vuelve muy densa,
el grado máximo está entre 34 y 27, y el grado mínimo varía entre 0 y 12.
¿Cómo generamos una red totalmente conectada? ¿Cuántas aristas tenemos para \(N\) nodos en general?
De los experimentos antes mencionados se percibe que conforme aumenta el tamaño de la red, cuando k es muy grande se tiene una red con
distribución aleatoria. De tal manera que podemos pasar de una red desconectada a una red aleatoria, esto es, la red pasa de una fase
desconectada a una fase aleatoria, aspecto que tiene que ver con el modelo.
Además, se dice que en una red aleatoria tienen valores de grado similares. Cuanto mayor es el tamaño de la red aleatoria, más parecidos
son los grados de los nodos. También se puede mostrar que la probabilidad de hallar un nodo con un grado muy alto o muy bajo es exponencialmente
pequeña. Lo que sí es determinante es que la estructura de la red aleatoria varía con el valor de \(p\).
En cuanto a su relación con las redes del mundo real, se afirma que sus propiedades no coinciden con él. Su construcción tiene un sentido más
analítico, tienen un coeficiente de agrupamiento bajo.
Erdös, P.; Rényi, A. (1959). "On Random Graphs. I.". Publicationes Mathematicae 6: 290–297
Erdös, P.; Rényi, A. (1960). "The Evolution of Random Graphs". Magyar Tud. Akad. Mat. Kutató Int. Közl. 5: 17–61.
Neuman, M. E. (2008). The Structure and Function of Complex Networks, SIAM REVIEW, Society for Industrial and Applied Mathematics, Vol. 45,No . 2, pp . 167–256.
Sayama, Hiroki (2015). Introduction to Modeling and Analysis of Complex Systems. Open SUNY textbooks, Milne Library, State University of New York at Geneseo, 485 pages, Print ISBN: 1942341083.
Solé, Ricard (2009). Redes Complejas, del genoma al internet. Grupo Editorial Planeta, 235 páginas, Barcelona, España.
Wilensky, U., Rand, W. (2008). NetLogo Random Network model. http://ccl.northwestern.edu/netlogo/models/RandomNetwork Center for Connected
Learning and Computer-Based Modeling, Northwestern Institute on Complex Systems, Northwestern University, Evanston, IL.Wilensky,
Wilensky, U. (1999). NetLogo. http://ccl.northwestern.edu/netlogo/ Center for Connected Learning and Computer-Based Modeling, Northwestern University, Evanston, IL..
Wilensky, U. & Rand, W. (2015). Introduction to Agent-Based Modeling: Modeling Natural, Social and Engineered Complex Systems with NetLogo. Cambridge, MA. MIT Press