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

View source: R/logLik.fitRMU.R

logLik.fitRMUR Documentation

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

Details

logLik.fitRMU Return Log Likelihood of a fit for fitRMU

Value

The Log Likelihood value for the fitted model with data

Author(s)

Marc Girondot

See Also

Other Fill gaps in RMU: CI.RMU(), fitRMU(), fitRMU_MHmcmc(), fitRMU_MHmcmc_p(), plot.fitRMU()

Examples

## Not run: 
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)

## End(Not run)

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