phenology (version 7.2)

logLik.ECFOCF: Return Log Likelihood of a fit done using fitCF

Description

Return Log Likelihood of a fit generated by fitCF.

Usage

# S3 method for ECFOCF
logLik(object, ...)

Arguments

object

A result file generated by fitCF

...

Not used

Value

The Log Likelihood value for the fitted model with data

Details

logLik.ECFOCF return Log Likelihood of a fit done using fitCF

See Also

Other Model of Clutch Frequency: ECFOCF_f, TableECFOCF, fitCF, lnLCF, plot.ECFOCF, plot.TableECFOCF

Examples

# NOT RUN {
library(phenology)
# Example
ECFOCF_2002 <- TableECFOCF(MarineTurtles_2002, date0=as.Date("2002-01-01"))
par <- c(mu = 2.6404831115214353, 
        size = 0.69362774786433479, 
        mu_season = 12.6404831115214353, 
        size_season = 1.69362774786433479, 
        a2=0)
fp <- rep(0, dim(ECFOCF_2002)[3])
names(fp) <- paste0("p.", formatC(1:(dim(ECFOCF_2002)[3]), width=2, flag="0"))
par <- c(par, fp[attributes(ECFOCF_2002)$table["begin"]:attributes(ECFOCF_2002)$table["end"]])
fixed.parameters <- c(a1=Inf, p=-Inf)

 lnLCF(x=par, data=ECFOCF_2002, fixed.parameters=fixed.parameters)
 
 o_mu1season1a2p <- fitCF(x=par, fixed.parameters=fixed.parameters, 
                          data=ECFOCF_2002)
 
logLik(o_mu1season1a2p)
AIC(o_mu1season1a2p)
# }