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