11 votos

La probabilidad de epidemia a nivel mundial

Considere la posibilidad de $\mathbb{Z}^2$ como un grafo, donde cada nodo tiene cuatro vecinos. 4 señales emitidas son de $(0,0)$ en cada una de las cuatro direcciones (1 por sentido) . Un nodo que recibe una señal (o más) en un intervalo de tiempo se volverá a emitir a lo largo de los 4 bordes de sus cuatro vecinos en el siguiente paso de tiempo. Un nodo que no recibe una señal en el anterior paso de tiempo no emiten una señal independientemente de lo anterior recibido una señal. Hay un $50\%$ de probabilidad de que una señal se pierde cuando se viaja a lo largo de una arista entre dos nodos vecinos. Un nodo que recibe más de 1 señal actúa de la misma como si sólo recibió 1. La emisión de una señal en cada una de las 4 direcciones son eventos independientes.

¿Cuál es la probabilidad de que la señal en algún momento llegan a $(10^5,10^5)$? Investigación: las Simulaciones muestran: Sí. ~90%

¿Cuál es la probabilidad de que la señal en algún momento llegan a $(x,y)$ si una señal que viaja a lo largo de un borde muere con una probabilidad de $0<p<1$?

¿Cuál es el mínimo de p, para el cual la probabilidad de que N aleatoria inicial de las células vivas de morir se acerca a 0, como N se acerca a infinito? Experimento muestra p cerca de 0.2872.

En $\mathbb{Z}^1$, $p_{min}=0.6445...$, cómo calcular esto?

En $\mathbb{Z}^3$, $p_{min}=0.1775...$, cómo calcular esto?

15voto

Tim Puntos 499

Esto puede ser demasiado simplista, y sólo obtener un $71.0\%$ resultado. Supongo que todo el nodo $(10^5,10^5)$ todos los nodos vecinos tienen la misma probabilidad de ser alcanzado, se $P$ para distinguirla de la $p$ en la pregunta.

Otra suposición: las transmisiones de obras $1/2$ el tiempo.

        *

   *    *    *
       (a,a)

        *

Por lo $(a,a)$ podría recibir de cualquiera de los cuatro vecinos.

En cada uno de los nodos vecinos yo uso la probabilidad condicional de que el nodo que recibe la señal antes de que el nodo centro de $(a,a)$. Supongo que es igual a $.75P$ desde el vecino nodo no recibe la señal de $(a,a)$.

A continuación, tengo el siguiente ecuación usando el binomio de probabilidad fórmula

$$\binom{4}{1} (.75 P) (1-.75P)^3(1-0.5)+ \binom{4}{2} (.75P)^2 (1-.75P)^2 (1-.25)+ \binom{4}{3} (.75)^3 (1-.75P) (1-.125)+ \binom{4}{3} (.75)^4 (1-.0625) = P$$

Solving, $P \aprox .703$.

Este enfoque falla cuando cerca del origen, creo.

12voto

Did Puntos 1

Esto no calcular un valor exacto, pero algunos de los rigurosos límites superiores de la probabilidad de $r$ que cada sitio en $\mathbb Z^2$ es alcanzado, finalmente.

Considerar la Galton-Watson proceso de ramificación donde la progenie de cada individuo tiene la distribución del número de señales transmitidas en tiempo $n+1$ por cualquier sitio llegó en el momento $n$. Por lo tanto la descendencia de cualquier nodo en la Galton-Watson árbol es $0$, $1$, $2$, $3$ y $4$ con sus respectivas probabilidades $1$, $4$, $6$, $4$ y $1$ dividido por $16$.

A continuación, la generación de $n$ del proceso sobreestima el número de sitios alcanzados en el tiempo $n$ ya que las señales recibidas de los diferentes vecinos a estar separados. Por lo tanto la probabilidad de $q$ que el árbol se extingue subestima la probabilidad de que la nube de sitios que reciben una señal de que se vacíe. Por otro lado, si la nube sobrevive para siempre, el camino de cualquier eternamente señal transmitida es una simple caminata aleatoria en $\mathbb Z^2$, y estos son recurrentes por lo tanto cada sitio es alcanzado por el rastro de la nube en la red, finalmente.

La probabilidad de extinción $q$ es clásicamente el más pequeño de la raíz en $[0,1]$ de la de un solo paso ecuación $$ 16q=1+4t+6q^2+4t^3+q^4, $$ que es $q=0.087378$. Por lo tanto cada sitio en la red es alcanzado por el proceso de transmisión con una probabilidad de $r\lt91.2622\%$.

Ahora, la probabilidad de $r_T$ que $T=(10^5,10^5)$ corresponde al hecho de que la generación de $n=2\cdot10^5$ de el árbol no está vacío, cuya probabilidad $1-q_n$ está muy cerca de a $1-q$, por lo tanto $r_T\lt1-q_n\approx1-q$.

Nota: se obtiene un simple obligado darse cuenta de que, si el sitio inicial transmite nada, no hay otro sitio es alcanzado, y esto ocurre con una probabilidad de $1/16=6.25\%$, por lo tanto $r\lt93.75\%$.

3voto

JiminyCricket Puntos 143

En la misma vena como Didier de la respuesta de la prestación de los límites de la probabilidad de extinción $q$, también podemos obtener los límites de la probabilidad de transmisión $p$ requerido para la probabilidad de $N$ señales de morir para ir a $0$$N\to\infty$, que es la probabilidad requerida para la probabilidad de extinción $q$ no $1$.

Didier de la ecuación de $q$ puede ser reescrita como $q=(1-(1-q)/2)^4$, que dice que la señal va extinto si todas las cuatro señales de morir, ya sea porque no lo hacen en todo el borde ($1/2$) o porque se extinguieron después ($1-q$). Generalizando $d$ dimensiones y la probabilidad de transmisión $p$, esto es

$$q=(1-p(1-q))^{2d}\;.$$

Esta ecuación siempre tiene una raíz en $1$. Por lo suficientemente grande $q$, también tiene una segunda raíz en $[0,1]$, y la probabilidad de extinción es dada por la raíz. El caso crítico en el que la probabilidad de extinción se convierte en $1$ se produce cuando estos dos raíces coinciden. Diferenciando con respecto a $q$ y la sustitución de $q=1$ los rendimientos de la condición de $1$ a ser una doble raíz:

$$ \begin{eqnarray} 1&=&2dp(1-p(1-q))\;, \\ 1&=&2dp\;, \\ p&=&\frac1{2d}\;. \end{eqnarray} $$

Este es un límite inferior, ya que en el caso de los confluentes señales de que hay menos señales para mantener el fuego. Sus experimentos parecen indicar que el obligado se vuelve mejor con el aumento de la dimensión, lo cual tiene sentido ya que las señales se vuelven menos propensos a chocar.

1voto

Daniel Schierbeck Puntos 962

Esta no es una respuesta, pero era demasiado largo para un comentario. Más como algunas observaciones adicionales, las preguntas relacionadas con el sobre el tema de si la generación de funciones de ayuda.

Mi elección particular de funciones de generación de modelo aditivo señales y por lo tanto podría ser utilizado como límites superiores para no aditivas.

Mi enfoque de abajo es probablemente no original. Bueno combinatoria interpretaciones probablemente ya existen. Esta es una de Galton-Watson estocástico de ramificación proceso.

El caso unidimensional ha de generación de función en el poder formal de la serie de anillo con indeterminates $x,x^{-1}$ (o su cociente con el ideal de $\left< x^{-1}x-1 \right>$) dada por $$ \eqalign{ p_t^*(x) &= p^t \left(x+x^{-1}\right)^t = p^t \sum_{s=0}^{t} {t \elegir s} x^{2-t} = p^t \sum_{s=0}^{t} {t \elegir s} x^{t-2s} \cr & = p^t \left( \left\lfloor\tfrac{t}2\right\rfloor- \left\lfloor\tfrac{t-1}2\right\rfloor \right) {t \elegir t/2} + p^t \sum_{s=0}^{\left\lfloor\tfrac{t-1}2\right\rfloor} {t \elegir s} \left(x^{t-2}+x^{2-t}\right) } $$ para el tiempo de $t$, es decir, el coeficiente de $x^m$ es la probabilidad de que la señal llegue a coordinar $x=m$ tiempo $t$. El primer término en la segunda línea, el término constante, que representa la probabilidad de de una revisión del origen, es positivo con coeficiente máximo de $t$, incluso, pero el $0$ $t$ impar. Las diferencias entre paréntesis anterior, $$ e(t)= \left\lfloor\tfrac{t}2\right\rfloor- \left\lfloor\tfrac{t-1}2\right\rfloor = \tfrac{t+1}2-2\left\lfloor \tfrac{t+1}2\right\rfloor =\frac{1+(-1)^{t+1}}{2} =(t+1)~\text{mod}~2 $$ indicala uniformidad. También podríamos tomar la última expresión y tirar el (estrictamente) negativo powered términos de $x^{2s-t}$, redefiniendo el coeficiente de $x^m$ a la probabilidad de la señal de partida en el origen en el tiempo de $0$ llega a $x=\pm m$ (uno o el otro) después de tiempo $t$, y llaman a esto la reducción de la generación de la función $$ p_t(x) = p^t \left[ e(t) {t \elegir t/2} + \sum_{s=0}^{\left\lfloor\tfrac{t-1}2\right\rfloor} {t \elegir s} x^{t-2s} \right]. $$ Nota en particular de que $e(t){t\choose t/2}$ es por lo tanto el número de caminos de longitud $t$ que comienzan y terminan en el origen. Suma más de $t\ge0$, obtenemos una función de la generación de $$ p(x)=\sum_{t=0}^{\infty} p_t(x) $$ para una cantidad de tiempo gastado en cada coordenada, es decir, el coeficiente de $x^m$ es la suma de las probabilidades que la señal está en uno de $x=\pm m$ sobre todos los tiempos $t\ge 0$. De particular interés para nosotros es el total estimado de tiempo que estuvo en el origen, $$ p(0)=\sum_{t=0}^\infty{2t\c}p^{2}. $$ Combinatoria nota: no hay doble contabilidad (no hay necesidad de P. I. E.), ya que cada término se cuenta sólo el terminal punto sobre todos los caminos de longitud $t$. Esto actúa como una "fuerza motriz" de un oscilador armónico amortiguado por una probabilidad de morir. Por Stirling aproximación y la raíz de la prueba, esto converge, para $p<\frac12$.

Reclamación$\left(\mathbb{Z}^1\right)$: Para $r=2p<1$, se espera que el número total de visitas al origen converge a $$p(0)=\frac{1}{\sqrt{1-r^2}}.$$

Prueba: esto se desprende de la general, teorema del binomio, como, por ejemplo, de Yuanes el año 2001 la proposición 4.3: $$ \left(-4\right)^n {-\tfrac12 \elegir n} = \left(-4\right)^n \frac{(-1)(-3)\cdots(1-2n)}{2^n\n!} = \frac{(2n)!}{(n!)^2}. $$

Ahora, para el caso bidimensional, se puede evitar el uso de la multinomial teorema de expansión, si nos gusta, con $$ \eqalign{ p_t^*(x,y) &= p^t \left(x+x^{-1}+y+y^{-1}\right)^t = \sum_{s=0}^{t} {t \elegir s} p_s^*(x) \, p_{t-s}^*(y) \cr &= p^t \sum_{s=0}^{ t } \sum_{u=0}^{ s } \sum_{v=0}^{t-s} { t \elegir s} { s \elegir u} {t-s\elegir v} x^{2u-s} y^{2v+s-t} } $$ o, tirar la negativa powered términos, $$ \eqalign{ p_t(x,y) &= \sum_{s=0}^{t} {t \elegir s} p_s(x) \, p_{t-s}(y) %% %% Esta fórmula podría simplificarse de forma análoga a p_t^*(x) anterior: %% %\cr %&= p^t % \sum_{s=0}^{ t } % \sum_{u=0}^{ s } % \sum_{v=0}^{t-s} % { t \elegir s} % { s \elegir u} % {t-s\elegir v} % x^{2u-s} y^{2v+s-t} } $$ y tomar el análogo de la suma de la $p(x,y)$. Tenga en cuenta que $p_t(0,0)=0$ por extraño $t$, y que (por $2t$ incluso) $$ p_{2}(0,0) =p^{2}\sum_{s=0}^t{2t\elegir 2s}{2\elegir s}{2t-2s\elegir t-s} =p^{2}\sum_{s=0}^t{t\elegir s}^2 =p^{2}{2t\elegir 2}^2 $$ que es $p^t$ veces el cuadrado de un central de coeficiente binomial; cf. secuencia A000984. (Las dos últimas fórmulas, con el cuadrado de la binimial los coeficientes de aquí se ofrecen sin la prueba.) Por lo tanto $$ p(0,0)=\sum_{t=0}^\infty{2t\elegir t}^2p^{2t}, $$ que convergen $p<\frac14$.

Reclamación$\left(\mathbb{Z}^2\right)$: Para $r=4p<1$, se espera que el número total de visitas al origen converge a $$p(0,0)={}_2F_1\left(\tfrac12;\tfrac12;1;r^2\right)=\frac1{M(1,\sqrt{1-r^2})}$$ converge a un caso especial de la hipergeométrica de Gauss serie ${}_2F_1$ que es una integral elíptica completa de primera especiey puede ser calculado con la extremadamente rápida convergencia (Gauss de 1799) por la aritmética-media geométrica $$M(a,b)=\int_0^{\frac\pi2}\frac{d\phi}{\sqrt{(a\cos\phi)^2+(b\sin\phi)^2}}.$$

El primer OP pregunta es tal vez relacionado con el tiempo de espera $$ c_{mn} = \mathbb{E} \left[ T_{mn} \right] = \frac1{m!\,n!} \a la izquierda. \left( \partial_x^{\,m} \partial_y^{\n} p(x,y) \right) \right|_{(0,0)} \qquad(m,n\ge0) $$ pasó en uno de los cuatro puntos de $(\pm m,\pm n)$ que es la suma de los $x^my^n$ sobre todos los términos $p_t(x,y)$ y puede ser calculada de manera determinista sin mayor simplificación en $O(m^3n^3)$ tiempo mediante la adición de los productos de los coeficientes binomiales. Tal vez hay maneras de simplificar esto la utilización de una combinación de binomio de identidades, la paridad y simetría, o con asymptotics, para obtener una aproximación razonable en un posible cantidad de tiempo.

La conocida fórmula con el mayor orden derivado de la utilizando la serie de Taylor truco para el coeficiente deseado es válido en la reducción de la generación de la función de $m,n\ge 0$.

Como una especulación salvaje, es la respuesta para el (primer y segundo OP pregunta entonces $$ \frac{c_{mn}}{1-p(0,0)} $$ o algo en una vena similar, de forma análoga a la el jugador de la ruina problema o el razonamiento de @joriki o @Didier?

Conjeturas para dimensiones superiores:

Conjetura: El número esperado de visitas al origen en $\mathbb{Z}^d$ es $$ \sum_{n=0}^\infty {2n \choose n}^d p^{2n}. $$

Conjetura: El número de rutas en $\mathbb{Z}^d$ de la longitud de la $2n$ que comienzan y terminan en el origen es $$ {2n \choose n}^d. $$ Nota: $d=2$ es Richard Stanley bijective prueba de problema #230.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X