Test the change of margins for multilogit and logit

The change uses the averages of the margins and std_err evaluation,
instead of evaluating at means.

R package erer does not produce the same result as Stata. It has
different behaviors.

Remove erer from the R script, using analytical derivative directly.