Extracting Profiled Restricted Log Likelihood

Share:

Description

Extracting maximized profiled restricted log likelihood (PREML) from a varComp object

Usage

1
2
## S3 method for class 'varComp'
logLik(object, ...)

Arguments

object

A varComp object

...

Not used

Details

The likelihood value is the profiled restricted log likelihood. The actual value depends on the residual contrast being chosen, but the location of maximum does not. Thus it is only comparable when the same residual contrast is used. This value is not to be compared to results reported using functions in other packages, e.g., nlme::lme or SAS. The degree of freedom reported is the number of variance components not on the boundary of parameter space, whereas the decision is based on object$control$boundary.eps.

Value

A numeric value of class logLik.

Author(s)

Long Qu

References

Qu L, Guennel T, Marshall SL. (2013) Linear Score Tests for Variance Components in Linear Mixed Models and Applications to Genetic Association Studies. Biometrics, Volume 69, Issue 4, pages 883–892.

See Also

varComp.

Examples

1
2
3
4
library(nlme)
data(Oxide)
vcf = varComp(Thickness~Source, Oxide, ~Lot/Wafer)
logLik(vcf)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.