Skip to contents

Log-Likelihood of a VLMC

Source: R/vlmc_likelihood.R, R/vlmc_likelihood_cpp.R
logLik.vlmc.Rd

This function evaluates the log-likelihood of a VLMC fitted on a discrete time series.

Usage

# S3 method for vlmc
logLik(object, initial = c("truncated", "specific", "extended"), ...)

# S3 method for vlmc_cpp
logLik(object, initial = c("truncated", "specific", "extended"), ...)

Arguments

object

the vlmc representation.

initial

specifies the likelihood function, more precisely the way the first few observations for which contexts cannot be calculated are integrated in the likelihood. Defaults to "truncated". See loglikelihood() for details.

...

additional parameters for logLik.

Value

an object of class logLik. This is a number, the log-likelihood of the (CO)VLMC with the following attributes:

  • df: the number of parameters used by the VLMC for this likelihood calculation

  • nobs: the number of observations included in this likelihood calculation

  • initial: the value of the initial parameter used to compute this likelihood

See also

loglikelihood()

Examples

pc <- powerconsumption[powerconsumption$week == 5, ]
breaks <- c(
  0,
  median(powerconsumption$active_power, na.rm = TRUE),
  max(powerconsumption$active_power, na.rm = TRUE)
)
labels <- c(0, 1)
dts <- cut(pc$active_power, breaks = breaks, labels = labels)
m_nocovariate <- vlmc(dts)
ll <- logLik(m_nocovariate)
ll
#> 'log Lik.' -207.9581 (df=69)
attributes(ll)
#> $nobs
#> [1] 945
#> 
#> $df
#> [1] 69
#> 
#> $initial
#> [1] "truncated"
#> 
#> $class
#> [1] "logLik"
#>