Differences

This shows you the differences between two versions of the page.

--- 진화생물학:한곳_짝짓기_경쟁 [2019/03/29 15:53] – [반수이배체] admin
+++ 진화생물학:한곳_짝짓기_경쟁 [2023/09/05 15:46] (current) – external edit 127.0.0.1
@@ Line 18: / Line 18: @@
    - $SR$  암컷  $\times$   $RR$  수컷: 이렇게 만나는 사건의 확률은
-$$ᅦP_{SR \times RR} = \left[\frac{(1-s)K/2}{(N-1)(1-r)K + (1-s)K}\right] \times \left[\frac{(N-1)rK+sK/2}{(N-1)rK + sK}\right]$$
+ $P_{SR \times RR} = \left[\frac{(1-s)K/2}{(N-1)(1-r)K + (1-s)K}\right] \times \left[\frac{(N-1)rK+sK/2}{(N-1)rK + sK}\right]$
 인데, 암컷의  $SR$  중에서  $S$ 가 나와야 하므로 다시 확률  $1/2$ 을 곱한다.
    - $RR$  암컷  $\times$   $SR$  수컷: 마찬가지로 이 사건의 확률은 아래와 같은데
@@ Line 45: / Line 45: @@
 여기에서는  $R$ 과  $S$  중  $S$ 가 열성인 돌연변이 유전자로서 암컷에 있을 때에 작용한다고 하자. 즉  $RR$ 이나  $RS$ 인 암컷은 수컷을 낳을 비율이  $r$ 이고  $SS$ 인 암컷은 수컷을 낳을 비율이  $s$ 이다.
-위에서처럼 하나의 패치에  $N$  마리의 암컷이 알을 낳을 텐데, 돌연변이가 희귀하기 때문에 이 중  $(N-1)$  마리는  $RR$   $R$ 이 만났던 경우이고 나머지 한 마리에 대해서만 돌연변이가 들어올 수도 있다고 한다. 이 마지막  $N$  번째 알은  $F$ 인 암컷과  $M$ 인 수컷이 만나는 다음과 같은 다섯 가지의 가능성이 있다 (암수가 만났더라도 알은 수정이 되지 않아 반수체 수컷이 나올 수도 있고 이것이 위에 적은  $r$ 이나  $s$ 로 결정된다):  $SS \times S$ ,  $SS \times R$ ,  $RS \times S$ ,  $RS \times R$ , 그리고  $RR \times S$ . 이 상태들의 수  $x_{F \times M}$ 이 어떻게 전이해 가는지를 추적하는 것이 관건이다.
+위에서처럼 하나의 패치에  $N$  마리의 암컷이 알을 낳을 텐데, 돌연변이가 희귀하기 때문에 이 중  $(N-1)$  마리는  $RR$ 과  $R$ 이 만났던 경우이고 나머지 한 마리에 대해서만 돌연변이가 들어올 수도 있다고 한다. 이 마지막  $N$  번째 알은  $F$ 인 암컷과  $M$ 인 수컷이 만나는 다음과 같은 다섯 가지의 가능성이 있다 (암수가 만났더라도 알은 수정이 되지 않아 반수체 수컷이 나올 수도 있고 이것이 위에 적은  $r$ 이나  $s$ 로 결정된다):  $SS \times S$ ,  $SS \times R$ ,  $RS \times S$ ,  $RS \times R$ , 그리고  $RR \times S$ . 이 상태들의 수  $x_{F \times M}$ 이 어떻게 전이해 가는지를 추적하는 것이 관건이다.
   *먼저  $SS \times S$ 가 하나 있을 때 몇 개의  $SS \times S$ 가 나올지 생각해보자. 암컷 한 마리당  $K$  마리의 자식을 낳는다고 하면, 이 패치에서 태어나는 암컷의 수는  $(N-1)(1-r)K + (1-s)K$ 이고 이 중  $(N-1)(1-r)K$ 는  $RR$ 인 암컷, 나머지  $(1-s)K$ 는  $SS$ 인 암컷이다. 마찬가지로 수컷의 수는  $(N-1)rK + sK$ 로서  $(N-1)rK$ 는  $R$ 인 수컷, 나머지  $sK$ 는  $S$ 인 수컷이다. 따라서 이들 중 다시  $SS \times S$ 가 나오는 수는 (암컷 중  $SS$ 의 수) 곱하기 (수컷 중  $S$ 의 비율 $\equiv p$ )로서 아래와 같다:
@@ Line 90: / Line 90: @@
 ======검증======
+위 내용의 검증을 위해 이 군집에 두 쌍만이 존재한다고 가정하고 계산을 수행해보자. 먼저 군집에 존재할 수 있는 두 쌍의 경우의 수는 총 21개이고, 각각의 비율을  $x_1,...,x_{21}$ 이라고 하자.
+$$(RRXR)(RRXR) = x_1,\quad (RRXS)(RRXR) = x_2,\quad (RSXR)(RRXR) = x_3,\quad (RSXS)(RRXR) = x_4\\
+(SSXR)(RRXR)=x_5,\quad (RSXR)(RRXS)=x_6,\quad  (RRXS)(RRXS)=x_7,\quad  (RSXR)(RSXR)=x_8\\
+(SSXS)(RRXR)=x_9,\quad  (RSXS)(RRXS)=x_{10},\quad (SSXR)(RRXS)=x_{11},\quad (RSXS)(RSXR)=x_{12}\\
+(SSXR)(RSXR)=x_{13},\quad (SSXS)(RRXS)=x_{14},\quad (SSXS)(RSXR)=x_{15},\quad (RSXS)(RSXS)=x_{16}\\
+(SSXR)(RSXS)=x_{17},\quad (SSXR)(SSXR)=x_{18},\quad (SSXS)(RSXS)=x_{19},\quad (SSXS)(SSXR)=x_{20}\\
+(SSXS)(SSXS)=x_{21}$$
+그리고 이들의 자녀 세대가 무작위 짝짓기를 한다고 가정하고, 자녀 세대에서 나올 수 있는 각 쌍의 수는 다음과 같다.
+$$
+N_{RRXR} = (1-r)K\left[2x_1+x_2+\frac{9}{8}x_3+ \frac{3}{4}x_4+\frac{3}{8}x_6+\frac{1}{2}x_8+\frac{1}{4}x_{12}+\frac{r}{r+s}\left(x_5+x_9+\frac{1}{4}x_{13}+\frac{1}{4}x_{15}\right)\right]\\
+N_{RRXS} = (1-r)K\left[\frac{3}{8}x_3+\frac{1}{4}x_4+\frac{1}{8}x_6+\frac{1}{2}x_8+\frac{1}{4}x_{12}+\frac{r}{r+s}\left(x_5+x_9\right)+\frac{2s+r}{4(s+r)}\left(x_{13}+x_{15}\right)\right]\\
+N_{RSXR} = (1-r)K\left[x_2+\frac{3}{8}x_3+\frac{3}{8}x4+\frac{9}{8}x_6+2x_7+\frac{1}{2}x_8+\frac{9}{8}x_{10}+\frac{1}{2}x_{12}+\frac{1}{2}x_{16}+\frac{r}{r+s}\left(x_{11}+x{14}+\frac{1}{4}x_{15}+\frac{1}{4}x_{17}+\frac{1}{4}x_{19}+\right)\right]\\+(1-s)K\frac{r}{r+s}\left[x_5+x_{11}+\frac{1}{2}x_{13}+\frac{1}{2}x_{17}\right]\\
+N_{RSXS} = (1-r)K\left[\frac{1}{8}x_3+\frac{1}{8}x_4+\frac{3}{8}x_6+\frac{1}{2}x_8+\frac{3}{8}x_{10}+\frac{1}{2}x_{12}+\frac{1}{2}x_{16}+\frac{s}{r+s}\left(x_{11}+x_{14}\right)+\frac{2s+r}{4(s+r)}\left(x_{13}+x_{15}+x_{17}+x_{19}\right)\right]\\+(1-s)K\left[2x_{18}+x_{20}+\frac{s}{r+s}\left(x_5+x_{11}\right)+\frac{2s+r}{2(s+r)}\left(x_{13}+x_{17}\right)\right]\\
+N_{SSXR} = (1-r)K\left[\frac{3}{8}x_4+\frac{3}{8}x_{10}+\frac{1}{4}x_{12}+\frac{1}{2}x_{16}+\frac{r}{4(s+r)}\left(x_{17}+x_{19}\right)\right]+(1-s)K\frac{r}{r+s}\left[x_9+x_{14}+\frac{1}{2}x_{15}+\frac{1}{2}x_{19}\right]\\
+N_{SSXS} = (1-r)K\left[\frac{1}{8}x_4+\frac{1}{8}x_{10}+\frac{1}{4}x_{12}+\frac{1}{2}x_{16}+\frac{2s+r}{4(s+r)}\left(x_{17}+x_{19}\right)\right]+(1-s)K\left[x_{10}+2x_{21}+\frac{s}{s+r}\left(x_9+x_{14}\right)+\frac{2s+r}{2(r+s)}\left(x_{15}+x_{19}\right)\right]
+$$
+이 때 전체 쌍의 수를  $N_{tot} = N_{RRXR}+N_{RRXS}+N_{RSXR}+N_{RSXS}+N_{SSXR}+N_{SSXS}$  이라고 하고 각 쌍의 비율을
+$$P_{RRXR} = \frac{N_{RRXR}}{N_{tot}},\enspace P_{RRXS} = \frac{N_{RRXS}}{N_{tot}},\enspace P_{RSXR} = \frac{N_{RSXR}}{N_{tot}},\enspace
+P_{RSXS} = \frac{N_{RSXS}}{N_{tot}},\enspace P_{SSXR} = \frac{N_{SSXR}}{N_{tot}},\enspace P_{SSXS} = \frac{N_{SSXS}}{N_{tot}}$$
+이라고 한다면 자녀 세대에서 특정 두 쌍이 존재할 확률을 다음과 같이 계산할 수 있다.
+$$
+x'_{1} = P_{RRXR}\times P_{RRXR},\quad x'_{2} = P_{RRXS}\times P_{RRXR},\quad x'_{3} = P_{RSXR}\times P_{RRXR},\quad x'_{4} = P_{RSXS}\times P_{RRXR}\\
+x'_{5} = P_{SSXR}\times P_{RRXR},\quad x'_{6} = P_{RSXR}\times P_{RRXS},\quad x'_{7} = P_{RRXS}\times P_{RRXS},\quad x'_{8} = P_{RSXR}\times P_{RSXR}\\
+x'_{9} = P_{SSXS}\times P_{RRXR},\quad x'_{10} = P_{RSXS}\times P_{RRXS},\quad x'_{11} = P_{SSXR}\times P_{RRXS},\quad x'_{12} = P_{RSXS}\times P_{RSXR}\\
+x'_{13} = P_{SSXR}\times P_{RSXR},\quad x'_{14} = P_{SSXS}\times P_{RRXS},\quad x'_{15} = P_{SSXS}\times P_{RSXR},\quad x'_{16} = P_{RSXS}\times P_{RSXS}\\
+x'_{17} = P_{SSXR}\times P_{RSXS},\quad x'_{18} = P_{SSXR}\times P_{SSXR},\quad x'_{19} = P_{SSXS}\times P_{RSXS},\quad x'_{20} = P_{SSXS}\times P_{SSXR}\\
+x'_{21} = P_{SSXS}\times P_{SSXS}
+$$
+이로써 자녀 세대의 쌍의 비율  ${x'}$ 를 부모 세대의 쌍의 비율  ${x}$ 들로 표현해내었다.
+지금 우리는 돌연변이가 희귀한 경우에 대해 관심이 있기 때문에  $x_1$ 을 제외한 나머지  $x_i$ 들이 매우 작다고 가정하고 다음과 같이  ${x'}$ 를  $x_i$ 에 대해 전개할 것이다.
+$$
+\begin{cases}
+x'_1 = f_1(x_1,...,x_{21}) \approx \left.\frac{\partial f_{1}}{\partial x_1}\right\vert_{(1,0,...,0)} (x_1 - 1) + \cdots + \left.\frac{\partial f_{1}}{\partial x_{21}}\right\vert_{(1,0,...,0)} x_{21}\\
+x'_2 = f_2(x_1,...,x_{21}) \approx \left.\frac{\partial f_{2}}{\partial x_1}\right\vert_{(1,0,...,0)} (x_1 - 1) + \cdots + \left.\frac{\partial f_{2}}{\partial x_{21}}\right\vert_{(1,0,...,0)} x_{21}\\
+\qquad\qquad\cdot\\
+\qquad\qquad\cdot\\
+\qquad\qquad\cdot\\
+x'_{21} = f_{21}(x_1,...,x_{21}) \approx \left.\frac{\partial f_{21}}{\partial x_1}\right\vert_{(1,0,...,0)} (x_1 - 1) + \cdots + \left.\frac{\partial f_{21}}{\partial x_{21}}\right\vert_{(1,0,...,0)} x_{21}\\
+\end{cases}
+$$
+이제 이 20x20 행렬의 가장 큰 고윳값을  $r$ 과  $s$ 에 따라 그리고, 고윳값이 1보다 작은 영역을 보면 다음을 얻는다.
+{{ :진화생물학:out.png?600 |}}
+이로써  $N=2$ 일 때 평형 성비는  $r=0.214$ 가 됨을 알 수 있다.
+그리고 한 쌍이  $(RRXR)$ 이라고 가정하고  $x'_2,x'_3,x'_4,x'_5,x'_9$ 를 가지고 5x5 행렬을 써보면 이 행렬이 윗 절의
+$$
+\begin{pmatrix}
+pu & 0 & \frac{1}{4N} & 0 & 0\\
+(1-p)u & 0 & \frac{2N-1}{4N} & 0 & 0\\
+& pu & \frac{1}{4N} & \frac{1}{4N} & 0\\
+& (1-p)u & \frac{2N-1}{4N} & \frac{2N-1}{4N} & 1\\
+(N-1)p & (N-1)p & \frac{N-1}{2N} & \frac{2N-1}{4N} & 0
+\end{pmatrix}
+$$
+와 같음을 확인할 수 있다.
 ======참고문헌======
+  * W. D. Hamilton, //Extraordinary Sex Ratios//, Science **156**, 477--488 (1967).
+  * P. D. Taylor and M. G. Bulmer, //Local Mate Competition and the Sex Ratio//, J. Theor. Biol. **86** 409--419 (1980).