# configuration file for program src/newAge.py # Configuration saved at at Tue Jul 27 16:11:00 2010 # label: Replicates # description: # Number of replicates to run the simulation. Between 2 and 1000 reps = '10' # label: Generations # description: # Time units (generations) of evolution. Max: 1000 time = '40' # label: Burn-in # description: # Lenght of burn-in stage. Between 2 and 1000. burnin = '20' # label: Save stats # description: # Save stats for posterior analysis of distribution saveDecays = 'False' # label: Years to save # description: # List of years to save population in a genepop-format file saveGens = '[0, 1, 2, 3, 6]' # label: Calculate F/pedigree # description: # Calculate F from pedigree. This calculation is expensive. doFPed = 'False' # label: number of Microsats # description: # Number of microsatellite-type loci (multi-allelic). Max: 200 nMSats = '10' # label: maximum of alleles # description: # Maximum number of alleles in microsatellite loci. max:100 maxAllele = '20' # label: Frequencies initalization # description: # Choose model for intialization of allelic frequencies in microsat loci. initFreq = 'DIRICHLET' # label: Frequencies datafile # description: # File name for user-defined frequencies. Used only if initFreq = FILE. fileFreq = '' # label: Number of SNPs # description: # Number of SNP-type loci (bi-allelic). Max: 200 nSNPs = '0' # label: Demographic model # description: # Single pop, Island or Stepping-Stone demoModel = 'SINGLE' # label: Initial N # description: # Initial population size (per deme). Max: 100000 initialN = '430' # label: Number of Demes # description: # Number of demes. Max:100 demes = '5' # label: Migration rate # description: # Migration rate migration = '0.001' # label: Maximum age (max 25) # description: # Maximum age of individuals. Note that 0 is the age of newborns maxAge = '21' # label: Adult age # description: # age of individuals when they become adults adultAge = '2' # label: Initial Age Structure # description: # List of probability of age of individuals. Relative (no need to sum up to 100/1) ASinit = '[19.448, 13.64, 10.574, 8.977, 7.622, 6.471, 5.494, 4.665, 3.961, 3.363, 2.855, 2.424, 2.058, 1.747, 1.484, 1.26, 1.069, 0.908, 0.771, 0.655, 0.556]' # label: Mating mode # description: # Mode of mating: RANDOM: random mating with replacement; MONOGFEMS1: monogamic females with one offspring per female; MONOG: monogamy matingMode = 'MONOGFEMS1' # label: Constant born size # description: # Force population with constant born size over time constSize = 'False' # label: Cohort size N0 # description: # Cohort size. Number of newborns generated if constant born size cohortN = '100' # label: Male survivorship # description: # List of age-specific survivorship in males. Last value must be 0 sMales = '[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.0]' # label: Female survivorship # description: # List of age-specific survivorship in females. Last value must be 0. sFemales = '[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.0]' # label: Male probability of mating # description: # List of age-specific probability of mates in males. They are relative, there is no need to sum up to 1. bMales = '[0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]' # label: Female fecundity (b) # description: # List of age-specific fecundity in females. If the population has not constant born size then mean the average number of offspring per indivdual or each age. bFemales = '[0, 0, 0.46, 0.46, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7]' # label: Alpha males # description: # Enable the alpha males mating system alphaMating = 'False' # label: Prop. offspring by alphas # description: # Proportion of offspring sired by alphas alphaOff = '0.0' # label: Minimum age of alphas # description: # Minimum age of an individual who can become an alpha alphaMinAge = '0' # label: Maximum age of alphas # description: # Maximum age of an individual who can become an alpha alphaMaxAge = '0' # label: Duration of alpha # description: # Duration of alpha trait alphaDuration = '0' # label: Prop. of alphas # description: # Proportion of males being alphas alphaProp = '0.0' # label: Has Carrying capacity? # description: # Environment imposes a population limit hasCarrying = 'False' # label: Carrying capacity # description: # Population limit capacity = '100' # label: Culling # description: # Allow culling culling = 'True' # label: Culling objective # description: # Population size after culling cullObjective = '400' # label: Culling target # description: # Target for culling: ALL: random culling; YOUNG: Random culling of non-adult individuals; ADULTS: random culling of adult individuals cullTarget = 'ALL' # label: Culling mode # description: # Mode for culling: INTERVAL: culling occurs every x years; THRESHOLD: culling occurs when population size is larger than a threshold cullMode = 'THRESHOLD' # label: Culling threshold # description: # Population size threshold to do culling if mode = THRESHOLD cullLevel = '700' # label: Culling Interval # description: # Interval of time between culling events if mode = INTERVAL cullInterval = '0' #The same options can be given by command line options (subject to minor changes) # --gui=False --reps='10' --time='40' --burnin='20' \ # --saveGens='[0, 1, 2, 3, 6]' --nMSats='10' --maxAllele='20' \ # --initFreq='DIRICHLET' --fileFreq='' --nSNPs='0' --demoModel='SINGLE' \ # --initialN='430' --demes='5' --migration='0.001' --maxAge='21' \ # --adultAge='2' --ASinit='[19.448, 13.64, 10.574, 8.977, 7.622, 6.471, \ # 5.494, 4.665, 3.961, 3.363, 2.855, 2.424, 2.058, 1.747, 1.484, 1.26, \ # 1.069, 0.908, 0.771, 0.655, 0.556]' --matingMode='MONOGFEMS1' \ # --cohortN='100' --sMales='[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.0]' --sFemales='[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.0]' --bMales='[0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, \ # 1, 1, 1, 1, 1, 1, 1, 1, 1]' --bFemales='[0, 0, 0.46, 0.46, 0.7, 0.7, \ # 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, \ # 0.7]' --alphaOff='0.0' --alphaMinAge='0' --alphaMaxAge='0' \ # --alphaDuration='0' --alphaProp='0.0' --capacity='100' --culling \ # --cullObjective='400' --cullTarget='ALL' --cullMode='THRESHOLD' \ # --cullLevel='700' --cullInterval='0' #Or a shorter version if default arguments are ignored # --gui=False --reps='10' --time='40' --burnin='20' \ # --saveGens='[0, 1, 2, 3, 6]' --nMSats='10' --maxAllele='20' \ # --initFreq='DIRICHLET' --fileFreq='' --nSNPs='0' --demoModel='SINGLE' \ # --initialN='430' --demes='5' --migration='0.001' --maxAge='21' \ # --adultAge='2' --ASinit='[19.448, 13.64, 10.574, 8.977, 7.622, 6.471, \ # 5.494, 4.665, 3.961, 3.363, 2.855, 2.424, 2.058, 1.747, 1.484, 1.26, \ # 1.069, 0.908, 0.771, 0.655, 0.556]' --matingMode='MONOGFEMS1' \ # --cohortN='100' --sMales='[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.0]' --sFemales='[0.76, 0.84, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, 0.92, \ # 0.92, 0.92, 0.92, 0.0]' --bMales='[0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, \ # 1, 1, 1, 1, 1, 1, 1, 1, 1]' --bFemales='[0, 0, 0.46, 0.46, 0.7, 0.7, \ # 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, \ # 0.7]' --alphaOff='0.0' --alphaMinAge='0' --alphaMaxAge='0' \ # --alphaDuration='0' --alphaProp='0.0' --capacity='100' --culling \ # --cullObjective='400' --cullTarget='ALL' --cullMode='THRESHOLD' \ # --cullLevel='700' --cullInterval='0'