Make likelihood() an S3 generic

[jamesmbaazam]

The broader question is whether it would not be clearer to define S3 methods for numeric, epichains and epichains_summary where

likelihood.epichains <- function(chains, ...) {
  likelihood(summary(chains, ...))
}

likelihood.epichains_summary <- function(chains, stat_max = Inf, ...) {
  stat_max_from_sim <- attr(chains, "stat_max")
  if (missing(stat_max)) {
    stat_max <- stat_max_from_sim
  } else if (stat_max != stat_max_from_sim) {
    warning(`stat_max` is different from the one used in the simulations given in `chains`")
  }
  likelihood(as.numeric(chains, stat_max, ...))
}

and the current likelihood() becomes likelihood.numeric()

Originally posted by @sbfnk in #213 (comment)

[jamesmbaazam]

This blog post by Hugo is relevant to this issue: https://epiverse-trace.github.io/posts/s3-generic/.