[R] New version of coxme / lmekin

Terry Therneau therneau at mayo.edu
Fri Dec 30 21:22:09 CET 2011


Version 2.2 of coxme has been posted to CRAN, Windows versions and
mirrors should appear in due course.  This is a major update with three
features of note:

1. A non-upwardly compatable change:

Extractor functions: beta= fixed effects, b=random effects
             nlme    lme4     coxme <2.2   coxme 2.2   lmekin 2.2
           ------------------------------------------------------
beta        fixef    fixef     fixef       fixef         fixef
b           ranef    ranef     ----        ranef         ranef
var(beta)   vcov     vcov      ----        vcov          vcov
var(b)      VarCorr  VarCorr   ranef       VarCorr       VarCorr

  Notice that "ranef" in the prior release was out of step with everyone
else, which became obvious when working on lmekin. Also, the constructs
"fit$coef$random" and "fit$coef$fixed" no longer work due to some
internal rearrangments.


2. The lmekin function is now part of the package.  This is the last
portion of the older "kinship" package to be updated; I will be talking
to the current maintainer of kinship about making it depreciated. This
version of lmekin is both more stable and more capable than the old
one.  
  Why use lmekin?  It can fit most of the models that lme or lmer can,
and uses the same formula notation as lmer.  However, if lmer can fit
your needs I see no good reason to use lmekin and lots of reasons not to
use it (plots, residuals, documentation, ...).  The lmekin routine can
also fit models with complex user-specified variance structures,
however, of which familial kinship correlation matrices are the
principle example.

3. More vignettes.

Terry Therneau



More information about the R-help mailing list