topic Re: How to calculate standard errors of marginal effects in Statistical Procedures

How to calculate standard errors of marginal effects

burtsm — Wed, 23 Sep 2015 23:50:49 GMT

I am running a probit regression and want to report marginal effects and the standard errors of the marginal effects. Here is my code for the regression and marginal effects:

proc qlim data=sa.Ipumsi_sa_reg;
    model CHDIED=SADC ROA AGE / DISCRETE(d=probit);
   where OSEAS=0 and 20 le AGE le 35;
   weight WTPER;
    output out=outme marginal;
run;
quit;

proc means data=outme n mean;
      var Meff_P2_SADC Meff_P2_ROA Meff_P2_AGE;
      title 'Average of the Individual Marginal Effects';
   run;
   quit;

How can I calculate the standard errors of the marginal effects?

Thank you.

Re: How to calculate standard errors of marginal effects

StatDave — Thu, 24 Sep 2015 15:24:57 GMT

You can get the estimated marginal effects and their standard errors by fitting the model in PROC NLMIXED and using the PREDICT statement as shown in this note on marginal effects. For example, these statements use QLIM and NLMIXED to fit the same probit model to the cancer remission data shown in the first example in the LOGISTIC documentation. The OUTNL output data set from PROC NLMIXED contains the estimated marginal effects (Pred) and the standard errors (StdErrPred).

proc qlim data=remiss;
model remiss=blast smear / discrete;
output out=outqlim marginal;
run;

proc nlmixed data=remiss;
lp=int+bblast*blast+bsmear*smear;
p=probnorm(lp);
model remiss ~ binary(p);
predict pdf('NORMAL',lp)*bblast out=outnl;
run;

Re: How to calculate standard errors of marginal effects

SlutskyFan — Tue, 07 Mar 2017 15:26:50 GMT

This example seems to give SE and CI for each value at which a ME is calculated. If we are interested in average marginal effects, other examples use PROC MEANS in combination with PROC QLIM to average all of the MEs to get the desired result. (https://support.sas.com/rnd/app/ets/examples/margeff/index.html) Would we get the correct standard errors and confidence intervals simiply using proc means as such:

PROC MEANS DATA=OUTNL MEAN STDERR CLM;

VAR PRED;

RUN;

Or is there another way to leverage NLMIXED in a such a way to produce average marginal effects and their standard errors? I'm thinking to get standard errors for marginal effects calculated this way some sort of loop/bootstrap approach might be necessary?

I know this is an old discussion but looking for a solution today and I ran across this string.

Thanks.

Re: How to calculate standard errors of marginal effects

StatDave — Thu, 18 Oct 2018 13:01:00 GMT

The average marginal effect and its standard error can be computed with the Margins macro. The standard error if determined using the delta method.