Incomplete check for sampling start time occurring before stationarity; the Watkins approximation in the pytransmission package gives negative times for small population sizes and per locus innovation rates that are small (not a bug, just a feature of using an approximation with small populations).  Commented out.