Find the expected time of fixation given an initial proportion.











up vote
0
down vote

favorite












Consider the Wright-Fisher model with selection. This means that for a given $lambda>-1$ and $forall n$ the probability of any individual to choose parent of type $A$ is $frac{(1+s)p}{1+sp}$ subject to the condition that the $p= frac{X_n}{N},$ where $X_n$ is the number of individuals with type $A$ allele at generation $n.$



Define the event $tau_N = inf{ngeq 0:X_n=0text{ or }X_n=H}.$ My goal is to find or estimate $E[tau_N|p(0)=x]$ where $p(0)$ denotes the initial proportion of $A$ allele and $xin [0,1].$



Say $xin (0,1)$ in order to avoid fixation at $t=0.$ Then, using conditional expectation we get that
$$E[tau_N|p(0)=x]=sum_{t=0}^{infty}tP[tau_N=t|p(0)=x]=
sum_{t=0}^{infty}tP[X_t=0|p(0)=x]+sum_{t=0}^{infty}tP[X_t=N|p(0)=x].
$$



I don't know how to proceed after this step. Any hints/suggestions will be much appreciated.










share|cite|improve this question




























    up vote
    0
    down vote

    favorite












    Consider the Wright-Fisher model with selection. This means that for a given $lambda>-1$ and $forall n$ the probability of any individual to choose parent of type $A$ is $frac{(1+s)p}{1+sp}$ subject to the condition that the $p= frac{X_n}{N},$ where $X_n$ is the number of individuals with type $A$ allele at generation $n.$



    Define the event $tau_N = inf{ngeq 0:X_n=0text{ or }X_n=H}.$ My goal is to find or estimate $E[tau_N|p(0)=x]$ where $p(0)$ denotes the initial proportion of $A$ allele and $xin [0,1].$



    Say $xin (0,1)$ in order to avoid fixation at $t=0.$ Then, using conditional expectation we get that
    $$E[tau_N|p(0)=x]=sum_{t=0}^{infty}tP[tau_N=t|p(0)=x]=
    sum_{t=0}^{infty}tP[X_t=0|p(0)=x]+sum_{t=0}^{infty}tP[X_t=N|p(0)=x].
    $$



    I don't know how to proceed after this step. Any hints/suggestions will be much appreciated.










    share|cite|improve this question


























      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      Consider the Wright-Fisher model with selection. This means that for a given $lambda>-1$ and $forall n$ the probability of any individual to choose parent of type $A$ is $frac{(1+s)p}{1+sp}$ subject to the condition that the $p= frac{X_n}{N},$ where $X_n$ is the number of individuals with type $A$ allele at generation $n.$



      Define the event $tau_N = inf{ngeq 0:X_n=0text{ or }X_n=H}.$ My goal is to find or estimate $E[tau_N|p(0)=x]$ where $p(0)$ denotes the initial proportion of $A$ allele and $xin [0,1].$



      Say $xin (0,1)$ in order to avoid fixation at $t=0.$ Then, using conditional expectation we get that
      $$E[tau_N|p(0)=x]=sum_{t=0}^{infty}tP[tau_N=t|p(0)=x]=
      sum_{t=0}^{infty}tP[X_t=0|p(0)=x]+sum_{t=0}^{infty}tP[X_t=N|p(0)=x].
      $$



      I don't know how to proceed after this step. Any hints/suggestions will be much appreciated.










      share|cite|improve this question















      Consider the Wright-Fisher model with selection. This means that for a given $lambda>-1$ and $forall n$ the probability of any individual to choose parent of type $A$ is $frac{(1+s)p}{1+sp}$ subject to the condition that the $p= frac{X_n}{N},$ where $X_n$ is the number of individuals with type $A$ allele at generation $n.$



      Define the event $tau_N = inf{ngeq 0:X_n=0text{ or }X_n=H}.$ My goal is to find or estimate $E[tau_N|p(0)=x]$ where $p(0)$ denotes the initial proportion of $A$ allele and $xin [0,1].$



      Say $xin (0,1)$ in order to avoid fixation at $t=0.$ Then, using conditional expectation we get that
      $$E[tau_N|p(0)=x]=sum_{t=0}^{infty}tP[tau_N=t|p(0)=x]=
      sum_{t=0}^{infty}tP[X_t=0|p(0)=x]+sum_{t=0}^{infty}tP[X_t=N|p(0)=x].
      $$



      I don't know how to proceed after this step. Any hints/suggestions will be much appreciated.







      probability






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Nov 14 at 15:39

























      asked Nov 12 at 17:22









      Hello_World

      3,69821630




      3,69821630



























          active

          oldest

          votes











          Your Answer





          StackExchange.ifUsing("editor", function () {
          return StackExchange.using("mathjaxEditing", function () {
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          });
          });
          }, "mathjax-editing");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "69"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });














           

          draft saved


          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2995583%2ffind-the-expected-time-of-fixation-given-an-initial-proportion%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown






























          active

          oldest

          votes













          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes
















           

          draft saved


          draft discarded



















































           


          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2995583%2ffind-the-expected-time-of-fixation-given-an-initial-proportion%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          How do I know what Microsoft account the skydrive app is syncing to?

          When does type information flow backwards in C++?

          Grease: Live!