13.2: Un Modelo de Diversificación Dependiente del Estado
- Page ID
- 54314
\( \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}\)Los modelos que consideraremos en este capítulo incluyen la evolución de rasgos y la diversificación de linajes asociados. En el caso más simple, podemos considerar un modelo donde el personaje tiene dos estados, 0 y 1, y las tasas de diversificación dependen de esos estados. Necesitamos modelar las transiciones entre estos estados, lo que podemos hacer de manera idéntica a lo que hicimos en el Capítulo 7 usando un modelo de Markov en tiempo continuo. Expresamos este modelo usando dos parámetros de velocidad, una tasa de avance q 01 y una tasa hacia atrás q 10.
Consideramos ahora la idea de que las tasas de diversificación pueden depender del estado del carácter. Suponemos que las especies con estado de carácter 0 tienen una cierta tasa de especiación (λ 0) y tasa de extinción (μ 0), y que las especies en 1 tienen tasas potencialmente diferentes tanto de especiación (λ 1) como de extinción (μ 1 ). Es decir, cuando el personaje evoluciona, afecta la tasa de especiación y/o extinción de los linajes. Así, tenemos un modelo de seis parámetros (Maddison et al. 2007). Suponemos que los linajes padres dan a luz a hijas con el mismo estado de carácter, es decir, que los estados de carácter no cambian en la especiación.
Es sencillo simular la evolución bajo nuestro modelo de diversificación dependiente del estado. Procedemos de la misma manera que lo hicimos para los modelos de nacimiento-muerte, dibujando tiempos de espera, pero estos tiempos de espera pueden ser tiempos de espera para el siguiente cambio de estado de personaje, especiación o evento de extinción. En particular, imagine que hay n linajes presentes en el tiempo t, y que k de estos linajes están en el estado 0 (y n − k están en el estado 1). El tiempo de espera para el siguiente evento seguirá una distribución exponencial con un parámetro de tasa de:
\[ρ = k(q_{01} + λ_0 + μ_0)+(n − k)(q_{10} + λ_1 + μ_1) \label{13.1}\]
Esta ecuación dice que la tasa total de eventos es la suma de los eventos que pueden suceder a linajes con estado 0 (cambio de estado a 1, especiación o extinción) y los eventos análogos que pueden suceder a linajes con estado 1. Una vez que tenemos un tiempo de espera, podemos asignar un tipo de evento dependiendo de las probabilidades. Por ejemplo, la probabilidad de que el evento sea un cambio de estado de carácter de 0 a 1 es:
\[p_{q_{01}} = (n ⋅ q_{01})/ρ\label{13.2}\]
Y la probabilidad de que el evento sea la extinción de un linaje con estado de carácter 1 es:
\[p_{μ_1} = \dfrac{(n − k)⋅μ_1}{ρ} \label{13.3}\]
Y así sucesivamente para los otros cuatro posibles eventos.
Una vez que hayamos elegido un evento de esta manera, podemos asignarlo aleatoriamente a uno de los linajes en el estado apropiado, siendo cada linaje igualmente probable que sea elegido. Luego procedemos hacia adelante en el tiempo hasta que tengamos un conjunto de datos con el tamaño deseado o profundidad de tiempo total.
Un ejemplo de simulación se muestra en la Figura 13.1. Como puede ver, bajo estos parámetros del modelo el impacto de los estados de carácter en la diversificación es fácilmente evidente. En la siguiente sección descubriremos cómo extraer esa información de nuestros datos.