R/logLik.R

In flexrsurv: Flexible Relative Survival Analysis

## log-likelihood for flexrsurv objects
logLik.flexrsurv <- function(object, ...)
{
	if(!missing(...)) warning("extra arguments discarded")
	val <- object$loglik
	## Note: weights are not implemented
	attr(val, "nobs") <- object$nobs
	attr(val, "df") <-  length(coef(object))
	class(val) <- "logLik"
	val
}

nobs.flexrsurv <- function(object, ...) object$nobs 


vcov.flexrsurv <- function (object, ...) {
	vmat <- object$var
	vname <- names(object$coefficients)
	dimnames(vmat) <- list(vname, vname)
	vmat
}

extractAIC.flexrsurv <- function (fit, scale = 0, k = 2, ...) 
{
	loglik <- logLik(fit)
	edf <-  attr(loglik, "nobs") - attr(loglik, "df")
	aic <- stats::AIC(fit)
	c(edf, aic + (k - 2) * edf)
}