The function is a wrapper around derive_vcov_mnl
and
derive_vcov_rpl
and calculates the variance-covariance matrix
of the specified model and design given the priors.
derive_vcov(design_env, model)
The variance covariance matrix. If the Fisher information matrix is singular, then return NULL