optim(par = as.numeric(c(0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1)),
+ fn = loglik_function,
+ dat = reported_data,
+ control = list(fnscale=-1))
# Where can I start to look Error in eval(substitute(expr), data, enclos = parent.frame()) : object 'S' not found