Other stats from the model
θ = 980321.63
L = sum of all the sites used in the site frequency spectrum = 11288846.4
u = mutation rate = 2.9*10-9 (from Heliconius melpomene )
Generation time =0.3