Define the Log Likelihood function for the observational model

This function is defined by the user and describes the observational likelihood of the data given the model.

Arguments

titre_val

The titre value from the data

titre_est

The model-estimate titre value.

pars

The fitted parameters needed to calculate the log likelihood. These are defined in the prior entry of the observationalModel list.

Value

A functoion that returns the log likelihood value for the COP model.

Details

Add information here.

Author

Your Name

Examples

# Example usage:
obsLogLikelihood = function(titre_val, titre_est, pars) {
   if (titre_val <= log10(40)) {
       ll <- pcauchy(log10(40), titre_est, pars[1], log.p = TRUE)
   } else {
       ll <- dcauchy(titre_val, titre_est, pars[1], log = TRUE)
   }
   ll
}