R/logLik.R

## Retunrs the log likelihood of a VLMC model

setMethod("logLik", "PSTf", function(object) {

	n <- nrow(object@data)

	message(" [>] model fitted to ", n, " sequence(s) - ", nobs(object), " symbols")

	pstsum <- summary(object)
	res <- object@logLik

	class(res) <- "logLik"
	attr(res, "df") <- pstsum@freepar

	return(res)
}
)
 
 

Try the PST package in your browser

Any scripts or data that you put into this service are public.

PST documentation built on May 2, 2019, 3:43 a.m.