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

View source: R/logLik.ECFOCF.R

logLik.ECFOCFR Documentation

Return Log Likelihood of a fit done using fitCF

Description

Return Log Likelihood of a fit generated by fitCF.

Usage

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

Arguments

object

A result file generated by fitCF

...

Not used

Details

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

Value

The Log Likelihood value for the fitted model with data

Author(s)

Marc Girondot

See Also

Other Model of Clutch Frequency: ECFOCF_f(), ECFOCF_full(), TableECFOCF(), fitCF(), fitCF_MHmcmc(), fitCF_MHmcmc_p(), generateCF(), 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["final"]])
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)

## End(Not run)

phenology documentation built on Sept. 11, 2024, 6:07 p.m.