topic Re: PROC CAUSALMED warning messages in Statistical Procedures

PROC CAUSALMED warning messages

mepidemiology — Fri, 01 Apr 2022 14:49:33 GMT

Hello!

I am trying to perform a mediation analysis using PROC CAUSALMED and the outcome is a binary variable. I've set dist=binomial and link=log, and I've tried dist=binary and link=log. Each time I receive these warnings:

WARNING: At least one element of the gradient for the outcome model is greater than 1e-3.
WARNING: The Hessian for the outcome model has been ridged with a maximum value of 524288.
Standard errors might not be accurate.

Does anyone know a) what these mean and b) how I can fix them? Thank you!!

Re: PROC CAUSALMED warning messages

StatDave — Fri, 01 Apr 2022 15:01:24 GMT

Both messages indicate numerical problems fitting the model. This is extremely common when fitting a binary response model using a link function other than the default logit link. The log link does not guarantee that the estimated values are in the valid range for a binomial mean - between 0 and 1. So, use the logit link and see if that fixes the problem.

Re: PROC CAUSALMED warning messages

SAS_Rob — Fri, 01 Apr 2022 15:01:56 GMT

Without seeing the LOG with the CAUSALMED code in it, it is hard to be specific, but in general the message means that there are convergence issues with the outcome model. The usual approach to solve the issue would be to simplify the model by either removing covariates from the COVAR statement or by removing the interaction between the mediator and treatment.

Re: PROC CAUSALMED warning messages

ballardw — Fri, 01 Apr 2022 15:02:50 GMT

You should include the entire code of the procedure call along with all of the messages from the log. IF you tried multiple sets of options then include them as well.

The folks that know more about this procedure may see some other option(s) that might be contributing to the issue. Plus then you don't get suggestions for options that you are already using.

Re: PROC CAUSALMED warning messages

mepidemiology — Tue, 05 Apr 2022 14:22:12 GMT

Thank you! I have tried this and while it does work, it only provides the Odds Ratio. I am trying to get the Risk Ratio (and link=log was suggested) because I'm trying to compare the output from R to SAS, but this is the problem I'm encountering.

Re: PROC CAUSALMED warning messages

SAS_Rob — Tue, 05 Apr 2022 14:27:54 GMT

When you use the LOG link then they are actually relative risks and not odds ratios. They are mislabeled in earlier releases, but the more recent releases reflects the proper labeling, but the label will say Odds Ratios until then.

Re: PROC CAUSALMED warning messages

llllllll1 — Sat, 07 May 2022 11:32:50 GMT

Hello, I'm handling the same problem recently. Since the outcome variable is not rare, the log link was suggested in SAS help.

I've read these solutions you suggested, but as you can see, I've no reason to reduce those covariates, and there is no interaction option initially, so what else can I do to treat this problem, or if I can ignore the warnings?

code is as follows:

proc casusalmed data=data all;

class X(ref=first) M(ref=last) Y(ref=first) C(ref=first);

model Y=X M /link=log;

mediator M=X;

covar c;

run;