09-23-2015 07:50 PM
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;
output out=outme marginal;
proc means data=outme n mean;
var Meff_P2_SADC Meff_P2_ROA Meff_P2_AGE;
title 'Average of the Individual Marginal Effects';
How can I calculate the standard errors of the marginal effects?
09-24-2015 11:07 AM - edited 09-24-2015 11:24 AM
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;
proc nlmixed data=remiss;
model remiss ~ binary(p);
predict pdf('NORMAL',lp)*bblast out=outnl;
03-07-2017 10:26 AM
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;
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.