chr18 78077248 1.5075956950023575e-08
chr19 59128983 1.8220141872466202e-08
chr20 63025520 1.7178269031631664e-08
chr21 48129895 1.3045214034879191e-08
chr22 51304566 1.4445022767788226e-08
chrX 155270560 1.164662223273842e-08
chrY 59373566 0.0
"""
_chromosomes = []
for line in _chromosome_data.splitlines():
name, length, mean_rr = line.split()[:3]
_chromosomes.append(
genomes.Chromosome(
name=name,
length=int(length),
default_mutation_rate=1e-8, # WRONG!,
default_recombination_rate=float(mean_rr)))
#: :class:`stdpopsim.Genome` definition for humans.
genome = genomes.Genome(species="homo_sapiens",
chromosomes=_chromosomes,
default_genetic_map=HapmapII_GRCh37.name)
#
# Experimental interface used to develop the CLI.
#
class TempChromosome(object):
"""
chr2L 23513712
chr2R 25286936
chr3L 28110227
chr3R 32079331
chr4 1348131
chrY 3667352
chrM 19524
"""
_chromosomes = []
for line in _chromosome_data.splitlines():
name, length = line.split()[:2]
_chromosomes.append(
genomes.Chromosome(
name=name,
length=int(length),
default_mutation_rate=8.4e-9, # WRONG!, underestimate used in S&S
default_recombination_rate=8.4e-9)
) # WRONG, underestimate used in S&S!
#: :class:`stdpopsim.Genome` definition for D. melanogaster. Chromosome length data is
#: based on `dm6 <https://www.ncbi.nlm.nih.gov/assembly/GCF_000001215.4/>`_.
genome = genomes.Genome(species="drosophila_melanogaster",
chromosomes=_chromosomes,
default_genetic_map=Comeron2012_dm6.name)
###########################################################
#
# Demographic models
#
###########################################################
_chromosome_data = """\
chr1 30427671
chr2 19698289
chr3 23459830
chr4 18585056
chr5 26975502
"""
# mutation rate from Ossowski 2010 Science
# recombination value from Huber et al 2014 MBE
# rho=200/Mb, assume Ne=124,000, rho=2*Ne*r
_chromosomes = []
for line in _chromosome_data.splitlines():
name, length = line.split()[:2]
_chromosomes.append(
genomes.Chromosome(name=name,
length=int(length),
default_mutation_rate=7e-9,
default_recombination_rate=8.1e-9))
genome = genomes.Genome(species="arabidopsis_thaliana",
chromosomes=_chromosomes,
default_genetic_map=Salome2012.name)
###########################################################
#
# Demographic models
#
###########################################################
class Durvasula2017MSMC(models.Model):
"""
"""
import msprime
import stdpopsim.models as models
import stdpopsim.genomes as genomes
###########################################################
#
# Genome definition
#
###########################################################
_chromosomes = []
_chromosomes.append(
genomes.Chromosome(name="chr",
length=4641652,
default_mutation_rate=1e-5 + 2e-4,
default_recombination_rate=0.0))
# mean_conversion_rate=8.9e-11 # not implemented yet!
# mean_conversion_length=542 # not implemented yet!
#: :class:`stdpopsim.Genome` definition for E. Coli.
# Chromosome length data is based on strain K-12.
genome = genomes.Genome(species="e_coli",
chromosomes=_chromosomes,
default_genetic_map=None)
###########################################################
#
# Demographic models
#
###########################################################
chr18 80373285
chr19 58617616
chr20 64444167
chr21 46709983
chr22 50818468
chrX 156040895
chrY 57227415
"""
_chromosomes = []
for line in _chromosome_data.splitlines():
name, length = line.split()[:2]
_chromosomes.append(
genomes.Chromosome(
name=name,
length=int(length),
mean_mutation_rate=1e-8, # WRONG!,
mean_recombination_rate=1e-8)) # WRONG!
#: :class:`stdpopsim.Genome` definition for humans. Chromosome length data is
#: based on `GRCh38 <https://www.ncbi.nlm.nih.gov/grc/human/data>`_.
genome = genomes.Genome(species="homo_sapiens",
chromosomes=_chromosomes,
default_genetic_map=HapmapII_GRCh37.name)
###########################################################
#
# Demographic models
#
###########################################################