phenology (version 4.2.4)

logLik.fitRMU: Return Log Likelihood of a fit generated by fitRMU

Description

Return Log Likelihood of a fit generated by fitRMU

Usage

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

Arguments

object

A result file generated by fitRMU

...

Not used

Value

• The Log Likelihood value for the fitted model with data

Details

logLik.fitRMU Return Log Likelihood of a fit for fitRMU

See Also

Other Fill gaps in RMU: fitRMU_MHmcmc_p; fitRMU_MHmcmc; fitRMU; plot.fitRMU

Examples

library(phenology)
RMU.name.AtlanticW <- data.frame(mean=c("Yalimapo.French.Guiana",
                                         "Galibi.Suriname",
                                         "Irakumpapy.French.Guiana"),
                                 se=c("se_Yalimapo.French.Guiana",
                                      "se_Galibi.Suriname",
                                      "se_Irakumpapy.French.Guiana"))
data.AtlanticW <- data.frame(Year=c(1990:2000),
      Yalimapo.French.Guiana=c(2076, 2765, 2890, 2678, NA,
                               6542, 5678, 1243, NA, 1566, 1566),
      se_Yalimapo.French.Guiana=c(123.2, 27.7, 62.5, 126, NA,
                                 230, 129, 167, NA, 145, 20),
      Galibi.Suriname=c(276, 275, 290, NA, 267,
                       542, 678, NA, 243, 156, 123),
      se_Galibi.Suriname=c(22.3, 34.2, 23.2, NA, 23.2,
                           4.3, 2.3, NA, 10.3, 10.1, 8.9),
      Irakumpapy.French.Guiana=c(1076, 1765, 1390, 1678, NA,
                               3542, 2678, 243, NA, 566, 566),
      se_Irakumpapy.French.Guiana=c(23.2, 29.7, 22.5, 226, NA,
                                 130, 29, 67, NA, 15, 20))

cst <- fitRMU(data=data.AtlanticW, RMU.name=RMU.name.AtlanticW,
               colname.year="Year", model.trend="Constant",
               model.SD="Zero")
logLik(cst)
AIC(cst)