topic Re: Trouble with the %CIF macro in SAS Procedures

Trouble with the %CIF macro

ChrisMac — Sat, 12 Sep 2015 14:53:27 GMT

Hi everyone,

I'm having a bit trouble with the %CIF macro. I can't change the title even though I'm doing it correct (I think).

This statement works fine:
%CIF (DATA=final, TIME=follow_up, STATUS=hypertension, EVENT=1, censored=0, GROUP=case);

When adding TITLE=xxxx, the title won't change, although the cumulative incidence figure looks like the above:

%CIF (DATA=final, TIME=follow_up, STATUS=hypertension, EVENT=1, censored=0, GROUP=case, TITLE=xxxx);

Hope, you guys can help!

Re: Trouble with the %CIF macro

Tom — Sat, 12 Sep 2015 15:11:54 GMT

I think that the TITLE= parameter only applies to a limited part of the output.

TITLE=value
    specifies the first title for the cumulative incidence function plot. Do not use quotation marks.

TITLE2=value
    specifies the second title for the cumulative incidence function plot. Do not use quotation marks.

But even for that it looks like the TITLE statements are using single quotes instead of double quotes. You should report to SAS.

title '&title';
title2 '&title2';

Re: Trouble with the %CIF macro

Reeza — Sun, 13 Sep 2015 23:57:45 GMT

Not sure if this is useful, but If you have SAS 9.4 you can use Proc PHREG for Fine and Gray's Competing Risks model in SAS now.

http://www.sas.com/content/dam/SAS/en_ca/User%20Group%20Presentations/Vancouver-User-Group/Gondara_CompetingRiskAnalysis_May2015.pdf

Re: Trouble with the %CIF macro

EastwoodDC — Wed, 16 Sep 2015 22:00:24 GMT

I just discovered the same problem, and it works just fine if you change those sinle quotes to double.

Previously I was using an older version of this macro I found online, and I didn't encounter the bug until I use the autocall version.

And I'm new here - where do I report this?

Re: Trouble with the %CIF macro

Reeza — Thu, 17 Sep 2015 01:20:27 GMT

Contact SAS tech support, open a track.

Re: Trouble with the %CIF macro

ChrisMac — Thu, 17 Sep 2015 11:49:25 GMT

Hi,

After running the following command - this time with double quotes - the title stil wont change. Am I doing it wrong?

%CIF (DATA=final, TIME=follow_up, STATUS=hypertension, EVENT=1, censored=0, GROUP=case, TITLE="Survival data");

Re: Trouble with the %CIF macro

Reeza — Thu, 17 Sep 2015 15:37:15 GMT

I believe you need to modify the actual macro code not the macro call. In the code change the single quotes to double quotes.

Re: Trouble with the %CIF macro

Tom — Fri, 18 Sep 2015 10:08:27 GMT

On the old site the link was in the page footer. On this new site click on the Contact link under the Connect column in the page footer. Then on that page click on the support link:

http://support.sas.com/ctx/supportform/createForm

Re: Trouble with the %CIF macro

EastwoodDC — Wed, 06 Apr 2016 17:56:10 GMT

I encountered a different problem with the %CIF macro.

The error occurs if the Cumulative Incidence Function (CIF) ever reaches 1.0, which can occur if there are no competing or censored events in one of the groups. The EXP() function generates an error because it may be given a value larger than IML can evaluate. This occurs in the following code when cifEST[cifloc,2] is very close to 1.0:

lowci[cifloc] = cifEST[cifloc, 2]##exp((-z*cifEST[cifloc, 3]) /(cifEST[cifloc,2]#LOG(cifEST[cifloc,2]))); 
upci[cifloc] = cifEST[cifloc, 2]##exp((z*cifEST[cifloc, 3]) /(cifEST[cifloc,2]#LOG(cifEST[cifloc,2])));

which occurs 3 times in the macro.

This error can be avoided by inserting a MIN() function to prevent an overly large number from being evaluated, as follows:

lowci[cifloc] = cifEST[cifloc, 2]##exp((-z*cifEST[cifloc, 3]) /min(9999999,cifEST[cifloc,2]#LOG(cifEST[cifloc,2]))); 
upci[cifloc] = cifEST[cifloc, 2]##exp((z*cifEST[cifloc, 3]) /min(9999999,cifEST[cifloc,2]#LOG(cifEST[cifloc,2])));

Divide by zero warning can also occur here, but those do not cause the macro to fail.