ranef {ggmix} | R Documentation |
Generic function for extracting the random effects. This is the same generic (and same name) defined in the nlme and lme4 package.
ranef(object, ...)
random.effects(object, ...)
## Default S3 method:
random.effects(object, ...)
## Default S3 method:
ranef(object, ...)
## S3 method for class 'ggmix_gic'
ranef(object, s = "lambda.min", ...)
object |
any fitted model object from which random effects estimates can
be extracted. Currently supports "ggmix_gic" objects outputted by the
|
... |
other parameters. currently ignored |
s |
Value(s) of the penalty parameter |
For objects of class "ggmix_gic", this function returns the subject-specific random effect value for the model which minimizes the GIC using the maximum a posteriori principle
a numeric vector of length equal to the number of observations of subject-specific random effects
data("admixed")
fit <- ggmix(x = admixed$xtrain, y = admixed$ytrain,
kinship = admixed$kin_train)
gicfit <- gic(fit)
# random effect at selected value of lambda
plot(ggmix::ranef(gicfit))
# random effects at specific values of lambda
head(ggmix::ranef(gicfit, s = c(0.1,0.2)))