*exported meanEBV
and added varEBV
to
complement meanP
/varP
and
meanG
/varG
*Changed all parameters of the CATTLE demographic model to exactly match Macleod et al. (2013) - specifically reducing the mutation rate from 2.5e-8 (from human literature) to 1.2e-8 (used in Macleod et al., 2013) and recombination rate from 1e-8 (generic) to 9.26e-9 (used in Macleod et al., 2013). These changes will reduce number of segregating sites to ~240K per chromosome for 100 samples and will run faster.
*changed misc slot in Pop class from a list organised as ind x nodes to to a list organised as nodes x ind (this simplified code and increased speed)
*removed setMisc
and getMisc
because the
new misc slot structure makes it easy to set and get misc components
with base R code
*added length
method for Pop class that returns number
of individuals (like nInd
)
*added length
method for MultiPop class that returns
number of populations
*fixed bug in quadrivalent pairing resulting in distribution of double reductions not respecting the centromere
*fixed bug in SimParam$restrSegSites
with excluding
sites at end of chromosome
*fix SimParam examples for CRAN
*deleted bad example code for setMisc
*changed examples to use a single thread for CRAN testing this change is not shown in the documentation
*renamed MegaPop
to MultiPop
*fixed bug in writePlink
to correctly export map
positions in cM
*fixed bug in writeRecords
due to removed reps slot in
pops
*added altAddTraitAD
for specifying traits with
dominance effects using dominance variance and inbreeding depression
*add miscPop slot to class Pop
*updated MaCS citation to https site
*Changed citation to use bibentry
instead of
citEntry
*fixed a bug in IBD tracking
*add setFounderHap
to SimParam for applying custom
haplotypes to founders
*added addSnpChipByName
to SimParam for defining SNP
chips by marker names
*changed C++ using sprintf
to use
snprintf
*fixed bug in calculation of genic variance
*fixed importHaplo
not passing ploidy to
newMapPop
*fixed bug with correlated error variances
*fixed column name bug with multiple traits in
setEBV
*fixed CTD caused by runMacs
when too many segSites are
requested
*fixed missing names in GV when using resetPop
*fixed bug in importTrait
*popVar
now deals with matrices having 1 row
*updated link to Gaynor, 2017
*added ability to exclude loci by name in
SimParam$restrSegSites
*pullMarkerGeno
and pullMarkerHaplo
now
work with a MapPop class
*added setMarkerHaplo
to manually change genotypes in a
Pop or MapPop
*added addSegSite
for manually adding segregating sites
to a MapPop class
*simParam$setCorE
has been deprecated in favor of a corE
argument in simParam$setVarE
*setPheno
now takes corE as an argument
*setPheno
now allows the user to set phenotypes for a
subset of traits
*add newEmptyPop
to create populations with zero
individuals
*removed reps slot from populations and heterogeneous residual variance GS models
*added h2, H2, and corE to setPhenoGCA
and
setPhenoProgTest
*the “EUROPEAN” species history was removed from runMacs
due to lengthy runtime
*added getPed
to quick extract a population’s
pedigree
*added getGenMap
to pull a genetic map in data.frame
format
*fixed bugs relating to importData
functions
*fixed writePlink
errors and no longer requires equal
length chromosomes
*added importGenMap
to format genetic maps for
AlphaSimR
*added importInbredGeno
and importHaplo
to
make it easier to create a simulation from external data
*added importSnpChip
, importTrait
to
SimParam
to make it easier to manually define traits
*added pullMarkerGeno
and pullMarkerHaplo
to make it easier to extract genotypes and haplotypes of specific loci
without defining a trait or SNP chip
*reduceGenome
, mergeGenome
and
doubleGenome
should really now work with pedigree and
recombination tracking
*added missing #ifdef _OPENMP to OCS.cpp
*removed use of PI variable in C++ code due to it being compiler specific
*added snpChip argument to pullIbdHaplo
for backwards
compatibility
*exposed internal mixed model solvers
*all selection functions now return a warning when there are not enough individuals
*fixed error in pullIbdHaplo
when chr isn’t NULL
*fixed an error with assigning 1 QTL and/or SNP
*changed geno slot from matrix to list to support future RcppArmadillo changes
*doubleGenome
and reduceGenome
now work
with IBD tracking
*fixed errors in implementation of Gamma Sprinkling model
*fixed formatting error in genetic maps created by runMacs that broke genotype extraction functions
*added h2 and H2 to setPhenoGCA
*pullGeno
and pullHaplo
functions now
report marker names from the genetic map
*removed lazyData field in DESCRIPTION
*AlphaSimR manuscript has been published in G3 (citation added)
*changed to a Gamma Sprinkling model for crossovers, default is still a Gamma model
*change default interference parameter (v) to 2.6 to be consistent with the Kosambi mapping function (was 1, consistent with the Haldane mapping function)
*new internal id (iid) that allows user to freely change id slot in populations
*runMacs2
now adjusts Ne for autopolyploids
*parent populations are now passed to finalizePop
*check added that throws an error when use of discontinued “gender” argument is detected
*added experimental MegaPop-class
*references to gender have been changed to the more appropriate terms sex or sexes
*added misc slot to populations
*added finalizePop
to SimParam
*added physical positions to getSnpMap
and
getQtlMap
*you can now use h2 and H2 to specify error variance in
setPheno
*SimParam$setVarE
now accepts a matrix for varE
*fixed a bug in editGenome
when making multiple
edits
*adding merging of centromere vector in cChr
*GxE traits now default to random sampling of p-values
*fixed a bug in restrSegSites
*fixed a bug in selection of segSites
*changed output of genParam
to match Bulmer, 1976
*nProgeny added to makeCross
and
makeCross2
*all SimParam
documentation is now in
?SimParam
*non-overlapping QTL and SNP is now the default
*new interface for restrSegSites
in
SimParam
*fixed subset by id for populations
*fixed major bug in newMapPop
*switched to a circular design for the balance option in
randCross
and randCross2
*added reduceGenome
and doubleGenome
for
changing ploidy levels
*added minSnpFreq to SimParam_addSnpChip for any reference population
*the c
function now merges individuals for MapPop
objects (was chromosomes before)
*the cChr
function new merges chromosomes for MapPop
objects
*fixed broken SimParam_addStructuredSnpChip
*removed broken pullMultipleSnpGeno
and
pullMultipleSnpHaplo
*fixed broken writePlink
*rework of setEBV
(breaks some scripts)
*genotype data now stored as bits (was bytes)
*implemented a gamma model for crossover interference
*added the mutate function to model random mutations
*added a vignette explaining the biological model for traits
*GS models now handle polyploids
*heterogenous error variance is now optional in GS models (default is homogeneous error)
*improved gene drop functionality of pedigreeCross
*added keepParents option to makeDH and self (indirectly extends
selectFam
and selectWithinFam
)
*added RRBLUP_SCA2
*set methods for the “show” function when applied to populations
*fixed a bug returning the first individual when selecting 0
*fixed error in recombination track when using
makeDH
*fixed error causing epistatic effects to mask GxE effects
*fixed an error with pullSegSiteGeno
and
pullSegSiteHaplo
with variable number of sites per
chromosome
*added traits with epistasis
*Max number of threads automatically detected
*added RRBLUP_D2
*added version tracking to SimParam
*removed trackHaploPop
(super-ceded by
pullIbdHaplo
)
*added fastRRBLUP
*fixed faulty double crossover logic
*fixed broken writePlink
*fixed broken pullIbdHaplo
*mergePops
no longer assumes diploidy
*added support for autopolyploids
*added RRBLUP_GCA2
*randCross2
can now “balance” crossing when not using
gender
*fixed recombination tracking bug in createDH2
*removed bug in setEBV
with append=TRUE
*fixed ambiguous overloading in optimize.cpp
*setPheno
(not setPhenoGCA
) passes the
number of reps to populations
*fixed bug in editGenomeTopQtl
*fixed bug in RRBLUP_D
*fixed bug in resetPop
*fixed bug in SimParam_rescaleTraits
*removed unimplemented SimParam_restrSnpSites and SimParam_restrQtlSites
*add error message for no traits in calcGCA
*added GxE traits with zero environmental variance
*faster trait scaling
*faster calculation of genetic values
*dsyevr now called via arma_fortran
*added OpenMP support
*parallelized cross2
*parallelized runMacs
*parallelized calculation of genetic values
*variance calculations now account for inbreeding
*fixes for male selection in selectOP
*add fixEff to setPhenoGCA
*added default runMacs
option to return all segSites
*added ability to specify separate male and female genetic maps
*pullGeno
and pullHaplo
functions can now
specify chromosomes
*added RRBLUP2
for special GS cases
*improved speed by replacing Rcpp random number generators
*changed available MaCS species
*GS functions now use populations directly
*added pullIbdHaplo
*added writePlink
*fixed population sub-setting checks to prevent invalid selections
*fixed slow calcGCA
*fixed error in addTraitAG
preventing multiple
traits
*fixed bug with mergePops
when merging ebv
*fixed bug in setVarE
when using H2 and multiple
traits
*selectFam
now handles half-sib families
*selectWithinFa
m now handles half-sib families
*Removed restriction on varE=NULL in setPhenoGCA
*Added NEWS file
*Added selectOP
to model selection in open pollinating
plants
*Added runMacs2
as a wrapper for
runMacs
*Fixed error when using H2 in SimParam_setVarE