PoonLab / twt

Generic framework for the coalescent simulation of pathogen trees within host trees in R
GNU Affero General Public License v3.0
3 stars 2 forks source link

Initial Comp at epidemic start time w/ unknown infection time; unable to infer its lineages' coalescent waiting times #46

Closed gtng92 closed 5 years ago

gtng92 commented 5 years ago

Drawing coalescent waiting times from linear population growth dynamics. Current implementation implies that the final piece remains at a constant population size with a slope of zero. First infected compartment in epidemic has no recorded infection time, I originally assumed the initial infection time to be arbitrarily large.

Instead, I'm going to implement it so the maximum time is right before it hits the saturation point (before the population size goes constant). Since there shouldn't be any transmissions possible after this point in time (since we are only doing births and no deaths).

ArtPoon commented 5 years ago

Rewording the problem: this is an issue at the root of the outer tree, where we are at the last infected individual carrying two or more lineages (one of which was just traced back through a transmission event). If individuals are immortal (no death rate specified), then we have way of asserting a time frame of that root individual's infection, so we can't say much about their within-host population dynamics that we need to simulate the coalescent process. So I propose to impose this constraint:

ArtPoon commented 5 years ago

We are proceeding with option 1 for now, until death rates are implemented.