Cost-Effective Acceptability Curve

I prepared this plot for:

Filion KB, Xie X, van der Avoort CJ, Dendukuri N, Brophy JM. Microvolt T-wave alternans and the selective use of implantable cardioverter defibrillators for primary prevention: a cost-effectiveness study. Int J Technol Assess Health Care 2009; 25(2):151-60 Link: http://www.ncbi.nlm.nih.gov/pubmed/19331707 (The abstract of this article also published in Can J Cardiol 2008; 24(Suppl E): 330E)

 

proc import

    datafile="D:\TAU project\2008_Mtwa\MTWA-ICD(Jan07).xls"

    out=bb

    DBMS=xls REPLACE;

    getnames=yes;

      RANGE="Results_new1$a1602:d1623";

run;

 

 

data dd;

    set bb;

      label Med= 'Medical therapy'

          spp='Selected Primary Prevention'

           pp='Primary Prevention';

run;

 

 

proc format;

   value WTP     20000='20K'

                 40000='40K'

                 60000='60K'

                 80000='80K'

                 100000='100K'

                 120000='120K'

                     140000='140K'

                 160000='160K'

                     180000='180K'

                 200000='200K';

run;

 

 

proc gplot data=dd;

      plot (med spp pp)*wtp/  overlay  legend=legend1 haxis=axis1 vaxis=axis2;

    symbol1  i=join ci=blue   v=star         l=7 h=1.5 w=2;

    symbol2  i=join ci=black  v=dot        l=1 h=1.5 w=2 ;

      symbol3  i=join ci=red    v=square     l=3 h=1.5 w=2;

     format  wtp wtp.;

    axis1 label=(h=2 f=Arial  'Threshold willingness to pay (CAD)') width=1 length=100

          order=(0 to 200000 by 20000) value=(f='Times New Roman' h=16pt) minor=none;

    axis2 label=(h=2  f=Arial  a=90 'Probability cost-effective') width=1 length=35

          order=(0 to 1 by 0.1 )  value=(f='Times New Roman' h=16pt) minor=none;

      legend1 down=3 position=(top right inside)label=(f='Times New Roman' h=16pt 'Strategy')

          frame value=(f='Times New Roman' h=16pt) mode=share;

run;

quit;