mvmeta (version 0.4.7)

logLik.mvmeta: Extract Log-Likelihood from mvmeta Objects

Description

This method function returns the log-likelihood for fitted univariate or multivariate meta-analytical models represented in objects of class "mvmeta".

Usage

"logLik"(object, ...)

Arguments

object

an object of class "mvmeta".

...

further arguments passed to or from other methods.

Value

A numeric scalar of class "logLik" with attributes, providing the (restricted) log likelihood of the model. Attributes correspond to the component df of mvmeta objects, namely the following scalars: nall (number of observations used for estimation, excluding missing values), nobs (equal to nall, minus the number of fixed-effects coefficients for REML models, fixed (number of estimated fixed-effects coefficients), random (number of estimated (co)variance terms).

See Also

See the default method logLik. See mvmeta-package for an overview of the package and modelling framework.

Examples

# RUN THE MODEL 
model <- mvmeta(cbind(PD,AL)~pubyear,S=berkey98[5:7],data=berkey98)

# LOG-LIKELIHOOD
ll <- logLik(model)
ll
attributes(ll)

# AIC and BIC
AIC(model)
BIC(model)