תהליך מורן

תהליך מורן או מודל מורן הוא תהליך סטוכסטי פשוט המשמש בביולוגיה לתיאור אוכלוסיות סופיות. התהליך נקרא על שמו של פטריק מורן, שהציע לראשונה את המודל בשנת 1958. <ref></ref> מודל מאפשר תיאור תהליכים שהם מגדילי מגוון כגון מוטציות וכן השפעות מפחיתות מגוון כגון סחף גנטי וברירה טבעית. תהליך מורן יכול לתאר את הדינמיקה ההסתברותית באוכלוסייה סופית בגודל קבוע <math>N</math> שבה שני אללים <math>A</math> ו-<math>B</math> מתחרים על הדומיננטיות. שני האללים מסוגלים לייצר עותקים של עצמם.

בכל שלב נבחר פרט אקראי (שהוא עם אלל <math>A</math> או אלל <math>B</math>) לרבייה (שכפול של עצמו) ופרט אקראי למוות; כך שגודל האוכלוסייה נשאר קבוע (ניתן לבחור את אותו פרט למוות ולרבייה באותו שלב.) כדי למדל ברירה טבעית, סוג אחד צריך להיות בעל כשירות גבוהה יותר ולכן יש סיכוי גבוה יותר שייבחר עבור רבייה.

== סחף נייטרלי ==
סחף נייטרלי הוא מצב שבו מוטציה נייטרלית יכולה להתפשט בכל אוכלוסייה, כך שבסופו של דבר האלל המקורי יאבד. במוטציה נייטרלית אין שום יתרון או חיסרון בכשירות של בעליה. המקרה הפשוט של תהליך מורן יכול לתאר תופעה זו.

תהליך מורן הוא שרשרת מרקוב זמן בדיד, עם מספר מצבים סופי, ולכן ניתן לתיאור באמצעות מטריצת מעבר. תהליך מורן מוגדר על מרחב מצבים <math>i=0,...,N</math>. הימצאות התהליך בזמן מסוים במצב <math>i</math> אומר שמספר ה-<math>A</math> (מספר נושאי אלל <math>A</math>) הוא <math>i</math>. הרכיב <math>P_{i,j}</math> של המטריצת המעבר שווה להסתברות לעבור ממצב <math>i</math>, בשלב זמן נתון למצב <math>j</math> בשלב הזמן הבא. מכיוון שמספר ה-<math>A</math> לפי המודל יכול להשתנות לכל היותר באחד בכל שלב זמן (שכפול או מוות), קיים מעבר רק בין מצב <math>i</math> למצב <math>i-1</math> ו-<math>i+1</math> (או הישארות במצב <math>i</math>.) לפיכך מטריצת המעבר של התהליך הסטוכסטי היא תלת אלכסונית והסתברויות המעבר עבור <math>0<i<N</math> הן:

: <math>\begin{align}

P_{i,i-1} &= \frac{N-i}{N} \frac{i}{N}\\
P_{i,i} &= 1- P_{i,i-1} - P_{i,i+1}\\
P_{i,i+1} &= \frac{i}{N} \frac{N-i}{N}\\

\end{align}</math>
: כמו כן, <math>P_{0,0}=1</math> ו-<math>P_{N,N}=1</math>, וכל שאר רכיבי המטריצה אפסים. כדי להבין את הסתברויות המעבר יש לחזור להגדרת התהליך הקובעת שבכל שלב ייבחר פרט אחד לשכפול ואחד למוות. ברגע שכל ה-<math>A</math> מתו, הם לא ישובו להיות חלק מהאוכלוסייה. זאת מכיוון שלפי הנחות המודל הדרך היחידה להיווצרות <math>A</math> היא באמצעות שכפול של <math>A</math> שכבר קיים ולכן <math>P_{0,0}=1</math>. מאותה סיבה אוכלוסיית ה-<math>A</math> תמיד תישאר <math>N</math> מהרגע שהם הגיעו למספר הזה והגיעו לרוויה ולכן <math>P_{N,N}=1</math> (במצב זה אין <math>B</math> ולכן אין שכפולי <math>B</math>.) המצבים <math>0</math> ו- <math>N</math> נקראים מצבים סופגים ואילו המצבים <math>1,...,N-1</math> נקראים חולפים. ניתן להבין את הסתברויות המעבר ממצב חולף כמכפלה שלההסתברות לבחור את האלל שהכמות שלו תגדל באחד בהסתברות לבחירה באלל מהסוג השני שהכמות שלו תקטן באחד. ברור שאם אותו אלל נבחר לרבייה ולמוות, אז הכמות של אף אחד מהאללים לא משתנה.

בסופו של דבר התהליך יגיע בהכרח (בהסתברות <math>1</math>) לאחד המצבים הסופגים ואז יישאר שם לנצח. במצבים החולפים, תנודות אקראיות יתרחשו אך בסופו של דבר אוכלוסיית ה-<math>A</math> תיכחד או תשתלט. זהו אחד ההבדלים החשובים תהליכים סטוכסטיים בהשוואה לתהליכים דטרמיניסטיים שהאחרונים לא יכולים למדל מאורעות אקראיים. ניתן לחשב את התוחלת ואת השונות של <math>X(t)</math> (המצב בזמן <math>t</math>) בהינתן המצב ההתחלתי <math>X(0)=i</math>:

: <math> \begin{align}
\operatorname{E}[X(t)\mid X(0) = i] &= i \\
\operatorname{Var}(X(t)\mid X(0) = i) &= \tfrac{2i}{N} \left(1-\tfrac{i}{N} \right ) \frac{1- \left(1-\frac{2}{N^2} \right )^t}{\frac{2}{N^2}}
\end{align}</math> 

בהינתן <math>X(0)=i</math> ההסתברות ש-<math>A</math> יגיע למצב רוויה נקראת הסתברות רוויה . עבור תהליך מורן המתואר הסתברות זו היא <math>p_i=i/N</math>.

מכיוון שלכל הפרטים יש את אותה הכשירות, יש להם גם את אותו הסיכוי להפוך לאב הקדמון של כל האוכלוסייה; לכן ההסתברות לכך היא <math>1/N</math>. תוחלת הזמן להגעה לספיגה בהינתן מצב התחלתי <math>X(0)=i</math> היא:

: <math>k_i = N \left[ \sum_{j=1}^{i} \frac{N-i}{N-j} + \sum_{j=i+1}^{N-1} \frac{i}{j} \right]</math>

עבור <math>N</math> גדול מקבלים

: <math>\lim_{N\to \infty} \frac{k_i}{N^2} = - \left[ (1-x_i) \ln(1-x_i) + x_i \ln(x_i) \right]</math>

== ברירה טבעית ==
אם לאלל אחד יש יתרון כשירות על פני האלל השני, סביר יותר שהוא ייבחר לשכפול. ניתן לשלב זאת במודל אם לפרטים עם אלל A יש כשירות <math>f_i</math> ולפרטים עם אלל B יש כשירות כאשר i הוא מספר ה- <math>A</math>. מטריצת המעבר של התהליך הסטוכסטי היא בצורת תלת אלכסונית והסתברויות המעבר עבור <math>0<i<N</math> הן:

: <math> \begin{align}

P_{i,i-1} &= \frac{g_i (N-i) }{f_i \cdot i + g_i (N-i)} \cdot \frac{i}{N}\\
P_{i,i} &= 1- P_{i,i-1} - P_{i,i+1}\\
P_{i,i+1} &= \frac{f_i \cdot i}{f_i \cdot i + g_i (N-i)} \cdot \frac{N-i}{N}\\

\end{align}</math>

כמו כן, <math>P_{0,0}=1</math> ו-<math>P_{N,N}=1</math>, וכל שאר רכיבי המטריצה אפסים. כדי להבין את הסתברויות המעבר יש לשים לב שהכשירות נכנסת לאותו מרכיב במכפלה שנותן את ההסתברות של כמות אלל מסוים לגדול באחד. לפיכך ההסתברות שפרט <math>A</math> נבחר לשכפול ותלויה בכשירותו של <math>A</math> ולכן שווה ל-

: <math>\frac{f_i \cdot i}{f_i \cdot i + g_i (N-i)}</math>.

ההסתברות <math>p_{i}</math> להגיע למצב רוויה בהינתן שהמצב ההתחלתי הוא i נתונה באמצעות נוסחת נסיגה:

: <math> p_i = \begin{cases} 0 & i=0\\ \beta_i p_{i-1}+(1-\alpha_i-\beta_i) p_i+\alpha_i p_{i+1} & 1 \leq i \leq N-1\\ 1 & i =N
\end{cases}</math>

ונוסחה סגורה ניתנת על ידי

: <math>p_i = \frac \qquad \text{(1)}</math>

כאשר <math>\gamma_i = g_i / f_i</math> .

המקרה הכללי הזה שבו הכשירות של <math>A</math> ו-<math>B</math> תלויים בכמות של כל סוג נחקר בתורת המשחקים האבולוציונית .

מודל פשוט יותר מתקבל אם מניחים יחס כשירות קבוע <math> r</math> . פרטים מסוג <math>A</math> מתרבים בקצב קבוע <math> r</math> ופרטים עם אלל <math>B</math> מתרבים בקצב <math>1</math>. לפיכך אם ל-<math>A</math> יש יתרון כשירות על <math>B</math> אז <math> r</math> יהיה גדול מאחד, ואם ל- <math>B</math> יש יתרון <math> r</math> יהיה קטן מאחד. לפיכך מטריצת המעבר של התהליך הסטוכסטי היא בצורת תלת אלכסונית והסתברויות המעבר הן

: <math> \begin{align}
P_{0,0}&=1\\
P_{i,i-1} &= \frac{N-i}{r \cdot i + N-i} \cdot \frac{i}{N}\\
P_{i,i} &= 1- P_{i,i-1} - P_{i,i+1}\\
P_{i,i+1} &= \frac{r \cdot i}{r \cdot i + N-i} \cdot \frac{N-i}{N}\\
P_{N,N}&=1.
\end{align}</math>

במקרה הזה <math>\gamma_i = 1/r</math> מהווה גורם קבוע לכל הרכב אוכלוסייה ולכן הסתברות הרוויה ממשוואה (1) היא:

: <math>p_i = \frac{1-r^{-i}} { 1-r^{-N} } \quad \Rightarrow \quad p_1 = \rho = \frac{1-r^{-1}} {1-r^{-N}} \qquad \text{(2)}</math>

כאשר הסתברות הרוויה של <math>A</math> בודד באוכלוסייה (כל השאר הם <math>B</math>) היא לעיתים קרובות מעניינת ומסומנת ב- .

במקרה של ברירה טבעית, ניתן גם לחשב את התוחלת והשונות לאחר יחידת זמן אחת מנקודת מוצא של מצב i.

: <math> \begin{align}
\operatorname{E}[ X(t+1) \mid X(t) = i ] &= p s \dfrac{1-p}{p s + 1} + i \\
\operatorname{Var}( X(t+1) \mid X(t)=i) &=p(1-p)\dfrac{ (s+1) + (p s + 1)^2 }{(p s +1)^2}
\end{align}</math>

כאשר <math>p = 1/n</math> ו- <math>s=r-1 </math>.

== הערות שוליים ==

== לקריאה נוספת ==
* 
* 

== קישורים חיצוניים ==
* 

קטגוריה:תהליכים סטוכסטיים