deltaSE() function to calculate approximate standard errors for functions of (co)variance parameters (e.g., h2, standard deviations of variances, or correlations)
Rcon arguments to
gremlin() for constraining parameters
siremodel, we could restrain the
grSf <- gremlin(WWG11 ~ sex, random= ~ sire, data = Mrode11, Gstart = list(matrix(0.38)), Gcon = list("F"), control = gremlinControl(lambda = FALSE))
Similar to above change (
Rcon), introduced steps to deal with parameters outside of the boundaries of their parameter space (e.g., variance < 0).
change version numbering to just 3 numbers (instead of 4)
Implement "step-halving" algorithm for AI updates
gremlinControl()function for advanced changes to the way gremlin runs
ai, and elsewhere (where relevant) in next version
Completely revised way models are built and called
grModis the model structure for which a log-likelihood can be calculated
gremlinRclass distinguishes from
gremlinclass in that
gremlinRobjects will only use
Rcode written by the package in order to run the model. Class
gremlinwill execute underlying c++ code written in the package.
Average Information algorithm has been vastly improved
ai()efficiently calculates the AI matrix without directly computing several matrix inverses (as previously coded)
lambda and alternative parameterizations now possible and executed by the same code
lambdaparameterization is the REML likelihood of the variance ratios after factoring out a residual variance from the Mixed Model Equations.
M) matrix from which the Cholesky factorization (and
tyPycalculations are made)
C) and obtain
C, now do a
solvewith Cholesky of
c++code) to calculate
tyPybased off Boldman and Van Vleck
logLikmethods as well
Improved algorithm that reduces computational resources and time! Also implemented c++ code in
gremlin(), while keeping
gremlinR() purely the R implementation (at least from the package writing standpoint).
Documentation has switched from filling out the
.Rd files manually to providing
documentation next to the function code in the
.R files using
Congratulations, its a gremlin!
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.