Saltar al contenido principal
LibreTexts Español

10.4: Simulación de árboles de nacimiento y muerte

  • Page ID
    53891
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    Podemos utilizar las propiedades estadísticas de los modelos de nacimiento y muerte para simular árboles filogenéticos a través del tiempo. Podríamos comenzar con un solo linaje en el tiempo 0. Sin embargo, el árbol filogenético a menudo comienza con el primer evento de especiación en el clado, por lo que también se puede comenzar la simulación con dos linajes en el tiempo 0 (esta diferencia se relaciona con la distinción entre las edades de copa y tallo de los clados; véase también el Capítulo 11).

    Para simular nuestro árbol, necesitamos dibujar tiempos de espera entre eventos de especiación y extinción, conectar nuevos linajes con el árbol y podar los linajes cuando se extinguen. También necesitamos un criterio de parada, que puede tener que ver con un número particular de taxones o un intervalo de tiempo fijo. Consideraremos esto último, y dejaremos árboles en crecimiento a un número fijo de taxones como ejercicio para el lector. Nuestro algoritmo de simulación es el siguiente.

    algoritmo de simulación

    Supongamos que tenemos un cierto número de linajes “vivos” en nuestro árbol (1 o 2 inicialmente), un tiempo actual (t c = 0 inicialmente) y un tiempo de parada t s t o p.

    1. Dibuja un tiempo de espera t i al siguiente evento de especiación o extinción. Los tiempos de espera se extraen de una distribución exponencial con el parámetro de tasa N a l i v e * (λ + μ) donde N a l i v e es el número actual de linajes vivos en el árbol.
    2. Verifique si la simulación termina antes del siguiente evento. Es decir, si t c + t i > t s t o p, finalizar la simulación.
    3. Decidir si el siguiente evento es un evento de especiación [con probabilidad λ/(λ + μ)] o un evento de extinción [con probabilidad μ/(λ + μ)]. Esto se puede hacer dibujando un número aleatorio uniforme u i del intervalo [0, 1] y asignando especiación al evento si u i < λ/(λ + μ) y extinción de lo contrario.
    4. Si (3) es un evento de especiación, entonces elija un linaje vivo aleatorio en el árbol. Adjunte una nueva rama al árbol en este punto y agrega un nuevo linaje vivo a la simulación. Regresar al paso 1.
    5. Si (3) es un evento de extinción, elija un linaje vivo aleatorio en el árbol. Ese linaje ya está muerto. Mientras aún haya al menos un linaje vivo en el árbol, regrese a (1); de lo contrario, todo tu clado se ha extinguido, y puedes detener la simulación.

    Este procedimiento devuelve un árbol filogenético que incluye linajes vivos y muertos. Uno puede podar cualquier taxa extinto para devolver un árbol de nacimiento y muerte de sobrevivientes, lo que está más en línea con lo que normalmente estudiamos usando especies existentes. También vale la pena señalar que clados enteros pueden —y muchas veces lo hacen— extinguirse bajo este protocolo antes de que uno alcance el tiempo t s t o p. Tenga en cuenta también que existe una forma mucho más eficiente de simular árboles (Stadler 2011).

    Podemos pensar en las predicciones filogenéticas de los modelos de nacimiento y muerte de dos maneras: considerando la topología de los árboles y considerando la longitud de las ramas de los árboles. Consideraré cada uno de estos dos aspectos de los árboles a continuación.


    This page titled 10.4: Simulación de árboles de nacimiento y muerte is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Luke J. Harmon via source content that was edited to the style and standards of the LibreTexts platform.