Computes the Fisher information matrix by Monte Carlo simulation. Uses the outer product of scores.
# S3 method for class 'likelihood_model'
fim(model, ...)Function that takes (theta, n_obs, n_samples, ...) and returns FIM matrix
This default requires the model to implement rdata and score methods.